@stdlib/blas-base-wasm-sdsdot 0.1.0

package/lib/module.js

@@ -0,0 +1,220 @@
+/**
+* @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.
+*/
+
+/* eslint-disable no-restricted-syntax, no-invalid-this */
+
+'use strict';
+
+// MODULES //
+
+var isWebAssemblyMemory = require( '@stdlib/assert-is-wasm-memory' );
+var setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );
+var inherits = require( '@stdlib/utils-inherit' );
+var WasmModule = require( '@stdlib/wasm-module-wrapper' );
+var format = require( '@stdlib/string-format' );
+var wasmBinary = require( './binary.js' );
+
+
+// MAIN //
+
+/**
+* BLAS routine WebAssembly module wrapper constructor.
+*
+* @constructor
+* @param {Object} memory - WebAssembly memory instance
+* @throws {TypeError} must provide a WebAssembly memory instance
+* @returns {Module} module instance
+*
+* @example
+* var Memory = require( '@stdlib/wasm-memory' );
+* var oneTo = require( '@stdlib/array-one-to' );
+* var ones = require( '@stdlib/array-ones' );
+* var zeros = require( '@stdlib/array-zeros' );
+* var bytesPerElement = require( '@stdlib/ndarray-base-bytes-per-element' );
+*
+* // Create a new memory instance with an initial size of 10 pages (640KiB) and a maximum size of 100 pages (6.4MiB):
+* var mem = new Memory({
+*     'initial': 10,
+*     'maximum': 100
+* });
+*
+* // Create a BLAS routine:
+* var sdsdot = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* sdsdot.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float32';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define pointers (i.e., byte offsets) for storing two vectors:
+* var xptr = 0;
+* var yptr = N * bytesPerElement( dtype );
+*
+* // Write vector values to module memory:
+* sdsdot.write( xptr, oneTo( N, dtype ) );
+* sdsdot.write( yptr, ones( N, dtype ) );
+*
+* // Perform computation:
+* var dot = sdsdot.main( N, 0.0, xptr, 1, yptr, 1 );
+* // returns 15.0
+*/
+function Module( memory ) {
+	if ( !( this instanceof Module ) ) {
+		return new Module( memory );
+	}
+	if ( !isWebAssemblyMemory( memory ) ) {
+		throw new TypeError( format( 'invalid argument. Must provide a WebAssembly memory instance. Value: `%s`.', memory ) );
+	}
+	// Call the parent constructor:
+	WasmModule.call( this, wasmBinary, memory, {
+		'env': {
+			'memory': memory
+		}
+	});
+
+	return this;
+}
+
+// Inherit from the parent constructor:
+inherits( Module, WasmModule );
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation.
+*
+* @name main
+* @memberof Module.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} scalar - scalar constant to add to dot product
+* @param {NonNegativeInteger} xptr - first input array pointer (i.e., byte offset)
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} yptr - second input array pointer (i.e., byte offset)
+* @param {integer} strideY - `y` stride length
+* @returns {number} dot product
+*
+* @example
+* var Memory = require( '@stdlib/wasm-memory' );
+* var oneTo = require( '@stdlib/array-one-to' );
+* var ones = require( '@stdlib/array-ones' );
+* var zeros = require( '@stdlib/array-zeros' );
+* var bytesPerElement = require( '@stdlib/ndarray-base-bytes-per-element' );
+*
+* // Create a new memory instance with an initial size of 10 pages (640KiB) and a maximum size of 100 pages (6.4MiB):
+* var mem = new Memory({
+*     'initial': 10,
+*     'maximum': 100
+* });
+*
+* // Create a BLAS routine:
+* var sdsdot = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* sdsdot.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float32';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define pointers (i.e., byte offsets) for storing two vectors:
+* var xptr = 0;
+* var yptr = N * bytesPerElement( dtype );
+*
+* // Write vector values to module memory:
+* sdsdot.write( xptr, oneTo( N, dtype ) );
+* sdsdot.write( yptr, ones( N, dtype ) );
+*
+* // Perform computation:
+* var dot = sdsdot.main( N, 0.0, xptr, 1, yptr, 1 );
+* // returns 15.0
+*/
+setReadOnly( Module.prototype, 'main', function sdsdot( N, scalar, xptr, strideX, yptr, strideY ) {
+	return this._instance.exports.c_sdsdot( N, scalar, xptr, strideX, yptr, strideY );
+});
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*
+* @name ndarray
+* @memberof Module.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} scalar - scalar constant to add to dot product
+* @param {NonNegativeInteger} xptr - first input array pointer (i.e., byte offset)
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting `x` index
+* @param {NonNegativeInteger} yptr - second input array pointer (i.e., byte offset)
+* @param {integer} strideY - `y` stride length
+* @param {NonNegativeInteger} offsetY - starting `y` index
+* @returns {number} dot product
+*
+* @example
+* var Memory = require( '@stdlib/wasm-memory' );
+* var oneTo = require( '@stdlib/array-one-to' );
+* var ones = require( '@stdlib/array-ones' );
+* var zeros = require( '@stdlib/array-zeros' );
+* var bytesPerElement = require( '@stdlib/ndarray-base-bytes-per-element' );
+*
+* // Create a new memory instance with an initial size of 10 pages (640KiB) and a maximum size of 100 pages (6.4MiB):
+* var mem = new Memory({
+*     'initial': 10,
+*     'maximum': 100
+* });
+*
+* // Create a BLAS routine:
+* var sdsdot = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* sdsdot.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float32';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define pointers (i.e., byte offsets) for storing two vectors:
+* var xptr = 0;
+* var yptr = N * bytesPerElement( dtype );
+*
+* // Write vector values to module memory:
+* sdsdot.write( xptr, oneTo( N, dtype ) );
+* sdsdot.write( yptr, ones( N, dtype ) );
+*
+* // Perform computation:
+* var sdsdot = sdsdot.ndarray( N, 0.0, xptr, 1, 0, yptr, 1, 0 );
+* // returns 15.0
+*/
+setReadOnly( Module.prototype, 'ndarray', function sdsdot( N, scalar, xptr, strideX, offsetX, yptr, strideY, offsetY ) {
+	return this._instance.exports.c_sdsdot_ndarray( N, scalar, xptr, strideX, offsetX, yptr, strideY, offsetY ); // eslint-disable-line max-len
+});
+
+
+// EXPORTS //
+
+module.exports = Module;

