vrsa_log2f: Array single precision base-2 logarithm

• C Prototype:
  • void vrsa_log2f (int n, float *x, float *y)

  • Inputs:

    • int	n	- the number of values in both the input and output arrays.
      float	*x	- pointer to the array of input values.
      float	*y	- pointer to the array of output values.

  • Outputs:
    • The base-2 logarithm of each x value, filled into the y array.
• Fortran Subroutine Interface:
  • SUBROUTINE VRSA_LOG2F(N,X,Y)
  • Inputs:
    • INTEGER N - the number of values in both the input and output arrays.
    • REAL X(N) - array of single precision input values.
  • Outputs:
    • REAL Y(N) - array of base-2 logarithms of input values.
• Notes:
  • vrsa_log2f computes the single precision base-2 logarithm for each element of an array of input arguments.
    This routine accepts an array of single precision input values, computes the base-2 log for each input value, and stores the result in the array pointed to by the y pointer input. It is the responsibility of the calling program to allocate/deallocate enough storage for the output array. This is a relaxed version of log2f, suitable for use with fastmath compiler flags or applications not requiring full error handling. Denormal inputs may produce unpredictable results. Special case inputs produce C99 return values. The routine is accurate to better than 1 ulp over the valid input range.
• Special case return values:

  • Input	Output
    +/- 0	-infinity
    negative	QNaN
    QNaN	same QNaN
    SNaN	same NaN converted to QNaN
    +infinity	+infinity
    -infinity	QNaN

• Performance:
  • 29 cycles per value for valid inputs, n = 24, longer for input values very close to 1.0.