Repository navigation
feat: add C implementation for blas/ext/base/ndarray/dcircshift
#16501
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
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 | ||||
|---|---|---|---|---|---|---|
|
|
@@ -108,6 +108,175 @@ console.log( ndarray2array( x ) ); | |||||
|
|
||||||
| <!-- /.examples --> | ||||||
|
|
||||||
| <!-- C interface documentation. --> | ||||||
|
|
||||||
| * * * | ||||||
|
|
||||||
| <section class="c"> | ||||||
|
|
||||||
| ## C APIs | ||||||
|
|
||||||
| <!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. --> | ||||||
|
|
||||||
| <section class="intro"> | ||||||
|
|
||||||
| </section> | ||||||
|
|
||||||
| <!-- /.intro --> | ||||||
|
|
||||||
| <!-- C usage documentation. --> | ||||||
|
|
||||||
| <section class="usage"> | ||||||
|
|
||||||
| ### Usage | ||||||
|
|
||||||
| ```c | ||||||
| #include "stdlib/blas/ext/base/ndarray/dcircshift.h" | ||||||
| ``` | ||||||
|
|
||||||
| #### stdlib_blas_ext_dcircshift( arrays ) | ||||||
|
|
||||||
| Circularly shifts the elements of a one-dimensional double-precision floating-point ndarray by a specified number of positions. | ||||||
|
|
||||||
| ```c | ||||||
| #include "stdlib/ndarray/ctor.h" | ||||||
| #include "stdlib/ndarray/dtypes.h" | ||||||
| #include "stdlib/ndarray/index_modes.h" | ||||||
| #include "stdlib/ndarray/orders.h" | ||||||
| #include "stdlib/ndarray/base/bytes_per_element.h" | ||||||
| #include <stdint.h> | ||||||
|
|
||||||
| // Create an input ndarray: | ||||||
| double dataX[] = { 1.0, 2.0, 3.0, 4.0, 5.0 }; | ||||||
| int64_t shapeX[] = { 5 }; | ||||||
| int64_t stridesX[] = { STDLIB_NDARRAY_FLOAT64_BYTES_PER_ELEMENT }; | ||||||
| int8_t submodes[] = { STDLIB_NDARRAY_INDEX_ERROR }; | ||||||
|
|
||||||
| struct ndarray *x = stdlib_ndarray_allocate( STDLIB_NDARRAY_FLOAT64, (uint8_t *)dataX, 1, shapeX, stridesX, 0, STDLIB_NDARRAY_ROW_MAJOR, STDLIB_NDARRAY_INDEX_ERROR, 1, submodes ); | ||||||
|
|
||||||
| // Create an ndarray containing the number of positions to shift: | ||||||
| const int64_t kdata[] = { 2 }; | ||||||
| int64_t kstrides[] = { 0 }; | ||||||
|
|
||||||
| struct ndarray *k = stdlib_ndarray_allocate( STDLIB_NDARRAY_INT64, (uint8_t *)kdata, 0, NULL, kstrides, 0, STDLIB_NDARRAY_ROW_MAJOR, STDLIB_NDARRAY_INDEX_ERROR, 1, submodes ); | ||||||
|
Member
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
|
||||||
|
|
||||||
| // Perform computation: | ||||||
| const struct ndarray *arrays[] = { x, k }; | ||||||
| stdlib_blas_ext_dcircshift( arrays ); | ||||||
|
|
||||||
| // Free allocated memory: | ||||||
| stdlib_ndarray_free( x ); | ||||||
| stdlib_ndarray_free( k ); | ||||||
| ``` | ||||||
|
|
||||||
| The function accepts the following arguments: | ||||||
|
|
||||||
| - **arrays**: `[in] struct ndarray**` list containing the following ndarrays: | ||||||
|
|
||||||
| - `[inout] struct ndarray*` a one-dimensional input ndarray. | ||||||
| - `[in] struct ndarray*` a zero-dimensional ndarray specifying the number of positions to shift. | ||||||
|
|
||||||
| ```c | ||||||
| void stdlib_blas_ext_dcircshift( const struct ndarray *arrays[] ); | ||||||
| ``` | ||||||
|
|
||||||
| </section> | ||||||
|
|
||||||
| <!-- /.usage --> | ||||||
|
|
||||||
| <!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. --> | ||||||
|
|
||||||
| <section class="notes"> | ||||||
|
|
||||||
| </section> | ||||||
|
|
||||||
| <!-- /.notes --> | ||||||
|
|
||||||
| <!-- C API usage examples. --> | ||||||
|
|
||||||
| <section class="examples"> | ||||||
|
|
||||||
| ### Examples | ||||||
|
|
||||||
| ```c | ||||||
| #include "stdlib/blas/ext/base/ndarray/dcircshift.h" | ||||||
| #include "stdlib/ndarray/ctor.h" | ||||||
| #include "stdlib/ndarray/dtypes.h" | ||||||
| #include "stdlib/ndarray/index_modes.h" | ||||||
| #include "stdlib/ndarray/orders.h" | ||||||
| #include "stdlib/ndarray/base/bytes_per_element.h" | ||||||
| #include <stdint.h> | ||||||
| #include <stdlib.h> | ||||||
| #include <stdio.h> | ||||||
|
|
||||||
| int main( void ) { | ||||||
| // Create a data buffer for the input vector: | ||||||
| double dataX[] = { 1.0, 2.0, 3.0, 4.0, 5.0 }; | ||||||
|
|
||||||
| // Specify the number of array dimensions: | ||||||
| const int64_t ndimsX = 1; | ||||||
|
|
||||||
| // Specify the array shape: | ||||||
| int64_t shapeX[] = { 5 }; | ||||||
|
|
||||||
| // Specify the array strides: | ||||||
| int64_t stridesX[] = { STDLIB_NDARRAY_FLOAT64_BYTES_PER_ELEMENT }; | ||||||
|
|
||||||
| // Specify the byte offset: | ||||||
| const int64_t offset = 0; | ||||||
|
|
||||||
| // Specify the array order: | ||||||
| const enum STDLIB_NDARRAY_ORDER order = STDLIB_NDARRAY_ROW_MAJOR; | ||||||
|
|
||||||
| // Specify the index mode: | ||||||
| const enum STDLIB_NDARRAY_INDEX_MODE imode = STDLIB_NDARRAY_INDEX_ERROR; | ||||||
|
|
||||||
| // Specify the subscript index modes: | ||||||
| int8_t submodes[] = { STDLIB_NDARRAY_INDEX_ERROR }; | ||||||
| const int64_t nsubmodes = 1; | ||||||
|
|
||||||
| // Create an input ndarray: | ||||||
| struct ndarray *x = stdlib_ndarray_allocate( STDLIB_NDARRAY_FLOAT64, (uint8_t *)dataX, ndimsX, shapeX, stridesX, offset, order, imode, nsubmodes, submodes ); | ||||||
|
|
||||||
| // Create a data buffer for an ndarray containing the number of positions to shift: | ||||||
| const int64_t kdata[] = { 2 }; | ||||||
|
Member
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
|
||||||
|
|
||||||
| // Specify the array strides for a zero-dimensional ndarray: | ||||||
| int64_t kstrides[] = { 0 }; | ||||||
|
|
||||||
| // Create an ndarray containing the shift amount: | ||||||
| struct ndarray *k = stdlib_ndarray_allocate( STDLIB_NDARRAY_INT64, (uint8_t *)kdata, 0, NULL, kstrides, 0, order, imode, nsubmodes, submodes ); | ||||||
|
Member
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
|
||||||
|
|
||||||
| if ( x == NULL || k == NULL ) { | ||||||
| fprintf( stderr, "Error allocating memory.\n" ); | ||||||
| exit( 1 ); | ||||||
| } | ||||||
|
|
||||||
| // Define a list of ndarrays: | ||||||
| const struct ndarray *arrays[] = { x, k }; | ||||||
|
|
||||||
| // Perform computation: | ||||||
| stdlib_blas_ext_dcircshift( arrays ); | ||||||
|
|
||||||
| // Print the result: | ||||||
| for ( int i = 0; i < 5; i++ ) { | ||||||
| printf( "x[ %i ] = %lf\n", i, dataX[ i ] ); | ||||||
| } | ||||||
|
|
||||||
| // Free allocated memory: | ||||||
| stdlib_ndarray_free( x ); | ||||||
| stdlib_ndarray_free( k ); | ||||||
| } | ||||||
| ``` | ||||||
|
|
||||||
| </section> | ||||||
|
|
||||||
| <!-- /.examples --> | ||||||
|
|
||||||
| </section> | ||||||
|
|
||||||
| <!-- /.c --> | ||||||
|
|
||||||
| <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. --> | ||||||
|
|
||||||
| <section class="related"> | ||||||
|
|
||||||
| Original file line number | Diff line number | Diff line change | ||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| @@ -0,0 +1,111 @@ | ||||||||||||||||||||||||||||||||||||
| /** | ||||||||||||||||||||||||||||||||||||
| * @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. | ||||||||||||||||||||||||||||||||||||
| */ | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| 'use strict'; | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| // MODULES // | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| var resolve = require( 'path' ).resolve; | ||||||||||||||||||||||||||||||||||||
| var bench = require( '@stdlib/bench' ); | ||||||||||||||||||||||||||||||||||||
| var uniform = require( '@stdlib/random/uniform' ); | ||||||||||||||||||||||||||||||||||||
| var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' ); | ||||||||||||||||||||||||||||||||||||
| var pow = require( '@stdlib/math/base/special/pow' ); | ||||||||||||||||||||||||||||||||||||
| var format = require( '@stdlib/string/format' ); | ||||||||||||||||||||||||||||||||||||
| var tryRequire = require( '@stdlib/utils/try-require' ); | ||||||||||||||||||||||||||||||||||||
| var pkg = require( './../package.json' ).name; | ||||||||||||||||||||||||||||||||||||
|
Comment on lines
+23
to
+30
Member
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
|
||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| // VARIABLES // | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| var dcircshift = tryRequire( resolve( __dirname, './../lib/native.js' ) ); | ||||||||||||||||||||||||||||||||||||
| var opts = { | ||||||||||||||||||||||||||||||||||||
| 'skip': ( dcircshift instanceof Error ) | ||||||||||||||||||||||||||||||||||||
| }; | ||||||||||||||||||||||||||||||||||||
| var options = { | ||||||||||||||||||||||||||||||||||||
| 'dtype': 'float64' | ||||||||||||||||||||||||||||||||||||
| }; | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| // FUNCTIONS // | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| /** | ||||||||||||||||||||||||||||||||||||
| * Creates a benchmark function. | ||||||||||||||||||||||||||||||||||||
| * | ||||||||||||||||||||||||||||||||||||
| * @private | ||||||||||||||||||||||||||||||||||||
| * @param {PositiveInteger} len - ndarray length | ||||||||||||||||||||||||||||||||||||
|
Member
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
|
||||||||||||||||||||||||||||||||||||
| * @returns {Function} benchmark function | ||||||||||||||||||||||||||||||||||||
| */ | ||||||||||||||||||||||||||||||||||||
| function createBenchmark( len ) { | ||||||||||||||||||||||||||||||||||||
| var x = uniform( [ len ], -100.0, 100.0, options ); | ||||||||||||||||||||||||||||||||||||
| var k = scalar2ndarray( 2, { | ||||||||||||||||||||||||||||||||||||
| 'dtype': 'int64' | ||||||||||||||||||||||||||||||||||||
| }); | ||||||||||||||||||||||||||||||||||||
|
Comment on lines
+54
to
+57
Member
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
|
||||||||||||||||||||||||||||||||||||
| return benchmark; | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| /** | ||||||||||||||||||||||||||||||||||||
| * Benchmark function. | ||||||||||||||||||||||||||||||||||||
| * | ||||||||||||||||||||||||||||||||||||
| * @private | ||||||||||||||||||||||||||||||||||||
| * @param {Benchmark} b - benchmark instance | ||||||||||||||||||||||||||||||||||||
| */ | ||||||||||||||||||||||||||||||||||||
| function benchmark( b ) { | ||||||||||||||||||||||||||||||||||||
| var out; | ||||||||||||||||||||||||||||||||||||
| var i; | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| b.tic(); | ||||||||||||||||||||||||||||||||||||
| for ( i = 0; i < b.iterations; i++ ) { | ||||||||||||||||||||||||||||||||||||
| out = dcircshift( [ x, k ] ); | ||||||||||||||||||||||||||||||||||||
| if ( typeof out !== 'object' ) { | ||||||||||||||||||||||||||||||||||||
| b.fail( 'should return an ndarray' ); | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
| b.toc(); | ||||||||||||||||||||||||||||||||||||
| if ( typeof out !== 'object' ) { | ||||||||||||||||||||||||||||||||||||
| b.fail( 'should return an ndarray' ); | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
|
Comment on lines
+78
to
+80
Member
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
|
||||||||||||||||||||||||||||||||||||
| b.pass( 'benchmark finished' ); | ||||||||||||||||||||||||||||||||||||
| b.end(); | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| // MAIN // | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| /** | ||||||||||||||||||||||||||||||||||||
| * Main execution sequence. | ||||||||||||||||||||||||||||||||||||
| * | ||||||||||||||||||||||||||||||||||||
| * @private | ||||||||||||||||||||||||||||||||||||
| */ | ||||||||||||||||||||||||||||||||||||
| function main() { | ||||||||||||||||||||||||||||||||||||
| var len; | ||||||||||||||||||||||||||||||||||||
| var min; | ||||||||||||||||||||||||||||||||||||
| var max; | ||||||||||||||||||||||||||||||||||||
| var f; | ||||||||||||||||||||||||||||||||||||
| var i; | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| min = 1; // 10^min | ||||||||||||||||||||||||||||||||||||
| max = 6; // 10^max | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| for ( i = min; i <= max; i++ ) { | ||||||||||||||||||||||||||||||||||||
| len = pow( 10, i ); | ||||||||||||||||||||||||||||||||||||
| f = createBenchmark( len ); | ||||||||||||||||||||||||||||||||||||
| bench( format( '%s::native:len=%d', pkg, len ), opts, f ); | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
| } | ||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||
| main(); | ||||||||||||||||||||||||||||||||||||
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.