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
94 changes: 83 additions & 11 deletions lib/node_modules/@stdlib/lapack/base/dlaswp/README.md
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@ Performs a series of row interchanges on an input matrix `A` using pivot indices
var Int32Array = require( '@stdlib/array/int32' );
var Float64Array = require( '@stdlib/array/float64' );

var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); // => [ [ 1.0, 2.0 ], [ 3.0, 4.0 ], [ 5.0, 6.0 ] ]
var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
var IPIV = new Int32Array( [ 2, 0, 1 ] );

dlaswp( 'row-major', 2, A, 2, 0, 2, IPIV, 1 );
Expand All @@ -62,7 +62,7 @@ The sign of the increment parameter `incx` determines the order in which pivots
var Int32Array = require( '@stdlib/array/int32' );
var Float64Array = require( '@stdlib/array/float64' );

var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); // => [ [ 1.0, 2.0 ], [ 3.0, 4.0 ], [ 5.0, 6.0 ] ]
var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
var IPIV = new Int32Array( [ 2, 0, 1 ] );

dlaswp( 'row-major', 2, A, 2, 0, 2, IPIV, -1 );
Expand All @@ -75,7 +75,7 @@ To perform strided access over `IPIV`, provide an `abs(incx)` value greater than
var Int32Array = require( '@stdlib/array/int32' );
var Float64Array = require( '@stdlib/array/float64' );

var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); // => [ [ 1.0, 2.0 ], [ 3.0, 4.0 ], [ 5.0, 6.0 ] ]
var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
var IPIV = new Int32Array( [ 2, 999, 0, 999, 1 ] );

dlaswp( 'row-major', 2, A, 2, 0, 2, IPIV, 2 );
Expand Down Expand Up @@ -110,7 +110,7 @@ Performs a series of row interchanges on the matrix `A` using pivot indices stor
var Int32Array = require( '@stdlib/array/int32' );
var Float64Array = require( '@stdlib/array/float64' );

var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] ); // => [ [ 1.0, 2.0 ], [ 3.0, 4.0 ], [ 5.0, 6.0 ] ]
var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
var IPIV = new Int32Array( [ 2, 0, 1 ] );

dlaswp.ndarray( 2, A, 2, 1, 0, 0, 2, 1, IPIV, 1, 0 );
Expand Down Expand Up @@ -219,21 +219,66 @@ console.log( ndarray2array( A, shape, strides, offset, order ) );
### Usage

```c
TODO
#include "stdlib/lapack/base/dlaswp.h"
```

#### TODO
#### c\_dlaswp( layout, N, \*A, LDA, k1, k2, \*IPIV, incx )

TODO.
Performs a series of row interchanges on a matrix `A` using pivot indices stored in `IPIV`.

```c
TODO
#include "stdlib/lapack/base/shared.h"

double A[ 3*2 ] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 };
const LAPACK_INT IPIV[ 3 ] = { 2, 0, 1 };

c_dlaswp( LAPACK_ROW_MAJOR, 2, A, 2, 0, 2, IPIV, 1 );
```

The function accepts the following arguments:

- **layout**: `[in] LAPACK_LAYOUT` storage layout.
- **N**: `[in] LAPACK_INT` number of columns in `A`.
- **A**: `[inout] double*` input matrix.
- **LDA**: `[in] LAPACK_INT` stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`).
- **k1**: `[in] LAPACK_INT` index of first row to interchange.
- **k2**: `[in] LAPACK_INT` index of last row to interchange.
- **IPIV**: `[in] LAPACK_INT*` vector of pivot indices.
- **incx**: `[in] LAPACK_INT` increment between successive values of `IPIV`.

```c
void c_dlaswp( const LAPACK_LAYOUT layout, const LAPACK_INT N, double *A, const LAPACK_INT LDA, const LAPACK_INT k1, const LAPACK_INT k2, const LAPACK_INT *IPIV, const LAPACK_INT incx );
```

TODO
#### c\_dlaswp\_ndarray( N, \*A, sa1, sa2, oa, k1, k2, inck, \*IPIV, si, oi )

Performs a series of row interchanges on a matrix `A` using pivot indices stored in `IPIV` and alternative indexing semantics.

```c
#include "stdlib/lapack/base/shared.h"

double A[ 3*2 ] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 };
const LAPACK_INT IPIV[ 3 ] = { 2, 0, 1 };

