@stdlib/blas-base-wasm-idamax 0.1.0

package/lib/module.js

@@ -0,0 +1,198 @@
+/**
+* @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.
+*/
+
+/* 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' );
+*
+* // 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 idamax = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* idamax.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float64';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define a pointer (i.e., byte offset) to the first vector element:
+* var xptr = 0;
+*
+* // Write vector values to module memory:
+* idamax.write( xptr, oneTo( N, dtype ) );
+*
+* // Perform computation:
+* var idx = idamax.main( N, xptr, 1 );
+* // returns 4
+*/
+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 );
+
+/**
+* Finds the index of the first element having the maximum absolute value.
+*
+* @name main
+* @memberof Module.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {NonNegativeInteger} xptr - input array pointer (i.e., byte offset)
+* @param {integer} strideX - `x` stride length
+* @returns {integer} index value
+*
+* @example
+* var Memory = require( '@stdlib/wasm-memory' );
+* var oneTo = require( '@stdlib/array-one-to' );
+*
+* // 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 idamax = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* idamax.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float64';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define a pointer (i.e., byte offset) to the first vector element:
+* var xptr = 0;
+*
+* // Write vector values to module memory:
+* idamax.write( xptr, oneTo( N, dtype ) );
+*
+* // Perform computation:
+* var idx = idamax.main( N, xptr, 1 );
+* // returns 4
+*/
+setReadOnly( Module.prototype, 'main', function idamax( N, xptr, strideX ) {
+	return this._instance.exports.c_idamax( N, xptr, strideX );
+});
+
+/**
+* Finds the index of the first element having the maximum absolute value using alternative indexing semantics.
+*
+* @name ndarray
+* @memberof Module.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {NonNegativeInteger} xptr - input array pointer (i.e., byte offset)
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @returns {integer} index value
+*
+* @example
+* var Memory = require( '@stdlib/wasm-memory' );
+* var oneTo = require( '@stdlib/array-one-to' );
+*
+* // 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 idamax = new Module( mem );
+* // returns <Module>
+*
+* // Initialize the routine:
+* idamax.initializeSync();
+*
+* // Define a vector data type:
+* var dtype = 'float64';
+*
+* // Specify a vector length:
+* var N = 5;
+*
+* // Define a pointer (i.e., byte offset) to the first vector element:
+* var xptr = 0;
+*
+* // Write vector values to module memory:
+* idamax.write( xptr, oneTo( N, dtype ) );
+*
+* // Perform computation:
+* var idx = idamax.ndarray( N, xptr, 1, 0 );
+* // returns 4
+*/
+setReadOnly( Module.prototype, 'ndarray', function idamax( N, xptr, strideX, offsetX ) {
+	return this._instance.exports.c_idamax_ndarray( N, xptr, strideX, offsetX );
+});
+
+
+// EXPORTS //
+
+module.exports = Module;

package/lib/routine.js

@@ -0,0 +1,166 @@
+/**
+* @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.
+*/
+
+/* 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 Float64Array = require( '@stdlib/array-float64' );
+*
+* // Create a new routine:
+* var idamax = new Routine();
+*
+* // Initialize the module:
+* idamax.initializeSync();
+*
+* // Define a strided array:
+* var x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, -5.0 ] );
+*
+* // Perform operation:
+* var idx = idamax.main( x.length, x, 1 );
+* // returns 4
+*
+* @example
+* var Float64Array = require( '@stdlib/array-float64' );
+*
+* // Create a new routine:
+* var idamax = new Routine();
+*
+* // Initialize the module:
+* idamax.initializeSync();
+*
+* // Define a strided array:
+* var x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, -5.0 ] );
+*
+* // Perform operation:
+* var idx = idamax.ndarray( x.length, x, 1, 0 );
+* // returns 4
+*/
+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 );
+
+/**
+* Finds the index of the first element having the maximum absolute value.
+*
+* @name main
+* @memberof Routine.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {Float64Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @returns {integer} index value
+*
+* @example
+* var Float64Array = require( '@stdlib/array-float64' );
+*
+* // Create a new routine:
+* var idamax = new Routine();
+*
+* // Initialize the module:
+* idamax.initializeSync();
+*
+* // Define a strided array:
+* var x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, -5.0 ] );
+*
+* // Perform operation:
+* var idx = idamax.main( x.length, x, 1 );
+* // returns 4
+*/
+setReadOnly( Routine.prototype, 'main', function idamax( N, x, strideX ) {
+	return this.ndarray( N, x, strideX, stride2offset( N, strideX ) );
+});
+
+/**
+* Finds the index of the first element having the maximum absolute value using alternative indexing semantics.
+*
+* @name ndarray
+* @memberof Routine.prototype
+* @readonly
+* @type {Function}
+* @param {PositiveInteger} N - number of indexed elements
+* @param {Float64Array} x - input array
+* @param {integer} strideX - `x` stride length
+* @param {NonNegativeInteger} offsetX - starting index for `x`
+* @returns {integer} index value
+*
+* @example
+* var Float64Array = require( '@stdlib/array-float64' );
+*
+* // Create a new routine:
+* var idamax = new Routine();
+*
+* // Initialize the module:
+* idamax.initializeSync();
+*
+* // Define a strided array:
+* var x = new Float64Array( [ 1.0, 2.0, -3.0, 4.0, -5.0 ] );
+*
+* // Perform operation:
+* var idx = idamax.ndarray( x.length, x, 1, 0 );
+* // returns 4
+*/
+setReadOnly( Routine.prototype, 'ndarray', function idamax( N, x, strideX, offsetX ) {
+	var ptrs;
+	var p0;
+
+	// Convert the input array to "pointers" in the module's memory:
+	ptrs = arrays2ptrs( this, [
+		strided2object( N, x, strideX, offsetX )
+	]);
+	p0 = ptrs[ 0 ];
+
+	// Perform computation by calling the corresponding parent method:
+	return Module.prototype.ndarray.call( this, N, p0.ptr, p0.stride, p0.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-idamax"
+      ]
+    }
+  ]
+}

