-
-
Notifications
You must be signed in to change notification settings - Fork 862
feat: add ndarray/base/nullary-strided1d
#7772
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: develop
Are you sure you want to change the base?
Changes from all commits
4bc9fed
070159e
1fcc2fd
d2680d3
4c5c84a
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
Original file line number | Diff line number | Diff line change | ||||
---|---|---|---|---|---|---|
@@ -0,0 +1,214 @@ | ||||||
<!-- | ||||||
|
||||||
@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. | ||||||
|
||||||
--> | ||||||
|
||||||
# nullaryStrided1d | ||||||
|
||||||
> Apply a one-dimensional strided array function to a list of specified dimensions in an input ndarray. | ||||||
|
||||||
<section class="intro"> | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.intro --> | ||||||
|
||||||
<section class="usage"> | ||||||
|
||||||
## Usage | ||||||
|
||||||
```javascript | ||||||
var nullaryStrided1d = require( '@stdlib/ndarray/base/nullary-strided1d' ); | ||||||
``` | ||||||
|
||||||
#### nullaryStrided1d( fcn, arrays, dims\[, options] ) | ||||||
|
||||||
Applies a one-dimensional strided array function to a list of specified dimensions in an input ndarray. | ||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||
|
||||||
<!-- eslint-disable max-len --> | ||||||
|
||||||
```javascript | ||||||
var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); | ||||||
var getStride = require( '@stdlib/ndarray/base/stride' ); | ||||||
var getOffset = require( '@stdlib/ndarray/base/offset' ); | ||||||
var getData = require( '@stdlib/ndarray/base/data-buffer' ); | ||||||
var numelDimension = require( '@stdlib/ndarray/base/numel-dimension' ); | ||||||
var ndarraylike2scalar = require( '@stdlib/ndarray/base/ndarraylike2scalar' ); | ||||||
var gsorthp = require( '@stdlib/blas/ext/base/gsorthp' ).ndarray; | ||||||
|
||||||
function wrapper( arrays ) { | ||||||
var x = arrays[ 0 ]; | ||||||
var o = arrays[ 1 ]; | ||||||
return gsorthp( numelDimension( x, 0 ), ndarraylike2scalar( o ), getData( x ), getStride( x, 0 ), getOffset( x ) ); | ||||||
} | ||||||
|
||||||
// Create data buffers: | ||||||
var xbuf = [ 12.0, 11.0, 10.0, 9.0, 8.0, 7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0 ]; | ||||||
|
||||||
// Define the array shapes: | ||||||
var xsh = [ 1, 3, 2, 2 ]; | ||||||
|
||||||
// Define the array strides: | ||||||
var sx = [ 12, 4, 2, 1 ]; | ||||||
|
||||||
// Define the index offsets: | ||||||
var ox = 0; | ||||||
|
||||||
// Create an input ndarray-like object: | ||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||
var x = { | ||||||
'dtype': 'generic', | ||||||
'data': xbuf, | ||||||
'shape': xsh, | ||||||
'strides': sx, | ||||||
'offset': ox, | ||||||
'order': 'row-major' | ||||||
}; | ||||||
|
||||||
// Create an ndarray-like object for the sort order: | ||||||
var sortOrder = { | ||||||
'dtype': 'generic', | ||||||
'data': [ 1.0 ], | ||||||
'shape': [ 1, 3 ], | ||||||
'strides': [ 0, 0 ], | ||||||
'offset': 0, | ||||||
'order': 'row-major' | ||||||
}; | ||||||
|
||||||
// Apply strided function: | ||||||
nullaryStrided1d( wrapper, [ x, sortOrder ], [ 2, 3 ] ); | ||||||
|
||||||
var arr = ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ); | ||||||
// returns [ [ [ [ 9.0, 10.0 ], [ 11.0, 12.0 ] ], [ [ 5.0, 6.0 ], [ 7.0, 8.0 ] ], [ [ 1.0, 2.0 ], [ 3.0, 4.0 ] ] ] ] | ||||||
``` | ||||||
|
||||||
The function accepts the following arguments: | ||||||
|
||||||
- **fcn**: function which will be applied to a one-dimensional input subarray. | ||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||
- **arrays**: array-like object containing one input ndarray followed by any additional ndarray arguments. | ||||||
- **dims**: list of dimensions to which to apply a strided array function. | ||||||
- **options**: function options which are passed through to `fcn` (_optional_). | ||||||
|
||||||
Each provided ndarray should be an object with the following properties: | ||||||
|
||||||
- **dtype**: data type. | ||||||
- **data**: data buffer. | ||||||
- **shape**: dimensions. | ||||||
- **strides**: stride lengths. | ||||||
- **offset**: index offset. | ||||||
- **order**: specifies whether an ndarray is row-major (C-style) or column major (Fortran-style). | ||||||
|
||||||
#### TODO: document factory method | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.usage --> | ||||||
|
||||||
<section class="notes"> | ||||||
|
||||||
## Notes | ||||||
|
||||||
- Any additional ndarray arguments are expected to have the same dimensions as the loop dimensions of the input ndarray. When calling the strided array function, any additional ndarray arguments are provided as zero-dimensional ndarray-like objects. | ||||||
|
||||||
- The strided array function is expected to have the following signature: | ||||||
|
||||||
```text | ||||||
fcn( arrays[, options] ) | ||||||
``` | ||||||
|
||||||
where | ||||||
|
||||||
- **arrays**: array containing a one-dimensional subarray of the input ndarray and any additional ndarray arguments as zero-dimensional ndarrays. | ||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||
- **options**: function options (_optional_). | ||||||
|
||||||
- The function iterates over ndarray elements according to the memory layout of the input ndarray. | ||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||
|
||||||
- For very high-dimensional ndarrays which are non-contiguous, one should consider copying the underlying data to contiguous memory before performing an operation in order to achieve better performance. | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.notes --> | ||||||
|
||||||
<section class="examples"> | ||||||
|
||||||
## Examples | ||||||
|
||||||
<!-- eslint-disable max-len --> | ||||||
|
||||||
<!-- eslint no-undef: "error" --> | ||||||
|
||||||
```javascript | ||||||
var discreteUniform = require( '@stdlib/random/array/discrete-uniform' ); | ||||||
var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); | ||||||
var numelDimension = require( '@stdlib/ndarray/base/numel-dimension' ); | ||||||
var getData = require( '@stdlib/ndarray/base/data-buffer' ); | ||||||
var getStride = require( '@stdlib/ndarray/base/stride' ); | ||||||
var getOffset = require( '@stdlib/ndarray/base/offset' ); | ||||||
var ndarraylike2scalar = require( '@stdlib/ndarray/base/ndarraylike2scalar' ); | ||||||
var gsorthp = require( '@stdlib/blas/ext/base/gsorthp' ).ndarray; | ||||||
var nullaryStrided1d = require( '@stdlib/ndarray/base/nullary-strided1d' ); | ||||||
|
||||||
function wrapper( arrays ) { | ||||||
var x = arrays[ 0 ]; | ||||||
var o = arrays[ 1 ]; | ||||||
return gsorthp( numelDimension( x, 0 ), ndarraylike2scalar( o ), getData( x ), getStride( x, 0 ), getOffset( x ) ); | ||||||
} | ||||||
|
||||||
var N = 10; | ||||||
var x = { | ||||||
'dtype': 'generic', | ||||||
'data': discreteUniform( N, -5, 5, { | ||||||
'dtype': 'generic' | ||||||
}), | ||||||
'shape': [ 1, 5, 2 ], | ||||||
'strides': [ 10, 2, 1 ], | ||||||
'offset': 0, | ||||||
'order': 'row-major' | ||||||
}; | ||||||
console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); | ||||||
|
||||||
var sortOrder = { | ||||||
'dtype': 'generic', | ||||||
'data': [ 1.0 ], | ||||||
'shape': [ 2 ], | ||||||
'strides': [ 0, 0 ], | ||||||
'offset': 0, | ||||||
'order': 'row-major' | ||||||
}; | ||||||
|
||||||
nullaryStrided1d( wrapper, [ x, sortOrder ], [ 0, 1 ] ); | ||||||
|
||||||
console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); | ||||||
``` | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.examples --> | ||||||
|
||||||
<!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. --> | ||||||
|
||||||
<section class="related"> | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.related --> | ||||||
|
||||||
<section class="links"> | ||||||
|
||||||
</section> | ||||||
|
||||||
<!-- /.links --> |
Original file line number | Diff line number | Diff line change | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
@@ -0,0 +1,92 @@ | ||||||||||
|
||||||||||
{{alias}}( fcn, arrays, dims[, options] ) | ||||||||||
Applies a one-dimensional strided array function to a list of specified | ||||||||||
dimensions in an input ndarray. | ||||||||||
Comment on lines
+3
to
+4
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||||||
|
||||||||||
Each provided "ndarray" should be an object with the following properties: | ||||||||||
|
||||||||||
- dtype: data type. | ||||||||||
- data: data buffer. | ||||||||||
- shape: dimensions. | ||||||||||
- strides: stride lengths. | ||||||||||
- offset: index offset. | ||||||||||
- order: specifies whether an ndarray is row-major (C-style) or column-major | ||||||||||
(Fortran-style). | ||||||||||
|
||||||||||
Any additional ndarray arguments are expected to have the same dimensions as | ||||||||||
the loop dimensions of the input ndarray. When calling the strided array | ||||||||||
function, any additional ndarray arguments are provided as zero-dimensional | ||||||||||
ndarray-like objects. | ||||||||||
|
||||||||||
Parameters | ||||||||||
---------- | ||||||||||
fcn: Function | ||||||||||
Function which will be applied to a one-dimensional input subarray. The | ||||||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||||||
function should have the following signature: | ||||||||||
|
||||||||||
fcn( arrays[, options] ) | ||||||||||
|
||||||||||
where | ||||||||||
|
||||||||||
- arrays: array containing a one-dimensional subarray of the input | ||||||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||||||
ndarray and any additional ndarray arguments as zero-dimensional | ||||||||||
ndarrays. | ||||||||||
- options: function options. | ||||||||||
|
||||||||||
arrays: ArrayLikeObject<ndarray> | ||||||||||
Array-like object containing an input ndarray followed by any additional | ||||||||||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more.
Suggested change
|
||||||||||
ndarray arguments. | ||||||||||
|
||||||||||
dims: Array<integer> | ||||||||||
List of dimensions to which to apply a strided array function. | ||||||||||
|
||||||||||
options: Object (optional) | ||||||||||
Function options. | ||||||||||
|
||||||||||
Examples | ||||||||||
-------- | ||||||||||
// Define ndarray data and meta data... | ||||||||||
> var xbuf = [ 4.0, 3.0, 2.0, 1.0 ]; | ||||||||||
> var dtype = 'generic'; | ||||||||||
> var shx = [ 2, 2 ]; | ||||||||||
> var sx = [ 2, 1 ]; | ||||||||||
> var ox = 0; | ||||||||||
> var order = 'row-major'; | ||||||||||
|
||||||||||
// Define a wrapper for an extended BLAS function... | ||||||||||
> var f = {{alias:@stdlib/blas/ext/base/gsorthp}}.ndarray; | ||||||||||
> function fcn( arrays ) { | ||||||||||
... var x = arrays[ 0 ]; | ||||||||||
... var o = arrays[ 1 ]; | ||||||||||
... var N = x.shape[ 0 ]; | ||||||||||
... var dx = x.data; | ||||||||||
... var sx = x.strides[ 0 ]; | ||||||||||
... var ox = x.offset; | ||||||||||
... var init = o.data[ o.offset ]; | ||||||||||
... return f( N, init, dx, sx, ox ); | ||||||||||
... }; | ||||||||||
|
||||||||||
// Using minimal ndarray-like objects... | ||||||||||
> var x = { | ||||||||||
... 'dtype': dtype, | ||||||||||
... 'data': xbuf, | ||||||||||
... 'shape': shx, | ||||||||||
... 'strides': sx, | ||||||||||
... 'offset': ox, | ||||||||||
... 'order': order | ||||||||||
... }; | ||||||||||
> var sortOrder = { | ||||||||||
... 'dtype': dtype, | ||||||||||
... 'data': [ 1.0 ], | ||||||||||
... 'shape': [], | ||||||||||
... 'strides': [ 0 ], | ||||||||||
... 'offset': 0, | ||||||||||
... 'order': order | ||||||||||
... }; | ||||||||||
> {{alias}}( fcn, [ x, sortOrder ], [ 0, 1 ] ); | ||||||||||
> x.data | ||||||||||
[ 1.0, 2.0, 3.0, 4.0 ] | ||||||||||
|
||||||||||
See Also | ||||||||||
-------- | ||||||||||
|
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,63 @@ | ||
/** | ||
* @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 max-len */ | ||
|
||
'use strict'; | ||
|
||
var discreteUniform = require( '@stdlib/random/array/discrete-uniform' ); | ||
var ndarray2array = require( '@stdlib/ndarray/base/to-array' ); | ||
var numelDimension = require( '@stdlib/ndarray/base/numel-dimension' ); | ||
var getData = require( '@stdlib/ndarray/base/data-buffer' ); | ||
var getStride = require( '@stdlib/ndarray/base/stride' ); | ||
var getOffset = require( '@stdlib/ndarray/base/offset' ); | ||
var ndarraylike2scalar = require( '@stdlib/ndarray/base/ndarraylike2scalar' ); | ||
var gsorthp = require( '@stdlib/blas/ext/base/gsorthp' ).ndarray; | ||
var nullaryStrided1d = require( './../lib' ); | ||
|
||
function wrapper( arrays ) { | ||
var x = arrays[ 0 ]; | ||
var o = arrays[ 1 ]; | ||
return gsorthp( numelDimension( x, 0 ), ndarraylike2scalar( o ), getData( x ), getStride( x, 0 ), getOffset( x ) ); | ||
} | ||
|
||
var N = 10; | ||
var x = { | ||
'dtype': 'generic', | ||
'data': discreteUniform( N, -5, 5, { | ||
'dtype': 'generic' | ||
}), | ||
'shape': [ 1, 5, 2 ], | ||
'strides': [ 10, 2, 1 ], | ||
'offset': 0, | ||
'order': 'row-major' | ||
}; | ||
console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); | ||
|
||
var sortOrder = { | ||
'dtype': 'generic', | ||
'data': [ 1.0 ], | ||
'shape': [ 2 ], | ||
'strides': [ 0 ], | ||
'offset': 0, | ||
'order': 'row-major' | ||
}; | ||
|
||
nullaryStrided1d( wrapper, [ x, sortOrder ], [ 0, 1 ] ); | ||
|
||
console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) ); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.