diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/README.md b/lib/node_modules/@stdlib/math/base/special/cflipsignf/README.md index 48eb8e377adf..20d11c6904a5 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/README.md +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/README.md @@ -131,19 +131,28 @@ for ( i = 0; i < 100; i++ ) { Returns a single-precision complex floating-point number with the same magnitude as `z` and the sign of `y*z`. ```c -#include +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/real.h" +#include "stdlib/complex/float32/imag.h" -float complex y = stdlib_base_cflipsignf( 2.0-1.0*I, -1.0 ); -// returns -2.0+1.0*I +stdlib_complex64_t z = stdlib_complex64( 2.5f, -1.5f ); + +stdlib_complex64_t out = stdlib_base_cflipsignf( z, -1.0f ); + +float re = stdlib_complex64_real( out ); +// returns -2.5f + +float im = stdlib_complex64_imag( out ); +// returns 1.5f ``` The function accepts the following arguments: -- **z**: `[in] float complex` input value. +- **z**: `[in] stdlib_complex64_t` input value. - **y**: `[in] float` number from which to derive the sign. ```c -float complex stdlib_base_cflipsignf( const float complex z, const float y ); +stdlib_complex64_t stdlib_base_cflipsignf( const stdlib_complex64_t z, const float y ); ``` @@ -166,19 +175,31 @@ float complex stdlib_base_cflipsignf( const float complex z, const float y ); ```c #include "stdlib/math/base/special/cflipsignf.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/reim.h" #include -#include int main( void ) { - const float complex x[] = { 3.14f+1.5f*I, -3.14f-1.5f*I, 0.0f+0.0f*I, 0.0f/0.0f+0.0f/0.0f*I }; - - float complex v; - float complex y; + const stdlib_complex64_t x[] = { + stdlib_complex64( 3.14f, 1.5f ), + stdlib_complex64( -3.14f, -1.5f ), + stdlib_complex64( 0.0f, 0.0f ), + stdlib_complex64( 0.0f/0.0f, 0.0f/0.0f ) + }; + + stdlib_complex64_t v; + stdlib_complex64_t y; + float re1; + float im1; + float re2; + float im2; int i; for ( i = 0; i < 4; i++ ) { v = x[ i ]; y = stdlib_base_cflipsignf( v, -1.0f ); - printf( "cflipsignf(%f + %fi, %f) = %f + %fi\n", crealf( v ), cimagf( v ), -1.0f, crealf( y ), cimagf( y ) ); + stdlib_complex64_reim( v, &re1, &im1 ); + stdlib_complex64_reim( y, &re2, &im2 ); + printf( "cflipsignf(%f + %fi, %f) = %f + %fi\n", re1, im1, -1.0f, re2, im2 ); } } ``` diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/benchmark.native.js new file mode 100644 index 000000000000..0acd10248fa9 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/benchmark.native.js @@ -0,0 +1,71 @@ +/** +* @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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var real = require( '@stdlib/complex/float32/real' ); +var imag = require( '@stdlib/complex/float32/imag' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var cflipsignf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( cflipsignf instanceof Error ) +}; + + +// MAIN // + +bench( pkg+'::native', opts, function benchmark( b ) { + var values; + var y; + var i; + var v; + + values = [ + new Complex64( uniform( -500.0, 500.0 ), uniform( -500.0, 500.0 ) ), + new Complex64( uniform( -500.0, 500.0 ), uniform( -500.0, 500.0 ) ) + ]; + v = uniform( 100, -500.0, 500.0, { + 'dtype': 'float32' + }); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = cflipsignf( values[ i%values.length ], v ); + if ( isnanf( real( y ) ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( imag( y ) ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/c/native/benchmark.c b/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/c/native/benchmark.c index bdab52de6d07..865b5a8d99b4 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/c/native/benchmark.c +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/benchmark/c/native/benchmark.c @@ -17,7 +17,8 @@ */ #include "stdlib/math/base/special/cflipsignf.h" -#include +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/reim.h" #include #include #include @@ -91,28 +92,32 @@ static float rand_float( void ) { * @return elapsed time in seconds */ static double benchmark( void ) { - float v[ 100 ]; - float complex x; - float complex y; double elapsed; + float v[ 100 ]; double t; + float re; + float im; int i; + stdlib_complex64_t x[ 100 ]; + stdlib_complex64_t y; + for ( i = 0; i < 100; i++ ) { v[ i ] = ( 1000.0f*rand_float() ) - 500.0f; + x[ i ] = stdlib_complex64( v[ i ], v[ i ] ); } t = tic(); for ( i = 0; i < ITERATIONS; i++ ) { - x = v[ i%100 ] + v[ i%100 ]*I; - y = stdlib_base_cflipsignf( x, -v[ i%100 ] ); - if ( crealf( y ) != crealf( y ) ) { + y = stdlib_base_cflipsignf( x[ i%100 ], -v[ i%100 ] ); + stdlib_complex64_reim( y, &re, &im ); + if ( re != re ) { printf( "unexpected result\n" ); break; } } elapsed = tic() - t; - if ( cimagf( y ) != cimagf( y ) ) { + if ( im != im ) { printf( "unexpected result\n" ); } return elapsed; diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/binding.gyp b/lib/node_modules/@stdlib/math/base/special/cflipsignf/binding.gyp new file mode 100644 index 000000000000..68a1ca11d160 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/binding.gyp @@ -0,0 +1,170 @@ +# @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. + +# A `.gyp` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/examples/c/example.c b/lib/node_modules/@stdlib/math/base/special/cflipsignf/examples/c/example.c index a49f358715eb..cdaaa72da438 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/examples/c/example.c +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/examples/c/example.c @@ -17,18 +17,30 @@ */ #include "stdlib/math/base/special/cflipsignf.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/reim.h" #include -#include int main( void ) { - const float complex x[] = { 3.14f+1.5f*I, -3.14f-1.5f*I, 0.0f+0.0f*I, 0.0f/0.0f+0.0f/0.0f*I }; + const stdlib_complex64_t x[] = { + stdlib_complex64( 3.14f, 1.5f ), + stdlib_complex64( -3.14f, -1.5f ), + stdlib_complex64( 0.0f, 0.0f ), + stdlib_complex64( 0.0f/0.0f, 0.0f/0.0f ) + }; - float complex v; - float complex y; + stdlib_complex64_t v; + stdlib_complex64_t y; + float re1; + float im1; + float re2; + float im2; int i; for ( i = 0; i < 4; i++ ) { v = x[ i ]; y = stdlib_base_cflipsignf( v, -1.0f ); - printf( "cflipsignf(%f + %fi, %f) = %f + %fi\n", crealf( v ), cimagf( v ), -1.0f, crealf( y ), cimagf( y ) ); + stdlib_complex64_reim( v, &re1, &im1 ); + stdlib_complex64_reim( y, &re2, &im2 ); + printf( "cflipsignf(%f + %fi, %f) = %f + %fi\n", re1, im1, -1.0f, re2, im2 ); } } diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/include.gypi b/lib/node_modules/@stdlib/math/base/special/cflipsignf/include.gypi new file mode 100644 index 000000000000..ecfaf82a3279 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/include.gypi @@ -0,0 +1,53 @@ +# @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. + +# A GYP include file for building a Node.js native add-on. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + ' +#include "stdlib/complex/float32/ctor.h" /* * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler. @@ -31,7 +31,7 @@ extern "C" { /** * Returns a single-precision complex floating-point number with the same magnitude as `z` and the sign of `y*z`. */ -float complex stdlib_base_cflipsignf( const float complex z, const float y ); +stdlib_complex64_t stdlib_base_cflipsignf( const stdlib_complex64_t z, const float y ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/lib/native.js b/lib/node_modules/@stdlib/math/base/special/cflipsignf/lib/native.js new file mode 100644 index 000000000000..2919e602e20c --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/lib/native.js @@ -0,0 +1,59 @@ +/** +* @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 Complex64 = require( '@stdlib/complex/float32/ctor' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Returns a single-precision complex floating-point number with the same magnitude as `z` and the sign of `y*z`. +* +* @private +* @param {Complex64} z - complex number +* @param {number} y - number from which to derive the sign +* @returns {Complex64} result +* +* @example +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var real = require( '@stdlib/complex/float32/real' ); +* var imag = require( '@stdlib/complex/float32/imag' ); +* +* var v = cflipsignf( new Complex64( -4.0, 5.5 ), -55.0 ); +* // returns +* +* var re = real( v ); +* // returns 4.0 +* +* var im = imag( v ); +* // returns -5.5 +*/ +function cflipsignf( z, y ) { + var v = addon( z, y ); + return new Complex64( v.re, v.im ); +} + + +// EXPORTS // + +module.exports = cflipsignf; diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/manifest.json b/lib/node_modules/@stdlib/math/base/special/cflipsignf/manifest.json index a0526daf69d2..fe738b2e5f09 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/manifest.json +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/manifest.json @@ -1,40 +1,78 @@ { - "options": {}, - "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": [ - { - "src": [ - "./src/main.c" - ], - "include": [ - "./include" - ], - "libraries": [], - "libpath": [], - "dependencies": [ + "options": { + "task": "build" + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/napi/binary", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/reim", "@stdlib/math/base/assert/is-negative-zerof" ] - } - ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/reim", + "@stdlib/math/base/assert/is-negative-zerof" + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/reim", + "@stdlib/math/base/assert/is-negative-zerof" + ] + } + ] } diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/package.json b/lib/node_modules/@stdlib/math/base/special/cflipsignf/package.json index 076b80e0d49b..bf0b786c8c11 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/package.json +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/package.json @@ -14,6 +14,7 @@ } ], "main": "./lib", + "gypfile": true, "directories": { "benchmark": "./benchmark", "doc": "./docs", diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/Makefile b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/Makefile new file mode 100644 index 000000000000..7733b6180cb4 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @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. +#/ + +# 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/math/base/special/cflipsignf/src/addon.c b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/addon.c new file mode 100644 index 000000000000..73b3ba0b196c --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/addon.c @@ -0,0 +1,22 @@ +/** +* @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. +*/ + +#include "stdlib/math/base/special/cflipsignf.h" +#include "stdlib/math/base/napi/binary.h" + +STDLIB_MATH_BASE_NAPI_MODULE_CF_C( stdlib_base_cflipsignf ) diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/main.c b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/main.c index ef8f618b4985..0b9eb79afa09 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/src/main.c @@ -18,7 +18,8 @@ #include "stdlib/math/base/special/cflipsignf.h" #include "stdlib/math/base/assert/is_negative_zerof.h" -#include +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/reim.h" /** * Returns a single-precision complex floating-point number with the same magnitude as `z` and the sign of `y*z`. @@ -28,15 +29,29 @@ * @return result * * @example -* float complex y = stdlib_base_cflipsignf( 3.0-2.0*I, -1.0 ); -* // returns -3.0+2.0*I +* #include "stdlib/complex/float32/ctor.h" +* #include "stdlib/complex/float32/real.h" +* #include "stdlib/complex/float32/imag.h" +* +* stdlib_complex64_t z = stdlib_complex64( 3.5f, -2.5f ); +* +* stdlib_complex64_t out = stdlib_base_cflipsignf( z, -1.0f ); +* +* float re = stdlib_complex64_real( out ); +* // returns -3.5f +* +* float im = stdlib_complex64_imag( out ); +* // returns 2.5f */ -float complex stdlib_base_cflipsignf( const float complex z, const float y ) { - float re = crealf( z ); - float im = cimagf( z ); +stdlib_complex64_t stdlib_base_cflipsignf( const stdlib_complex64_t z, const float y ) { + float re; + float im; + + stdlib_complex64_reim( z, &re, &im ); + if ( y < 0 || stdlib_base_is_negative_zerof( y ) ) { re = -re; im = -im; } - return re + im*I; + return stdlib_complex64( re, im ); } diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.js b/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.js index 8f3379df368b..ced5ccdcbb07 100644 --- a/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.js +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.js @@ -25,7 +25,7 @@ var PINF = require( '@stdlib/constants/float32/pinf' ); var NINF = require( '@stdlib/constants/float32/ninf' ); var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); var EPS = require( '@stdlib/constants/float32/eps' ); -var abs = require( '@stdlib/math/base/special/abs' ); +var absf = require( '@stdlib/math/base/special/absf' ); var isNegativeZerof = require( '@stdlib/math/base/assert/is-negative-zerof' ); var isPositiveZerof = require( '@stdlib/math/base/assert/is-positive-zerof' ); var Complex64 = require( '@stdlib/complex/float32/ctor' ); @@ -76,12 +76,12 @@ tape( 'the function evaluates the flipsign function', function test( t ) { t.strictEqual( are, ere[ i ], 're: '+re[ i ]+'. Expected: '+ere[ i ] ); t.strictEqual( aim, eim[ i ], 'im: '+im[ i ]+'. Expected: '+eim[ i ] ); } else { - delta = abs( are - ere[ i ] ); - tol = EPS * abs( ere[ i ] ); + delta = absf( are - ere[ i ] ); + tol = EPS * absf( ere[ i ] ); t.ok( delta <= tol, 'within tolerance. re: '+re[ i ]+'. im: '+im[ i ]+'. Actual re: '+are+'. Expected re: '+ere[ i ]+'. delta: '+delta+'. tol: '+tol+'.' ); - delta = abs( aim - eim[ i ] ); - tol = EPS * abs( eim[ i ] ); + delta = absf( aim - eim[ i ] ); + tol = EPS * absf( eim[ i ] ); t.ok( delta <= tol, 'within tolerance. re: '+re[ i ]+'. im: '+im[ i ]+'. Actual im: '+aim+'. Expected im: '+eim[ i ]+'. delta: '+delta+'. tol: '+tol+'.' ); } } diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.native.js b/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.native.js new file mode 100644 index 000000000000..96b5bad56fec --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/cflipsignf/test/test.native.js @@ -0,0 +1,228 @@ +/** +* @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 PINF = require( '@stdlib/constants/float32/pinf' ); +var NINF = require( '@stdlib/constants/float32/ninf' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var EPS = require( '@stdlib/constants/float32/eps' ); +var absf = require( '@stdlib/math/base/special/absf' ); +var isNegativeZerof = require( '@stdlib/math/base/assert/is-negative-zerof' ); +var isPositiveZerof = require( '@stdlib/math/base/assert/is-positive-zerof' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var real = require( '@stdlib/complex/float32/real' ); +var imag = require( '@stdlib/complex/float32/imag' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var cflipsignf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( cflipsignf instanceof Error ) +}; + + +// FIXTURES // + +var data = require( './fixtures/julia/data.json' ); + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof cflipsignf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function evaluates the flipsign function', opts, function test( t ) { + var delta; + var ans; + var ere; + var eim; + var are; + var aim; + var tol; + var re; + var im; + var y; + var z; + var i; + + re = data.re; + im = data.im; + y = data.y; + ere = data.expected_re; + eim = data.expected_im; + + for ( i = 0; i < re.length; i++ ) { + z = new Complex64( re[ i ], im[ i ] ); + ans = cflipsignf( z, y[ i ] ); + are = real( ans ); + aim = imag( ans ); + if ( are === ere[ i ] && aim === eim[ i ] ) { + t.equal( are, ere[ i ], 're: '+re[ i ]+'. Expected: '+ere[ i ] ); + t.equal( aim, eim[ i ], 'im: '+im[ i ]+'. Expected: '+eim[ i ] ); + } else { + delta = absf( are - ere[ i ] ); + tol = EPS * absf( ere[ i ] ); + t.ok( delta <= tol, 'within tolerance. re: '+re[ i ]+'. im: '+im[ i ]+'. Actual re: '+are+'. Expected re: '+ere[ i ]+'. delta: '+delta+'. tol: '+tol+'.' ); + + delta = absf( aim - eim[ i ] ); + tol = EPS * absf( eim[ i ] ); + t.ok( delta <= tol, 'within tolerance. re: '+re[ i ]+'. im: '+im[ i ]+'. Actual im: '+aim+'. Expected im: '+eim[ i ]+'. delta: '+delta+'. tol: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function does not flip the sign of real and imaginary input arguments if `y` is `+infinity`', opts, function test( t ) { + var actual; + var z; + + // Tested against Julia... + z = new Complex64( -4.0, 5.0 ); + actual = cflipsignf( z, PINF ); + + t.strictEqual( real( actual ), -4.0, 'returns expected value' ); + t.strictEqual( imag( actual ), 5.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function flips the sign of real and imaginary input arguments if `y` is `-infinity`', opts, function test( t ) { + var actual; + var z; + + // Tested against Julia... + z = new Complex64( -4.0, 5.0 ); + actual = cflipsignf( z, NINF ); + + t.strictEqual( real( actual ), 4.0, 'returns expected value' ); + t.strictEqual( imag( actual ), -5.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function flips the sign of real and imaginary input arguments if `y` is `-0`', opts, function test( t ) { + var actual; + var z; + + // Tested against Julia... + z = new Complex64( -4.0, 5.0 ); + actual = cflipsignf( z, -0.0 ); + + t.strictEqual( real( actual ), 4.0, 'returns expected value' ); + t.strictEqual( imag( actual ), -5.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function does not flip the sign of real and imaginary input arguments if `y` is `+0`', opts, function test( t ) { + var actual; + var z; + + // Tested against Julia... + z = new Complex64( -4.0, 5.0 ); + actual = cflipsignf( z, +0.0 ); + + t.strictEqual( real( actual ), -4.0, 'returns expected value' ); + t.strictEqual( imag( actual ), 5.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a `NaN` if provided a `NaN` independent of `y`', opts, function test( t ) { + var v; + + v = cflipsignf( new Complex64( NaN, NaN ), 2.0 ); + t.strictEqual( isnanf( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isnanf( imag( v ) ), true, 'returns expected value' ); + + v = cflipsignf( new Complex64( NaN, NaN ), -2.0 ); + t.strictEqual( isnanf( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isnanf( imag( v ) ), true, 'returns expected value' ); + + v = cflipsignf( new Complex64( NaN, NaN ), NaN ); + t.strictEqual( isnanf( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isnanf( imag( v ) ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns `+0` if provided `+0` and a positive `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( +0.0, +0.0 ), 4.0 ); + t.strictEqual( isPositiveZerof( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isPositiveZerof( imag( v ) ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `-0` if provided `-0` and a positive `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( -0.0, -0.0 ), 4.0 ); + t.strictEqual( isNegativeZerof( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isNegativeZerof( imag( v ) ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `+0` if provided `-0` and a negative `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( -0.0, -0.0 ), -4.0 ); + t.strictEqual( isPositiveZerof( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isPositiveZerof( imag( v ) ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `-0` if provided `+0` and a negative `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( +0.0, +0.0 ), -4.0 ); + t.strictEqual( isNegativeZerof( real( v ) ), true, 'returns expected value' ); + t.strictEqual( isNegativeZerof( imag( v ) ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `+infinity` if provided `+infinity` and a positive `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( PINF, PINF ), 4.0 ); + t.strictEqual( real( v ), PINF, 'returns expected value' ); + t.strictEqual( imag( v ), PINF, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `-infinity` if provided `-infinity` and a positive `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( NINF, NINF ), 4.0 ); + t.strictEqual( real( v ), NINF, 'returns expected value' ); + t.strictEqual( imag( v ), NINF, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `+infinity` if provided `-infinity` and a negative `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( NINF, NINF ), -4.0 ); + t.strictEqual( real( v ), PINF, 'returns expected value' ); + t.strictEqual( imag( v ), PINF, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns `-infinity` if provided `+infinity` and a negative `y`', opts, function test( t ) { + var v = cflipsignf( new Complex64( PINF, PINF ), -4.0 ); + t.strictEqual( real( v ), NINF, 'returns expected value' ); + t.strictEqual( imag( v ), NINF, 'returns expected value' ); + t.end(); +});