package/package.json

@@ -0,0 +1,96 @@
+{
+  "name": "@stdlib/blas-base-wasm-idamax",
+  "version": "0.1.0",
+  "description": "Find the index of the first element having the maximum absolute value.",
+  "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-idamax.git"
+  },
+  "bugs": {
+    "url": "https://github.com/stdlib-js/stdlib/issues"
+  },
+  "dependencies": {
+    "@stdlib/assert-is-wasm-memory": "^0.1.0",
+    "@stdlib/blas-base-idamax": "^0.1.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",
+    "idamax",
+    "maximum",
+    "abs",
+    "absolute",
+    "find",
+    "index",
+    "linear",
+    "algebra",
+    "subroutines",
+    "vector",
+    "array",
+    "ndarray",
+    "float64",
+    "double",
+    "float64array",
+    "webassembly",
+    "wasm"
+  ],
+  "funding": {
+    "type": "opencollective",
+    "url": "https://opencollective.com/stdlib"
+  }
+}

package/src/exports.json

@@ -0,0 +1,4 @@
+[
+  "_c_idamax",
+  "_c_idamax_ndarray"
+]

package/src/main.wat

@@ -0,0 +1,107 @@
+;; @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.
+
+(module
+  (type (;0;) (func))
+  (type (;1;) (func (param i32 i32 i32) (result i32)))
+  (type (;2;) (func (param i32 i32 i32 i32) (result i32)))
+  (import "env" "memory" (memory (;0;) 0))
+  (func (;0;) (type 0)
+    nop)
+  (func (;1;) (type 1) (param i32 i32 i32) (result i32)
+    local.get 0
+    local.get 1
+    local.get 2
+    i32.const 1
+    local.get 0
+    i32.sub
+    local.get 2
+    i32.mul
+    i32.const 0
+    local.get 2
+    i32.const 0
+    i32.le_s
+    select
+    call 2)
+  (func (;2;) (type 2) (param i32 i32 i32 i32) (result i32)
+    (local i32 i32 i32 f64 f64)
+    block  ;; label = @1
+      local.get 0
+      i32.const 0
+      i32.le_s
+      if  ;; label = @2
+        i32.const -1
+        local.set 4
+        br 1 (;@1;)
+      end
+      local.get 0
+      i32.const 1
+      i32.eq
+      if  ;; label = @2
+        i32.const 0
+        return
+      end
+      local.get 1
+      local.get 3
+      i32.const 3
+      i32.shl
+      i32.add
+      f64.load
+      f64.abs
+      local.set 7
+      i32.const 1
+      local.set 5
+      loop  ;; label = @2
+        local.get 0
+        local.get 5
+        i32.eq
+        br_if 1 (;@1;)
+        local.get 1
+        local.get 2
+        local.get 3
+        i32.add
+        local.tee 3
+        i32.const 3
+        i32.shl
+        i32.add
+        f64.load
+        f64.abs
+        local.tee 8
+        local.get 7
+        local.get 7
+        local.get 8
+        f64.lt
+        local.tee 6
+        select
+        local.set 7
+        local.get 5
+        local.get 4
+        local.get 6
+        select
+        local.set 4
+        local.get 5
+        i32.const 1
+        i32.add
+        local.set 5
+        br 0 (;@2;)
+      end
+      unreachable
+    end
+    local.get 4)
+  (export "__wasm_call_ctors" (func 0))
+  (export "__wasm_apply_data_relocs" (func 0))
+  (export "c_idamax" (func 1))
+  (export "c_idamax_ndarray" (func 2)))