From 0538961650558fa4de2def08141e5132ae260862 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Sun, 25 May 2025 10:39:51 +0530 Subject: [PATCH 01/18] feat: add C implementation for blas/base/zdscal --- .../blas/base/zdscal/src/zdscal_ndarray.c | 47 +++++++++++++++++++ 1 file changed, 47 insertions(+) create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c new file mode 100644 index 000000000000..1eaf30d1661e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -0,0 +1,47 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float64/base/scale.h" + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +* @param offsetX starting index for X +*/ +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + double *x = (double *)X; + CBLAS_INT ix; + CBLAS_INT sx; + CBLAS_INT i; + double tmp; + + if ( N <= 0 || alpha == 1.0 ) { + return; + } + sx = strideX * 2; + ix = offsetX * 2; + for ( i = 0; i < N; i++ ) { + } + return; +} From 0dc7cda814df856231eef8323e1689407bb4e77f Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Sun, 25 May 2025 23:12:34 +0530 Subject: [PATCH 02/18] chore: add implementation --- 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 status: passed - task: lint_package_json status: passed - task: lint_repl_help status: na - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - 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: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/base/zdscal/README.md | 130 +++++++ .../base/zdscal/benchmark/benchmark.native.js | 113 ++++++ .../benchmark/benchmark.ndarray.native.js | 113 ++++++ .../blas/base/zdscal/benchmark/c/Makefile | 146 ++++++++ .../zdscal/benchmark/c/benchmark.length.c | 192 +++++++++++ .../base/zdscal/benchmark/fortran/Makefile | 141 ++++++++ .../benchmark/fortran/benchmark.length.f | 210 ++++++++++++ .../@stdlib/blas/base/zdscal/binding.gyp | 265 ++++++++++++++ .../blas/base/zdscal/examples/c/Makefile | 146 ++++++++ .../blas/base/zdscal/examples/c/example.c | 48 +++ .../@stdlib/blas/base/zdscal/include.gypi | 70 ++++ .../zdscal/include/stdlib/blas/base/zdscal.h | 49 +++ .../include/stdlib/blas/base/zdscal_cblas.h | 44 +++ .../include/stdlib/blas/base/zdscal_fortran.h | 43 +++ .../@stdlib/blas/base/zdscal/lib/native.js | 35 ++ .../blas/base/zdscal/lib/ndarray.native.js | 57 ++++ .../blas/base/zdscal/lib/zdscal.native.js | 57 ++++ .../@stdlib/blas/base/zdscal/manifest.json | 322 ++++++++++++++++++ .../@stdlib/blas/base/zdscal/package.json | 3 + .../@stdlib/blas/base/zdscal/src/Makefile | 70 ++++ .../@stdlib/blas/base/zdscal/src/addon.c | 65 ++++ .../@stdlib/blas/base/zdscal/src/zdscal.c | 35 ++ .../@stdlib/blas/base/zdscal/src/zdscal.f | 80 +++++ .../blas/base/zdscal/src/zdscal_cblas.c | 50 +++ .../@stdlib/blas/base/zdscal/src/zdscal_f.c | 51 +++ .../blas/base/zdscal/src/zdscal_ndarray.c | 11 +- .../base/zdscal/test/test.ndarray.native.js | 276 +++++++++++++++ .../base/zdscal/test/test.zdscal.native.js | 226 ++++++++++++ 28 files changed, 3043 insertions(+), 5 deletions(-) create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/include.gypi create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/manifest.json create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/base/zdscal/test/test.zdscal.native.js diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/README.md b/lib/node_modules/@stdlib/blas/base/zdscal/README.md index 49d5fc189255..8bc72da40065 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/README.md +++ b/lib/node_modules/@stdlib/blas/base/zdscal/README.md @@ -149,6 +149,136 @@ console.log( zx.toString() ); + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/base/zdscal.h" +``` + +#### c_zdscal( N, alpha, \*X, strideX ) + +Scales a double-precision complex floating-point vector by a double-precision floating-point constant. + +```c +#include "stdlib/complex/float64/ctor.h" + +stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ) }; // interleaved real and imaginary components + +c_zdscal( 3, 2.0, x, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **alpha**: `[in] double` constant. +- **x**: `[inout] stdlib_complex128_t*` input array. +- **strideX**: `[in] CBLAS_INT` index increment for `x`. + +```c +void c_zdscal( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ); +``` + +#### c_zdscal_ndarray( N, alpha, \*X, strideX, offsetX ) + +Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. + +```c +stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ) }; // interleaved real and imaginary components + +c_zdscal_ndarray( 3, 2.0, x, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **alpha**: `[in] double` constant. +- **x**: `[inout] void*` input array. +- **strideX**: `[in] CBLAS_INT` index increment for `x`. +- **offsetX**: `[in] CBLAS_INT` starting index for `x`. + +```c +void c_zdscal_ndarray( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/blas/base/zdscal.h" +#include "stdlib/complex/float64/ctor.h" +#include "stdlib/complex/float64/real.h" +#include "stdlib/complex/float64/imag.h" +#include + +int main( void ) { + // Create strided arrays: + stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ), stdlib_complex128( 7.0, 8.0 ) }; + + // Specify the number of elements: + const int N = 4; + + // Specify stride lengths: + const int strideX = 1; + + c_zdscal( N, 2.0, (void *)x, strideX ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %lf + %lfj\n", i, stdlib_complex128_real( x[ i ] ), stdlib_complex128_imag( x[ i ] ) ); + } + + c_zdscal_ndarray( N, 2.0, (void *)x, strideX, 0 ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %lf + %lfj\n", i, stdlib_complex128_real( x[ i ] ), stdlib_complex128_imag( x[ i ] ) ); + } +} +``` + +
+ + + +
+ + +
diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js new file mode 100644 index 000000000000..5d7c2adce3cb --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js @@ -0,0 +1,113 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex128Array = require( '@stdlib/array/complex128' ); +var reinterpret = require( '@stdlib/strided/base/reinterpret-complex128' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var zdscal = tryRequire( resolve( __dirname, './../lib/zdscal.native.js' ) ); +var opts = { + 'skip': ( zdscal instanceof Error ) +}; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + + x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var viewX; + var i; + + viewX = reinterpret( x, 0 ); + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + zdscal( x.length, 2.0, x, 1 ); + if ( isnan( viewX[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( viewX[ 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( pkg+'::native:len='+len, opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..3abc62abe781 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js @@ -0,0 +1,113 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex128Array = require( '@stdlib/array/complex128' ); +var reinterpret = require( '@stdlib/strided/base/reinterpret-complex128' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var zdscal = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( zdscal instanceof Error ) +}; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + + x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var viewX; + var i; + + viewX = reinterpret( x, 0 ); + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + zdscal( x.length, 2.0, x, 1, 0 ); + if ( isnan( viewX[ i%(len*2) ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( viewX[ 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( pkg+'::native:ndarray:len='+len, opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile new file mode 100644 index 000000000000..9f97140e7cb0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2024 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/base/zdscal/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c new file mode 100644 index 000000000000..716a9ff3dee6 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c @@ -0,0 +1,192 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/complex/float64/ctor.h" +#include "stdlib/complex/float64/real.h" +#include "stdlib/complex/float64/imag.h" +#include +#include +#include +#include +#include + +#define NAME "zdscal" +#define ITERATIONS 10000000 +#define REPEATS 3 +#define MIN 1 +#define MAX 5 + +/** +* 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 stdlib_complex128_t rand_complex_128( void ) { + double real = (double)rand() / ((double)RAND_MAX + 1.0); + double imag = (double)rand() / ((double)RAND_MAX + 1.0); + return stdlib_complex128(real, imag); +} + +/** +* 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_complex128_t x[ len ]; + double elapsed; + double t; + int i; + + for ( i = 0; i < len; i++ ) { + x[ i ] = stdlib_complex128( ( stdlib_complex128_real( rand_complex_128() ) *10000.0 ) - 5000.0, ( stdlib_complex128_imag( rand_complex_128() ) *10000.0 ) - 5000.0 ) ; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + c_zdscal( len, 2.0, x, 1 ); + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + printf( "should not return NaN\n" ); + } + 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_complex128_t x[ len ]; + double elapsed; + double t; + int i; + + for ( i = 0; i < len; i++ ) { + x[ i ] = stdlib_complex128( ( stdlib_complex128_real( rand_complex_128() ) *10000.0 ) - 5000.0, ( stdlib_complex128_imag( rand_complex_128() ) *10000.0 ) - 5000.0 ) ; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + c_zdscal_ndarray( len, 2.0, x, 1, 0 ); + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + printf( "should not return NaN\n" ); + } + 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/base/zdscal/benchmark/fortran/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile new file mode 100644 index 000000000000..22405312ef4f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile @@ -0,0 +1,141 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2024 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 Fortran source files: +ifdef FORTRAN_COMPILER + FC := $(FORTRAN_COMPILER) +else + FC := gfortran +endif + +# Define the command-line options when compiling Fortran files: +FFLAGS ?= \ + -std=f95 \ + -ffree-form \ + -O3 \ + -Wall \ + -Wextra \ + -Wno-compare-reals \ + -Wimplicit-interface \ + -fno-underscoring \ + -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 ?= + +# List of Fortran source files: +SOURCE_FILES ?= ../../src/zdrot.f + +# List of Fortran targets: +f_targets := benchmark.length.out + + +# RULES # + +#/ +# Compiles Fortran source files. +# +# @param {string} SOURCE_FILES - list of Fortran source files +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop`) +# @param {string} [FORTRAN_COMPILER] - Fortran compiler +# @param {string} [FFLAGS] - Fortran compiler flags +# @param {(string|void)} [fPIC] - compiler flag indicating whether to generate position independent code +# +# @example +# make +# +# @example +# make all +#/ +all: $(f_targets) + +.PHONY: all + +#/ +# Compiles Fortran source files. +# +# @private +# @param {string} SOURCE_FILES - list of Fortran source files +# @param {(string|void)} INCLUDE - list of includes (e.g., `-I /foo/bar -I /beep/boop`) +# @param {string} FC - Fortran compiler +# @param {string} FFLAGS - Fortran compiler flags +# @param {(string|void)} fPIC - compiler flag indicating whether to generate position independent code +#/ +$(f_targets): %.out: %.f + $(QUIET) $(FC) $(FFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< + +#/ +# Runs compiled benchmarks. +# +# @example +# make run +#/ +run: $(f_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/base/zdscal/benchmark/fortran/benchmark.length.f b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f new file mode 100644 index 000000000000..382fb5bf2fce --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f @@ -0,0 +1,210 @@ +!> +! @license Apache-2.0 +! +! Copyright (c) 2024 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. +!< + +program bench + implicit none + ! .. + ! Local constants: + character(5), parameter :: name = 'zdscal' ! if changed, be sure to adjust length + integer, parameter :: iterations = 10000000 + integer, parameter :: repeats = 3 + integer, parameter :: min = 1 + integer, parameter :: max = 6 + ! .. + ! Run the benchmarks: + call main() + ! .. + ! Functions: +contains + ! .. + ! Prints the TAP version. + ! .. + subroutine print_version() + print '(A)', 'TAP version 13' + end subroutine print_version + ! .. + ! Prints the TAP summary. + ! + ! @param {integer} total - total number of tests + ! @param {integer} passing - total number of passing tests + ! .. + subroutine print_summary( total, passing ) + ! .. + ! Scalar arguments: + integer, intent(in) :: total, passing + ! .. + ! Local variables: + character(len=999) :: str, tmp + ! .. + ! Intrinsic functions: + intrinsic adjustl, trim + ! .. + print '(A)', '#' + ! .. + write (str, '(I15)') total ! TAP plan + tmp = adjustl( str ) + print '(A,A)', '1..', trim( tmp ) + ! .. + print '(A,A)', '# total ', trim( tmp ) + ! .. + write (str, '(I15)') passing + tmp = adjustl( str ) + print '(A,A)', '# pass ', trim( tmp ) + ! .. + print '(A)', '#' + print '(A)', '# ok' + end subroutine print_summary + ! .. + ! Prints benchmarks results. + ! + ! @param {integer} iterations - number of iterations + ! @param {double} elapsed - elapsed time in seconds + ! .. + subroutine print_results( iterations, elapsed ) + ! .. + ! Scalar arguments: + double precision, intent(in) :: elapsed + integer, intent(in) :: iterations + ! .. + ! Local variables: + double precision :: rate + character(len=999) :: str, tmp + ! .. + ! Intrinsic functions: + intrinsic dble, adjustl, trim + ! .. + rate = dble( iterations ) / elapsed + ! .. + print '(A)', ' ---' + ! .. + write (str, '(I15)') iterations + tmp = adjustl( str ) + print '(A,A)', ' iterations: ', trim( tmp ) + ! .. + write (str, '(f0.9)') elapsed + tmp = adjustl( str ) + print '(A,A)', ' elapsed: ', trim( tmp ) + ! .. + write( str, '(f0.9)') rate + tmp = adjustl( str ) + print '(A,A)', ' rate: ', trim( tmp ) + ! .. + print '(A)', ' ...' + end subroutine print_results + ! .. + ! Runs a benchmark. + ! + ! @param {integer} iterations - number of iterations + ! @param {integer} len - array length + ! @return {double} elapsed time in seconds + ! .. + double precision function benchmark( iterations, len ) + ! .. + ! External functions: + interface + subroutine zdscal( N, alpha, x, strideX ) + complex(kind=kind(0.0d0)) :: x(*) + integer :: strideX, N + double precision :: alpha + end subroutine zdscal + end interface + ! .. + ! Scalar arguments: + integer, intent(in) :: iterations, len + ! .. + ! Local scalars: + double precision :: elapsed, r1, r2 + real :: t1, t2 + integer :: i + ! .. + ! Local arrays: + complex(kind=kind(0.0d0)), allocatable :: x(:) + ! .. + ! Intrinsic functions: + intrinsic random_number, cpu_time, cmplx + ! .. + ! Allocate arrays: + allocate( x(len) ) + ! .. + do i = 1, len + call random_number( r1 ) + call random_number( r2 ) + x( i ) = cmplx( (r1*10000.0d0)-5000.0d0, (r2*10000.0d0)-5000.0d0, kind=kind(0.0d0) ) + end do + ! .. + call cpu_time( t1 ) + ! .. + do i = 1, iterations + call zdscal( len, 2.0d0, x, 1 ); + if ( x( 1 ) /= x( 1 ) ) then + print '(A)', 'should not return NaN' + exit + end if + end do + ! .. + call cpu_time( t2 ) + ! .. + elapsed = t2 - t1 + ! .. + if ( x( 1 ) /= x( 1 ) ) then + print '(A)', 'should not return NaN' + end if + ! .. + ! Deallocate arrays: + deallocate( x ) + ! .. + benchmark = elapsed + return + end function benchmark + ! .. + ! Main execution sequence. + ! .. + subroutine main() + ! .. + ! Local variables: + integer :: count, iter, len, i, j + double precision :: elapsed + character(len=999) :: str, tmp + ! .. + ! Intrinsic functions: + intrinsic adjustl, trim + ! .. + call print_version() + count = 0 + do i = min, max + len = 10**i + iter = iterations / 10**(i-1) + do j = 1, repeats + count = count + 1 + ! .. + write (str, '(I15)') len + tmp = adjustl( str ) + print '(A,A,A,A)', '# fortran::', name, ':len=', trim( tmp ) + ! .. + elapsed = benchmark( iter, len ) + ! .. + call print_results( iter, elapsed ) + ! .. + write (str, '(I15)') count + tmp = adjustl( str ) + print '(A,A,A)', 'ok ', trim( tmp ), ' benchmark finished' + end do + end do + call print_summary( count, count ) + end subroutine main +end program bench diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp b/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp new file mode 100644 index 000000000000..02a2799da097 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp @@ -0,0 +1,265 @@ +# @license Apache-2.0 +# +# Copyright (c) 2024 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', + + # Fortran compiler (to override -Dfortran_compiler=): + 'fortran_compiler%': 'gfortran', + + # Fortran compiler flags: + 'fflags': [ + # Specify the Fortran standard to which a program is expected to conform: + '-std=f95', + + # Indicate that the layout is free-form source code: + '-ffree-form', + + # Aggressive optimization: + '-O3', + + # Enable commonly used warning options: + '-Wall', + + # Warn if source code contains problematic language features: + '-Wextra', + + # Warn if a procedure is called without an explicit interface: + '-Wimplicit-interface', + + # Do not transform names of entities specified in Fortran source files by appending underscores (i.e., don't mangle names, thus allowing easier usage in C wrappers): + '-fno-underscoring', + + # Warn if source code contains Fortran 95 extensions and C-language constructs: + '-pedantic', + + # Compile but do not link (output is an object file): + '-c', + ], + + # 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 + + # Define custom build actions for particular inputs: + 'rules': [ + { + # Define a rule for processing Fortran files: + 'extension': 'f', + + # Define the pathnames to be used as inputs when performing processing: + 'inputs': [ + # Full path of the current input: + '<(RULE_INPUT_PATH)' + ], + + # Define the outputs produced during processing: + 'outputs': [ + # Store an output object file in a directory for placing intermediate results (only accessible within a single target): + '<(INTERMEDIATE_DIR)/<(RULE_INPUT_ROOT).<(obj)' + ], + + # Define the rule for compiling Fortran based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + + # Rule to compile Fortran on Windows: + { + 'rule_name': 'compile_fortran_windows', + 'message': 'Compiling Fortran file <(RULE_INPUT_PATH) on Windows...', + + 'process_outputs_as_sources': 0, + + # Define the command-line invocation: + 'action': [ + '<(fortran_compiler)', + '<@(fflags)', + '<@(_inputs)', + '-o', + '<@(_outputs)', + ], + }, + + # Rule to compile Fortran on non-Windows: + { + 'rule_name': 'compile_fortran_linux', + 'message': 'Compiling Fortran file <(RULE_INPUT_PATH) on Linux...', + + 'process_outputs_as_sources': 1, + + # Define the command-line invocation: + 'action': [ + '<(fortran_compiler)', + '<@(fflags)', + '-fPIC', # generate platform-independent code + '<@(_inputs)', + '-o', + '<@(_outputs)', + ], + } + ], # end condition (OS=="win") + ], # end conditions + }, # end rule (extension=="f") + ], # end rules + }, # 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/base/zdscal/examples/c/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile new file mode 100644 index 000000000000..6aed70daf167 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2024 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/base/zdscal/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c new file mode 100644 index 000000000000..a55b78aebc99 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c @@ -0,0 +1,48 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/complex/float64/ctor.h" +#include "stdlib/complex/float64/real.h" +#include "stdlib/complex/float64/imag.h" +#include + +int main( void ) { + // Create strided arrays: + stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ), stdlib_complex128( 7.0, 8.0 ) }; + + // Specify the number of elements: + const int N = 4; + + // Specify stride lengths: + const int strideX = 1; + + c_zdscal( N, 2.0, (void *)x, strideX ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %lf + %lfj\n", i, stdlib_complex128_real( x[ i ] ), stdlib_complex128_imag( x[ i ] ) ); + } + + c_zdscal_ndarray( N, 2.0, (void *)x, strideX, 0 ); + + // Print the result: + for ( int i = 0; i < N; i++ ) { + printf( "x[ %i ] = %lf + %lfj\n", i, stdlib_complex128_real( x[ i ] ), stdlib_complex128_imag( x[ i ] ) ); + } +} diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi b/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi new file mode 100644 index 000000000000..497aeca15320 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi @@ -0,0 +1,70 @@ +# @license Apache-2.0 +# +# Copyright (c) 2024 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. +# +# Note that nesting variables is required due to how GYP processes a configuration. Any variables defined within a nested 'variables' section is defined in the outer scope. Thus, conditions in the outer variable scope are free to use these variables without running into "variable undefined" errors. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +# +# Variable nesting hacks: +# +# [3]: https://chromium.googlesource.com/external/skia/gyp/+/master/common_variables.gypi +# [4]: https://src.chromium.org/viewvc/chrome/trunk/src/build/common.gypi?revision=127004 +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + 'variables': { + # Host BLAS library (to override -Dblas=): + 'blas%': '', + + # Path to BLAS library (to override -Dblas_dir=): + 'blas_dir%': '', + }, # end variables + + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + '<@(blas_dir)', + ' [ 2.0, 4.0, 6.0, 8.0, 10.0, 12.0 ] +*/ +function zdscal( N, alpha, x, strideX, offsetX ) { + var viewX = reinterpret( x, 0 ); + addon.ndarray( N, alpha, viewX, strideX, offsetX ); + return x; +} + + +// EXPORTS // + +module.exports = zdscal; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js new file mode 100644 index 000000000000..df219d310835 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js @@ -0,0 +1,57 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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-complex128' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} alpha - constant +* @param {Complex128Array} x - input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting `x` index +* @returns {Complex128Array} input array +* +* @example +* var Complex128Array = require( '@stdlib/array/complex128' ); +* var Complex128 = require( '@stdlib/complex/float64/ctor' ); +* +* var x = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); +* +* zdscal( 3, 2.0, x, 1, 0 ); +* // x => [ 2.0, 4.0, 6.0, 8.0, 10.0, 12.0 ] +*/ +function zdscal( N, alpha, x, strideX ) { + var viewX = reinterpret( x, 0 ); + addon( N, alpha, viewX, strideX ); + return x; +} + + +// EXPORTS // + +module.exports = zdscal; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json new file mode 100644 index 000000000000..4490c55137c8 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json @@ -0,0 +1,322 @@ +{ + "options": { + "task": "build", + "os": "linux", + "blas": "", + "wasm": false + }, + "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", + "os": "win", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/napi/argv-double" + ] + }, + { + "task": "benchmark", + "os": "win", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + { + "task": "examples", + "os": "win", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/napi/argv-double" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/napi/argv-double" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/zdscal_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared" + ] + }, + + { + "task": "build", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/napi/argv-double", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" + ] + }, + { + "task": "benchmark", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/math/base/assert/is-nan", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", + "@stdlib/complex/float64/base/scale" + ] + }, + { + "task": "examples", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", + "@stdlib/complex/float64/base/scale" + ] + }, + + { + "task": "build", + "os": "", + "blas": "", + "wasm": true, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/package.json b/lib/node_modules/@stdlib/blas/base/zdscal/package.json index 2107fd81e8c7..8e53524825b0 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/package.json +++ b/lib/node_modules/@stdlib/blas/base/zdscal/package.json @@ -15,11 +15,14 @@ ], "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", diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile new file mode 100644 index 000000000000..bcf18aa46655 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2024 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/base/zdscal/src/addon.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c new file mode 100644 index 000000000000..f3b7e8d66ed0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c @@ -0,0 +1,65 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_double.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_strided_complex128array.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, 4 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_DOUBLE( env, alpha, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 ); + stdlib_complex128_t *x = (stdlib_complex128_t *)X; + API_SUFFIX(c_zdscal)( N, alpha, 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, 5 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); + STDLIB_NAPI_ARGV_DOUBLE( env, alpha, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 ); + stdlib_complex128_t *x = (stdlib_complex128_t *)X; + API_SUFFIX(c_zdscal_ndarray)( N, alpha, x, strideX, offsetX ); + return NULL; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c new file mode 100644 index 000000000000..7b79f98f46e8 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/strided/base/stride2offset.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float64/ctor.h" + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +*/ +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { + CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + API_SUFFIX(c_zdscal_ndarray)( N, alpha, X, strideX, ox ); +} diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f new file mode 100644 index 000000000000..3a1e05b6bb9c --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f @@ -0,0 +1,80 @@ +!> +! @license Apache-2.0 +! +! Copyright (c) 2024 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. +!< + +!> Scales a double-precision complex floating-point vector by a double-precision floating-point constant. +! +! ## Notes +! +! * Modified version of reference BLAS level1 routine (version 3.9.0). Updated to "free form" Fortran 95. +! +! ## Authors +! +! * Univ. of Tennessee +! * Univ. of California Berkeley +! * Univ. of Colorado Denver +! * NAG Ltd. +! +! ## License +! +! From : +! +! > The reference BLAS is a freely-available software package. It is available from netlib via anonymous ftp and the World Wide Web. Thus, it can be included in commercial software packages (and has been). We only ask that proper credit be given to the authors. +! > +! > Like all software, it is copyrighted. It is not trademarked, but we do ask the following: +! > +! > * If you modify the source for these routines we ask that you change the name of the routine and comment the changes made to the original. +! > +! > * We will gladly answer any questions regarding the software. If a modification is done, however, it is the responsibility of the person who modified the routine to provide support. +! +! @param {integer} N - number of indexed elements +! @param {double} alpha - constant +! @param {Array>} x - input array +! @param {integer} strideX - `x` stride length +!< +subroutine zdscal( N, alpha, x, strideX ) + implicit none + ! .. + ! Scalar arguments: + integer :: strideX, N + double precision :: alpha + ! .. + ! Array arguments: + complex(kind=kind(0.0d0)) :: x(*) + ! .. + ! Local scalars: + integer :: ix, i, nstrideX + ! .. + ! Intrinsic functions: + intrinsic dble, dcmplx, dimag + ! .. + if ( N <= 0 .OR. strideX <= 0 .OR. alpha == 1.0 ) then + return + end if + ! .. + if ( strideX == 1 ) then + do i = 1, N + x( i ) = dcmplx(alpha*dble(x(i)),alpha*dimag(x(i))) + end do + else + nstrideX = N * strideX; + do i = 1, nstrideX, strideX + x( i ) = dcmplx(alpha*dble(x(i)),alpha*dimag(x(i))) + end do + end if + return +end subroutine zdscal diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c new file mode 100644 index 000000000000..6cf9c64ade52 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -0,0 +1,50 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/blas/base/zdscal_cblas.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float64/ctor.h" +#include "stdlib/strided/base/min_view_buffer_index.h" + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +*/ +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { + API_SUFFIX(cblas_zdscal)( N, alpha, X, strideX ); +} + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +* @param offsetX starting index for X +*/ +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + stdlib_complex128_t *x = (stdlib_complex128_t *)X; + x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer + API_SUFFIX(cblas_zdscal)( N, alpha, x, strideX ); +} diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c new file mode 100644 index 000000000000..37546e53dc22 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -0,0 +1,51 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2024 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/base/zdscal.h" +#include "stdlib/blas/base/zdscal_fortran.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float64/ctor.h" +#include "stdlib/strided/base/min_view_buffer_index.h" + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +*/ +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { + zdscal( &N, alpha, X, &strideX ); +} + +/** +* Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. +* +* @param N number of indexed elements +* @param alpha constant +* @param X input array +* @param strideX X stride length +* @param offsetX starting index for X +*/ +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + stdlib_complex128_t *x = (stdlib_complex128_t *)X; + + x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); + zdscal( &N, alpha, x, &strideX ); +} diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c index 1eaf30d1661e..18b421b8bf8e 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -18,6 +18,7 @@ #include "stdlib/blas/base/zdscal.h" #include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float64/ctor.h" #include "stdlib/complex/float64/base/scale.h" /** @@ -29,19 +30,19 @@ * @param strideX X stride length * @param offsetX starting index for X */ -void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { - double *x = (double *)X; +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { CBLAS_INT ix; CBLAS_INT sx; CBLAS_INT i; - double tmp; if ( N <= 0 || alpha == 1.0 ) { return; } - sx = strideX * 2; - ix = offsetX * 2; + sx = strideX; + ix = offsetX; for ( i = 0; i < N; i++ ) { + X[ ix ] = stdlib_base_complex128_scale( alpha, X[ ix ] ); + ix += sx; } return; } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js new file mode 100644 index 000000000000..b4526d5c5670 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js @@ -0,0 +1,276 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 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 isSameComplex128Array = require( '@stdlib/assert/is-same-complex128array' ); +var Complex128Array = require( '@stdlib/array/complex128' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var zdscal = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( zdscal instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof zdscal, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 5', opts, function test( t ) { + t.strictEqual( zdscal.length, 5, 'arity of 5' ); + t.end(); +}); + +tape( 'the function scales elements from `zx` by `da`', opts, function test( t ) { + var expected; + var zx; + + zx = new Complex128Array([ + 0.3, // 1 + 0.1, // 1 + 0.5, // 2 + 0.0, // 2 + 0.0, // 3 + 0.5, // 3 + 0.0, // 4 + 0.2, // 4 + 2.0, + 3.0, + 2.0, + 3.0 + ]); + + zdscal( 4, 2.0, zx, 1, 0 ); + + expected = new Complex128Array([ + 0.6, // 1 + 0.2, // 1 + 1.0, // 2 + 0.0, // 2 + 0.0, // 3 + 1.0, // 3 + 0.0, // 4 + 0.4, // 4 + 2.0, + 3.0, + 2.0, + 3.0 + ]); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `zx` stride', opts, function test( t ) { + var expected; + var zx; + + zx = new Complex128Array([ + 0.1, // 1 + 0.1, // 1 + 3.0, + 6.0, + -0.6, // 2 + 0.1, // 2 + 4.0, + 7.0, + 0.1, // 3 + -0.3, // 3 + 7.0, + 2.0 + ]); + + zdscal( 3, 2.0, zx, 2, 0 ); + + expected = new Complex128Array([ + 0.2, // 1 + 0.2, // 1 + 3.0, + 6.0, + -1.2, // 2 + 0.2, // 2 + 4.0, + 7.0, + 0.2, // 3 + -0.6, // 3 + 7.0, + 2.0 + ]); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `zx` stride', opts, function test( t ) { + var expected; + var zx; + + zx = new Complex128Array([ + 0.1, // 3 + 0.1, // 3 + 3.0, + 6.0, + -0.6, // 2 + 0.1, // 2 + 4.0, + 7.0, + 0.1, // 1 + -0.3, // 1 + 7.0, + 2.0 + ]); + + zdscal( 3, 2.0, zx, -2, 4 ); + + expected = new Complex128Array([ + 0.2, // 3 + 0.2, // 3 + 3.0, + 6.0, + -1.2, // 2 + 0.2, // 2 + 4.0, + 7.0, + 0.2, // 1 + -0.6, // 1 + 7.0, + 2.0 + ]); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `zx` offset', opts, function test( t ) { + var expected; + var zx; + + zx = new Complex128Array([ + 0.1, + 0.1, + 3.0, + 6.0, + -0.6, + 0.1, + 4.0, // 1 + 6.0, // 1 + 0.1, // 2 + -0.3, // 2 + 7.0, // 3 + 2.0 // 3 + ]); + + zdscal( 3, 2.0, zx, 1, 3 ); + + expected = new Complex128Array([ + 0.1, + 0.1, + 3.0, + 6.0, + -0.6, + 0.1, + 8.0, // 1 + 12.0, // 1 + 0.2, // 2 + -0.6, // 2 + 14.0, // 3 + 4.0 // 3 + ]); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns a reference to the input array', opts, function test( t ) { + var out; + var zx; + + zx = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + out = zdscal( 4, 2.0, zx, 1, 0 ); + + t.strictEqual( out, zx, '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 zx; + + zx = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + + expected = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); // eslint-disable-line max-len + + zdscal( -1, 2.0, zx, 1, 0 ); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + + zdscal( 0, 2.0, zx, 1, 0 ); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports complex access patterns', opts, function test( t ) { + var expected; + var zx; + + zx = new Complex128Array([ + 0.1, + 0.1, + 3.0, + 6.0, + -0.6, + 0.1, + 4.0, // 2 + 6.0, // 2 + 0.1, + -0.3, + 7.0, + 2.0, + 2.0, // 1 + 3.0 // 1 + ]); + + zdscal( 2, 2.0, zx, -3, 6 ); + + expected = new Complex128Array([ + 0.1, + 0.1, + 3.0, + 6.0, + -0.6, + 0.1, + 8.0, // 2 + 12.0, // 2 + 0.1, + -0.3, + 7.0, + 2.0, + 4.0, // 1 + 6.0 // 1 + ]); + t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/test/test.zdscal.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.zdscal.native.js new file mode 100644 index 000000000000..41bbe6f74e02 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.zdscal.native.js @@ -0,0 +1,226 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 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 isSameComplex128Array = require( '@stdlib/assert/is-same-complex128array' ); +var Complex128Array = require( '@stdlib/array/complex128' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var zdscal = tryRequire( resolve( __dirname, './../lib/zdscal.native.js' ) ); +var opts = { + 'skip': ( zdscal instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof zdscal, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 4', opts, function test( t ) { + t.strictEqual( zdscal.length, 4, 'arity of 4' ); + t.end(); +}); + +tape( 'the function scales elements from `x` by `alpha`', opts, function test( t ) { + var expected; + var x; + + x = new Complex128Array([ + 0.3, // 1 + 0.1, // 1 + 0.5, // 2 + 0.0, // 2 + 0.0, // 3 + 0.5, // 3 + 0.0, // 4 + 0.2 // 4 + ]); + + zdscal( 4, 2.0, x, 1 ); + + expected = new Complex128Array([ + 0.6, // 1 + 0.2, // 1 + 1.0, // 2 + 0.0, // 2 + 0.0, // 3 + 1.0, // 3 + 0.0, // 4 + 0.4 // 4 + ]); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex128Array([ + 0.1, // 1 + 0.1, // 1 + 3.0, + 6.0, + -0.6, // 2 + 0.1, // 2 + 4.0, + 7.0, + 0.1, // 3 + -0.3, // 3 + 7.0, + 2.0 + ]); + + zdscal( 3, 2.0, x, 2 ); + + expected = new Complex128Array([ + 0.2, // 1 + 0.2, // 1 + 3.0, + 6.0, + -1.2, // 2 + 0.2, // 2 + 4.0, + 7.0, + 0.2, // 3 + -0.6, // 3 + 7.0, + 2.0 + ]); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `x` stride', opts, function test( t ) { + var expected; + var x; + + x = new Complex128Array([ + 0.1, // 3 + 0.1, // 3 + 3.0, + 6.0, + -0.6, // 2 + 0.1, // 2 + 4.0, + 7.0, + 0.1, // 1 + -0.3, // 1 + 7.0, + 2.0 + ]); + + zdscal( 3, 2.0, x, -2 ); + + expected = new Complex128Array([ + 0.2, // 3 + 0.2, // 3 + 3.0, + 6.0, + -1.2, // 2 + 0.2, // 2 + 4.0, + 7.0, + 0.2, // 1 + -0.6, // 1 + 7.0, + 2.0 + ]); + t.strictEqual( isSameComplex128Array( 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 Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + out = zdscal( 4, 2.0, 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 Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + + expected = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); // eslint-disable-line max-len + + zdscal( -1, 2.0, x, 1 ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); + + zdscal( 0, 2.0, x, 1 ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports view offsets', opts, function test( t ) { + var expected; + var x0; + var x1; + + // Initial array: + x0 = new Complex128Array([ + 0.1, + -0.3, + 8.0, // 1 + 9.0, // 1 + 0.5, // 2 + -0.1, // 2 + 2.0, // 3 + 5.0, // 3 + 2.0, + 3.0 + ]); + + // Create offset view: + x1 = new Complex128Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // begin at 2nd element + + zdscal( 3, 2.0, x1, 1 ); + + expected = new Complex128Array([ + 0.1, + -0.3, + 16.0, // 1 + 18.0, // 1 + 1.0, // 2 + -0.2, // 2 + 4.0, // 3 + 10.0, // 3 + 2.0, + 3.0 + ]); + t.strictEqual( isSameComplex128Array( x0, expected ), true, 'returns expected value' ); + t.end(); +}); From 586ab55d467ccc2d6a56f409f420dd083814e879 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Sun, 25 May 2025 23:19:55 +0530 Subject: [PATCH 03/18] chore: add implementation --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/base/zdscal/benchmark/fortran/Makefile | 2 +- .../blas/base/zdscal/benchmark/fortran/benchmark.length.f | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f | 6 +++--- 3 files changed, 5 insertions(+), 5 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile index 22405312ef4f..b06f29ef1f9a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile @@ -78,7 +78,7 @@ endif INCLUDE ?= # List of Fortran source files: -SOURCE_FILES ?= ../../src/zdrot.f +SOURCE_FILES ?= ../../src/zdscal.f # List of Fortran targets: f_targets := benchmark.length.out diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f index 382fb5bf2fce..b95e07c976e9 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f @@ -20,7 +20,7 @@ program bench implicit none ! .. ! Local constants: - character(5), parameter :: name = 'zdscal' ! if changed, be sure to adjust length + character(6), parameter :: name = 'zdscal' ! if changed, be sure to adjust length integer, parameter :: iterations = 10000000 integer, parameter :: repeats = 3 integer, parameter :: min = 1 diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f index 3a1e05b6bb9c..2bc980b14df9 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f @@ -60,7 +60,7 @@ subroutine zdscal( N, alpha, x, strideX ) integer :: ix, i, nstrideX ! .. ! Intrinsic functions: - intrinsic dble, dcmplx, dimag + intrinsic dble, cmplx, aimag ! .. if ( N <= 0 .OR. strideX <= 0 .OR. alpha == 1.0 ) then return @@ -68,12 +68,12 @@ subroutine zdscal( N, alpha, x, strideX ) ! .. if ( strideX == 1 ) then do i = 1, N - x( i ) = dcmplx(alpha*dble(x(i)),alpha*dimag(x(i))) + x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i))) end do else nstrideX = N * strideX; do i = 1, nstrideX, strideX - x( i ) = dcmplx(alpha*dble(x(i)),alpha*dimag(x(i))) + x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i))) end do end if return From c5280c8f8b4945c3035990f75175cb4986c34ef5 Mon Sep 17 00:00:00 2001 From: stdlib-bot <82920195+stdlib-bot@users.noreply.github.com> Date: Sun, 25 May 2025 17:52:41 +0000 Subject: [PATCH 04/18] chore: update copyright years --- .../@stdlib/blas/base/zdscal/benchmark/benchmark.native.js | 2 +- .../blas/base/zdscal/benchmark/benchmark.ndarray.native.js | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile | 2 +- .../@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c | 2 +- .../@stdlib/blas/base/zdscal/benchmark/fortran/Makefile | 2 +- .../blas/base/zdscal/benchmark/fortran/benchmark.length.f | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/include.gypi | 2 +- .../@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h | 2 +- .../blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h | 2 +- .../blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c | 2 +- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c | 2 +- 23 files changed, 23 insertions(+), 23 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js index 5d7c2adce3cb..71694db0edde 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js index 3abc62abe781..bf25815c4ebe 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile index 9f97140e7cb0..cce2c865d7ad 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c index 716a9ff3dee6..eaf11c32d6f5 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile index b06f29ef1f9a..7b2a53414298 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f index b95e07c976e9..699a75f41c46 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/fortran/benchmark.length.f @@ -1,7 +1,7 @@ !> ! @license Apache-2.0 ! -! Copyright (c) 2024 The Stdlib Authors. +! Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp b/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp index 02a2799da097..08de71a2020e 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp +++ b/lib/node_modules/@stdlib/blas/base/zdscal/binding.gyp @@ -1,6 +1,6 @@ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile index 6aed70daf167..25ced822f96a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile +++ b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c index a55b78aebc99..eb27698a3606 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi b/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi index 497aeca15320..4217944b5d20 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include.gypi @@ -1,6 +1,6 @@ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h index 1d8269b00143..36a4cd8fa62f 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h index ea203d5c8b1f..0a7a8779960d 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h index 528b5f3fed65..fc9497d7b1bf 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_fortran.h @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js index ade5d91a44c1..4badc2e583b8 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js index e66e7429e00b..6046cfe9f480 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js index df219d310835..e36a0e48dd2b 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile b/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile index bcf18aa46655..7733b6180cb4 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2024 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c index f3b7e8d66ed0..127887792c0e 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c index 7b79f98f46e8..673e0a9a16f9 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f index 2bc980b14df9..e93c1c950ee2 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f @@ -1,7 +1,7 @@ !> ! @license Apache-2.0 ! -! Copyright (c) 2024 The Stdlib Authors. +! Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c index 6cf9c64ade52..b18784da9a06 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index 37546e53dc22..c375974d7ea6 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c index 18b421b8bf8e..90af1f2484bc 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2024 The Stdlib Authors. +* Copyright (c) 2025 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. From 11629ecb953c918fe5a148aa2318a0fe5c04d891 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Tue, 27 May 2025 21:37:28 +0530 Subject: [PATCH 05/18] chore: update implementation --- 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 status: passed - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - 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: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/base/zdscal/README.md | 21 +- .../base/zdscal/benchmark/benchmark.native.js | 4 +- .../benchmark/benchmark.ndarray.native.js | 4 +- .../zdscal/benchmark/c/benchmark.length.c | 19 +- .../blas/base/zdscal/examples/c/example.c | 7 +- .../zdscal/include/stdlib/blas/base/zdscal.h | 4 +- .../include/stdlib/blas/base/zdscal_cblas.h | 2 +- .../@stdlib/blas/base/zdscal/manifest.json | 199 +++++++++++++++--- .../@stdlib/blas/base/zdscal/src/addon.c | 6 +- .../@stdlib/blas/base/zdscal/src/zdscal.c | 3 +- .../blas/base/zdscal/src/zdscal_cblas.c | 4 +- .../@stdlib/blas/base/zdscal/src/zdscal_f.c | 1 - .../blas/base/zdscal/src/zdscal_ndarray.c | 5 +- 13 files changed, 220 insertions(+), 59 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/README.md b/lib/node_modules/@stdlib/blas/base/zdscal/README.md index 8bc72da40065..bafcd446aa95 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/README.md +++ b/lib/node_modules/@stdlib/blas/base/zdscal/README.md @@ -182,7 +182,11 @@ Scales a double-precision complex floating-point vector by a double-precision fl ```c #include "stdlib/complex/float64/ctor.h" -stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ) }; // interleaved real and imaginary components +stdlib_complex128_t x[] = { + stdlib_complex128( 1.0, 2.0 ), + stdlib_complex128( 3.0, 4.0 ), + stdlib_complex128( 5.0, 6.0 ) +}; c_zdscal( 3, 2.0, x, 1 ); ``` @@ -203,7 +207,13 @@ void c_zdscal( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, co Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. ```c -stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ) }; // interleaved real and imaginary components +#include "stdlib/complex/float64/ctor.h" + +stdlib_complex128_t x[] = { + stdlib_complex128( 1.0, 2.0 ), + stdlib_complex128( 3.0, 4.0 ), + stdlib_complex128( 5.0, 6.0 ) +}; c_zdscal_ndarray( 3, 2.0, x, 1, 0 ); ``` @@ -247,7 +257,12 @@ void c_zdscal_ndarray( const CBLAS_INT N, const double alpha, stdlib_complex128_ int main( void ) { // Create strided arrays: - stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ), stdlib_complex128( 7.0, 8.0 ) }; + stdlib_complex128_t x[] = { + stdlib_complex128( 1.0, 2.0 ), + stdlib_complex128( 3.0, 4.0 ), + stdlib_complex128( 5.0, 6.0 ), + stdlib_complex128( 7.0, 8.0 ) + }; // Specify the number of elements: const int N = 4; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js index 71694db0edde..20df78899a50 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.native.js @@ -52,9 +52,7 @@ var options = { * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x; - - x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); + var x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); return benchmark; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js index bf25815c4ebe..b03e31c087ea 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/benchmark.ndarray.native.js @@ -52,9 +52,7 @@ var options = { * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x; - - x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); + var x = new Complex128Array( uniform( len*2, -100.0, 100.0, options ) ); return benchmark; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c index eaf11c32d6f5..0494dfdd0407 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/benchmark/c/benchmark.length.c @@ -86,10 +86,9 @@ static double tic( void ) { * * @return random number */ -static stdlib_complex128_t rand_complex_128( void ) { - double real = (double)rand() / ((double)RAND_MAX + 1.0); - double imag = (double)rand() / ((double)RAND_MAX + 1.0); - return stdlib_complex128(real, imag); +static double rand_double( void ) { + int r = rand(); + return (double)r / ( (double)RAND_MAX + 1.0 ); } /** @@ -106,18 +105,18 @@ static double benchmark1( int iterations, int len ) { int i; for ( i = 0; i < len; i++ ) { - x[ i ] = stdlib_complex128( ( stdlib_complex128_real( rand_complex_128() ) *10000.0 ) - 5000.0, ( stdlib_complex128_imag( rand_complex_128() ) *10000.0 ) - 5000.0 ) ; + x[ i ] = stdlib_complex128( (rand_double()*10000.0)-5000.0, (rand_double()*10000.0)-5000.0 ); } t = tic(); for ( i = 0; i < iterations; i++ ) { c_zdscal( len, 2.0, x, 1 ); - if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ i%len ] ) ) ) { printf( "should not return NaN\n" ); break; } } elapsed = tic() - t; - if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + if ( stdlib_base_is_nan( stdlib_complex128_imag( x[ i%len ] ) ) ) { printf( "should not return NaN\n" ); } return elapsed; @@ -137,18 +136,18 @@ static double benchmark2( int iterations, int len ) { int i; for ( i = 0; i < len; i++ ) { - x[ i ] = stdlib_complex128( ( stdlib_complex128_real( rand_complex_128() ) *10000.0 ) - 5000.0, ( stdlib_complex128_imag( rand_complex_128() ) *10000.0 ) - 5000.0 ) ; + x[ i ] = stdlib_complex128( (rand_double()*10000.0)-5000.0, (rand_double()*10000.0)-5000.0 ); } t = tic(); for ( i = 0; i < iterations; i++ ) { c_zdscal_ndarray( len, 2.0, x, 1, 0 ); - if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ i%len ] ) ) ) { printf( "should not return NaN\n" ); break; } } elapsed = tic() - t; - if ( stdlib_base_is_nan( stdlib_complex128_real( x[ 0 ] ) ) ) { + if ( stdlib_base_is_nan( stdlib_complex128_real( x[ i%len ] ) ) ) { printf( "should not return NaN\n" ); } return elapsed; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c index eb27698a3606..ebd10457f2a3 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/examples/c/example.c @@ -24,7 +24,12 @@ int main( void ) { // Create strided arrays: - stdlib_complex128_t x[] = { stdlib_complex128( 1.0, 2.0 ), stdlib_complex128( 3.0, 4.0 ), stdlib_complex128( 5.0, 6.0 ), stdlib_complex128( 7.0, 8.0 ) }; + stdlib_complex128_t x[] = { + stdlib_complex128( 1.0, 2.0 ), + stdlib_complex128( 3.0, 4.0 ), + stdlib_complex128( 5.0, 6.0 ), + stdlib_complex128( 7.0, 8.0 ) + }; // Specify the number of elements: const int N = 4; diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h index 36a4cd8fa62f..80598f26325c 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal.h @@ -35,12 +35,12 @@ extern "C" { /** * Scales a double-precision complex floating-point vector by a double-precision floating-point constant. */ -void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ); +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ); /** * Scales a double-precision complex floating-point vector by a double-precision floating-point constant using alternative indexing semantics. */ -void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h index 0a7a8779960d..8a7ab3118617 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h +++ b/lib/node_modules/@stdlib/blas/base/zdscal/include/stdlib/blas/base/zdscal_cblas.h @@ -35,7 +35,7 @@ extern "C" { /** * Scales a double-precision complex floating-point vector by a double-precision floating-point constant. */ -void API_SUFFIX(cblas_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ); +void API_SUFFIX(cblas_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json index 4490c55137c8..b49f78266f8a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json +++ b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json @@ -30,7 +30,75 @@ "confs": [ { "task": "build", - "os": "win", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.f", + "./src/zdscal_f.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" + ] + }, + { + "task": "benchmark", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" + ] + }, + { + "task": "examples", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" + ] + }, + + { + "task": "build", + "os": "linux", "blas": "openblas", "wasm": false, "src": [ @@ -45,17 +113,19 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared", "@stdlib/napi/export", "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", "@stdlib/napi/argv-int64", "@stdlib/napi/argv-strided-complex128array", - "@stdlib/napi/argv-double" + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { "task": "benchmark", - "os": "win", + "os": "linux", "blas": "openblas", "wasm": false, "src": [ @@ -70,12 +140,14 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { "task": "examples", - "os": "win", + "os": "linux", "blas": "openblas", "wasm": false, "src": [ @@ -90,7 +162,77 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.f", + "./src/zdscal_f.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex128array", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/zdscal.c", + "./src/zdscal_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/base/scale" ] }, @@ -110,12 +252,14 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared", "@stdlib/napi/export", "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", "@stdlib/napi/argv-int64", "@stdlib/napi/argv-strided-complex128array", - "@stdlib/napi/argv-double" + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { @@ -134,7 +278,9 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { @@ -153,7 +299,9 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, @@ -174,12 +322,14 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared", "@stdlib/napi/export", "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", "@stdlib/napi/argv-int64", "@stdlib/napi/argv-strided-complex128array", - "@stdlib/napi/argv-double" + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { @@ -199,7 +349,9 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { @@ -219,13 +371,15 @@ ], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared" + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float64/ctor" ] }, { "task": "build", - "os": "linux", + "os": "win", "blas": "", "wasm": false, "src": [ @@ -238,12 +392,12 @@ "libraries": [], "libpath": [], "dependencies": [ - "@stdlib/blas/base/shared", "@stdlib/napi/export", "@stdlib/napi/argv", + "@stdlib/blas/base/shared", + "@stdlib/napi/argv-double", "@stdlib/napi/argv-int64", "@stdlib/napi/argv-strided-complex128array", - "@stdlib/napi/argv-double", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/base/scale" @@ -251,7 +405,7 @@ }, { "task": "benchmark", - "os": "linux", + "os": "win", "blas": "", "wasm": false, "src": [ @@ -265,17 +419,14 @@ "libpath": [], "dependencies": [ "@stdlib/blas/base/shared", - "@stdlib/math/base/assert/is-nan", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", - "@stdlib/complex/float64/real", - "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, { "task": "examples", - "os": "linux", + "os": "win", "blas": "", "wasm": false, "src": [ @@ -291,8 +442,6 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", - "@stdlib/complex/float64/real", - "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c index 127887792c0e..d17f7847b585 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/addon.c @@ -38,8 +38,7 @@ static napi_value addon( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); STDLIB_NAPI_ARGV_DOUBLE( env, alpha, argv, 1 ); STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 ); - stdlib_complex128_t *x = (stdlib_complex128_t *)X; - API_SUFFIX(c_zdscal)( N, alpha, x, strideX ); + API_SUFFIX(c_zdscal)( N, alpha, (void *)X, strideX ); return NULL; } @@ -57,8 +56,7 @@ static napi_value addon_method( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); STDLIB_NAPI_ARGV_DOUBLE( env, alpha, argv, 1 ); STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 ); - stdlib_complex128_t *x = (stdlib_complex128_t *)X; - API_SUFFIX(c_zdscal_ndarray)( N, alpha, x, strideX, offsetX ); + API_SUFFIX(c_zdscal_ndarray)( N, alpha, (void *)X, strideX, offsetX ); return NULL; } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c index 673e0a9a16f9..bb11514b96d5 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.c @@ -19,7 +19,6 @@ #include "stdlib/blas/base/zdscal.h" #include "stdlib/strided/base/stride2offset.h" #include "stdlib/blas/base/shared.h" -#include "stdlib/complex/float64/ctor.h" /** * Scales a double-precision complex floating-point vector by a double-precision floating-point constant. @@ -29,7 +28,7 @@ * @param X input array * @param strideX X stride length */ -void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); API_SUFFIX(c_zdscal_ndarray)( N, alpha, X, strideX, ox ); } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c index b18784da9a06..b5eb669a0338 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -30,7 +30,7 @@ * @param X input array * @param strideX X stride length */ -void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { API_SUFFIX(cblas_zdscal)( N, alpha, X, strideX ); } @@ -43,7 +43,7 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex * @param strideX X stride length * @param offsetX starting index for X */ -void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer API_SUFFIX(cblas_zdscal)( N, alpha, x, strideX ); diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index c375974d7ea6..8a71d7aaa946 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -45,7 +45,6 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex */ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; - x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); zdscal( &N, alpha, x, &strideX ); } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c index 90af1f2484bc..7c1a33c44df2 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -30,7 +30,8 @@ * @param strideX X stride length * @param offsetX starting index for X */ -void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + stdlib_complex128_t *x = (stdlib_complex128_t *)X; CBLAS_INT ix; CBLAS_INT sx; CBLAS_INT i; @@ -41,7 +42,7 @@ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib sx = strideX; ix = offsetX; for ( i = 0; i < N; i++ ) { - X[ ix ] = stdlib_base_complex128_scale( alpha, X[ ix ] ); + x[ ix ] = stdlib_base_complex128_scale( alpha, x[ ix ] ); ix += sx; } return; From be3bded914a1dfae6b77949986e8ac93088efab5 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Tue, 27 May 2025 23:51:05 +0530 Subject: [PATCH 06/18] chore: update implementation --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/base/zdscal/manifest.json | 36 +++++++++++++++---- 1 file changed, 30 insertions(+), 6 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json index b49f78266f8a..6397d8de8d98 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json +++ b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json @@ -71,6 +71,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, @@ -92,6 +94,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, @@ -142,7 +146,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, { @@ -164,7 +170,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, @@ -211,6 +219,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, @@ -232,6 +242,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, @@ -280,7 +292,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, { @@ -301,7 +315,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, @@ -351,7 +367,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, { @@ -373,7 +391,9 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", - "@stdlib/complex/float64/ctor" + "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag" ] }, @@ -421,6 +441,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, @@ -442,6 +464,8 @@ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", "@stdlib/complex/float64/ctor", + "@stdlib/complex/float64/real", + "@stdlib/complex/float64/imag", "@stdlib/complex/float64/base/scale" ] }, From 9a513b5d087f92bdf60803a0fdcf45179c30c3b7 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Tue, 27 May 2025 23:53:20 +0530 Subject: [PATCH 07/18] chore: update implementation --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- lib/node_modules/@stdlib/blas/base/zdscal/manifest.json | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json index 6397d8de8d98..44cbd65825c5 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json +++ b/lib/node_modules/@stdlib/blas/base/zdscal/manifest.json @@ -70,6 +70,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag", @@ -146,6 +147,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag" @@ -218,6 +220,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag", @@ -292,6 +295,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag" @@ -367,6 +371,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag" @@ -440,6 +445,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", "@stdlib/complex/float64/ctor", "@stdlib/complex/float64/real", "@stdlib/complex/float64/imag", From d485742395fddccf2b1194b5ea7ac5e3329e9d53 Mon Sep 17 00:00:00 2001 From: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com> Date: Tue, 27 May 2025 23:59:02 +0530 Subject: [PATCH 08/18] chore: add void datatype Co-authored-by: Athan Signed-off-by: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com> --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c index b5eb669a0338..74b5ffb07e26 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -46,5 +46,5 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer - API_SUFFIX(cblas_zdscal)( N, alpha, x, strideX ); + API_SUFFIX(cblas_zdscal)( N, alpha, (void *)x, strideX ); } From f3bcc1c6fb5d13c427be5c861360500cf0891f63 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 00:00:53 +0530 Subject: [PATCH 09/18] chore: add void datatype --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index 8a71d7aaa946..fa1c5fff7dc1 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -30,7 +30,7 @@ * @param X input array * @param strideX X stride length */ -void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX ) { +void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { zdscal( &N, alpha, X, &strideX ); } @@ -43,8 +43,8 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, stdlib_complex * @param strideX X stride length * @param offsetX starting index for X */ -void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, stdlib_complex128_t *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { +void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); - zdscal( &N, alpha, x, &strideX ); + zdscal( &N, alpha, (void *)x, &strideX ); } From 3128670bd0b550ce834a9bb1a3650e6899f4d07b Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 00:16:51 +0530 Subject: [PATCH 10/18] chore: update fortran --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f index e93c1c950ee2..95bfc5eea6e2 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f @@ -57,7 +57,7 @@ subroutine zdscal( N, alpha, x, strideX ) complex(kind=kind(0.0d0)) :: x(*) ! .. ! Local scalars: - integer :: ix, i, nstrideX + integer :: i, nstrideX ! .. ! Intrinsic functions: intrinsic dble, cmplx, aimag @@ -68,12 +68,12 @@ subroutine zdscal( N, alpha, x, strideX ) ! .. if ( strideX == 1 ) then do i = 1, N - x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i))) + x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i)), kind=kind(0.0d0)) end do else nstrideX = N * strideX; do i = 1, nstrideX, strideX - x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i))) + x( i ) = cmplx(alpha*dble(x(i)),alpha*aimag(x(i)), kind=kind(0.0d0)) end do end if return From 04a6ba874cfce1ca87019ffc3bd46dd1525d9033 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 00:21:02 +0530 Subject: [PATCH 11/18] chore: update pointer --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index fa1c5fff7dc1..40527ce6936a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -31,7 +31,7 @@ * @param strideX X stride length */ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { - zdscal( &N, alpha, X, &strideX ); + zdscal( &N, &alpha, X, &strideX ); } /** @@ -46,5 +46,5 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); - zdscal( &N, alpha, (void *)x, &strideX ); + zdscal( &N, &alpha, (void *)x, &strideX ); } From 48d6ebb58815acecae890e627aaab5ea675a4067 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 00:24:06 +0530 Subject: [PATCH 12/18] chore: update precision --- 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 status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - 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: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f index 95bfc5eea6e2..990e222afccd 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal.f @@ -62,7 +62,7 @@ subroutine zdscal( N, alpha, x, strideX ) ! Intrinsic functions: intrinsic dble, cmplx, aimag ! .. - if ( N <= 0 .OR. strideX <= 0 .OR. alpha == 1.0 ) then + if ( N <= 0 .OR. strideX <= 0 .OR. alpha == 1.0d0 ) then return end if ! .. From 6d81a329a71c69d5f3e5cd358dc8165dbb21a6c6 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 10:42:53 +0530 Subject: [PATCH 13/18] chore: update implementation --- .../blas/base/zdscal/src/zdscal_cblas.c | 13 +++++++-- .../blas/base/zdscal/src/zdscal_ndarray.c | 28 +++++++++---------- 2 files changed, 25 insertions(+), 16 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c index 74b5ffb07e26..5aa41e144f57 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -31,7 +31,11 @@ * @param strideX X stride length */ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { - API_SUFFIX(cblas_zdscal)( N, alpha, X, strideX ); + CBLAS_INT sx = strideX; + if ( sx < 0 ) { + sx = -sx; + } + API_SUFFIX(cblas_zdscal)( N, alpha, X, sx ); } /** @@ -45,6 +49,11 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const */ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; + CBLAS_INT sx = strideX; + x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer - API_SUFFIX(cblas_zdscal)( N, alpha, (void *)x, strideX ); + if ( sx < 0 ) { + sx = -sx; + } + API_SUFFIX(cblas_zdscal)( N, alpha, (void *)x, sx ); } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c index 7c1a33c44df2..71746f72325f 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -31,19 +31,19 @@ * @param offsetX starting index for X */ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { - stdlib_complex128_t *x = (stdlib_complex128_t *)X; - CBLAS_INT ix; - CBLAS_INT sx; - CBLAS_INT i; + stdlib_complex128_t z; + int64_t is1; + int64_t i; - if ( N <= 0 || alpha == 1.0 ) { - return; - } - sx = strideX; - ix = offsetX; - for ( i = 0; i < N; i++ ) { - x[ ix ] = stdlib_base_complex128_scale( alpha, x[ ix ] ); - ix += sx; - } - return; + if ( N <= 0 || alpha == 1.0 ) { + return; + } + stdlib_complex128_t *ip1 = (stdlib_complex128_t *)X; + is1 = (int64_t)strideX; + ip1 += offsetX; + for ( i = 0; i < N; i++, ip1 += is1 ) { + z = *ip1; + *ip1 = stdlib_base_complex128_scale( alpha, z ); + } + return; } From 5c0bf7301a3415f6063155c89e63f7660b9ae803 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 11:20:52 +0530 Subject: [PATCH 14/18] chore: update indentation --- .../blas/base/zdscal/lib/ndarray.native.js | 19 ++++++ .../blas/base/zdscal/src/zdscal_cblas.c | 4 +- .../blas/base/zdscal/src/zdscal_ndarray.c | 28 ++++---- .../base/zdscal/test/test.ndarray.native.js | 68 +++++++++---------- 4 files changed, 69 insertions(+), 50 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js index 6046cfe9f480..d0555ec2675a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js @@ -50,7 +50,26 @@ function zdscal( N, alpha, x, strideX, offsetX ) { addon.ndarray( N, alpha, viewX, strideX, offsetX ); return x; } +var Complex128Array = require( '@stdlib/array/complex128' ); +var x = new Complex128Array( [ 0.1, // 3 + 0.1, // 3 + 3.0, + 6.0, + -0.6, // 2 + 0.1, // 2 + 4.0, + 7.0, + 0.1, // 1 + -0.3, // 1 + 7.0, + 2.0 ] ); + +console.log(zdscal(3, 2.0, x, -2, 4)); +for (let i = 0; i < x.length; i++) { + const v = x.get(i); // returns { re, im } + console.log(`x[${i}] = ${v.re} + ${v.im}i`); + } // EXPORTS // diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c index 5aa41e144f57..83dd18f6b6de 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_cblas.c @@ -53,7 +53,7 @@ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void * x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer if ( sx < 0 ) { - sx = -sx; - } + sx = -sx; + } API_SUFFIX(cblas_zdscal)( N, alpha, (void *)x, sx ); } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c index 71746f72325f..7c1a33c44df2 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_ndarray.c @@ -31,19 +31,19 @@ * @param offsetX starting index for X */ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { - stdlib_complex128_t z; - int64_t is1; - int64_t i; + stdlib_complex128_t *x = (stdlib_complex128_t *)X; + CBLAS_INT ix; + CBLAS_INT sx; + CBLAS_INT i; - if ( N <= 0 || alpha == 1.0 ) { - return; - } - stdlib_complex128_t *ip1 = (stdlib_complex128_t *)X; - is1 = (int64_t)strideX; - ip1 += offsetX; - for ( i = 0; i < N; i++, ip1 += is1 ) { - z = *ip1; - *ip1 = stdlib_base_complex128_scale( alpha, z ); - } - return; + if ( N <= 0 || alpha == 1.0 ) { + return; + } + sx = strideX; + ix = offsetX; + for ( i = 0; i < N; i++ ) { + x[ ix ] = stdlib_base_complex128_scale( alpha, x[ ix ] ); + ix += sx; + } + return; } diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js index b4526d5c5670..fb8c25592580 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/test/test.ndarray.native.js @@ -48,11 +48,11 @@ tape( 'the function has an arity of 5', opts, function test( t ) { t.end(); }); -tape( 'the function scales elements from `zx` by `da`', opts, function test( t ) { +tape( 'the function scales elements from `x` by `alpha`', opts, function test( t ) { var expected; - var zx; + var x; - zx = new Complex128Array([ + x = new Complex128Array([ 0.3, // 1 0.1, // 1 0.5, // 2 @@ -67,7 +67,7 @@ tape( 'the function scales elements from `zx` by `da`', opts, function test( t ) 3.0 ]); - zdscal( 4, 2.0, zx, 1, 0 ); + zdscal( 4, 2.0, x, 1, 0 ); expected = new Complex128Array([ 0.6, // 1 @@ -83,15 +83,15 @@ tape( 'the function scales elements from `zx` by `da`', opts, function test( t ) 2.0, 3.0 ]); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); t.end(); }); -tape( 'the function supports a `zx` stride', opts, function test( t ) { +tape( 'the function supports a `x` stride', opts, function test( t ) { var expected; - var zx; + var x; - zx = new Complex128Array([ + x = new Complex128Array([ 0.1, // 1 0.1, // 1 3.0, @@ -106,7 +106,7 @@ tape( 'the function supports a `zx` stride', opts, function test( t ) { 2.0 ]); - zdscal( 3, 2.0, zx, 2, 0 ); + zdscal( 3, 2.0, x, 2, 0 ); expected = new Complex128Array([ 0.2, // 1 @@ -122,15 +122,15 @@ tape( 'the function supports a `zx` stride', opts, function test( t ) { 7.0, 2.0 ]); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); t.end(); }); -tape( 'the function supports a negative `zx` stride', opts, function test( t ) { +tape( 'the function supports a negative `x` stride', opts, function test( t ) { var expected; - var zx; + var x; - zx = new Complex128Array([ + x = new Complex128Array([ 0.1, // 3 0.1, // 3 3.0, @@ -145,7 +145,7 @@ tape( 'the function supports a negative `zx` stride', opts, function test( t ) { 2.0 ]); - zdscal( 3, 2.0, zx, -2, 4 ); + zdscal( 3, 2.0, x, -2, 4 ); expected = new Complex128Array([ 0.2, // 3 @@ -161,15 +161,15 @@ tape( 'the function supports a negative `zx` stride', opts, function test( t ) { 7.0, 2.0 ]); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); t.end(); }); -tape( 'the function supports a `zx` offset', opts, function test( t ) { +tape( 'the function supports a `x` offset', opts, function test( t ) { var expected; - var zx; + var x; - zx = new Complex128Array([ + x = new Complex128Array([ 0.1, 0.1, 3.0, @@ -184,7 +184,7 @@ tape( 'the function supports a `zx` offset', opts, function test( t ) { 2.0 // 3 ]); - zdscal( 3, 2.0, zx, 1, 3 ); + zdscal( 3, 2.0, x, 1, 3 ); expected = new Complex128Array([ 0.1, @@ -200,43 +200,43 @@ tape( 'the function supports a `zx` offset', opts, function test( t ) { 14.0, // 3 4.0 // 3 ]); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.strictEqual( isSameComplex128Array( 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 zx; + var x; - zx = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); - out = zdscal( 4, 2.0, zx, 1, 0 ); + x = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + out = zdscal( 4, 2.0, x, 1, 0 ); - t.strictEqual( out, zx, 'same reference' ); + 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 zx; + var x; - zx = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); + x = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); expected = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 ] ); // eslint-disable-line max-len - zdscal( -1, 2.0, zx, 1, 0 ); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + zdscal( -1, 2.0, x, 1, 0 ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); - zdscal( 0, 2.0, zx, 1, 0 ); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + zdscal( 0, 2.0, x, 1, 0 ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); t.end(); }); tape( 'the function supports complex access patterns', opts, function test( t ) { var expected; - var zx; + var x; - zx = new Complex128Array([ + x = new Complex128Array([ 0.1, 0.1, 3.0, @@ -253,7 +253,7 @@ tape( 'the function supports complex access patterns', opts, function test( t ) 3.0 // 1 ]); - zdscal( 2, 2.0, zx, -3, 6 ); + zdscal( 2, 2.0, x, -3, 6 ); expected = new Complex128Array([ 0.1, @@ -271,6 +271,6 @@ tape( 'the function supports complex access patterns', opts, function test( t ) 4.0, // 1 6.0 // 1 ]); - t.strictEqual( isSameComplex128Array( zx, expected ), true, 'returns expected value' ); + t.strictEqual( isSameComplex128Array( x, expected ), true, 'returns expected value' ); t.end(); }); From 12d01095ea17f2c66c410a62c70f1f7e3ba87124 Mon Sep 17 00:00:00 2001 From: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com> Date: Wed, 28 May 2025 11:25:26 +0530 Subject: [PATCH 15/18] chore: update implementation Signed-off-by: Shabareesh Shetty <139731143+ShabiShett07@users.noreply.github.com> --- .../blas/base/zdscal/lib/ndarray.native.js | 19 ------------------- 1 file changed, 19 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js index d0555ec2675a..6046cfe9f480 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js @@ -50,26 +50,7 @@ function zdscal( N, alpha, x, strideX, offsetX ) { addon.ndarray( N, alpha, viewX, strideX, offsetX ); return x; } -var Complex128Array = require( '@stdlib/array/complex128' ); -var x = new Complex128Array( [ 0.1, // 3 - 0.1, // 3 - 3.0, - 6.0, - -0.6, // 2 - 0.1, // 2 - 4.0, - 7.0, - 0.1, // 1 - -0.3, // 1 - 7.0, - 2.0 ] ); - -console.log(zdscal(3, 2.0, x, -2, 4)); -for (let i = 0; i < x.length; i++) { - const v = x.get(i); // returns { re, im } - console.log(`x[${i}] = ${v.re} + ${v.im}i`); - } // EXPORTS // From e424d87e0bb0f74eefa34a412223f497c2150b54 Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 12:01:55 +0530 Subject: [PATCH 16/18] chore: update fortan wrapper --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index 40527ce6936a..c28d8c1ffcb2 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -31,6 +31,10 @@ * @param strideX X stride length */ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX ) { + CBLAS_INT sx = strideX; + if ( sx < 0 ) { + sx = -sx; + } zdscal( &N, &alpha, X, &strideX ); } @@ -45,6 +49,11 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const */ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { stdlib_complex128_t *x = (stdlib_complex128_t *)X; + CBLAS_INT sx = strideX; + x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); + if ( sx < 0 ) { + sx = -sx; + } zdscal( &N, &alpha, (void *)x, &strideX ); } From 4c6be0336c46b38521568631cb37163d263b190a Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 12:12:19 +0530 Subject: [PATCH 17/18] chore: update fortan wrapper --- lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c index c28d8c1ffcb2..4adb8380f59a 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c +++ b/lib/node_modules/@stdlib/blas/base/zdscal/src/zdscal_f.c @@ -35,7 +35,7 @@ void API_SUFFIX(c_zdscal)( const CBLAS_INT N, const double alpha, void *X, const if ( sx < 0 ) { sx = -sx; } - zdscal( &N, &alpha, X, &strideX ); + zdscal( &N, &alpha, X, &sx ); } /** @@ -55,5 +55,5 @@ void API_SUFFIX(c_zdscal_ndarray)( const CBLAS_INT N, const double alpha, void * if ( sx < 0 ) { sx = -sx; } - zdscal( &N, &alpha, (void *)x, &strideX ); + zdscal( &N, &alpha, (void *)x, &sx ); } From a2a37b8deceb932d2e1ba7c7a077abd0fd7075cb Mon Sep 17 00:00:00 2001 From: ShabiShett07 Date: Wed, 28 May 2025 12:26:22 +0530 Subject: [PATCH 18/18] remove: unused packages --- lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.js | 1 - .../@stdlib/blas/base/zdscal/lib/ndarray.native.js | 1 - lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js | 3 +-- 3 files changed, 1 insertion(+), 4 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.js index be9eab9e337d..986aaf400b9c 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.js @@ -38,7 +38,6 @@ var scale = require( '@stdlib/complex/float64/base/scale' ).strided; * * @example * var Complex128Array = require( '@stdlib/array/complex128' ); -* var Complex128 = require( '@stdlib/complex/float64/ctor' ); * * var zx = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js index 6046cfe9f480..2cb554ed7a1c 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/ndarray.native.js @@ -38,7 +38,6 @@ var addon = require( './../src/addon.node' ); * * @example * var Complex128Array = require( '@stdlib/array/complex128' ); -* var Complex128 = require( '@stdlib/complex/float64/ctor' ); * * var x = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * diff --git a/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js index e36a0e48dd2b..f50bef504da6 100644 --- a/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js +++ b/lib/node_modules/@stdlib/blas/base/zdscal/lib/zdscal.native.js @@ -38,11 +38,10 @@ var addon = require( './../src/addon.node' ); * * @example * var Complex128Array = require( '@stdlib/array/complex128' ); -* var Complex128 = require( '@stdlib/complex/float64/ctor' ); * * var x = new Complex128Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); * -* zdscal( 3, 2.0, x, 1, 0 ); +* zdscal( 3, 2.0, x, 1 ); * // x => [ 2.0, 4.0, 6.0, 8.0, 10.0, 12.0 ] */ function zdscal( N, alpha, x, strideX ) {