Press n or j to go to the next uncovered block, b, p or k for the previous block.
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 | 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 1x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 7x 1x 1x 1x 1x 1x | /**
* @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 getShape = require( '@stdlib/ndarray/base/shape' );
var getStrides = require( '@stdlib/ndarray/base/strides' );
var getStride = require( '@stdlib/ndarray/base/stride' );
var getOffset = require( '@stdlib/ndarray/base/offset' );
var getData = require( '@stdlib/ndarray/base/data-buffer' );
var ndarraylike2scalar = require( '@stdlib/ndarray/base/ndarraylike2scalar' );
var resolveTriangleStr = require( '@stdlib/blas/base/matrix-triangle-resolve-str' );
var resolveTransposeStr = require( '@stdlib/blas/base/transpose-operation-resolve-str' );
var resolveDiagonalStr = require( '@stdlib/blas/base/diagonal-type-resolve-str' );
var strided = require( '@stdlib/blas/base/dtrsv' ).ndarray;
// MAIN //
/**
* Solves one of the systems of equations `A*x = b` or `A^T*x = b`, where `b` and `x` are one-dimensional ndarrays and `A` is an `N` by `N` unit, or non-unit, upper or lower triangular matrix.
*
* ## Notes
*
* - The function expects the following ndarrays:
*
* - a one-dimensional input/output ndarray corresponding to `x`.
* - a two-dimensional input ndarray corresponding to `A`.
* - a zero-dimensional ndarray specifying whether the upper or lower triangular part of `A` should be referenced.
* - a zero-dimensional ndarray specifying whether `A` should be transposed, conjugate-transposed, or not transposed.
* - a zero-dimensional ndarray specifying whether `A` has a unit or non-unit diagonal.
*
* @param {ArrayLikeObject<Object>} arrays - array-like object containing ndarrays
* @returns {Object} output ndarray
*
* @example
* var Float64Matrix = require( '@stdlib/ndarray/matrix/float64' );
* var Float64Vector = require( '@stdlib/ndarray/vector/float64' );
* var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
* var resolveTriangle = require( '@stdlib/blas/base/matrix-triangle-resolve-enum' );
* var resolveTranspose = require( '@stdlib/blas/base/transpose-operation-resolve-enum' );
* var resolveDiagonal = require( '@stdlib/blas/base/diagonal-type-resolve-enum' );
*
* var x = new Float64Vector( [ 1.0, 2.0, 3.0 ] );
* var A = new Float64Matrix( [ [ 1.0, 2.0, 3.0 ], [ 0.0, 4.0, 5.0 ], [ 0.0, 0.0, 6.0 ] ] );
*
* var uplo = scalar2ndarray( resolveTriangle( 'upper' ), {
* 'dtype': 'int32'
* });
* var trans = scalar2ndarray( resolveTranspose( 'no-transpose' ), {
* 'dtype': 'int32'
* });
* var diag = scalar2ndarray( resolveDiagonal( 'non-unit' ), {
* 'dtype': 'int32'
* });
*
* var out = dtrsv( [ x, A, uplo, trans, diag ] );
* // returns <ndarray>[ -0.25, -0.125, 0.5 ]
*
* var bool = ( out === x );
* // returns true
*/
function dtrsv( arrays ) {
var trans;
var diag;
var uplo;
var sh;
var st;
var x;
var A;
x = arrays[ 0 ];
A = arrays[ 1 ];
uplo = resolveTriangleStr( ndarraylike2scalar( arrays[ 2 ] ) );
trans = resolveTransposeStr( ndarraylike2scalar( arrays[ 3 ] ) );
diag = resolveDiagonalStr( ndarraylike2scalar( arrays[ 4 ] ) );
sh = getShape( A, false );
st = getStrides( A, false );
strided( uplo, trans, diag, sh[ 0 ], getData( A ), st[ 0 ], st[ 1 ], getOffset( A ), getData( x ), getStride( x, 0 ), getOffset( x ) );
return x;
}
// EXPORTS //
module.exports = dtrsv;
|