kernelSin
Compute the sine of a double-precision floating-point number on
[-π/4, π/4]
.
Usage
var kernelSin = require( '@stdlib/math/base/special/kernel-sin' );
kernelSin( x, y )
Computes the sine of a double-precision floating-point number on [-π/4, π/4]
.
var v = kernelSin( 0.0, 0.0 );
// returns ~0.0
v = kernelSin( 3.141592653589793/6.0, 0.0 );
// returns ~0.5
v = kernelSin( 0.619, 9.279e-18 );
// returns ~0.58
v = kernelSin( NaN, 0.0 );
// returns NaN
v = kernelSin( 3.0, NaN );
// returns NaN
v = kernelSin( NaN, NaN );
// returns NaN
Notes
For increased accuracy, the number for which the sine should be evaluated can be supplied as a double-double number (i.e., a non-evaluated sum of two double-precision floating-point numbers
x
andy
).As components of a double-double number, the two double-precision floating-point numbers
x
andy
must satisfywhere
ulp
stands for units in the last place.
Examples
var linspace = require( '@stdlib/array/base/linspace' );
var PI = require( '@stdlib/constants/float64/pi' );
var kernelSin = require( '@stdlib/math/base/special/kernel-sin' );
var x = linspace( -PI/4.0, PI/4.0, 100 );
var i;
for ( i = 0; i < x.length; i++ ) {
console.log( 'sine(%d) = %d', x[ i ], kernelSin( x[ i ], 0.0 ) );
}
C APIs
Usage
#include "stdlib/math/base/special/kernel_sin.h"
stdlib_base_kernel_sin( x, y )
Computes the sine of a double-precision floating-point number on [-π/4, π/4]
.
double out = stdlib_base_kernel_sin( 0.0, 0.0 );
// returns ~0.0
out = stdlib_base_kernel_sin( 0.619, 9.279e-18 );
// returns ~0.58
The function accepts the following arguments:
- x:
[in] double
input value (in radians, assumed to be bounded by~pi/4
in magnitude). - y:
[in] double
tail ofx
.
double stdlib_base_kernel_sin( const double x, const double y );
Notes
- For increased accuracy, the number for which the sine should be evaluated can be supplied as a double-double number (i.e., a non-evaluated sum of two double-precision floating-point numbers
x
andy
).
Examples
#include "stdlib/math/base/special/kernel_sin.h"
#include <stdio.h>
int main( void ) {
const double x[] = { -0.7853981633974483, -0.6108652381980153, -0.4363323129985824, -0.26179938779914946, -0.08726646259971649, 0.08726646259971649, 0.26179938779914935, 0.43633231299858233, 0.6108652381980153, 0.7853981633974483 };
double out;
int i;
for ( i = 0; i < 10; i++ ) {
out = stdlib_base_kernel_sin( x[ i ], 0.0 );
printf ( "x: %lf, y: %lf, out: %lf\n", x[ i ], 0.0, out );
}
}