package/lib/routine.js

@@ -0,0 +1,180 @@
+/**
+* @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.
+*/
+
+/* eslint-disable no-restricted-syntax, no-invalid-this */
+
+'use strict';
+
+// MODULES //
+
+var setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );
+var inherits = require( '@stdlib/utils-inherit' );
+var stride2offset = require( '@stdlib/strided-base-stride2offset' );
+var Memory = require( '@stdlib/wasm-memory' );
+var arrays2ptrs = require( '@stdlib/wasm-base-arrays2ptrs' );
+var strided2object = require( '@stdlib/wasm-base-strided2object' );
+var Module = require( './module.js' );
+
+
+// MAIN //
+
+/**
+* Routine constructor.
+*
+* @private
+* @constructor
+* @returns {Routine} routine instance
+*
+* @example
+* var Float32Array = require( '@stdlib/array-float32' );
+*
+* // Create a new routine:
+* var sdsdot = new Routine();
+*
+* // Initialize the module:
+* sdsdot.initializeSync();
+*
+* // Define strided arrays:
+* var x = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] );
+* var y = new Float32Array( [ 1.0, 1.0, 1.0, 1.0, 1.0 ] );
+*
+* // Perform operation:
+* var dot = sdsdot.main( x.length, 0.0, x, 1, y, 1 );
+* // returns 15.0
+*
+* @example
+* var Float32Array = require( '@stdlib/array-float32' );
+*
+* // Create a new routine:
+* var sdsdot = new Routine();
+*
+* // Initialize the module:
+* sdsdot.initializeSync();
+*
+* // Define strided arrays:
+* var x = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] );
+* var y = new Float32Array( [ 1.0, 1.0, 1.0, 1.0, 1.0 ] );
+*
+* // Perform operation:
+* var dot = sdsdot.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 );
+* // returns 15.0
+*/
+function Routine() {
+	if ( !( this instanceof Routine ) ) {
+		return new Routine();
+	}
+	Module.call( this, new Memory({
+		'initial': 0
+	}));
+	return this;
+}
+
+// Inherit from the parent constructor:
+inherits( Routine, Module );
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation.
+*
+* @name main
+* @memberof Routine.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} scalar - scalar constant to add to dot product
+* @param {Float32Array} x - first input array
+* @param {integer} strideX - `x` stride length
+* @param {Float32Array} y - second input array
+* @param {integer} strideY - `y` stride length
+* @returns {number} dot product
+*
+* @example
+* var Float32Array = require( '@stdlib/array-float32' );
+*
+* // Create a new routine:
+* var sdsdot = new Routine();
+*
+* // Initialize the module:
+* sdsdot.initializeSync();
+*
+* // Define strided arrays:
+* var x = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] );
+* var y = new Float32Array( [ 1.0, 1.0, 1.0, 1.0, 1.0 ] );
+*
+* // Perform operation:
+* var dot = sdsdot.main( x.length, 0.0, x, 1, y, 1 );
+* // returns 15.0
+*/
+setReadOnly( Routine.prototype, 'main', function sdsdot( N, scalar, x, strideX, y, strideY ) {
+	return this.ndarray( N, scalar, x, strideX, stride2offset( N, strideX ), y, strideY, stride2offset( N, strideY ) ); // eslint-disable-line max-len
+});
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*
+* @name ndarray
+* @memberof Routine.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {number} scalar - scalar constant to add to dot product
+* @param {Float32Array} x - first input array
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting `x` index
+* @param {Float32Array} y - second input array
+* @param {integer} strideY - `y` stride length
+* @param {NonNegativeInteger} offsetY - starting `y` index
+* @returns {number} dot product
+*
+* @example
+* var Float32Array = require( '@stdlib/array-float32' );
+*
+* // Create a new routine:
+* var sdsdot = new Routine();
+*
+* // Initialize the module:
+* sdsdot.initializeSync();
+*
+* // Define strided arrays:
+* var x = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] );
+* var y = new Float32Array( [ 1.0, 1.0, 1.0, 1.0, 1.0 ] );
+*
+* // Perform operation:
+* var dot = sdsdot.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 );
+* // returns 15.0
+*/
+setReadOnly( Routine.prototype, 'ndarray', function sdsdot( N, scalar, x, strideX, offsetX, y, strideY, offsetY ) {
+	var ptrs;
+	var p0;
+	var p1;
+
+	// Convert the input arrays to "pointers" in the module's memory:
+	ptrs = arrays2ptrs( this, [
+		strided2object( N, x, strideX, offsetX ),
+		strided2object( N, y, strideY, offsetY )
+	]);
+	p0 = ptrs[ 0 ];
+	p1 = ptrs[ 1 ];
+
+	// Perform computation by calling the corresponding parent method:
+	return Module.prototype.ndarray.call( this, N, scalar, p0.ptr, p0.stride, p0.offset, p1.ptr, p1.stride, p1.offset ); // eslint-disable-line max-len
+});
+
+
+// EXPORTS //
+
+module.exports = Routine;