c_dlaswp_ndarray( 2, A, 2, 1, 0, 0, 2, 1, IPIV, 1, 0 );
```

The function accepts the following arguments:

- **N**: `[in] LAPACK_INT` number of columns in `A`.
- **A**: `[inout] double*` input matrix.
- **sa1**: `[in] LAPACK_INT` stride of the first dimension of `A`.
- **sa2**: `[in] LAPACK_INT` stride of the second dimension of `A`.
- **oa**: `[in] LAPACK_INT` starting index for `A`.
- **k1**: `[in] LAPACK_INT` index of first row to interchange.
- **k2**: `[in] LAPACK_INT` index of last row to interchange.
- **inck**: `[in] LAPACK_INT` direction in which to apply pivots (-1 to apply pivots in reverse order; otherwise, apply in provided order).
- **IPIV**: `[in] LAPACK_INT*` vector of pivot indices.
- **si**: `[in] LAPACK_INT` index increment for `IPIV`.
- **oi**: `[in] LAPACK_INT` starting index for `IPIV`.

```c
TODO
void c_dlaswp_ndarray( const LAPACK_INT N, double *A, const LAPACK_INT strideA1, const LAPACK_INT strideA2, const LAPACK_INT offsetA, const LAPACK_INT k1, const LAPACK_INT k2, const LAPACK_INT inck, const LAPACK_INT *IPIV, const LAPACK_INT strideIPIV, const LAPACK_INT offsetIPIV );
```

</section>
Expand All @@ -255,7 +300,34 @@ TODO
### Examples

```c
TODO
#include "stdlib/lapack/base/dlaswp.h"
#include "stdlib/lapack/base/shared.h"
#include <stdio.h>

int main( void ) {
// Define a 3x3 matrix stored in row-major order:
double A[ 3*3 ] = {
1.0, 2.0, 3.0,
4.0, 5.0, 6.0,
7.0, 8.0, 9.0
};

// Define a vector of pivot indices:
const LAPACK_INT IPIV[ 3 ] = { 2, 0, 1 };

// Specify the number of columns:
const int N = 3;

// Perform row interchanges:
c_dlaswp( LAPACK_ROW_MAJOR, N, A, N, 0, 2, IPIV, 1 );

// Print the result:
for ( int i = 0; i < 3; i++ ) {
for ( int j = 0; j < N; j++ ) {
printf( "A[ %i, %i ] = %lf\n", i, j, A[ (i*N)+j ] );
}
}
}
```

</section>
Expand Down
146 changes: 146 additions & 0 deletions lib/node_modules/@stdlib/lapack/base/dlaswp/benchmark/c/Makefile
Original file line number Diff line number Diff line change
@@ -0,0 +1,146 @@
#/
# @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.
#/

# VARIABLES #

ifndef VERBOSE
QUIET := @
else
QUIET :=
endif

# Determine the OS ([1][1], [2][2]).
#
# [1]: https://en.wikipedia.org/wiki/Uname#Examples
# [2]: http://stackoverflow.com/a/27776822/2225624
OS ?= $(shell uname)
ifneq (, $(findstring MINGW,$(OS)))
OS := WINNT
else
ifneq (, $(findstring MSYS,$(OS)))
OS := WINNT
else
ifneq (, $(findstring CYGWIN,$(OS)))
OS := WINNT
else
ifneq (, $(findstring Windows_NT,$(OS)))
OS := WINNT
endif
endif
endif
endif

# Define the program used for compiling C source files:
ifdef C_COMPILER
CC := $(C_COMPILER)
else
CC := gcc
endif

# Define the command-line options when compiling C files:
CFLAGS ?= \
-std=c99 \
-O3 \
-Wall \
-pedantic

# Determine whether to generate position independent code ([1][1], [2][2]).
#
# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options
# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option
ifeq ($(OS), WINNT)
fPIC ?=
else
fPIC ?= -fPIC
endif

# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`):
INCLUDE ?=

# List of source files:
SOURCE_FILES ?=

# List of libraries (e.g., `-lopenblas -lpthread`):
LIBRARIES ?=

# List of library paths (e.g., `-L /foo/bar -L /beep/boop`):
LIBPATH ?=

# List of C targets:
c_targets := benchmark.size.out


# RULES #

#/
# Compiles source files.
#
# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`)
# @param {string} [CFLAGS] - C compiler options
# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`)
# @param {string} [SOURCE_FILES] - list of source files
# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`)
# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`)
#
# @example
# make
#
# @example
# make all
#/
all: $(c_targets)

.PHONY: all

#/
# Compiles C source files.
#
# @private
# @param {string} CC - C compiler (e.g., `gcc`)
# @param {string} CFLAGS - C compiler options
# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`)
# @param {string} SOURCE_FILES - list of source files
# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`)
# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`)
#/
$(c_targets): %.out: %.c
$(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES)

#/
# Runs compiled benchmarks.
#
# @example
# make run
#/
run: $(c_targets)
$(QUIET) ./$<

.PHONY: run

#/
# Removes generated files.
#
# @example
# make clean
#/
clean:
$(QUIET) -rm -f *.o *.out

.PHONY: clean
Loading