Fill a double-precision floating-point strided array with linearly spaced values over a specified interval.
var dlinspace = require( '@stdlib/blas/ext/base/dlinspace' );Fills a double-precision floating-point strided array with linearly spaced values over a specified interval.
var Float64Array = require( '@stdlib/array/float64' );
var x = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
dlinspace( x.length, 0.0, 7.0, true, x, 1 );
// x => <Float64Array>[ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0 ]The function has the following parameters:
- N: number of indexed elements.
- start: start of interval.
- stop: end of interval.
- endpoint: boolean indicating whether to include the
stopvalue when writing values to the input array. Iftrue, the input array is filled with evenly spaced values over the closed interval[start, stop]. Iffalse, the input array is filled with evenly spaced values over the half-open interval[start, stop). - x: input
Float64Array. - strideX: stride length.
The N and stride parameters determine which elements in the strided array are accessed at runtime. For example, to fill every other element:
var Float64Array = require( '@stdlib/array/float64' );
var x = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
dlinspace( 4, 1.0, 5.0, false, x, 2 );
// x => <Float64Array>[ 1.0, 0.0, 2.0, 0.0, 3.0, 0.0, 4.0, 0.0 ]Note that indexing is relative to the first index. To introduce an offset, use typed array views.
var Float64Array = require( '@stdlib/array/float64' );
// Initial array...
var x0 = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
// Create an offset view...
var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
// Fill every other element...
dlinspace( 3, 1.0, 3.0, true, x1, 2 );
// x0 => <Float64Array>[ 0.0, 1.0, 0.0, 2.0, 0.0, 3.0 ]Fills a double-precision floating-point strided array with linearly spaced values over a specified interval using alternative indexing semantics.
var Float64Array = require( '@stdlib/array/float64' );
var x = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
dlinspace.ndarray( x.length, 0.0, 7.0, true, x, 1, 0 );
// x => <Float64Array>[ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0 ]The function has the following additional parameters:
- offsetX: starting index.
While typed array views mandate a view offset based on the underlying buffer, the offset parameter supports indexing semantics based on a starting index. For example, to access only the last three elements:
var Float64Array = require( '@stdlib/array/float64' );
var x = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
dlinspace.ndarray( 3, 1.0, 3.0, true, x, 1, x.length-3 );
// x => <Float64Array>[ 0.0, 0.0, 0.0, 1.0, 2.0, 3.0 ]-
Let
Mbe the number of generated values (which is eitherNorN+1depending on whetherendpointistrueorfalse, respectively). The spacing between values is thus given byΔ = (stop-start)/(M-1) -
When the number of generated values is greater than
1andendpointistrue, the set of values written to a provided input array is guaranteed to include thestartandstopvalues. Beware, however, that values betweenstartandstopare subject to floating-point rounding errors. Hence,var Float64Array = require( '@stdlib/array/float64' ); var x = new Float64Array( [ 0.0, 0.0, 0.0 ] ); dlinspace( 3, 0.0, 1.0, true, x, 1 ); // x => <Float64Array>[ 0.0, ~0.5, 1.0 ]
where
x[1]is only guaranteed to be approximately equal to0.5. -
When
N = 1andendpointisfalse, only thestartvalue is written to a provided input array. WhenN = 1andendpointistrue, only thestopvalue is written to a provided input array. -
If
start < stop, values are written to a provided input array in ascending order; otherwise, values are written in descending order. -
If
N <= 0, both functions returnxunchanged.
var discreteUniform = require( '@stdlib/random/array/discrete-uniform' );
var dlinspace = require( '@stdlib/blas/ext/base/dlinspace' );
var x = discreteUniform( 10, -100, 100, {
'dtype': 'float64'
});
console.log( x );
dlinspace( x.length, 0.0, 10.0, true, x, 1 );
console.log( x );#include "stdlib/blas/ext/base/dlinspace.h"Fills a double-precision floating-point strided array with linearly spaced values over a specified interval.
#include <stdbool.h>
double x[] = { 0.0, 0.0, 0.0, 0.0 };
stdlib_strided_dlinspace( 4, 1.0, 5.0, true, x, 1 );The function accepts the following arguments:
- N:
[in] CBLAS_INTnumber of indexed elements. - start:
[in] doublestart of interval. - stop:
[in] doubleend of interval. - endpoint:
[in] boolboolean indicating whether to include thestopvalue when writing values to the input array. Iftrue, the input array is filled with evenly spaced values over the closed interval[start, stop]. Iffalse, the input array is filled with evenly spaced values over the half-open interval[start, stop). - X:
[out] double*input array. - strideX:
[in] CBLAS_INTstride length.
void stdlib_strided_dlinspace( const CBLAS_INT N, const double start, const double stop, const bool endpoint, double *X, const CBLAS_INT strideX );Fills a double-precision floating-point strided array with linearly spaced values over a specified interval using alternative indexing semantics.
#include <stdbool.h>
double x[] = { 0.0, 0.0, 0.0, 0.0 };
stdlib_strided_dlinspace_ndarray( 4, 1.0, 5.0, true, x, 1, 0 );The function accepts the following arguments:
- N:
[in] CBLAS_INTnumber of indexed elements. - start:
[in] doublestart of interval. - stop:
[in] doubleend of interval. - endpoint:
[in] boolboolean indicating whether to include thestopvalue when writing values to the input array. Iftrue, the input array is filled with evenly spaced values over the closed interval[start, stop]. Iffalse, the input array is filled with evenly spaced values over the half-open interval[start, stop). - X:
[out] double*input array. - strideX:
[in] CBLAS_INTstride length. - offsetX:
[in] CBLAS_INTstarting index.
void stdlib_strided_dlinspace_ndarray( const CBLAS_INT N, const double start, const double stop, const bool endpoint, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX );#include "stdlib/blas/ext/base/dlinspace.h"
#include <stdio.h>
#include <stdbool.h>
int main( void ) {
// Create a strided array:
double x[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
// Specify the number of indexed elements:
const int N = 8;
// Specify a stride:
const int strideX = 1;
// Fill the array:
stdlib_strided_dlinspace( N, 0.0, 10.0, true, x, strideX );
// Print the result:
for ( int i = 0; i < 8; i++ ) {
printf( "x[ %i ] = %lf\n", i, x[ i ] );
}
}