43#ifndef __Panzer_BlockedTpetraLinearObjContainer_hpp__
44#define __Panzer_BlockedTpetraLinearObjContainer_hpp__
46#include "PanzerDiscFE_config.hpp"
49#include "Tpetra_Vector.hpp"
50#include "Tpetra_CrsMatrix.hpp"
52#include "Teuchos_RCP.hpp"
58#include "Thyra_PhysicallyBlockedLinearOpBase.hpp"
59#include "Thyra_ProductVectorBase.hpp"
60#include "Thyra_TpetraThyraWrappers.hpp"
62#include <unordered_map>
69template <
typename ScalarT,
typename LocalOrdinalT,
typename GlobalOrdinalT,
typename NodeT=panzer::TpetraNodeType>
76 typedef Tpetra::Map<LocalOrdinalT,GlobalOrdinalT,NodeT>
MapType;
92 inline void set_x(
const Teuchos::RCP<VectorType> & in) {
set_x_th(in); }
98 inline void set_f(
const Teuchos::RCP<VectorType> & in) {
set_f_th(in); }
109 void set_x_th(
const Teuchos::RCP<VectorType> & in) {
x = in; }
115 void set_f_th(
const Teuchos::RCP<VectorType> & in) {
f = in; }
118 void set_A_th(
const Teuchos::RCP<CrsMatrixType> & in) {
A = in; }
119 Teuchos::RCP<CrsMatrixType>
get_A_th()
const {
return A; }
126 Teuchos::RCP<CrsMatrixType>
A;
void set_dxdt_th(const Teuchos::RCP< VectorType > &in)
Teuchos::RCP< VectorType > get_dxdt_th() const
virtual void initialize()
bool checkCompatibility() const
Make sure row and column spaces match up.
Teuchos::RCP< VectorType > x
void set_A_th(const Teuchos::RCP< CrsMatrixType > &in)
Teuchos::RCP< VectorType > get_x_th() const
Teuchos::RCP< CrsMatrixType > get_A_th() const
Teuchos::RCP< CrsMatrixType > get_A() const
Teuchos::RCP< VectorType > get_f() const
Teuchos::RCP< CrsMatrixType > A
void setMapsForBlocks(const std::vector< Teuchos::RCP< const MapType > > &blockMaps)
Teuchos::RCP< const MapType > getMapForBlock(std::size_t i) const
Teuchos::RCP< VectorType > f
void set_A(const Teuchos::RCP< CrsMatrixType > &in)
void set_f(const Teuchos::RCP< VectorType > &in)
void initializeMatrix(ScalarT value)
Put a particular scalar in the matrix.
void set_x(const Teuchos::RCP< VectorType > &in)
Thyra::VectorBase< ScalarT > VectorType
Teuchos::RCP< VectorType > get_x() const
void set_x_th(const Teuchos::RCP< VectorType > &in)
Teuchos::RCP< VectorType > get_dxdt() const
Tpetra::Map< LocalOrdinalT, GlobalOrdinalT, NodeT > MapType
std::vector< Teuchos::RCP< const MapType > > blockMaps_
Teuchos::RCP< VectorType > dxdt
Teuchos::RCP< VectorType > get_f_th() const
void set_f_th(const Teuchos::RCP< VectorType > &in)
void set_dxdt(const Teuchos::RCP< VectorType > &in)
Thyra::LinearOpBase< ScalarT > CrsMatrixType