mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
284683b00d
and non-zeros. Note that this routine is for allocation only. Caller must create and manage sparsity structure if this is being used in a global assembly process. This routine finds utility in MsMFEM BF construction.
61 lines
904 B
C
61 lines
904 B
C
#ifndef SPARSE_SYS_H_INCLUDED
|
|
#define SPARSE_SYS_H_INCLUDED
|
|
|
|
#include <stddef.h>
|
|
|
|
#ifdef __cplusplus
|
|
extern "C" {
|
|
#endif
|
|
|
|
#if defined(MATLAB_MEX_FILE) && MATLAB_MEX_FILE
|
|
#include <mex.h>
|
|
#undef MAT_SIZE_T
|
|
#define MAT_SIZE_T mwSignedIndex
|
|
#endif
|
|
|
|
#ifndef MAT_SIZE_T
|
|
#define MAT_SIZE_T int
|
|
#endif
|
|
|
|
|
|
struct CSRMatrix
|
|
{
|
|
size_t m;
|
|
size_t n;
|
|
size_t nnz;
|
|
|
|
MAT_SIZE_T *ia;
|
|
MAT_SIZE_T *ja;
|
|
|
|
double *sa;
|
|
};
|
|
|
|
|
|
|
|
struct CSRMatrix *
|
|
csrmatrix_new_count_nnz(size_t m);
|
|
|
|
struct CSRMatrix *
|
|
csrmatrix_new_known_nnz(size_t m, size_t nnz);
|
|
|
|
size_t
|
|
csrmatrix_new_elms_pushback(struct CSRMatrix *A);
|
|
|
|
size_t
|
|
csrmatrix_elm_index(size_t i, MAT_SIZE_T j, const struct CSRMatrix *A);
|
|
|
|
void
|
|
csrmatrix_sortrows(struct CSRMatrix *A);
|
|
|
|
void
|
|
csrmatrix_delete(struct CSRMatrix *A);
|
|
|
|
void
|
|
csrmatrix_zero(struct CSRMatrix *A);
|
|
|
|
#ifdef __cplusplus
|
|
}
|
|
#endif
|
|
|
|
#endif /* SPARSE_SYS_H_INCLUDED */
|