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 |
#define STR(x) #x
#define STRING(x) STR(x) #include <flens/lapack/interface/include/config.h> namespace flens { namespace lapack { extern "C" { //-- dgetrs -------------------------------------------------------------------- void LAPACK_DECL(dgetrs)(const char *TRANS, const INTEGER *N, const INTEGER *NRHS, const DOUBLE *A, const INTEGER *LDA, const INTEGER *IPIV, DOUBLE *B, const INTEGER *LDB, INTEGER *INFO) { // // Test the input parameters so that we pass LAPACK error checks // *INFO = 0; if (*TRANS!='N' && *TRANS!='T' && *TRANS!='C') { *INFO = -1; } else if (*N<0) { *INFO = -2; } else if (*NRHS<0) { *INFO = -3; } else if (*LDA<std::max(INTEGER(1), *N)) { *INFO = -5; } else if (*LDB<std::max(INTEGER(1), *N)) { *INFO = -8; } if (*INFO!=0) { *INFO = -(*INFO); LAPACK_ERROR("DGETRS", INFO); *INFO = -(*INFO); return; } // // Call FLENS implementation // Transpose trans = convertTo<Transpose>(*TRANS); DConstGeMatrixView _A = DConstFSView(*N, *N, A, *LDA); IConstDenseVectorView _IPIV = IConstArrayView(*N, IPIV, 1); DGeMatrixView _B = DFSView(*N, *NRHS, B, *LDB); trs(trans, _A, _IPIV, _B); } //-- zgetrs -------------------------------------------------------------------- void LAPACK_DECL(zgetrs)(const char *TRANS, const INTEGER *N, const INTEGER *NRHS, const DOUBLE_COMPLEX *A, const INTEGER *LDA, const INTEGER *IPIV, DOUBLE_COMPLEX *B, const INTEGER *LDB, INTEGER *INFO) { // // Test the input parameters so that we pass LAPACK error checks // *INFO = 0; if (*TRANS!='N' && *TRANS!='T' && *TRANS!='C') { *INFO = -1; } else if (*N<0) { *INFO = -2; } else if (*NRHS<0) { *INFO = -3; } else if (*LDA<std::max(INTEGER(1), *N)) { *INFO = -5; } else if (*LDB<std::max(INTEGER(1), *N)) { *INFO = -8; } if (*INFO!=0) { *INFO = -(*INFO); LAPACK_ERROR("ZGETRS", INFO); *INFO = -(*INFO); return; } // // Call FLENS implementation // const auto zA = reinterpret_cast<const CXX_DOUBLE_COMPLEX *>(A); auto zB = reinterpret_cast<CXX_DOUBLE_COMPLEX *>(B); Transpose trans = convertTo<Transpose>(*TRANS); ZConstGeMatrixView _A = ZConstFSView(*N, *N, zA, *LDA); IConstDenseVectorView _IPIV = IConstArrayView(*N, IPIV, 1); ZGeMatrixView _B = ZFSView(*N, *NRHS, zB, *LDB); trs(trans, _A, _IPIV, _B); } } // extern "C" } } // namespace lapack, flens |