package/manifest.json

@@ -0,0 +1,36 @@
+{
+  "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": [],
+      "include": [],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/blas-base-sdsdot"
+      ]
+    }
+  ]
+}

package/package.json

@@ -0,0 +1,92 @@
+{
+  "name": "@stdlib/blas-base-wasm-sdsdot",
+  "version": "0.1.0",
+  "description": "Calculate the dot product of two single-precision floating-point vectors with extended accumulation.",
+  "license": "Apache-2.0",
+  "author": {
+    "name": "The Stdlib Authors",
+    "url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
+  },
+  "contributors": [
+    {
+      "name": "The Stdlib Authors",
+      "url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
+    }
+  ],
+  "main": "./lib",
+  "browser": {
+    "./lib/binary.js": "./lib/binary.browser.js"
+  },
+  "directories": {
+    "doc": "./docs",
+    "lib": "./lib",
+    "scripts": "./scripts",
+    "src": "./src",
+    "dist": "./dist"
+  },
+  "types": "./docs/types",
+  "scripts": {},
+  "homepage": "https://stdlib.io",
+  "repository": {
+    "type": "git",
+    "url": "git://github.com/stdlib-js/blas-base-wasm-sdsdot.git"
+  },
+  "bugs": {
+    "url": "https://github.com/stdlib-js/stdlib/issues"
+  },
+  "dependencies": {
+    "@stdlib/assert-is-wasm-memory": "^0.1.0",
+    "@stdlib/blas-base-sdsdot": "^0.3.0",
+    "@stdlib/fs-read-wasm": "^0.2.2",
+    "@stdlib/strided-base-stride2offset": "^0.1.0",
+    "@stdlib/string-format": "^0.2.2",
+    "@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.2",
+    "@stdlib/utils-inherit": "^0.2.2",
+    "@stdlib/wasm-base-arrays2ptrs": "github:stdlib-js/wasm-base-arrays2ptrs#main",
+    "@stdlib/wasm-base-strided2object": "github:stdlib-js/wasm-base-strided2object#main",
+    "@stdlib/wasm-memory": "^0.1.0",
+    "@stdlib/wasm-module-wrapper": "github:stdlib-js/wasm-module-wrapper#main",
+    "@stdlib/error-tools-fmtprodmsg": "^0.2.2"
+  },
+  "devDependencies": {},
+  "engines": {
+    "node": ">=0.10.0",
+    "npm": ">2.7.0"
+  },
+  "os": [
+    "aix",
+    "darwin",
+    "freebsd",
+    "linux",
+    "macos",
+    "openbsd",
+    "sunos",
+    "win32",
+    "windows"
+  ],
+  "keywords": [
+    "stdlib",
+    "stdmath",
+    "mathematics",
+    "math",
+    "blas",
+    "level 1",
+    "sdsdot",
+    "linear",
+    "algebra",
+    "subroutines",
+    "dot",
+    "vector",
+    "array",
+    "ndarray",
+    "float32",
+    "float",
+    "float32array",
+    "webassembly",
+    "wasm"
+  ],
+  "funding": {
+    "type": "opencollective",
+    "url": "https://opencollective.com/stdlib"
+  }
+}

package/src/exports.json

@@ -0,0 +1,4 @@
+[
+  "_c_sdsdot",
+  "_c_sdsdot_ndarray"
+]