Press n or j to go to the next uncovered block, b, p or k for the previous block.
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 | 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 4x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 19543x 21x 21x 21x 21x 19522x 19543x 19522x 19522x 19522x 19522x 19522x 19522x 19522x 19522x 19522x 19543x 4x 4x 4x 4x 4x | /**
* @license Apache-2.0
*
* Copyright (c) 2026 The Stdlib Authors.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
'use strict';
// MODULES //
var isnan = require( '@stdlib/math/base/assert/is-nan' );
var isInfinite = require( '@stdlib/math/base/assert/is-infinite' );
var f32 = require( '@stdlib/number/float64/base/to-float32' );
// MAIN //
/**
* Divides two single-precision complex floating-point numbers and assigns results to a provided output array.
*
* @param {number} re1 - real component of the first complex number
* @param {number} im1 - imaginary component of the first complex number
* @param {number} re2 - real component of the second complex number
* @param {number} im2 - imaginary component of the second complex number
* @param {Collection} out - output array
* @param {integer} strideOut - stride length
* @param {NonNegativeInteger} offsetOut - starting index
* @returns {Collection} output array
*
* @example
* var Float32Array = require( '@stdlib/array/float32' );
*
* var out = assign( -13.0, -1.0, -2.0, 1.0, new Float32Array( 2 ), 1, 0 );
* // returns <Float32Array>[ 5.0, 3.0 ]
*/
function assign( re1, im1, re2, im2, out, strideOut, offsetOut ) {
var rePart;
var imPart;
var mag;
var a;
var b;
var c;
var d;
// Widen to float64 for intermediate calculations (Julia's approach for Float32)
a = re1;
b = im1;
c = re2;
d = im2;
if ( isnan( a ) || isnan( b ) || isnan( c ) || isnan( d ) ) {
out[ offsetOut ] = NaN;
out[ offsetOut + strideOut ] = NaN;
return out;
}
if ( isInfinite( c ) || isInfinite( d ) ) {
if ( isInfinite( a ) || isInfinite( b ) ) {
out[ offsetOut ] = NaN;
out[ offsetOut + strideOut ] = NaN;
return out;
}
out[ offsetOut ] = 0.0;
out[ offsetOut + strideOut ] = 0.0;
return out;
}
mag = 1.0 / ( ( c * c ) + ( d * d ) );
rePart = ( ( a * c ) + ( b * d ) ) * mag;
imPart = ( ( b * c ) - ( a * d ) ) * mag;
out[ offsetOut ] = f32( rePart );
out[ offsetOut + strideOut ] = f32( imPart );
return out;
}
// EXPORTS //
module.exports = assign;
|