11 #ifndef __Panzer_BlockedTpetraLinearObjContainer_hpp__
12 #define __Panzer_BlockedTpetraLinearObjContainer_hpp__
14 #include "PanzerDiscFE_config.hpp"
17 #include "Tpetra_Vector.hpp"
18 #include "Tpetra_CrsMatrix.hpp"
26 #include "Thyra_PhysicallyBlockedLinearOpBase.hpp"
27 #include "Thyra_ProductVectorBase.hpp"
28 #include "Thyra_TpetraThyraWrappers.hpp"
30 #include <unordered_map>
37 template <
typename ScalarT,
typename LocalOrdinalT,
typename GlobalOrdinalT,
typename NodeT=panzer::TpetraNodeType>
38 class BlockedTpetraLinearObjContainer :
public LinearObjContainer
39 ,
public ThyraObjContainer<ScalarT> {
44 typedef Tpetra::Map<LocalOrdinalT,GlobalOrdinalT,NodeT>
MapType;
void set_x_th(const Teuchos::RCP< VectorType > &in)
Teuchos::RCP< VectorType > get_dxdt() const
bool checkCompatibility() const
Make sure row and column spaces match up.
std::vector< Teuchos::RCP< const MapType > > blockMaps_
Teuchos::RCP< VectorType > get_x() const
void set_f_th(const Teuchos::RCP< VectorType > &in)
void set_dxdt(const Teuchos::RCP< VectorType > &in)
void initializeMatrix(ScalarT value)
Put a particular scalar in the matrix.
Tpetra::Map< LocalOrdinalT, GlobalOrdinalT, NodeT > MapType
void set_dxdt_th(const Teuchos::RCP< VectorType > &in)
Thyra::VectorBase< ScalarT > VectorType
virtual void initialize()
Teuchos::RCP< VectorType > get_x_th() const
Teuchos::RCP< VectorType > dxdt
Teuchos::RCP< const MapType > getMapForBlock(std::size_t i) const
Teuchos::RCP< CrsMatrixType > A
Teuchos::RCP< VectorType > get_dxdt_th() const
void set_A_th(const Teuchos::RCP< CrsMatrixType > &in)
Thyra::LinearOpBase< ScalarT > CrsMatrixType
void set_x(const Teuchos::RCP< VectorType > &in)
Teuchos::RCP< VectorType > f
Teuchos::RCP< VectorType > x
Teuchos::RCP< VectorType > get_f() const
Teuchos::RCP< CrsMatrixType > get_A_th() const
Teuchos::RCP< CrsMatrixType > get_A() const
Teuchos::RCP< VectorType > get_f_th() const
void set_f(const Teuchos::RCP< VectorType > &in)
void setMapsForBlocks(const std::vector< Teuchos::RCP< const MapType > > &blockMaps)
void set_A(const Teuchos::RCP< CrsMatrixType > &in)