MOOCHO (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DenseLinAlgLAPack.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ***********************************************************************
3 //
4 // Moocho: Multi-functional Object-Oriented arCHitecture for Optimization
5 // Copyright (2003) Sandia Corporation
6 //
7 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8 // license for use of this work by or on behalf of the U.S. Government.
9 //
10 // Redistribution and use in source and binary forms, with or without
11 // modification, are permitted provided that the following conditions are
12 // met:
13 //
14 // 1. Redistributions of source code must retain the above copyright
15 // notice, this list of conditions and the following disclaimer.
16 //
17 // 2. Redistributions in binary form must reproduce the above copyright
18 // notice, this list of conditions and the following disclaimer in the
19 // documentation and/or other materials provided with the distribution.
20 //
21 // 3. Neither the name of the Corporation nor the names of the
22 // contributors may be used to endorse or promote products derived from
23 // this software without specific prior written permission.
24 //
25 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36 //
37 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov)
38 //
39 // ***********************************************************************
40 // @HEADER
41 
42 #ifndef LIN_ALG_LA_PACK_H
43 #define LIN_ALG_LA_PACK_H
44 
46 
47 namespace DenseLinAlgLAPack {
48 
60 void potrf( DMatrixSliceTriEle* A );
61 
83 void geqrf( DMatrixSlice* A, DVectorSlice* tau, DVectorSlice* work );
84 
104 void ormrq(
106  , const DMatrixSlice& A, const DVectorSlice& tau
107  , DMatrixSlice* C, DVectorSlice* work
108  );
109 
129 void sytrf( DMatrixSliceTriEle* A, FortranTypes::f_int ipiv[], DVectorSlice* work );
130 
148 void sytrs(
149  const DMatrixSliceTriEle& A, FortranTypes::f_int ipiv[]
150  ,DMatrixSlice* B, DVectorSlice* work
151  );
152 
171 
190 void getrs(
191  const DMatrixSlice& LU, const FortranTypes::f_int ipiv[], BLAS_Cpp::Transp transp
192  ,DMatrixSlice* B
193  );
194 
195 } // end namespace DenseLinAlgLAPack
196 
197 #endif // LIN_ALG_LA_PACK_H
FortranTypes::f_int f_int
void sytrf(DMatrixSliceTriEle *A, FortranTypes::f_int ipiv[], DVectorSlice *work)
Calls xSYTRF to compute the P*A*P' = L'*D*L factorization of a symmetric possibly indefinite matrix...
void getrs(const DMatrixSlice &LU, const FortranTypes::f_int ipiv[], BLAS_Cpp::Transp transp, DMatrixSlice *B)
Calls xGETRS to solve linear systems with the factors of P'*A = L*U generated by xGETRF.
void geqrf(DMatrixSlice *A, DVectorSlice *tau, DVectorSlice *work)
Calls xGEQRF to compute the QR factorization of a matrix A.
void getrf(DMatrixSlice *A, FortranTypes::f_int ipiv[], FortranTypes::f_int *rank)
Calls xGETRF to compute the P'*A = L*U factorization of a general rectuangular matrix.
Transposed.
void potrf(DMatrixSliceTriEle *A)
Calls xPOTRF to compute the cholesky factorization of a symmetric positive definte matrix...
DenseLinAlgPack::VectorSliceTmpl< value_type > DVectorSlice
void ormrq(BLAS_Cpp::Side side, BLAS_Cpp::Transp trans, const DMatrixSlice &A, const DVectorSlice &tau, DMatrixSlice *C, DVectorSlice *work)
Calls xORMRQ to compute a matrix matrix op(C)=op(Q)*op(C) where Q is stored in A and tau computed fro...
Transp
TRANS.
void sytrs(const DMatrixSliceTriEle &A, FortranTypes::f_int ipiv[], DMatrixSlice *B, DVectorSlice *work)
Calls xSYTRS(...) to compute the solution of the factorized system A * X = B where A was factorized b...