Go to the documentation of this file.
19 #if !defined escript_SystemMatrixException_20040608_H
20 #define escript_SystemMatrixException_20040608_H
35 class SystemMatrixException :
public EsysException
#define MATRIX_FORMAT_BLK1
Definition: Paso.h:63
void SparseMatrix_MatrixVector_CSC_OFFSET1(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:134
int mm_read_mtx_crd_size(std::istream &f, int *M, int *N, int *nz)
Definition: mmio.cpp:182
static dim_t N
Definition: SparseMatrix_saveHB.cpp:50
static SystemMatrix_ptr loadMM_toCSC(const char *filename)
Definition: SystemMatrix_loadMM.cpp:213
escript::JMPI mpi_info
Definition: SystemMatrix.h:334
bool is_balanced
Definition: SystemMatrix.h:348
boost::shared_ptr< Distribution > Distribution_ptr
Definition: Distribution.h:37
escript::Distribution_ptr col_distribution
Definition: SystemMatrix.h:333
#define MATRIX_FORMAT_CSC
Definition: Paso.h:62
void RHS_loadMM_toCSR(const char *filename, double *b, dim_t size)
Definition: SystemMatrix_loadMM.cpp:312
void print() const
Definition: SystemMatrix_debug.cpp:119
void mergeMainAndCouple_CSR_OFFSET0(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:78
void mergeMainAndCouple_CSC_OFFSET1(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:298
boost::shared_ptr< SystemMatrix > SystemMatrix_ptr
Definition: SystemMatrix.h:53
boost::shared_ptr< SystemMatrixPattern > SystemMatrixPattern_ptr
Definition: SystemMatrixPattern.h:52
void mergeMainAndCouple_CSR_OFFSET0_Block(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:190
int mm_read_mtx_array_size(std::istream &f, int *M, int *N)
Definition: mmio.cpp:211
#define mm_is_general(typecode)
Definition: mmio.h:62
SystemMatrixException(const std::string &str)
Definition: SystemMatrixException.h:63
void copyRemoteCoupleBlock(bool recreatePattern)
Definition: SystemMatrix_copyRemoteCoupleBlock.cpp:52
Definition: FunctionSpace.h:47
JMPI makeInfo(MPI_Comm comm, bool owncom)
Definition: EsysMPI.cpp:41
void SparseMatrix_MatrixVector_CSR_OFFSET1(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:206
#define MATRIX_FORMAT_OFFSET1
Definition: Paso.h:64
double * finishCollect() const
Definition: SystemMatrix.h:172
void startCollect(const double *in) const
Definition: SystemMatrix.h:167
Definition: Distribution.h:44
#define MPI_INT
Definition: EsysMPI.h:45
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:441
#define MPI_COMM_WORLD
Definition: EsysMPI.h:47
static dim_t nz
Definition: SparseMatrix_saveHB.cpp:50
int comparIndex(const void *index1, const void *index2)
this int-comparison function is used by qsort/bsearch in various places
Definition: PasoUtil.cpp:51
index_t dim_t
Definition: DataTypes.h:90
boost::shared_ptr< JMPI_ > JMPI
Definition: EsysMPI.h:71
SystemMatrixType type
Definition: SystemMatrix.h:322
Coupler_ptr< real_t > col_coupler
Definition: SystemMatrix.h:336
SparseMatrix_ptr mainBlock
main block
Definition: SystemMatrix.h:340
boost::shared_ptr< Coupler< T > > Coupler_ptr
Definition: Coupler.h:55
char MM_typecode[4]
Definition: mmio.h:34
void swap(index_t *r, index_t *c, double *v, int left, int right)
Definition: SparseMatrix.cpp:69
index_t * global_id
stores the global ids for all cols in col_coupleBlock
Definition: SystemMatrix.h:358
Coupler_ptr< real_t > row_coupler
Definition: SystemMatrix.h:337
SparseMatrix_ptr row_coupleBlock
coupling to neighbouring processors (col - row)
Definition: SystemMatrix.h:344
dim_t col_block_size
Definition: SystemMatrix.h:329
virtual ~SystemMatrixException()
Destructor.
Definition: SystemMatrixException.h:66
void extendedRowsForST(dim_t *degree_ST, index_t *offset_ST, index_t *ST)
Definition: SystemMatrix_extendedRows.cpp:54
boost::shared_ptr< SharedComponents > SharedComponents_ptr
Definition: SharedComponents.h:48
EsysException(const std::string &message)
Constructor which creates an Exception with the given message.
Definition: EsysException.h:63
dim_t row_block_size
Definition: SystemMatrix.h:328
SystemMatrixPattern_ptr pattern
Definition: SystemMatrix.h:323
#define MPI_DOUBLE
Definition: EsysMPI.h:46
void fillWithGlobalCoordinates(double f1)
Definition: SystemMatrix_debug.cpp:56
PasoException exception class.
Definition: PasoException.h:45
dim_t getNumCols() const
Definition: SystemMatrix.h:202
SparseMatrix_ptr remote_coupleBlock
coupling of rows-cols on neighbouring processors (may not be valid)
Definition: SystemMatrix.h:346
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:85
#define mm_is_real(typecode)
Definition: mmio.h:57
Definition: AbstractContinuousDomain.cpp:23
int mm_read_banner(std::istream &f, MM_typecode *matcode)
Definition: mmio.cpp:104
void mergeMainAndCouple(index_t **p_ptr, index_t **p_idx, double **p_val) const
Definition: SystemMatrix_mergeMainAndCouple.cpp:62
dim_t getNumRows() const
Definition: SystemMatrix.h:197
void SparseMatrix_MatrixVector_CSR_OFFSET0(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:294
static dim_t M
Definition: SparseMatrix_saveHB.cpp:50
boost::shared_ptr< Connector > Connector_ptr
Definition: Coupler.h:50
boost::shared_ptr< Pattern > Pattern_ptr
Definition: Pattern.h:50
void SparseMatrix_MatrixVector_CSC_OFFSET0(const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
Definition: SparseMatrix_MatrixVector.cpp:58
SystemMatrix()
default constructor - throws exception.
Definition: SystemMatrix.cpp:54
void MatrixVector(double alpha, const double *in, double beta, double *out) const
Definition: SystemMatrix_MatrixVector.cpp:47
static SystemMatrix_ptr loadMM_toCSR(const char *filename)
Definition: SystemMatrix_loadMM.cpp:113
dim_t block_size
Definition: SystemMatrix.h:330
Definition: BiCGStab.cpp:26
void MatrixVector_CSR_OFFSET0(double alpha, const double *in, double beta, double *out) const
Definition: SystemMatrix_MatrixVector.cpp:76
void q_sort(index_t *row, index_t *col, double *val, int begin, int end, int N)
Definition: SparseMatrix.cpp:87
#define MATRIX_FORMAT_DEFAULT
Definition: Paso.h:61
#define MATRIX_FORMAT_DIAGONAL_BLOCK
Definition: Paso.h:65
escript::Distribution_ptr row_distribution
Definition: SystemMatrix.h:332
#define mm_is_array(typecode)
Definition: mmio.h:54
#define mm_is_sparse(typecode)
Definition: mmio.h:51
SparseMatrix_ptr col_coupleBlock
coupling to neighbouring processors (row - col)
Definition: SystemMatrix.h:342