AbstractLinAlgPack: C++ Interfaces For Vectors, Matrices And Related Linear Algebra Objects  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
AbstractLinAlgPack_MatrixSymNonsingSerial.hpp
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 SLAP_MATRIX_SYM_NONSINGULAR_SERIAL_H
43 #define SLAP_MATRIX_SYM_NONSINGULAR_SERIAL_H
44 
45 #include "AbstractLinAlgPack_MatrixNonsingSerial.hpp"
46 #include "AbstractLinAlgPack_MatrixSymNonsing.hpp"
47 
48 namespace AbstractLinAlgPack {
49 
57  : virtual public MatrixNonsingSerial
58  , virtual public AbstractLinAlgPack::MatrixSymNonsing // doxygen needs full name
59 {
60 public:
61 
64 
67 
74  virtual void M_StMtInvMtM(
75  DMatrixSliceSym* sym_gms_lhs, value_type alpha
76  ,const MatrixOpSerial& mwo, BLAS_Cpp::Transp mwo_trans
78  ) const;
79 
81 
84 
85  void M_StMtInvMtM(
86  MatrixSymOp* sym_lhs, value_type alpha
87  ,const MatrixOp& mwo, BLAS_Cpp::Transp mwo_trans
89  ) const;
90 
92 
93 }; // end class MatrixSymNonsingSerial
94 
101 
102 inline
105  DMatrixSliceSym* sym_gms_lhs, value_type alpha
106  ,const MatrixOpSerial& mwo, BLAS_Cpp::Transp mwo_trans
107  ,const MatrixSymNonsingSerial& mswons
109  )
110 {
111  mswons.M_StMtInvMtM(sym_gms_lhs,alpha,mwo,mwo_trans,mwo_rhs);
112 }
113 
115 
116 } // end namespace AbstractLinAlgPack
117 
118 #endif // SLAP_MATRIX_SYM_NONSINGULAR_SERIAL_H
friend void M_StMtInvMtM(MatrixSymOp *sym_gms_lhs, value_type alpha, const MatrixOp &mwo, BLAS_Cpp::Transp mwo_trans, const MatrixSymNonsing &mswof, EMatrixDummyArg mwo_rhs)
Abstract base class for all polymorphic symmetrix nonsingular matrices that can be used to solve for ...
Abstract base class for all serial polymorphic symmetrix nonsingular matrices that can be used to sol...
Base class for all matrices implemented in a shared memory address space.
Interface adding operations specific for a symmetric matrix {abstract}.
Abstract base class for all AbstractLinAlgPack::MatrixNonsing objects implemented in shared memory sp...
void M_StMtInvMtM(DMatrixSliceSym *sym_gms_lhs, value_type alpha, const MatrixOpSerial &mwo, BLAS_Cpp::Transp mwo_trans, const MatrixSymNonsingSerial &mswons, MatrixSymNonsingSerial::EMatrixDummyArg mwo_rhs)
sym_gms_lhs = alpha * op(mwo) * inv(mswof) * op(mwo)'
Base class for all matrices that support basic matrix operations.
virtual void M_StMtInvMtM(DMatrixSliceSym *sym_gms_lhs, value_type alpha, const MatrixOpSerial &mwo, BLAS_Cpp::Transp mwo_trans, EMatrixDummyArg) const
sym_gms_lhs = alpha * op(mwo) * inv(M) * op(mwo)'.
Transp