Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view

Large diffs are not rendered by default.

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
@@ -0,0 +1,119 @@
/**
* @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.
*/

/* eslint-disable max-len */

'use strict';

var Float64Array = require( '@stdlib/array/float64' );
var ndarray2array = require( '@stdlib/ndarray/base/to-array' );
var gwxpy = require( '@stdlib/blas/ext/base/ndarray/gwxpy' );
var strategy = require( '@stdlib/ndarray/base/kernels/generic/unary-strided1d/strategy' );
var kernel = require( './../lib' );

// Create data buffers:
var xbuf = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0 ] );
var ybuf = new Float64Array( [ 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0 ] );
var zbuf = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );

// Define the array shapes:
var xsh = [ 1, 1, 1, 1, 3, 2, 2 ];
var ysh = [ 1, 1, 1, 1, 3, 2, 2 ];
var zsh = [ 1, 1, 1, 1, 3, 2, 2 ];

// Define the array strides:
var sx = [ 12, 12, 12, 12, 4, 2, 1 ];
var sy = [ 12, 12, 12, 12, 4, 2, 1 ];
var sz = [ 12, 12, 12, 12, 4, 2, 1 ];

// Define the index offsets:
var ox = 0;
var oy = 0;
var oz = 0;

// Create the input ndarray descriptors:
var x = {
'dtype': 'float64',
'data': xbuf,
'shape': xsh,
'strides': sx,
'offset': ox,
'order': 'row-major'
};

var y = {
'dtype': 'float64',
'data': ybuf,
'shape': ysh,
'strides': sy,
'offset': oy,
'order': 'row-major'
};

// Create an output ndarray descriptor:
var z = {
'dtype': 'float64',
'data': zbuf,
'shape': zsh,
'strides': sz,
'offset': oz,
'order': 'row-major'
};

// Initialize ndarray descriptors representing subarray views:
var views = [
{
'dtype': x.dtype,
'data': x.data,
'shape': [ 2, 2 ],
'strides': [ 2, 1 ],
'offset': x.offset,
'order': x.order
},
{
'dtype': y.dtype,
'data': y.data,
'shape': [ 2, 2 ],
'strides': [ 2, 1 ],
'offset': y.offset,
'order': y.order
},
{
'dtype': z.dtype,
'data': z.data,
'shape': [ 2, 2 ],
'strides': [ 2, 1 ],
'offset': z.offset,
'order': z.order
}
];

// Resolve input/output strategies when iterating over subarray views:
var strategyX = strategy( views[ 0 ] );
var strategyY = strategy( views[ 1 ] );
var strategyZ = strategy( views[ 2 ] );

// Resolve a kernel:
var f = kernel( 5 );

// Apply strided function:
f( gwxpy, [ x, y, z ], views, [ 1, 1, 1, 1, 3 ], [ 12, 12, 12, 12, 4 ], [ 12, 12, 12, 12, 4 ], [ 12, 12, 12, 12, 4 ], strategyX, strategyY, strategyZ, {} );

console.log( ndarray2array( x.data, x.shape, x.strides, x.offset, x.order ) );
console.log( ndarray2array( y.data, y.shape, y.strides, y.offset, y.order ) );
console.log( ndarray2array( z.data, z.shape, z.strides, z.offset, z.order ) );
Loading