UNPKG

@stdlib/math-strided-special-dmskrsqrt

Version:

Compute the reciprocal square root for each element in a double-precision floating-point strided array according to a strided mask array.

63 lines (53 loc) 2.11 kB
/** * @license Apache-2.0 * * Copyright (c) 2021 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 dmskmap = require( '@stdlib/strided-base-dmskmap' ).ndarray; var rsqrt = require( '@stdlib/math-base-special-rsqrt' ); // MAIN // /** * Computes the reciprocal square root for each element in a double-precision floating-point strided array `x` according to a strided mask array and assigns the results to elements in a double-precision floating-point strided array `y`. * * @param {NonNegativeInteger} N - number of indexed elements * @param {Float64Array} x - input array * @param {integer} sx - `x` stride length * @param {NonNegativeInteger} ox - starting `x` index * @param {Uint8Array} m - mask array * @param {integer} sm - `m` stride length * @param {NonNegativeInteger} om - starting `m` index * @param {Float64Array} y - destination array * @param {integer} sy - `y` stride length * @param {NonNegativeInteger} oy - starting `y` index * @returns {Float64Array} `y` * * @example * var Float64Array = require( '@stdlib/array-float64' ); * var Uint8Array = require( '@stdlib/array-uint8' ); * * var x = new Float64Array( [ 0.0, 4.0, 9.0, 12.0, 24.0 ] ); * var m = new Uint8Array( [ 0, 0, 1, 0, 1 ] ); * var y = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0 ] ); * * dmskrsqrt( x.length, x, 1, 0, m, 1, 0, y, 1, 0 ); * // y => <Float64Array>[ Infinity, 0.5, 0.0, ~0.289, 0.0 ] */ function dmskrsqrt( N, x, sx, ox, m, sm, om, y, sy, oy ) { return dmskmap( N, x, sx, ox, m, sm, om, y, sy, oy, rsqrt ); } // EXPORTS // module.exports = dmskrsqrt;