trsm.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_DENSE_TRSM_H
#define MYRAMATH_DENSE_TRSM_H

#include <myramath/MYRAMATH_EXPORT.h>
#include <myramath/utility/Number.h>

namespace myra {

// Forward declarations.
template<class Number> class MatrixRange;
template<class Number> class CMatrixRange;
template<class Number> class CLowerMatrixRange;

//    side (input): solve from the 'L'eft (solve A*X=B) or 'R'ight (solve X*A=B)
//    uplo (input): does A's data exist in its 'U'pper or 'L'ower triangle?
//    op (input): any operation to apply to A? Can 'T'ranspose it, 'H'ermitian it, 'C'onjugate it, or 'N'othing
//    A (input): a square matrix range, only specified triangle is accessed
//    B (input/output): forcing data B, overwritten with solution data X
//    diag (input): override A to have 'U'nit diagonal, or 'N'ot? Defaults to 'N'
//    alpha (input): scale B prior to solution? Defaults to 1
MYRAMATH_EXPORT void trsm_inplace(char side, char uplo, char op, const CMatrixRange<NumberS>& A, const MatrixRange<NumberS>& B, char diag='N', NumberS alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char uplo, char op, const CMatrixRange<NumberD>& A, const MatrixRange<NumberD>& B, char diag='N', NumberD alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char uplo, char op, const CMatrixRange<NumberC>& A, const MatrixRange<NumberC>& B, char diag='N', NumberC alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char uplo, char op, const CMatrixRange<NumberZ>& A, const MatrixRange<NumberZ>& B, char diag='N', NumberZ alpha=1);

//    side (input): solve from the 'L'eft (solve A*X=B) or 'R'ight (solve X*A=B)
//    op (input): any operation to apply to A? Can 'T'ranspose it, 'H'ermitian it, 'C'onjugate it, or 'N'othing
//    A (input): the LowerMatrix being solved by
//    B (input/output): forcing data B, overwritten with solution data X
//    diag (input): override A to have 'U'nit diagonal, or 'N'ot? Defaults to 'N'
//    alpha (input): scale B prior to solution? Defaults to 1
MYRAMATH_EXPORT void trsm_inplace(char side, char op, const CLowerMatrixRange<NumberS>& A, const MatrixRange<NumberS>& B, char diag='N', NumberS alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char op, const CLowerMatrixRange<NumberD>& A, const MatrixRange<NumberD>& B, char diag='N', NumberD alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char op, const CLowerMatrixRange<NumberC>& A, const MatrixRange<NumberC>& B, char diag='N', NumberC alpha=1);
MYRAMATH_EXPORT void trsm_inplace(char side, char op, const CLowerMatrixRange<NumberZ>& A, const MatrixRange<NumberZ>& B, char diag='N', NumberZ alpha=1);

} // namespace

#endif