18 #ifndef __PASO_TRANSPORT_H__ 19 #define __PASO_TRANSPORT_H__ 21 #define DT_FACTOR_MAX 100000. 40 void solve(
double* u,
double dt,
double* u0,
double* q,
Options* options);
80 return mass_matrix->mainBlock->borrowMainDiagonalPointer();
114 #endif // __PASO_TRANSPORT_H__ boost::shared_ptr< TransportProblem > TransportProblem_ptr
Definition: Transport.h:28
#define PASO_DLL_API
Definition: Paso.h:41
static index_t getTypeId(index_t solver, index_t preconditioner, index_t package, bool symmetry, const esysUtils::JMPI &mpi_info)
Definition: Transport.h:83
TransportProblem(SystemMatrixPattern_ptr pattern, int block_size)
Definition: Transport.cpp:35
SystemMatrix_ptr transport_matrix
Definition: Transport.h:90
double * constraint_mask
Definition: Transport.h:99
Definition: Transport.h:33
Connector_ptr borrowConnector() const
Definition: Transport.h:73
void reset()
Definition: Transport.cpp:82
dim_t getBlockSize() const
Definition: Transport.h:48
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:38
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:38
esysUtils::JMPI mpi_info
Definition: Transport.h:109
double dt_max_T
safe time step size for transport part
Definition: Transport.h:98
double * reactive_matrix
Definition: Transport.h:106
SystemMatrix_ptr borrowTransportMatrix() const
Definition: Transport.h:53
double * main_diagonal_low_order_transport_matrix
Definition: Transport.h:101
index_t * borrowMainDiagonalPointer() const
Definition: Transport.h:78
bool valid_matrices
Definition: Transport.h:94
boost::shared_ptr< const TransportProblem > const_TransportProblem_ptr
Definition: Transport.h:30
SystemMatrix_ptr mass_matrix
Definition: Transport.h:91
dim_t getTotalNumRows() const
Definition: Transport.h:68
#define MATRIX_FORMAT_DEFAULT
Definition: Paso.h:53
double * lumped_mass_matrix
Definition: Transport.h:105
double * borrowLumpedMassMatrix() const
Definition: Transport.h:63
int index_t
Definition: types.h:24
double * main_diagonal_mass_matrix
Definition: Transport.h:107
SystemMatrix_ptr borrowMassMatrix() const
Definition: Transport.h:58
void setUpConstraint(const double *q)
Definition: Transport.cpp:92
SystemMatrix_ptr iteration_matrix
Definition: Transport.h:92
boost::shared_ptr< Connector > Connector_ptr
Definition: Coupler.h:35
~TransportProblem()
Definition: Transport.cpp:73
index_t dim_t
Definition: types.h:27
#define MATRIX_FORMAT_BLK1
Definition: Paso.h:55
double getSafeTimeStepSize()
Definition: Transport_solve.cpp:199
void insertConstraint(const double *r, double *source)
Definition: Transport.cpp:111
boost::shared_ptr< JMPI_ > JMPI
Definition: Esys_MPI.h:79
double dt_max_R
safe time step size for reactive part
Definition: Transport.h:96
void solve(double *u, double dt, double *u0, double *q, Options *options)
Definition: Transport_solve.cpp:49