From 7e3e2917e2ed5d3d4df09dacccfdeee7fb31b072 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Wed, 24 Jun 2026 14:08:00 +0500 Subject: [PATCH 1/2] feat: add blas/ext/base/clogspace --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown_pkg_readmes status: passed - task: lint_markdown_docs status: na - task: lint_markdown status: na - task: lint_package_json status: passed - task: lint_repl_help status: passed - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: passed - task: lint_javascript_tests status: passed - task: lint_javascript_benchmarks status: passed - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: missing_dependencies - task: lint_c_benchmarks status: missing_dependencies - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: passed - task: lint_license_headers status: passed --- --- .../@stdlib/blas/ext/base/clogspace/README.md | 360 ++++++++++++++++++ .../ext/base/clogspace/benchmark/benchmark.js | 112 ++++++ .../clogspace/benchmark/benchmark.native.js | 117 ++++++ .../clogspace/benchmark/benchmark.ndarray.js | 112 ++++++ .../benchmark/benchmark.ndarray.native.js | 117 ++++++ .../ext/base/clogspace/benchmark/c/Makefile | 146 +++++++ .../clogspace/benchmark/c/benchmark.length.c | 205 ++++++++++ .../blas/ext/base/clogspace/binding.gyp | 170 +++++++++ .../blas/ext/base/clogspace/docs/repl.txt | 148 +++++++ .../ext/base/clogspace/docs/types/index.d.ts | 149 ++++++++ .../ext/base/clogspace/docs/types/test.ts | 317 +++++++++++++++ .../ext/base/clogspace/examples/c/Makefile | 146 +++++++ .../ext/base/clogspace/examples/c/example.c | 41 ++ .../blas/ext/base/clogspace/examples/index.js | 32 ++ .../blas/ext/base/clogspace/include.gypi | 53 +++ .../include/stdlib/blas/ext/base/clogspace.h | 47 +++ .../blas/ext/base/clogspace/lib/clogspace.js | 72 ++++ .../base/clogspace/lib/clogspace.native.js | 74 ++++ .../blas/ext/base/clogspace/lib/index.js | 76 ++++ .../blas/ext/base/clogspace/lib/main.js | 35 ++ .../blas/ext/base/clogspace/lib/native.js | 35 ++ .../blas/ext/base/clogspace/lib/ndarray.js | 151 ++++++++ .../ext/base/clogspace/lib/ndarray.native.js | 75 ++++ .../blas/ext/base/clogspace/manifest.json | 102 +++++ .../blas/ext/base/clogspace/package.json | 82 ++++ .../blas/ext/base/clogspace/src/Makefile | 70 ++++ .../blas/ext/base/clogspace/src/addon.c | 72 ++++ .../blas/ext/base/clogspace/src/main.c | 129 +++++++ .../ext/base/clogspace/test/test.clogspace.js | 266 +++++++++++++ .../clogspace/test/test.clogspace.native.js | 275 +++++++++++++ .../blas/ext/base/clogspace/test/test.js | 82 ++++ .../ext/base/clogspace/test/test.ndarray.js | 263 +++++++++++++ .../clogspace/test/test.ndarray.native.js | 272 +++++++++++++ 33 files changed, 4403 insertions(+) create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/benchmark.length.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/binding.gyp create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/index.d.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/test.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/include.gypi create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/include/stdlib/blas/ext/base/clogspace.h create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/main.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/src/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/src/addon.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.native.js diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md b/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md new file mode 100644 index 000000000000..db775bf7e9d7 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md @@ -0,0 +1,360 @@ + + +# clogspace + +> Fill a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. + +
+ +## Usage + +```javascript +var clogspace = require( '@stdlib/blas/ext/base/clogspace' ); +``` + +#### clogspace( N, base, start, stop, endpoint, x, strideX ) + +Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 2.0, 0.0 ); + +clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +// x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +``` + +The function has the following parameters: + +- **N**: number of indexed elements. +- **base**: base of the logarithmic scale. +- **start**: exponent of the starting value, where the starting value is given by `base^start`. +- **stop**: exponent of the final value, where the final value is given by `base^stop`. +- **endpoint**: boolean indicating whether to include the `base^stop` value when writing values to the input array. If `true`, the input array is filled with logarithmically spaced values over the closed interval `[base^start, base^stop]`. If `false`, the input array is filled with logarithmically spaced values over the half-open interval `[base^start, base^stop)`. +- **x**: input [`Complex64Array`][@stdlib/array/complex64]. +- **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: + + + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 3.0, 0.0 ); + +clogspace( 4, 10.0, strt, stp, true, x, 2 ); +// x => [ 1.0, 0.0, 0.0, 0.0, 10.0, 0.0, 0.0, 0.0, 100.0, 0.0, 0.0, 0.0, 1000.0, 0.0, 0.0, 0.0 ] +``` + +Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views. + + + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +// Initial array... +var x0 = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +// Create an offset view... +var x1 = new Complex64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element + +// Fill every other element... +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 2.0, 0.0 ); + +clogspace( 3, 10.0, strt, stp, true, x1, 2 ); +// x0 => [ 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 10.0, 0.0, 0.0, 0.0, 100.0, 0.0 ] +``` + +#### clogspace.ndarray( N, base, start, stop, endpoint, x, strideX, offsetX ) + +Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 2.0, 0.0 ); + +clogspace.ndarray( x.length, 10.0, strt, stp, true, x, 1, 0 ); +// x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +``` + +The function has the following additional parameters: + +- **offsetX**: starting index. + +While [`typed array`][mdn-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: + + + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 2.0, 0.0 ); + +clogspace.ndarray( 3, 10.0, strt, stp, true, x, 1, x.length-3 ); +// x => [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +``` + +
+ + + +
+ +## Notes + +- Let `M` be the number of generated values (which is either `N` or `N+1` depending on whether `endpoint` is `true` or `false`, respectively). The complex increment between exponents is thus given by + + ```text + Δ = (stop-start)/(M-1) + ``` + + and the generated values are equal to `base^(start+Δ*i)` for `i = 0, 1, ..., M-1`, where exponentiation is performed in the complex sense. + +- When the number of generated values is greater than `1` and `endpoint` is `true`, the set of values written to a provided input array is guaranteed to include the `base^start` and `base^stop` values. Beware, however, that values between `base^start` and `base^stop` are subject to floating-point rounding errors. Hence, + + ```javascript + var Complex64Array = require( '@stdlib/array/complex64' ); + var Complex64 = require( '@stdlib/complex/float32/ctor' ); + + var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + + var strt = new Complex64( 0.0, 0.0 ); + var stp = new Complex64( 1.0, 0.0 ); + + clogspace( 3, 10.0, strt, stp, true, x, 1 ); + // x => [ 1.0, 0.0, ~3.162, 0.0, 10.0, 0.0 ] + ``` + + where `x[1]` is only guaranteed to be approximately equal to the square root of `10`. + +- When `N = 1` and `endpoint` is `false`, only the `base^start` value is written to a provided input array. When `N = 1` and `endpoint` is `true`, only the `base^stop` value is written to a provided input array. + +- Generated values follow a logarithmic spiral in the complex plane, where imaginary components of `start` and `stop` control phase rotation. If `Re(start) < Re(stop)`, the magnitudes of the generated values are in ascending order; otherwise, the magnitudes are in descending order. + +- If `N <= 0`, both functions return `x` unchanged. + +
+ + + +
+ +## Examples + + + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var logEach = require( '@stdlib/console/log-each' ); +var clogspace = require( '@stdlib/blas/ext/base/clogspace' ); + +var x = new Complex64Array( 6 ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 5.0, 0.0 ); + +clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +logEach( '%s', x ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/ext/base/clogspace.h" +``` + +#### stdlib_strided_clogspace( N, base, start, stop, endpoint, \*X, strideX ) + +Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. + +```c +#include "stdlib/complex/float32/ctor.h" +#include + +float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + +stdlib_strided_clogspace( 4, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 3.0f, 0.0f ), true, (stdlib_complex64_t *)x, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **base**: `[in] float` base of the logarithmic scale. +- **start**: `[in] stdlib_complex64_t` exponent of the starting value. +- **stop**: `[in] stdlib_complex64_t` exponent of the final value. +- **endpoint**: `[in] bool` boolean indicating whether to include the `base^stop` value when writing values to the input array. If `true`, the input array is filled with logarithmically spaced values over the closed interval `[base^start, base^stop]`. If `false`, the input array is filled with logarithmically spaced values over the half-open interval `[base^start, base^stop)`. +- **X**: `[out] stdlib_complex64_t*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. + +```c +void stdlib_strided_clogspace( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX ); +``` + + + +#### stdlib_strided_clogspace_ndarray( N, base, start, stop, endpoint, \*X, strideX, offsetX ) + +Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. + +```c +#include "stdlib/complex/float32/ctor.h" +#include + +float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + +stdlib_strided_clogspace_ndarray( 4, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 3.0f, 0.0f ), true, (stdlib_complex64_t *)x, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **base**: `[in] float` base of the logarithmic scale. +- **start**: `[in] stdlib_complex64_t` exponent of the starting value. +- **stop**: `[in] stdlib_complex64_t` exponent of the final value. +- **endpoint**: `[in] bool` boolean indicating whether to include the `base^stop` value when writing values to the input array. If `true`, the input array is filled with logarithmically spaced values over the closed interval `[base^start, base^stop]`. If `false`, the input array is filled with logarithmically spaced values over the half-open interval `[base^start, base^stop)`. +- **X**: `[out] stdlib_complex64_t*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. +- **offsetX**: `[in] CBLAS_INT` starting index. + +```c +void stdlib_strided_clogspace_ndarray( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/blas/ext/base/clogspace.h" +#include "stdlib/complex/float32/ctor.h" +#include +#include + +int main( void ) { + // Create a strided array of interleaved real and imaginary components: + float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + + // Specify the number of indexed elements: + const int N = 6; + + // Specify a stride: + const int strideX = 1; + + // Fill the array: + stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 5.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ i*2+1 ] ); + } +} +``` + +
+ + + +
+ + + + + +
+ +
+ + + + + +
+ +[@stdlib/array/complex64]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/complex64 + +[@stdlib/complex/float32/ctor]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/complex/float32/ctor + +[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray + +
+ + diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.js new file mode 100644 index 000000000000..60d07562f06e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.js @@ -0,0 +1,112 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var clogspace = require( './../lib/clogspace.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Create a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var strt; + var xbuf; + var stp; + var x; + + xbuf = uniform( len*2, -100.0, 100.0, options ); + x = new Complex64Array( xbuf.buffer ); + strt = new Complex64( 0.0, 0.0 ); + stp = new Complex64( 9.0, 0.0 ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + clogspace( len, 10.0, strt, stp, true, x, 1 ); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.native.js new file mode 100644 index 000000000000..693ae88906f3 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.native.js @@ -0,0 +1,117 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var clogspace = tryRequire( resolve( __dirname, './../lib/clogspace.native.js' ) ); +var opts = { + 'skip': ( clogspace instanceof Error ) +}; +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Create a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var strt; + var xbuf; + var stp; + var x; + + xbuf = uniform( len*2, -100.0, 100.0, options ); + x = new Complex64Array( xbuf.buffer ); + strt = new Complex64( 0.0, 0.0 ); + stp = new Complex64( 9.0, 0.0 ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + clogspace( len, 10.0, strt, stp, true, x, 1 ); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..3d5ad53bb8e0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.js @@ -0,0 +1,112 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var clogspace = require( './../lib/ndarray.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Create a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var strt; + var xbuf; + var stp; + var x; + + xbuf = uniform( len*2, -100.0, 100.0, options ); + x = new Complex64Array( xbuf.buffer ); + strt = new Complex64( 0.0, 0.0 ); + stp = new Complex64( 9.0, 0.0 ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + clogspace( len, 10.0, strt, stp, true, x, 1, 0 ); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:ndarray:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..681680a1d31c --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/benchmark.ndarray.native.js @@ -0,0 +1,117 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var clogspace = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( clogspace instanceof Error ) +}; +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Create a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var strt; + var xbuf; + var stp; + var x; + + xbuf = uniform( len*2, -100.0, 100.0, options ); + x = new Complex64Array( xbuf.buffer ); + strt = new Complex64( 0.0, 0.0 ); + stp = new Complex64( 9.0, 0.0 ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + clogspace( len, 10.0, strt, stp, true, x, 1, 0 ); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( xbuf[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:ndarray:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/Makefile new file mode 100644 index 000000000000..0756dc7da20a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := benchmark.length.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled benchmarks. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/benchmark.length.c new file mode 100644 index 000000000000..038c18a482e3 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/benchmark/c/benchmark.length.c @@ -0,0 +1,205 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/clogspace.h" +#include "stdlib/complex/float32/ctor.h" +#include +#include +#include +#include +#include +#include + +#define NAME "clogspace" +#define ITERATIONS 10000000 +#define REPEATS 3 +#define MIN 1 +#define MAX 6 + +/** +* Prints the TAP version. +*/ +static void print_version( void ) { + printf( "TAP version 13\n" ); +} + +/** +* Prints the TAP summary. +* +* @param total total number of tests +* @param passing total number of passing tests +*/ +static void print_summary( int total, int passing ) { + printf( "#\n" ); + printf( "1..%d\n", total ); // TAP plan + printf( "# total %d\n", total ); + printf( "# pass %d\n", passing ); + printf( "#\n" ); + printf( "# ok\n" ); +} + +/** +* Prints benchmarks results. +* +* @param iterations number of iterations +* @param elapsed elapsed time in seconds +*/ +static void print_results( int iterations, double elapsed ) { + double rate = (double)iterations / elapsed; + printf( " ---\n" ); + printf( " iterations: %d\n", iterations ); + printf( " elapsed: %0.9f\n", elapsed ); + printf( " rate: %0.9f\n", rate ); + printf( " ...\n" ); +} + +/** +* Returns a clock time. +* +* @return clock time +*/ +static double tic( void ) { + struct timeval now; + gettimeofday( &now, NULL ); + return (double)now.tv_sec + (double)now.tv_usec/1.0e6; +} + +/** +* Generates a random number on the interval [0,1). +* +* @return random number +*/ +static float rand_float( void ) { + int r = rand(); + return (float)r / ( (float)RAND_MAX + 1.0f ); +} + +/** +* Runs a benchmark. +* +* @param iterations number of iterations +* @param len array length +* @return elapsed time in seconds +*/ +static double benchmark1( int iterations, int len ) { + stdlib_complex64_t start; + stdlib_complex64_t stop; + double elapsed; + double t; + float *x; + int i; + + x = (float *) malloc( len * 2 * sizeof( float ) ); + for ( i = 0; i < len*2; i += 2 ) { + x[ i ] = ( rand_float()*2.0f ) - 1.0f; + x[ i+1 ] = ( rand_float()*2.0f ) - 1.0f; + } + start = stdlib_complex64( 0.0f, 0.0f ); + stop = stdlib_complex64( 9.0f, 0.0f ); + t = tic(); + for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar + stdlib_strided_clogspace( len, 10.0f, start, stop, true, (stdlib_complex64_t *)x, 1 ); + if ( x[ 0 ] != x[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( x[ 0 ] != x[ 0 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + return elapsed; +} + +/** +* Runs a benchmark. +* +* @param iterations number of iterations +* @param len array length +* @return elapsed time in seconds +*/ +static double benchmark2( int iterations, int len ) { + stdlib_complex64_t start; + stdlib_complex64_t stop; + double elapsed; + double t; + float *x; + int i; + + x = (float *) malloc( len * 2 * sizeof( float ) ); + for ( i = 0; i < len*2; i += 2 ) { + x[ i ] = ( rand_float()*2.0f ) - 1.0f; + x[ i+1 ] = ( rand_float()*2.0f ) - 1.0f; + } + start = stdlib_complex64( 0.0f, 0.0f ); + stop = stdlib_complex64( 9.0f, 0.0f ); + t = tic(); + for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar + stdlib_strided_clogspace_ndarray( len, 10.0f, start, stop, true, (stdlib_complex64_t *)x, 1, 0 ); + if ( x[ 0 ] != x[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( x[ 0 ] != x[ 0 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + return elapsed; +} + +/** +* Main execution sequence. +*/ +int main( void ) { + double elapsed; + int count; + int iter; + int len; + int i; + int j; + + // Use the current time to seed the random number generator: + srand( time( NULL ) ); + + print_version(); + count = 0; + for ( i = MIN; i <= MAX; i++ ) { + len = pow( 10, i ); + iter = ITERATIONS / pow( 10, i-1 ); + for ( j = 0; j < REPEATS; j++ ) { + count += 1; + printf( "# c::%s:len=%d\n", NAME, len ); + elapsed = benchmark1( iter, len ); + print_results( iter, elapsed ); + printf( "ok %d benchmark finished\n", count ); + } + for ( j = 0; j < REPEATS; j++ ) { + count += 1; + printf( "# c::%s:ndarray:len=%d\n", NAME, len ); + elapsed = benchmark2( iter, len ); + print_results( iter, elapsed ); + printf( "ok %d benchmark finished\n", count ); + } + } + print_summary( count, count ); +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/binding.gyp b/lib/node_modules/@stdlib/blas/ext/base/clogspace/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/binding.gyp @@ -0,0 +1,170 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +# A `.gyp` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt new file mode 100644 index 000000000000..809d5fc7dc58 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt @@ -0,0 +1,148 @@ + +{{alias}}( N, base, start, stop, endpoint, x, strideX ) + Fills a single-precision complex floating-point strided array with + logarithmically spaced values over a specified interval. + + If `N = 1` and `endpoint` is `true`, the set of values written to an input + array only includes `base^stop`, but not `base^start`; otherwise, when + `N = 1` and `endpoint` is `false`, the set of values written to an input + array only includes `base^start`, but not `base^stop`. + + Generated values follow a logarithmic spiral in the complex plane, where + imaginary components of `start` and `stop` control phase rotation. If the + real part of `start` is less than that of `stop`, the magnitudes of the + generated values will be in ascending order; otherwise, the magnitudes will + be in descending order. + + When `N >= 2` and `endpoint` is `true`, the set of values written to an + input array is guaranteed to include the `base^start` and `base^stop` + values. Beware, however, that values between `base^start` and `base^stop` + are subject to floating-point rounding errors. + + The `N` and stride parameters determine which elements in the strided array + are accessed at runtime. + + Indexing is relative to the first index. To introduce an offset, use typed + array views. + + If `N <= 0`, the function returns `x` unchanged. + + Parameters + ---------- + N: integer + Number of indexed elements. + + base: number + Base of the logarithmic scale. + + start: Complex64 + Exponent of the starting value. + + stop: Complex64 + Exponent of the final value. + + endpoint: boolean + Boolean indicating whether to include the `base^stop` value when + writing values to the input array. If `false`, the function generates + `N+1` logarithmically spaced values over the interval `[base^start, + base^stop]` and only writes `N` values to the input array, thus + excluding `base^stop` from the input array. Accordingly, for a fixed + `N`, the spacing between adjacent values written to an input array + changes depending on the value of `endpoint`. + + x: Complex64Array + Input array. + + strideX: integer + Stride length. + + Returns + ------- + x: Complex64Array + Input array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > var strt = new {{alias:@stdlib/complex/float32/ctor}}( 0.0, 0.0 ); + > var stp = new {{alias:@stdlib/complex/float32/ctor}}( 2.0, 0.0 ); + > {{alias}}( x.length, 10.0, strt, stp, true, x, 1 ) + [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + + // Using `N` and stride parameters: + > x = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > {{alias}}( 2, 10.0, strt, stp, true, x, 2 ) + [ 1.0, 0.0, 0.0, 0.0, 100.0, 0.0 ] + + // Using view offsets: + > var x0 = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > var x1 = new {{alias:@stdlib/array/complex64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); + > {{alias}}( 3, 10.0, strt, stp, true, x1, 1 ) + [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + > x0 + [ 0.0, 0.0, 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + + +{{alias}}.ndarray( N, base, start, stop, endpoint, x, strideX, offsetX ) + Fills a single-precision complex floating-point strided array with + logarithmically spaced values over a specified interval using alternative + indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the offset parameter supports indexing semantics based on a starting + index. + + Parameters + ---------- + N: integer + Number of indexed elements. + + base: number + Base of the logarithmic scale. + + start: Complex64 + Exponent of the starting value. + + stop: Complex64 + Exponent of the final value. + + endpoint: boolean + Boolean indicating whether to include the `base^stop` value when + writing values to the input array. If `false`, the function generates + `N+1` logarithmically spaced values over the interval `[base^start, + base^stop]` and only writes `N` values to the input array, thus + excluding `base^stop` from the input array. Accordingly, for a fixed + `N`, the spacing between adjacent values written to an input array + changes depending on the value of `endpoint`. + + x: Complex64Array + Input array. + + strideX: integer + Stride length. + + offsetX: integer + Starting index. + + Returns + ------- + x: Complex64Array + Input array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > var strt = new {{alias:@stdlib/complex/float32/ctor}}( 0.0, 0.0 ); + > var stp = new {{alias:@stdlib/complex/float32/ctor}}( 2.0, 0.0 ); + > {{alias}}.ndarray( x.length, 10.0, strt, stp, true, x, 1, 0 ) + [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + + // Using an index offset: + > x = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > {{alias}}.ndarray( 3, 10.0, strt, stp, true, x, 1, 1 ) + [ 0.0, 0.0, 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + + See Also + -------- diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/index.d.ts new file mode 100644 index 000000000000..ff6f7bc65ddb --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/index.d.ts @@ -0,0 +1,149 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +// TypeScript Version: 4.1 + +/// + +import { Complex64Array } from '@stdlib/types/array'; +import { Complex64 } from '@stdlib/types/complex'; + +/** +* Interface describing `clogspace`. +*/ +interface Routine { + /** + * Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. + * + * @param N - number of indexed elements + * @param base - base of the logarithmic scale + * @param start - exponent of the starting value + * @param stop - exponent of the final value + * @param endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array + * @param x - input array + * @param strideX - stride length + * @returns input array + * + * @example + * var Complex64Array = require( '@stdlib/array/complex64' ); + * var Complex64 = require( '@stdlib/complex/float32/ctor' ); + * + * var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + * + * var strt = new Complex64( 0.0, 0.0 ); + * var stp = new Complex64( 2.0, 0.0 ); + * + * clogspace( x.length, 10.0, strt, stp, true, x, 1 ); + * // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + */ + ( N: number, base: number, start: Complex64, stop: Complex64, endpoint: boolean, x: Complex64Array, strideX: number ): Complex64Array; + + /** + * Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. + * + * @param N - number of indexed elements + * @param base - base of the logarithmic scale + * @param start - exponent of the starting value + * @param stop - exponent of the final value + * @param endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array + * @param x - input array + * @param strideX - stride length + * @param offsetX - starting index + * @returns input array + * + * @example + * var Complex64Array = require( '@stdlib/array/complex64' ); + * var Complex64 = require( '@stdlib/complex/float32/ctor' ); + * + * var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + * + * var strt = new Complex64( 0.0, 0.0 ); + * var stp = new Complex64( 2.0, 0.0 ); + * + * clogspace.ndarray( x.length, 10.0, strt, stp, true, x, 1, 0 ); + * // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + */ + ndarray( N: number, base: number, start: Complex64, stop: Complex64, endpoint: boolean, x: Complex64Array, strideX: number, offsetX: number ): Complex64Array; +} + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +* +* @param N - number of indexed elements +* @param base - base of the logarithmic scale +* @param start - exponent of the starting value +* @param stop - exponent of the final value +* @param endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param x - input array +* @param strideX - stride length +* @returns input array +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, false, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace.ndarray( x.length, 10.0, strt, stp, true, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace.ndarray( x.length, 10.0, strt, stp, false, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ +declare var clogspace: Routine; + + +// EXPORTS // + +export = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/test.ts new file mode 100644 index 000000000000..c27b27280d5d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/types/test.ts @@ -0,0 +1,317 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable space-in-parens */ + +import Complex64Array = require( '@stdlib/array/complex64' ); +import Complex64 = require( '@stdlib/complex/float32/ctor' ); +import clogspace = require( './index' ); + + +// TESTS // + +// The function returns a Complex64Array... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, 10.0, start, stop, true, x, 1 ); // $ExpectType Complex64Array +} + +// The compiler throws an error if the function is provided a first argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( '10', 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( true, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( false, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( null, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( undefined, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( [], 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( {}, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace( ( x: number ): number => x, 10.0, start, stop, true, x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, '10', start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, true, start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, false, start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, null, start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, undefined, start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, [], start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, {}, start, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, ( x: number ): number => x, start, stop, true, x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not complex-like... +{ + const x = new Complex64Array( 10 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, 10.0, 10, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, '10', stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, true, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, false, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, null, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, undefined, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, [ '1' ], stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, {}, stop, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, ( x: number ): number => x, stop, true, x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not complex-like... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + + clogspace( x.length, 10.0, start, 10, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, '10', true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, true, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, false, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, null, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, undefined, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, [ '1' ], true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, {}, true, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, ( x: number ): number => x, true, x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fifth argument which is not a boolean... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, 10.0, start, stop, '10', x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, 99.99, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, null, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, undefined, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, [], x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, {}, x, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, ( x: number ): number => x, x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a sixth argument which is not a Complex64Array... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, 10.0, start, stop, true, 10, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, '10', 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, true, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, false, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, null, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, undefined, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, [], 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, {}, 1 ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, ( x: number ): number => x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a seventh argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace( x.length, 10.0, start, stop, true, x, '10' ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, true ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, false ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, null ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, undefined ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, [] ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, {} ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace(); // $ExpectError + clogspace( x.length ); // $ExpectError + clogspace( x.length, 10.0 ); // $ExpectError + clogspace( x.length, 10.0, start ); // $ExpectError + clogspace( x.length, 10.0, start, stop ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x ); // $ExpectError + clogspace( x.length, 10.0, start, stop, true, x, 1, {} ); // $ExpectError +} + +// Attached to main export is an `ndarray` method which returns a Complex64Array... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, 0 ); // $ExpectType Complex64Array +} + +// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( '10', 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( true, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( false, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( null, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( undefined, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( [], 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( {}, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( ( x: number ): number => x, 10.0, start, stop, true, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a second argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, '10', start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, true, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, false, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, null, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, undefined, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, [], start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, {}, start, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, ( x: number ): number => x, start, stop, true, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a third argument which is not complex-like... +{ + const x = new Complex64Array( 10 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, 10, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, '10', stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, true, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, false, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, null, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, undefined, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, [ '1' ], stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, {}, stop, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, ( x: number ): number => x, stop, true, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fourth argument which is not complex-like... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, 10, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, '10', true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, true, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, false, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, null, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, undefined, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, [ '1' ], true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, {}, true, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, ( x: number ): number => x, true, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a boolean... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, stop, '10', x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, 99.99, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, null, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, undefined, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, [], x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, {}, x, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, ( x: number ): number => x, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a Complex64Array... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, stop, true, 10, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, '10', 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, true, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, false, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, null, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, undefined, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, [], 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, {}, 1, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, stop, true, x, '10', 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, true, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, false, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, null, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, undefined, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, [], 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, {}, 0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an eighth argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, '10' ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, true ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, false ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, null ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, undefined ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, [] ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, {} ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments... +{ + const x = new Complex64Array( 10 ); + const start = new Complex64( 0.0, 0.0 ); + const stop = new Complex64( 10.0, 0.0 ); + + clogspace.ndarray(); // $ExpectError + clogspace.ndarray( x.length ); // $ExpectError + clogspace.ndarray( x.length, 10.0 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1 ); // $ExpectError + clogspace.ndarray( x.length, 10.0, start, stop, true, x, 1, 0, {} ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := example.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled examples. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c new file mode 100644 index 000000000000..29d07bfa80ac --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c @@ -0,0 +1,41 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/clogspace.h" +#include "stdlib/complex/float32/ctor.h" +#include +#include + +int main( void ) { + // Create a strided array of interleaved real and imaginary components: + float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + + // Specify the number of indexed elements: + const int N = 6; + + // Specify a stride: + const int strideX = 1; + + // Fill the array: + stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 5.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ i*2+1 ] ); + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/index.js new file mode 100644 index 000000000000..9d79f15e72e0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/index.js @@ -0,0 +1,32 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var logEach = require( '@stdlib/console/log-each' ); +var clogspace = require( './../lib' ); + +var x = new Complex64Array( 6 ); + +var strt = new Complex64( 0.0, 0.0 ); +var stp = new Complex64( 5.0, 0.0 ); + +clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +logEach( '%s', x ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/include.gypi b/lib/node_modules/@stdlib/blas/ext/base/clogspace/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/include.gypi @@ -0,0 +1,53 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +# A GYP include file for building a Node.js native add-on. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + ' + +/* +* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler. +*/ +#ifdef __cplusplus +extern "C" { +#endif + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +*/ +void API_SUFFIX(stdlib_strided_clogspace)( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX ); + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. +*/ +void API_SUFFIX(stdlib_strided_clogspace_ndarray)( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); + +#ifdef __cplusplus +} +#endif + +#endif // !STDLIB_BLAS_EXT_BASE_CLOGSPACE_H diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.js new file mode 100644 index 000000000000..684b420cda6b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.js @@ -0,0 +1,72 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var stride2offset = require( '@stdlib/strided/base/stride2offset' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} base - base of the logarithmic scale +* @param {Complex64} start - exponent of the starting value +* @param {Complex64} stop - exponent of the final value +* @param {boolean} endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param {Complex64Array} x - input array +* @param {integer} strideX - stride length +* @returns {Complex64Array} input array +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, false, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ +function clogspace( N, base, start, stop, endpoint, x, strideX ) { + return ndarray( N, base, start, stop, endpoint, x, strideX, stride2offset( N, strideX ) ); // eslint-disable-line max-len +} + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.native.js new file mode 100644 index 000000000000..1f1c4f8e4b55 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/clogspace.native.js @@ -0,0 +1,74 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} base - base of the logarithmic scale +* @param {Complex64} start - exponent of the starting value +* @param {Complex64} stop - exponent of the final value +* @param {boolean} endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param {Complex64Array} x - input array +* @param {integer} strideX - stride length +* @returns {Complex64Array} input array +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, false, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ +function clogspace( N, base, start, stop, endpoint, x, strideX ) { + var view = reinterpret( x, 0 ); + addon( N, base, start, stop, endpoint, view, strideX ); + return x; +} + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/index.js new file mode 100644 index 000000000000..80823613c5f4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/index.js @@ -0,0 +1,76 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +/** +* Fill a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +* +* @module @stdlib/blas/ext/base/clogspace +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var clogspace = require( '@stdlib/blas/ext/base/clogspace' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var clogspace = require( '@stdlib/blas/ext/base/clogspace' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace.ndarray( x.length, 10.0, strt, stp, false, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ + +// MODULES // + +var join = require( 'path' ).join; +var tryRequire = require( '@stdlib/utils/try-require' ); +var isError = require( '@stdlib/assert/is-error' ); +var main = require( './main.js' ); + + +// MAIN // + +var clogspace; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + clogspace = main; +} else { + clogspace = tmp; +} + + +// EXPORTS // + +module.exports = clogspace; + +// exports: { "ndarray": "clogspace.ndarray" } diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/main.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/main.js new file mode 100644 index 000000000000..271582398229 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/main.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var clogspace = require( './clogspace.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( clogspace, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/native.js new file mode 100644 index 000000000000..309f2d985920 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/native.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var clogspace = require( './clogspace.native.js' ); +var ndarray = require( './ndarray.native.js' ); + + +// MAIN // + +setReadOnly( clogspace, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js new file mode 100644 index 000000000000..5e1f2692753d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js @@ -0,0 +1,151 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var realf = require( '@stdlib/complex/float32/real' ); +var imagf = require( '@stdlib/complex/float32/imag' ); +var powf = require( '@stdlib/math/base/special/powf' ); +var lnf = require( '@stdlib/math/base/special/lnf' ); +var cosf = require( '@stdlib/math/base/special/cosf' ); +var sinf = require( '@stdlib/math/base/special/sinf' ); + + +// MAIN // + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} base - base of the logarithmic scale +* @param {Complex64} start - exponent of the starting value +* @param {Complex64} stop - exponent of the final value +* @param {boolean} endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param {Complex64Array} x - input array +* @param {integer} strideX - stride length +* @param {NonNegativeInteger} offsetX - starting index +* @returns {Complex64Array} input array +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, false, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ +function clogspace( N, base, start, stop, endpoint, x, strideX, offsetX ) { + var startRe; + var startIm; + var stopRe; + var stopIm; + var expRe; + var expIm; + var scale; + var view; + var lnb; + var dre; + var dim; + var ix; + var i; + + if ( N <= 0 ) { + return x; + } + + // Decompose the exponent bounds into their real and imaginary components: + startRe = realf( start ); + startIm = imagf( start ); + stopRe = realf( stop ); + stopIm = imagf( stop ); + + // Compute ln(base) for use in base^z = exp(z * ln(base)): + lnb = lnf( base ); + + // Reinterpret the complex input array as a real-valued array: + view = reinterpret( x, 0 ); + + // Adjust the offset according to the real-valued array (2 floats per element): + ix = offsetX * 2; + + // Set the first value: + if ( N === 1 ) { + if ( endpoint ) { + scale = powf( base, stopRe ); + view[ ix ] = scale * cosf( stopIm * lnb ); + view[ ix+1 ] = scale * sinf( stopIm * lnb ); + } else { + scale = powf( base, startRe ); + view[ ix ] = scale * cosf( startIm * lnb ); + view[ ix+1 ] = scale * sinf( startIm * lnb ); + } + return x; + } + scale = powf( base, startRe ); + view[ ix ] = scale * cosf( startIm * lnb ); + view[ ix+1 ] = scale * sinf( startIm * lnb ); + ix += strideX * 2; + + // Calculate the complex increment: + if ( endpoint ) { + N -= 1; + } + dre = ( stopRe - startRe ) / N; + dim = ( stopIm - startIm ) / N; + + // Generate logarithmically spaced values: + for ( i = 1; i < N; i++ ) { + expRe = startRe + ( dre * i ); + expIm = startIm + ( dim * i ); + scale = powf( base, expRe ); + view[ ix ] = scale * cosf( expIm * lnb ); + view[ ix+1 ] = scale * sinf( expIm * lnb ); + ix += strideX * 2; + } + // Check whether to include the `base^stop` value: + if ( endpoint ) { + scale = powf( base, stopRe ); + view[ ix ] = scale * cosf( stopIm * lnb ); + view[ ix+1 ] = scale * sinf( stopIm * lnb ); + } + return x; +} + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.native.js new file mode 100644 index 000000000000..9f1bbf34f11d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.native.js @@ -0,0 +1,75 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} base - base of the logarithmic scale +* @param {Complex64} start - exponent of the starting value +* @param {Complex64} stop - exponent of the final value +* @param {boolean} endpoint - boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param {Complex64Array} x - input array +* @param {integer} strideX - stride length +* @param {NonNegativeInteger} offsetX - starting index +* @returns {Complex64Array} input array +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, true, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); +* +* var strt = new Complex64( 0.0, 0.0 ); +* var stp = new Complex64( 2.0, 0.0 ); +* +* clogspace( x.length, 10.0, strt, stp, false, x, 1, 0 ); +* // x => [ 1.0, 0.0, 10.0, 0.0 ] +*/ +function clogspace( N, base, start, stop, endpoint, x, strideX, offsetX ) { + var view = reinterpret( x, 0 ); + addon.ndarray( N, base, start, stop, endpoint, view, strideX, offsetX ); + return x; +} + + +// EXPORTS // + +module.exports = clogspace; diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json new file mode 100644 index 000000000000..ce26c7c8960b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json @@ -0,0 +1,102 @@ +{ + "options": { + "task": "build" + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-float", + "@stdlib/napi/argv-bool", + "@stdlib/napi/argv-complex64", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/special/powf", + "@stdlib/math/base/special/lnf", + "@stdlib/math/base/special/cosf", + "@stdlib/math/base/special/sinf" + ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/special/powf", + "@stdlib/math/base/special/lnf", + "@stdlib/math/base/special/cosf", + "@stdlib/math/base/special/sinf" + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/special/powf", + "@stdlib/math/base/special/lnf", + "@stdlib/math/base/special/cosf", + "@stdlib/math/base/special/sinf" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json b/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json new file mode 100644 index 000000000000..774b0177edf0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json @@ -0,0 +1,82 @@ +{ + "name": "@stdlib/blas/ext/base/clogspace", + "version": "0.0.0", + "description": "Fill a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "browser": "./lib/main.js", + "gypfile": true, + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "include": "./include", + "lib": "./lib", + "src": "./src", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "mathematics", + "math", + "blas", + "extended", + "fill", + "assign", + "set", + "logspace", + "logarithmic", + "log", + "sequence", + "seq", + "strided", + "array", + "ndarray", + "complex", + "complex64", + "single-precision", + "float32", + "float32array" + ], + "__stdlib__": { + "wasm": false + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/Makefile b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/addon.c new file mode 100644 index 000000000000..0320a059fb95 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/addon.c @@ -0,0 +1,72 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/clogspace.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_float.h" +#include "stdlib/napi/argv_bool.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_complex64.h" +#include "stdlib/napi/argv_strided_complex64array.h" +#include + +/** +* Receives JavaScript callback invocation data. +* +* @param env environment under which the function is invoked +* @param info callback data +* @return Node-API value +*/ +static napi_value addon( napi_env env, napi_callback_info info ) { + STDLIB_NAPI_ARGV( env, info, argv, argc, 7 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_FLOAT( env, base, argv, 1 ); + STDLIB_NAPI_ARGV_COMPLEX64( env, start, argv, 2 ); + STDLIB_NAPI_ARGV_COMPLEX64( env, stop, argv, 3 ); + STDLIB_NAPI_ARGV_BOOL( env, endpoint, argv, 4 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 6 ); + STDLIB_NAPI_ARGV_STRIDED_COMPLEX64ARRAY( env, X, N, strideX, argv, 5 ); + API_SUFFIX(stdlib_strided_clogspace)( N, base, start, stop, endpoint, (stdlib_complex64_t *)X, strideX ); + return NULL; +} + +/** +* Receives JavaScript callback invocation data. +* +* @param env environment under which the function is invoked +* @param info callback data +* @return Node-API value +*/ +static napi_value addon_method( napi_env env, napi_callback_info info ) { + STDLIB_NAPI_ARGV( env, info, argv, argc, 8 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_FLOAT( env, base, argv, 1 ); + STDLIB_NAPI_ARGV_COMPLEX64( env, start, argv, 2 ); + STDLIB_NAPI_ARGV_COMPLEX64( env, stop, argv, 3 ); + STDLIB_NAPI_ARGV_BOOL( env, endpoint, argv, 4 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 6 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 7 ); + STDLIB_NAPI_ARGV_STRIDED_COMPLEX64ARRAY( env, X, N, strideX, argv, 5 ); + API_SUFFIX(stdlib_strided_clogspace_ndarray)( N, base, start, stop, endpoint, (stdlib_complex64_t *)X, strideX, offsetX ); + return NULL; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c new file mode 100644 index 000000000000..3ea3956fcfe6 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c @@ -0,0 +1,129 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/clogspace.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/real.h" +#include "stdlib/complex/float32/imag.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" +#include "stdlib/math/base/special/powf.h" +#include "stdlib/math/base/special/lnf.h" +#include "stdlib/math/base/special/cosf.h" +#include "stdlib/math/base/special/sinf.h" +#include + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval. +* +* @param N number of indexed elements +* @param base base of the logarithmic scale +* @param start exponent of the starting value +* @param stop exponent of the final value +* @param endpoint boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param X input array +* @param strideX stride length +*/ +void API_SUFFIX(stdlib_strided_clogspace)( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX ) { + CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + API_SUFFIX(stdlib_strided_clogspace_ndarray)( N, base, start, stop, endpoint, X, strideX, ox ); +} + +/** +* Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. +* +* @param N number of indexed elements +* @param base base of the logarithmic scale +* @param start exponent of the starting value +* @param stop exponent of the final value +* @param endpoint boolean indicating whether to include the `base^stop` value when writing values to the input array +* @param X input array +* @param strideX stride length +* @param offsetX starting index +*/ +void API_SUFFIX(stdlib_strided_clogspace_ndarray)( const CBLAS_INT N, const float base, const stdlib_complex64_t start, const stdlib_complex64_t stop, const bool endpoint, stdlib_complex64_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + float start_re; + float start_im; + float stop_re; + float stop_im; + float exp_re; + float exp_im; + float scale; + float lnb; + float dre; + float dim; + CBLAS_INT ix; + CBLAS_INT M; + CBLAS_INT i; + + if ( N <= 0 ) { + return; + } + + // Decompose the exponent bounds into their real and imaginary components: + start_re = stdlib_complex64_real( start ); + start_im = stdlib_complex64_imag( start ); + stop_re = stdlib_complex64_real( stop ); + stop_im = stdlib_complex64_imag( stop ); + + // Compute ln(base) for use in base^z = exp(z * ln(base)): + lnb = stdlib_base_lnf( base ); + + ix = offsetX; + + // Set the first value: + if ( N == 1 ) { + if ( endpoint ) { + scale = stdlib_base_powf( base, stop_re ); + X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( stop_im * lnb ), scale * stdlib_base_sinf( stop_im * lnb ) ); + } else { + scale = stdlib_base_powf( base, start_re ); + X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( start_im * lnb ), scale * stdlib_base_sinf( start_im * lnb ) ); + } + return; + } + + // Write the first value: + scale = stdlib_base_powf( base, start_re ); + X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( start_im * lnb ), scale * stdlib_base_sinf( start_im * lnb ) ); + ix += strideX; + + // Calculate the complex increment: + if ( endpoint ) { + M = N - 1; + } else { + M = N; + } + dre = ( stop_re - start_re ) / (float)M; + dim = ( stop_im - start_im ) / (float)M; + + // Generate logarithmically spaced values: + for ( i = 1; i < M; i++ ) { + exp_re = start_re + ( dre * (float)i ); + exp_im = start_im + ( dim * (float)i ); + scale = stdlib_base_powf( base, exp_re ); + X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( exp_im * lnb ), scale * stdlib_base_sinf( exp_im * lnb ) ); + ix += strideX; + } + // Check whether to include the `base^stop` value: + if ( endpoint ) { + scale = stdlib_base_powf( base, stop_re ); + X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( stop_im * lnb ), scale * stdlib_base_sinf( stop_im * lnb ) ); + } + return; +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.js new file mode 100644 index 000000000000..d78dc5e61f3b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.js @@ -0,0 +1,266 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var isSameComplex64Array = require( '@stdlib/assert/is-same-complex64array' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var lnf = require( '@stdlib/math/base/special/lnf' ); +var cosf = require( '@stdlib/math/base/special/cosf' ); +var sinf = require( '@stdlib/math/base/special/sinf' ); +var E = require( '@stdlib/constants/float32/e' ); +var clogspace = require( './../lib/clogspace.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof clogspace, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', function test( t ) { + t.strictEqual( clogspace.length, 7, 'has expected arity' ); + t.end(); +}); + +tape( 'the function fills a strided array', function test( t ) { + var expected; + var lnb; + var x; + + x = new Complex64Array([ + 4.0, + 2.0, + -3.0, + 5.0, + -1.0, + 2.0, + -5.0, + 6.0, + 3.0, + -7.0, + 0.0, + 1.0, + 8.0, + -2.0, + 4.0, + 3.0 + ]); + expected = new Complex64Array([ + 1.0, + 0.0, + 2.0, + 0.0, + 4.0, + 0.0, + 8.0, + 0.0, + 16.0, + 0.0, + 32.0, + 0.0, + 64.0, + 0.0, + 128.0, + 0.0 + ]); + + clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 7.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 2.0, 0.0, 4.0, 0.0 ] ); + + clogspace( x.length, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // Non-zero imaginary exponents: + lnb = lnf( E ); + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0, cosf( lnb ), sinf( lnb ) ] ); + clogspace( 2, E, new Complex64( 0.0, 0.0 ), new Complex64( 0.0, 1.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter equal to `1`, the function writes a single value to the input array', function test( t ) { + var expected; + var x; + + // When `endpoint` is `true`, write `base^stop`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 8.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // When `endpoint` is `false`, write `base^start`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the input array', function test( t ) { + var out; + var x; + + x = new Complex64Array([ + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0 + ]); + out = clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 5.0, 0.0 ), false, x, 1 ); + + t.strictEqual( out, x, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns the input array unchanged', function test( t ) { + var expected; + var x; + + x = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + expected = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + + clogspace( 0, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + clogspace( -4, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports specifying a stride', function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 0 + 1.0, // 0 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 2 + 3.0 // 2 + ]); + expected = new Complex64Array([ + 2.0, // 0 + 0.0, // 0 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 8.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 4.0, 0.0 ), false, x, 2 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports specifying a negative stride', function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 2 + 1.0, // 2 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 0 + 3.0 // 0 + ]); + expected = new Complex64Array([ + 8.0, // 2 + 0.0, // 2 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 2.0, // 0 + 0.0 // 0 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, -2 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', function test( t ) { + var expected; + var x0; + var x1; + + x0 = new Complex64Array([ + 1.0, + -2.0, + 3.0, // 0 + -4.0, // 0 + 6.0, + -8.0, + 10.0, // 1 + -12.0 // 1 + ]); + expected = new Complex64Array([ + 1.0, + -2.0, + 1.0, // 0 + 0.0, // 0 + 6.0, + -8.0, + 10.0, // 1 + 0.0 // 1 + ]); + + x1 = new Complex64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); + + clogspace( 2, 10.0, new Complex64( 0.0, 0.0 ), new Complex64( 1.0, 0.0 ), true, x1, 2 ); + t.strictEqual( isSameComplex64Array( x0, expected ), true, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.native.js new file mode 100644 index 000000000000..25d5c32baff4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.clogspace.native.js @@ -0,0 +1,275 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var isSameComplex64Array = require( '@stdlib/assert/is-same-complex64array' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var lnf = require( '@stdlib/math/base/special/lnf' ); +var cosf = require( '@stdlib/math/base/special/cosf' ); +var sinf = require( '@stdlib/math/base/special/sinf' ); +var E = require( '@stdlib/constants/float32/e' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var clogspace = tryRequire( resolve( __dirname, './../lib/clogspace.native.js' ) ); +var opts = { + 'skip': ( clogspace instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof clogspace, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', opts, function test( t ) { + t.strictEqual( clogspace.length, 7, 'has expected arity' ); + t.end(); +}); + +tape( 'the function fills a strided array', opts, function test( t ) { + var expected; + var lnb; + var x; + + x = new Complex64Array([ + 4.0, + 2.0, + -3.0, + 5.0, + -1.0, + 2.0, + -5.0, + 6.0, + 3.0, + -7.0, + 0.0, + 1.0, + 8.0, + -2.0, + 4.0, + 3.0 + ]); + expected = new Complex64Array([ + 1.0, + 0.0, + 2.0, + 0.0, + 4.0, + 0.0, + 8.0, + 0.0, + 16.0, + 0.0, + 32.0, + 0.0, + 64.0, + 0.0, + 128.0, + 0.0 + ]); + + clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 7.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 2.0, 0.0, 4.0, 0.0 ] ); + + clogspace( x.length, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // Non-zero imaginary exponents: + lnb = lnf( E ); + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0, cosf( lnb ), sinf( lnb ) ] ); + clogspace( 2, E, new Complex64( 0.0, 0.0 ), new Complex64( 0.0, 1.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter equal to `1`, the function writes a single value to the input array', opts, function test( t ) { + var expected; + var x; + + // When `endpoint` is `true`, write `base^stop`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 8.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // When `endpoint` is `false`, write `base^start`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the input array', opts, function test( t ) { + var out; + var x; + + x = new Complex64Array([ + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0 + ]); + out = clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 5.0, 0.0 ), false, x, 1 ); + + t.strictEqual( out, x, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns the input array unchanged', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + expected = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + + clogspace( 0, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + clogspace( -4, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports specifying a stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 0 + 1.0, // 0 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 2 + 3.0 // 2 + ]); + expected = new Complex64Array([ + 2.0, // 0 + 0.0, // 0 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 8.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 4.0, 0.0 ), false, x, 2 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports specifying a negative stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 2 + 1.0, // 2 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 0 + 3.0 // 0 + ]); + expected = new Complex64Array([ + 8.0, // 2 + 0.0, // 2 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 2.0, // 0 + 0.0 // 0 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, -2 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', opts, function test( t ) { + var expected; + var x0; + var x1; + + x0 = new Complex64Array([ + 1.0, + -2.0, + 3.0, // 0 + -4.0, // 0 + 6.0, + -8.0, + 10.0, // 1 + -12.0 // 1 + ]); + expected = new Complex64Array([ + 1.0, + -2.0, + 1.0, // 0 + 0.0, // 0 + 6.0, + -8.0, + 10.0, // 1 + 0.0 // 1 + ]); + + x1 = new Complex64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); + + clogspace( 2, 10.0, new Complex64( 0.0, 0.0 ), new Complex64( 1.0, 0.0 ), true, x1, 2 ); + t.strictEqual( isSameComplex64Array( x0, expected ), true, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.js new file mode 100644 index 000000000000..ace252227cdc --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.js @@ -0,0 +1,82 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var proxyquire = require( 'proxyquire' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var clogspace = require( './../lib' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof clogspace, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'attached to the main export is a method providing an ndarray interface', function test( t ) { + t.strictEqual( typeof clogspace.ndarray, 'function', 'method is a function' ); + t.end(); +}); + +tape( 'if a native implementation is available, the main export is the native implementation', opts, function test( t ) { + var clogspace = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( clogspace, mock, 'returns native implementation' ); + t.end(); + + function tryRequire() { + return mock; + } + + function mock() { + // Mock... + } +}); + +tape( 'if a native implementation is not available, the main export is a JavaScript implementation', opts, function test( t ) { + var clogspace; + var main; + + main = require( './../lib/clogspace.js' ); + + clogspace = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( clogspace, main, 'returns JavaScript implementation' ); + t.end(); + + function tryRequire() { + return new Error( 'Cannot find module' ); + } +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.js new file mode 100644 index 000000000000..a5201a94ca51 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.js @@ -0,0 +1,263 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var isSameComplex64Array = require( '@stdlib/assert/is-same-complex64array' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var lnf = require( '@stdlib/math/base/special/lnf' ); +var cosf = require( '@stdlib/math/base/special/cosf' ); +var sinf = require( '@stdlib/math/base/special/sinf' ); +var E = require( '@stdlib/constants/float32/e' ); +var clogspace = require( './../lib/ndarray.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof clogspace, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 8', function test( t ) { + t.strictEqual( clogspace.length, 8, 'has expected arity' ); + t.end(); +}); + +tape( 'the function fills a strided array', function test( t ) { + var expected; + var lnb; + var x; + + x = new Complex64Array([ + 4.0, + 2.0, + -3.0, + 5.0, + -1.0, + 2.0, + -5.0, + 6.0, + 3.0, + -7.0, + 0.0, + 1.0, + 8.0, + -2.0, + 4.0, + 3.0 + ]); + expected = new Complex64Array([ + 1.0, + 0.0, + 2.0, + 0.0, + 4.0, + 0.0, + 8.0, + 0.0, + 16.0, + 0.0, + 32.0, + 0.0, + 64.0, + 0.0, + 128.0, + 0.0 + ]); + + clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 7.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 2.0, 0.0, 4.0, 0.0 ] ); + + clogspace( x.length, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // Non-zero imaginary exponents: + lnb = lnf( E ); + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0, cosf( lnb ), sinf( lnb ) ] ); + clogspace( 2, E, new Complex64( 0.0, 0.0 ), new Complex64( 0.0, 1.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter equal to `1`, the function writes a single value to the input array', function test( t ) { + var expected; + var x; + + // When `endpoint` is `true`, write `base^stop`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 8.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // When `endpoint` is `false`, write `base^start`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the input array', function test( t ) { + var out; + var x; + + x = new Complex64Array([ + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0 + ]); + out = clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 5.0, 0.0 ), false, x, 1, 0 ); + + t.strictEqual( out, x, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns the input array unchanged', function test( t ) { + var expected; + var x; + + x = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + expected = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + + clogspace( 0, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + clogspace( -4, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports specifying a stride', function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 0 + 1.0, // 0 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 2 + 3.0 // 2 + ]); + expected = new Complex64Array([ + 2.0, // 0 + 0.0, // 0 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 8.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 4.0, 0.0 ), false, x, 2, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports specifying a negative stride', function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 2 + 1.0, // 2 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 0 + 3.0 // 0 + ]); + expected = new Complex64Array([ + 8.0, // 2 + 0.0, // 2 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 2.0, // 0 + 0.0 // 0 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, -2, 4 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an index offset', function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 1.0, + -2.0, + 0.0, // 0 + 0.0, // 0 + 0.0, // 1 + 0.0, // 1 + 0.0, // 2 + 0.0 // 2 + ]); + expected = new Complex64Array([ + 1.0, + -2.0, + 1.0, // 0 + 0.0, // 0 + 10.0, // 1 + 0.0, // 1 + 100.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 10.0, new Complex64( 0.0, 0.0 ), new Complex64( 2.0, 0.0 ), true, x, 1, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.native.js new file mode 100644 index 000000000000..8f20b13c88df --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/test/test.ndarray.native.js @@ -0,0 +1,272 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var isSameComplex64Array = require( '@stdlib/assert/is-same-complex64array' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var lnf = require( '@stdlib/math/base/special/lnf' ); +var cosf = require( '@stdlib/math/base/special/cosf' ); +var sinf = require( '@stdlib/math/base/special/sinf' ); +var E = require( '@stdlib/constants/float32/e' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var clogspace = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( clogspace instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof clogspace, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 8', opts, function test( t ) { + t.strictEqual( clogspace.length, 8, 'has expected arity' ); + t.end(); +}); + +tape( 'the function fills a strided array', opts, function test( t ) { + var expected; + var lnb; + var x; + + x = new Complex64Array([ + 4.0, + 2.0, + -3.0, + 5.0, + -1.0, + 2.0, + -5.0, + 6.0, + 3.0, + -7.0, + 0.0, + 1.0, + 8.0, + -2.0, + 4.0, + 3.0 + ]); + expected = new Complex64Array([ + 1.0, + 0.0, + 2.0, + 0.0, + 4.0, + 0.0, + 8.0, + 0.0, + 16.0, + 0.0, + 32.0, + 0.0, + 64.0, + 0.0, + 128.0, + 0.0 + ]); + + clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 7.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 2.0, 0.0, 4.0, 0.0 ] ); + + clogspace( x.length, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // Non-zero imaginary exponents: + lnb = lnf( E ); + x = new Complex64Array( [ 0.0, 0.0, 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0, cosf( lnb ), sinf( lnb ) ] ); + clogspace( 2, E, new Complex64( 0.0, 0.0 ), new Complex64( 0.0, 1.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter equal to `1`, the function writes a single value to the input array', opts, function test( t ) { + var expected; + var x; + + // When `endpoint` is `true`, write `base^stop`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 8.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + // When `endpoint` is `false`, write `base^start`: + x = new Complex64Array( [ 0.0, 0.0 ] ); + expected = new Complex64Array( [ 1.0, 0.0 ] ); + + clogspace( 1, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 3.0, 0.0 ), false, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the input array', opts, function test( t ) { + var out; + var x; + + x = new Complex64Array([ + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0, + 0.0 + ]); + out = clogspace( x.length, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 5.0, 0.0 ), false, x, 1, 0 ); + + t.strictEqual( out, x, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns the input array unchanged', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + expected = new Complex64Array( [ 3.0, -4.0, 1.0, -2.0, 5.0, 1.0 ] ); + + clogspace( 0, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + clogspace( -4, 2.0, new Complex64( 0.0, 0.0 ), new Complex64( 10.0, 0.0 ), true, x, 1, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports specifying a stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 0 + 1.0, // 0 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 2 + 3.0 // 2 + ]); + expected = new Complex64Array([ + 2.0, // 0 + 0.0, // 0 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 8.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 4.0, 0.0 ), false, x, 2, 0 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports specifying a negative stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 2.0, // 2 + 1.0, // 2 + -3.0, + -2.0, + -5.0, // 1 + 4.0, // 1 + 7.0, + -6.0, + 6.0, // 0 + 3.0 // 0 + ]); + expected = new Complex64Array([ + 8.0, // 2 + 0.0, // 2 + -3.0, + -2.0, + 4.0, // 1 + 0.0, // 1 + 7.0, + -6.0, + 2.0, // 0 + 0.0 // 0 + ]); + + clogspace( 3, 2.0, new Complex64( 1.0, 0.0 ), new Complex64( 3.0, 0.0 ), true, x, -2, 4 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an index offset', opts, function test( t ) { + var expected; + var x; + + x = new Complex64Array([ + 1.0, + -2.0, + 0.0, // 0 + 0.0, // 0 + 0.0, // 1 + 0.0, // 1 + 0.0, // 2 + 0.0 // 2 + ]); + expected = new Complex64Array([ + 1.0, + -2.0, + 1.0, // 0 + 0.0, // 0 + 10.0, // 1 + 0.0, // 1 + 100.0, // 2 + 0.0 // 2 + ]); + + clogspace( 3, 10.0, new Complex64( 0.0, 0.0 ), new Complex64( 2.0, 0.0 ), true, x, 1, 1 ); + t.strictEqual( isSameComplex64Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); From 2eb9a8a626dcad70327eea05296f07ae19b5eba3 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Tue, 28 Jul 2026 17:20:51 +0500 Subject: [PATCH 2/2] fix: apply suggestions from code review --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown_pkg_readmes status: passed - task: lint_markdown_docs status: na - task: lint_markdown status: na - task: lint_package_json status: passed - task: lint_repl_help status: passed - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: missing_dependencies - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/ext/base/clogspace/README.md | 8 +- .../blas/ext/base/clogspace/docs/repl.txt | 8 +- .../ext/base/clogspace/examples/c/example.c | 8 +- .../blas/ext/base/clogspace/lib/ndarray.js | 76 ++++++++++++++---- .../blas/ext/base/clogspace/manifest.json | 9 +-- .../blas/ext/base/clogspace/package.json | 3 +- .../blas/ext/base/clogspace/src/main.c | 77 ++++++++++++++----- 7 files changed, 136 insertions(+), 53 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md b/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md index db775bf7e9d7..f26f76a407ff 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/README.md @@ -311,20 +311,20 @@ void stdlib_strided_clogspace_ndarray( const CBLAS_INT N, const float base, cons int main( void ) { // Create a strided array of interleaved real and imaginary components: - float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; // Specify the number of indexed elements: - const int N = 6; + const int N = 8; // Specify a stride: const int strideX = 1; // Fill the array: - stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 5.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); + stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 2.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); // Print the result: for ( int i = 0; i < N; i++ ) { - printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ i*2+1 ] ); + printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ (i*2)+1 ] ); } } ``` diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt index 809d5fc7dc58..f3e8642f3b97 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/docs/repl.txt @@ -76,12 +76,12 @@ [ 1.0, 0.0, 0.0, 0.0, 100.0, 0.0 ] // Using view offsets: - > var x0 = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + > var x0 = new {{alias:@stdlib/array/complex64}}( [ 0.0, 0.0, 0.0, 0.0 ] ); > var x1 = new {{alias:@stdlib/array/complex64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); - > {{alias}}( 3, 10.0, strt, stp, true, x1, 1 ) - [ 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + > {{alias}}( 1, 10.0, strt, stp, true, x1, 1 ) + [ 100.0, 0.0 ] > x0 - [ 0.0, 0.0, 1.0, 0.0, 10.0, 0.0, 100.0, 0.0 ] + [ 0.0, 0.0, 100.0, 0.0 ] {{alias}}.ndarray( N, base, start, stop, endpoint, x, strideX, offsetX ) diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c index 29d07bfa80ac..dfe8221a3ac0 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/examples/c/example.c @@ -23,19 +23,19 @@ int main( void ) { // Create a strided array of interleaved real and imaginary components: - float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; + float x[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f }; // Specify the number of indexed elements: - const int N = 6; + const int N = 8; // Specify a stride: const int strideX = 1; // Fill the array: - stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 5.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); + stdlib_strided_clogspace( N, 10.0f, stdlib_complex64( 0.0f, 0.0f ), stdlib_complex64( 2.0f, 0.0f ), true, (stdlib_complex64_t *)x, strideX ); // Print the result: for ( int i = 0; i < N; i++ ) { - printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ i*2+1 ] ); + printf( "x[ %i ] = %f + %fi\n", i, x[ i*2 ], x[ (i*2)+1 ] ); } } diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js index 5e1f2692753d..df1309ecf546 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/lib/ndarray.js @@ -25,8 +25,13 @@ var realf = require( '@stdlib/complex/float32/real' ); var imagf = require( '@stdlib/complex/float32/imag' ); var powf = require( '@stdlib/math/base/special/powf' ); var lnf = require( '@stdlib/math/base/special/lnf' ); -var cosf = require( '@stdlib/math/base/special/cosf' ); -var sinf = require( '@stdlib/math/base/special/sinf' ); +var sincosf = require( '@stdlib/math/base/special/sincosf' ).assign; + + +// VARIABLES // + +// Workspace for storing computed sines and cosines: +var WORKSPACE = [ 0.0, 0.0 ]; // MAIN // @@ -34,6 +39,40 @@ var sinf = require( '@stdlib/math/base/special/sinf' ); /** * Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. * +* ## Method +* +* 1. If \\( b \\) is a positive real number, exponentiation with base \\( b \\) and complex exponent \\( z \\) is defined by means of the exponential function with a complex argument as +* +* ```tex +* b^z = e^{z \ln b} +* ``` +* +* where \\( \ln b \\) denotes the natural logarithm of \\( b \\). +* +* 2. Euler's formula states that, for any real number \\( y \\), one has +* +* ```tex +* e^{iy} = \cos y + i \sin y +* ``` +* +* 3. Consequently, Euler's formula allows expressing the polar form of \\( b^z \\) in terms of the real and imaginary components of \\( z = x + iy \\). +* +* ```tex +* b^{x + iy} = b^x ( \cos(y \ln b) + i \sin(y \ln b) ) +* ``` +* +* which follows from +* +* ```tex +* \begin{align*} +* b^{x + iy} &= b^x b^{iy} \\ +* &= b^x e^{iy \ln b} \\ +* &= b^x ( \cos(y \ln b) + i \sin(y \ln b) ) \\ +* \end{align*} +* ``` +* +* 4. Using Euler's formulation thus provides an update equation for generating each logarithmically spaced element in which we compute the sine and cosine of the imaginary component of the exponent and multiply by a real-valued scalar. +* * @param {PositiveInteger} N - number of indexed elements * @param {number} base - base of the logarithmic scale * @param {Complex64} start - exponent of the starting value @@ -80,6 +119,7 @@ function clogspace( N, base, start, stop, endpoint, x, strideX, offsetX ) { var lnb; var dre; var dim; + var sx; var ix; var i; @@ -99,26 +139,30 @@ function clogspace( N, base, start, stop, endpoint, x, strideX, offsetX ) { // Reinterpret the complex input array as a real-valued array: view = reinterpret( x, 0 ); - // Adjust the offset according to the real-valued array (2 floats per element): + // Adjust the stride and offset according to the real-valued array (2 floats per element): + sx = strideX * 2; ix = offsetX * 2; // Set the first value: if ( N === 1 ) { if ( endpoint ) { scale = powf( base, stopRe ); - view[ ix ] = scale * cosf( stopIm * lnb ); - view[ ix+1 ] = scale * sinf( stopIm * lnb ); + sincosf( stopIm * lnb, WORKSPACE, 1, 0 ); + view[ ix ] = scale * WORKSPACE[ 1 ]; + view[ ix+1 ] = scale * WORKSPACE[ 0 ]; } else { scale = powf( base, startRe ); - view[ ix ] = scale * cosf( startIm * lnb ); - view[ ix+1 ] = scale * sinf( startIm * lnb ); + sincosf( startIm * lnb, WORKSPACE, 1, 0 ); + view[ ix ] = scale * WORKSPACE[ 1 ]; + view[ ix+1 ] = scale * WORKSPACE[ 0 ]; } return x; } scale = powf( base, startRe ); - view[ ix ] = scale * cosf( startIm * lnb ); - view[ ix+1 ] = scale * sinf( startIm * lnb ); - ix += strideX * 2; + sincosf( startIm * lnb, WORKSPACE, 1, 0 ); + view[ ix ] = scale * WORKSPACE[ 1 ]; + view[ ix+1 ] = scale * WORKSPACE[ 0 ]; + ix += sx; // Calculate the complex increment: if ( endpoint ) { @@ -132,15 +176,17 @@ function clogspace( N, base, start, stop, endpoint, x, strideX, offsetX ) { expRe = startRe + ( dre * i ); expIm = startIm + ( dim * i ); scale = powf( base, expRe ); - view[ ix ] = scale * cosf( expIm * lnb ); - view[ ix+1 ] = scale * sinf( expIm * lnb ); - ix += strideX * 2; + sincosf( expIm * lnb, WORKSPACE, 1, 0 ); + view[ ix ] = scale * WORKSPACE[ 1 ]; + view[ ix+1 ] = scale * WORKSPACE[ 0 ]; + ix += sx; } // Check whether to include the `base^stop` value: if ( endpoint ) { scale = powf( base, stopRe ); - view[ ix ] = scale * cosf( stopIm * lnb ); - view[ ix+1 ] = scale * sinf( stopIm * lnb ); + sincosf( stopIm * lnb, WORKSPACE, 1, 0 ); + view[ ix ] = scale * WORKSPACE[ 1 ]; + view[ ix+1 ] = scale * WORKSPACE[ 0 ]; } return x; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json index ce26c7c8960b..d9c31cfdb1c8 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/manifest.json @@ -50,8 +50,7 @@ "@stdlib/strided/base/stride2offset", "@stdlib/math/base/special/powf", "@stdlib/math/base/special/lnf", - "@stdlib/math/base/special/cosf", - "@stdlib/math/base/special/sinf" + "@stdlib/math/base/special/sincosf" ] }, { @@ -72,8 +71,7 @@ "@stdlib/strided/base/stride2offset", "@stdlib/math/base/special/powf", "@stdlib/math/base/special/lnf", - "@stdlib/math/base/special/cosf", - "@stdlib/math/base/special/sinf" + "@stdlib/math/base/special/sincosf" ] }, { @@ -94,8 +92,7 @@ "@stdlib/strided/base/stride2offset", "@stdlib/math/base/special/powf", "@stdlib/math/base/special/lnf", - "@stdlib/math/base/special/cosf", - "@stdlib/math/base/special/sinf" + "@stdlib/math/base/special/sincosf" ] } ] diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json b/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json index 774b0177edf0..ffb7bcc29961 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/package.json @@ -74,7 +74,8 @@ "complex64", "single-precision", "float32", - "float32array" + "float32array", + "complex64array" ], "__stdlib__": { "wasm": false diff --git a/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c index 3ea3956fcfe6..b6a8a16efe4c 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/clogspace/src/main.c @@ -24,8 +24,7 @@ #include "stdlib/strided/base/stride2offset.h" #include "stdlib/math/base/special/powf.h" #include "stdlib/math/base/special/lnf.h" -#include "stdlib/math/base/special/cosf.h" -#include "stdlib/math/base/special/sinf.h" +#include "stdlib/math/base/special/sincosf.h" #include /** @@ -47,6 +46,40 @@ void API_SUFFIX(stdlib_strided_clogspace)( const CBLAS_INT N, const float base, /** * Fills a single-precision complex floating-point strided array with logarithmically spaced values over a specified interval using alternative indexing semantics. * +* ## Method +* +* 1. If \\( b \\) is a positive real number, exponentiation with base \\( b \\) and complex exponent \\( z \\) is defined by means of the exponential function with a complex argument as +* +* ```tex +* b^z = e^{z \ln b} +* ``` +* +* where \\( \ln b \\) denotes the natural logarithm of \\( b \\). +* +* 2. Euler's formula states that, for any real number \\( y \\), one has +* +* ```tex +* e^{iy} = \cos y + i \sin y +* ``` +* +* 3. Consequently, Euler's formula allows expressing the polar form of \\( b^z \\) in terms of the real and imaginary components of \\( z = x + iy \\). +* +* ```tex +* b^{x + iy} = b^x ( \cos(y \ln b) + i \sin(y \ln b) ) +* ``` +* +* which follows from +* +* ```tex +* \begin{align*} +* b^{x + iy} &= b^x b^{iy} \\ +* &= b^x e^{iy \ln b} \\ +* &= b^x ( \cos(y \ln b) + i \sin(y \ln b) ) \\ +* \end{align*} +* ``` +* +* 4. Using Euler's formulation thus provides an update equation for generating each logarithmically spaced element in which we compute the sine and cosine of the imaginary component of the exponent and multiply by a real-valued scalar. +* * @param N number of indexed elements * @param base base of the logarithmic scale * @param start exponent of the starting value @@ -61,15 +94,17 @@ void API_SUFFIX(stdlib_strided_clogspace_ndarray)( const CBLAS_INT N, const floa float start_im; float stop_re; float stop_im; + CBLAS_INT ix; float exp_re; float exp_im; float scale; float lnb; float dre; float dim; - CBLAS_INT ix; - CBLAS_INT M; - CBLAS_INT i; + float dc; + float ds; + float M; + float i; if ( N <= 0 ) { return; @@ -84,46 +119,50 @@ void API_SUFFIX(stdlib_strided_clogspace_ndarray)( const CBLAS_INT N, const floa // Compute ln(base) for use in base^z = exp(z * ln(base)): lnb = stdlib_base_lnf( base ); - ix = offsetX; - // Set the first value: + ix = offsetX; if ( N == 1 ) { if ( endpoint ) { scale = stdlib_base_powf( base, stop_re ); - X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( stop_im * lnb ), scale * stdlib_base_sinf( stop_im * lnb ) ); + stdlib_base_sincosf( stop_im * lnb, &ds, &dc ); + X[ ix ] = stdlib_complex64( scale * dc, scale * ds ); } else { scale = stdlib_base_powf( base, start_re ); - X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( start_im * lnb ), scale * stdlib_base_sinf( start_im * lnb ) ); + stdlib_base_sincosf( start_im * lnb, &ds, &dc ); + X[ ix ] = stdlib_complex64( scale * dc, scale * ds ); } return; } // Write the first value: scale = stdlib_base_powf( base, start_re ); - X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( start_im * lnb ), scale * stdlib_base_sinf( start_im * lnb ) ); + stdlib_base_sincosf( start_im * lnb, &ds, &dc ); + X[ ix ] = stdlib_complex64( scale * dc, scale * ds ); ix += strideX; // Calculate the complex increment: if ( endpoint ) { - M = N - 1; + M = (float)( N - 1 ); } else { - M = N; + M = (float)N; } - dre = ( stop_re - start_re ) / (float)M; - dim = ( stop_im - start_im ) / (float)M; + dre = ( stop_re - start_re ) / M; + dim = ( stop_im - start_im ) / M; // Generate logarithmically spaced values: - for ( i = 1; i < M; i++ ) { - exp_re = start_re + ( dre * (float)i ); - exp_im = start_im + ( dim * (float)i ); + for ( i = 1.0f; i < M; i += 1.0f ) { + exp_re = start_re + ( dre * i ); + exp_im = start_im + ( dim * i ); scale = stdlib_base_powf( base, exp_re ); - X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( exp_im * lnb ), scale * stdlib_base_sinf( exp_im * lnb ) ); + stdlib_base_sincosf( exp_im * lnb, &ds, &dc ); + X[ ix ] = stdlib_complex64( scale * dc, scale * ds ); ix += strideX; } // Check whether to include the `base^stop` value: if ( endpoint ) { scale = stdlib_base_powf( base, stop_re ); - X[ ix ] = stdlib_complex64( scale * stdlib_base_cosf( stop_im * lnb ), scale * stdlib_base_sinf( stop_im * lnb ) ); + stdlib_base_sincosf( stop_im * lnb, &ds, &dc ); + X[ ix ] = stdlib_complex64( scale * dc, scale * ds ); } return; }