Go to the documentation of this file.
18 #if !defined escript_SystemMatrixException_20040608_H
19 #define escript_SystemMatrixException_20040608_H
34 class SystemMatrixException :
public EsysException
#define MATRIX_FORMAT_BLK1
Definition: Paso.h:62
void SparseMatrix_MatrixVector_CSC_OFFSET1(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:132
int mm_read_mtx_crd_size(std::istream &f, int *M, int *N, int *nz)
Definition: mmio.cpp:181
static dim_t N
Definition: SparseMatrix_saveHB.cpp:48
static SystemMatrix_ptr loadMM_toCSC(const char *filename)
Definition: SystemMatrix_loadMM.cpp:211
escript::JMPI mpi_info
Definition: SystemMatrix.h:332
bool is_balanced
Definition: SystemMatrix.h:346
boost::shared_ptr< Distribution > Distribution_ptr
Definition: Distribution.h:35
escript::Distribution_ptr col_distribution
Definition: SystemMatrix.h:331
#define MATRIX_FORMAT_CSC
Definition: Paso.h:61
void RHS_loadMM_toCSR(const char *filename, double *b, dim_t size)
Definition: SystemMatrix_loadMM.cpp:310
void print() const
Definition: SystemMatrix_debug.cpp:117
void mergeMainAndCouple_CSR_OFFSET0(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:76
void mergeMainAndCouple_CSC_OFFSET1(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:296
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:51
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:50
void mergeMainAndCouple_CSR_OFFSET0_Block(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:188
int mm_read_mtx_array_size(std::istream &f, int *M, int *N)
Definition: mmio.cpp:210
#define mm_is_general(typecode)
Definition: mmio.h:61
SystemMatrixException(const std::string &str)
Definition: SystemMatrixException.h:60
void copyRemoteCoupleBlock(bool recreatePattern)
Definition: SystemMatrix_copyRemoteCoupleBlock.cpp:50
Definition: FunctionSpace.h:45
JMPI makeInfo(MPI_Comm comm, bool owncom)
Definition: EsysMPI.cpp:39
void SparseMatrix_MatrixVector_CSR_OFFSET1(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:204
#define MATRIX_FORMAT_OFFSET1
Definition: Paso.h:63
double * finishCollect() const
Definition: SystemMatrix.h:170
void startCollect(const double *in) const
Definition: SystemMatrix.h:165
Definition: Distribution.h:42
#define MPI_INT
Definition: EsysMPI.h:44
void SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:439
#define MPI_COMM_WORLD
Definition: EsysMPI.h:46
static dim_t nz
Definition: SparseMatrix_saveHB.cpp:48
int comparIndex(const void *index1, const void *index2)
this int-comparison function is used by qsort/bsearch in various places
Definition: PasoUtil.cpp:48
index_t dim_t
Definition: DataTypes.h:87
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:70
SystemMatrixType type
Definition: SystemMatrix.h:320
Coupler_ptr< real_t > col_coupler
Definition: SystemMatrix.h:334
SparseMatrix_ptr mainBlock
main block
Definition: SystemMatrix.h:338
boost::shared_ptr< Coupler< T > > Coupler_ptr
Definition: Coupler.h:53
char MM_typecode[4]
Definition: mmio.h:33
void swap(index_t *r, index_t *c, double *v, int left, int right)
Definition: SparseMatrix.cpp:67
index_t * global_id
stores the global ids for all cols in col_coupleBlock
Definition: SystemMatrix.h:356
Coupler_ptr< real_t > row_coupler
Definition: SystemMatrix.h:335
SparseMatrix_ptr row_coupleBlock
coupling to neighbouring processors (col - row)
Definition: SystemMatrix.h:342
dim_t col_block_size
Definition: SystemMatrix.h:327
virtual ~SystemMatrixException()
Destructor.
Definition: SystemMatrixException.h:63
void extendedRowsForST(dim_t *degree_ST, index_t *offset_ST, index_t *ST)
Definition: SystemMatrix_extendedRows.cpp:52
boost::shared_ptr< SharedComponents > SharedComponents_ptr
Definition: SharedComponents.h:46
EsysException(const std::string &message)
Constructor which creates an Exception with the given message.
Definition: EsysException.h:60
dim_t row_block_size
Definition: SystemMatrix.h:326
SystemMatrixPattern_ptr pattern
Definition: SystemMatrix.h:321
#define MPI_DOUBLE
Definition: EsysMPI.h:45
void fillWithGlobalCoordinates(double f1)
Definition: SystemMatrix_debug.cpp:54
PasoException exception class.
Definition: PasoException.h:43
dim_t getNumCols() const
Definition: SystemMatrix.h:200
SparseMatrix_ptr remote_coupleBlock
coupling of rows-cols on neighbouring processors (may not be valid)
Definition: SystemMatrix.h:344
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:82
#define mm_is_real(typecode)
Definition: mmio.h:56
Definition: AbstractContinuousDomain.cpp:22
int mm_read_banner(std::istream &f, MM_typecode *matcode)
Definition: mmio.cpp:103
void mergeMainAndCouple(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:60
dim_t getNumRows() const
Definition: SystemMatrix.h:195
void SparseMatrix_MatrixVector_CSR_OFFSET0(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:292
static dim_t M
Definition: SparseMatrix_saveHB.cpp:48
boost::shared_ptr< Connector > Connector_ptr
Definition: Coupler.h:48
boost::shared_ptr< Pattern > Pattern_ptr
Definition: Pattern.h:48
void SparseMatrix_MatrixVector_CSC_OFFSET0(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:56
SystemMatrix()
default constructor - throws exception.
Definition: SystemMatrix.cpp:52
void MatrixVector(double alpha, const double *in, double beta, double *out) const
Definition: SystemMatrix_MatrixVector.cpp:45
static SystemMatrix_ptr loadMM_toCSR(const char *filename)
Definition: SystemMatrix_loadMM.cpp:111
dim_t block_size
Definition: SystemMatrix.h:328
Definition: BiCGStab.cpp:25
void MatrixVector_CSR_OFFSET0(double alpha, const double *in, double beta, double *out) const
Definition: SystemMatrix_MatrixVector.cpp:74
void q_sort(index_t *row, index_t *col, double *val, int begin, int end, int N)
Definition: SparseMatrix.cpp:85
#define MATRIX_FORMAT_DEFAULT
Definition: Paso.h:60
#define MATRIX_FORMAT_DIAGONAL_BLOCK
Definition: Paso.h:64
escript::Distribution_ptr row_distribution
Definition: SystemMatrix.h:330
#define mm_is_array(typecode)
Definition: mmio.h:53
#define mm_is_sparse(typecode)
Definition: mmio.h:50
SparseMatrix_ptr col_coupleBlock
coupling to neighbouring processors (row - col)
Definition: SystemMatrix.h:340