Go to the documentation of this file.
28 #ifndef __PASO_SPARSEMATRIX_H__
29 #define __PASO_SPARSEMATRIX_H__
42 struct SparseMatrix : boost::enable_shared_from_this<SparseMatrix>
46 bool patternIsUnrolled);
61 const double* b)
const;
69 const index_t* new_col_index)
const;
77 void saveMM(
const char* filename)
const;
81 return pattern->borrowMainDiagonalPointer();
86 return pattern->borrowColoringPointer();
133 const double* mask_col,
134 double main_diagonal_value);
137 const double* mask_col,
138 double main_diagonal_value);
141 double main_diagonal_value);
144 double main_diagonal_value);
147 double main_diagonal_value);
149 void nullifyRows_CSR(
const double* mask_row,
double main_diagonal_value);
183 const double beta,
double* out);
188 const double beta,
double* out);
193 const double beta,
double* out);
198 const double beta,
double* out);
203 const double beta,
double* out);
214 #endif // __PASO_SPARSEMATRIX_H__
#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
void UMFPACK_free(SparseMatrix *A)
frees any UMFPACK related data from the matrix
Definition: UMFPACK.cpp:46
dim_t getTotalNumRows() const
Definition: SparseMatrix.h:110
int mm_read_mtx_crd_size(std::istream &f, int *M, int *N, int *nz)
Definition: mmio.cpp:181
void nullifyRowsAndCols_CSR_BLK1(const double *mask_row, const double *mask_col, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:69
static dim_t N
Definition: SparseMatrix_saveHB.cpp:48
dim_t numCols
Definition: SparseMatrix.h:175
void addAbsRow_CSR_OFFSET0(double *array) const
Definition: SparseMatrix.cpp:354
int mm_write_mtx_crd_size(std::ostream &f, int M, int N, int nz)
Definition: mmio.cpp:173
dim_t block_size
Definition: SparseMatrix.h:173
void copyBlockFromMainDiagonal(double *out) const
Definition: SparseMatrix.cpp:422
dim_t row_block_size
Definition: SparseMatrix.h:171
#define MATRIX_FORMAT_CSC
Definition: Paso.h:61
static SparseMatrix_ptr loadMM_toCSR(const char *filename)
Definition: SparseMatrix.cpp:213
dim_t getTotalNumCols() const
Definition: SparseMatrix.h:115
Pattern_ptr pattern
Definition: SparseMatrix.h:176
dim_t getNumColors() const
Definition: SparseMatrix.h:100
void * solver_p
pointer to data needed by a solver
Definition: SparseMatrix.h:186
#define mm_set_real(typecode)
Definition: mmio.h:77
void BlockOps_invM_2(double *invA, const double *A, int *failed)
Definition: BlockOps.h:109
void nullifyRowsAndCols_CSR(const double *mask_row, const double *mask_col, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:111
SparseMatrix_ptr getBlock(int blockid) const
Definition: SparseMatrix_getSubmatrix.cpp:89
dim_t maxDeg() const
Definition: SparseMatrix.h:105
SparseMatrix(SparseMatrixType type, Pattern_ptr pattern, dim_t rowBlockSize, dim_t colBlockSize, bool patternIsUnrolled)
Definition: SparseMatrix.cpp:118
#define mm_is_general(typecode)
Definition: mmio.h:61
void setValues(double value)
Definition: SparseMatrix.cpp:487
void BlockOps_invM_N(dim_t N, double *mat, index_t *pivot, int *failed)
LU factorization of NxN matrix mat with partial pivoting.
Definition: BlockOps.h:158
void invMain(double *inv_diag, index_t *pivot) const
Definition: SparseMatrix.cpp:502
#define PASO_MKL
Definition: Options.h:48
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
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
void BlockOps_solveAll(dim_t n_block, dim_t n, double *D, index_t *pivot, double *x)
Definition: BlockOps.h:229
void copy(dim_t N, double *out, const double *in)
out = in
Definition: PasoUtil.h:110
void Preconditioner_LocalSmoother_free(Preconditioner_LocalSmoother *in)
Definition: Smoother.cpp:53
void copyBlockToMainDiagonal(const double *in)
Definition: SparseMatrix.cpp:410
SparseMatrixType type
Definition: SparseMatrix.h:170
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
void saveMM(const char *filename) const
Definition: SparseMatrix.cpp:291
void maxAbsRow_CSR_OFFSET0(double *array) const
Definition: SparseMatrix.cpp:373
boost::shared_ptr< const SparseMatrix > const_SparseMatrix_ptr
Definition: SparseMatrix.h:48
void BlockOps_invM_3(double *invA, const double *A, int *failed)
Definition: BlockOps.h:127
void nullifyRows_CSR_BLK1(const double *mask_row, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:135
double getSparsity() const
Definition: SparseMatrix.h:135
dim_t getNumCols() const
Definition: SparseMatrix.h:125
#define mm_set_matrix(typecode)
Definition: mmio.h:70
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
void nullifyRowsAndCols_CSC_BLK1(const double *mask_row, const double *mask_col, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:51
void MKL_free(SparseMatrix *A)
Definition: MKL.cpp:47
int SparseMatrixType
Definition: SparseMatrix.h:50
void nullifyRows_CSR(const double *mask_row, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:152
void copyToMainDiagonal(const double *in)
Definition: SparseMatrix.cpp:448
void saveHB_CSC(const char *filename) const
Definition: SparseMatrix_saveHB.cpp:190
~SparseMatrix()
Definition: SparseMatrix.cpp:195
#define PASO_PASO
Definition: Options.h:54
SparseMatrix_ptr SparseMatrix_MatrixMatrixTranspose(const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
Definition: SparseMatrix_MatrixMatrixTranspose.cpp:70
dim_t len
Definition: SparseMatrix.h:177
#define PASO_UMFPACK
Definition: Options.h:49
void applyBlockMatrix(double *block_diag, index_t *pivot, double *x, const double *b) const
Definition: SparseMatrix.cpp:552
void copyFromMainDiagonal(double *out) const
Definition: SparseMatrix.cpp:434
index_t solver_package
package controlling the solver pointer
Definition: SparseMatrix.h:183
boost::shared_ptr< SparseMatrix > SparseMatrix_ptr
Definition: SparseMatrix.h:46
void addRow_CSR_OFFSET0(double *array) const
Definition: SparseMatrix.cpp:392
void BlockOps_Cpy_N(dim_t N, double *R, const double *V)
Definition: BlockOps.h:50
PasoException exception class.
Definition: PasoException.h:43
double * val
this is used for classical CSR or CSC
Definition: SparseMatrix.h:180
#define mm_initialize_typecode(typecode)
Definition: mmio.h:90
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
int mm_read_banner(std::istream &f, MM_typecode *matcode)
Definition: mmio.cpp:103
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
SparseMatrix_ptr getTranspose() const
Definition: SparseMatrix.cpp:561
SparseMatrix_ptr getSubmatrix(dim_t n_row_sub, dim_t n_col_sub, const index_t *row_list, const index_t *new_col_index) const
Definition: SparseMatrix_getSubmatrix.cpp:53
void nullifyRowsAndCols_CSC(const double *mask_row, const double *mask_col, double main_diagonal_value)
Definition: SparseMatrix_nullifyRowsAndCols.cpp:87
double getSize() const
Definition: SparseMatrix.h:130
int mm_write_banner(std::ostream &f, MM_typecode matcode)
Definition: mmio.cpp:355
static dim_t M
Definition: SparseMatrix_saveHB.cpp:48
dim_t numRows
Definition: SparseMatrix.h:174
dim_t getNumRows() const
Definition: SparseMatrix.h:120
boost::shared_ptr< Pattern > Pattern_ptr
Definition: Pattern.h:48
index_t * borrowMainDiagonalPointer() const
Definition: SparseMatrix.h:90
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
#define mm_set_coordinate(typecode)
Definition: mmio.h:71
Definition: escriptcore/src/IndexList.h:27
SparseMatrix_ptr SparseMatrix_MatrixMatrix(const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
Definition: SparseMatrix_MatrixMatrix.cpp:61
Definition: BiCGStab.cpp:25
#define PASO_SMOOTHER
Definition: Options.h:72
Definition: SparseMatrix.h:53
SparseMatrix_ptr unroll(SparseMatrixType type) const
Definition: SparseMatrix.cpp:637
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
index_t * borrowColoringPointer() const
Definition: SparseMatrix.h:95
#define MATRIX_FORMAT_DIAGONAL_BLOCK
Definition: Paso.h:64
void applyDiagonal_CSR_OFFSET0(const double *left, const double *right)
Definition: SparseMatrix.cpp:462
dim_t col_block_size
Definition: SparseMatrix.h:172
#define mm_is_sparse(typecode)
Definition: mmio.h:50