mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
edd470dd7d
There is no need for each file to contain a separate, though trivial, implementation of this feature.
86 lines
1.8 KiB
C
86 lines
1.8 KiB
C
/*
|
|
Copyright 2010 SINTEF ICT, Applied Mathematics.
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#ifndef OPM_SPARSE_SYS_HEADER_INCLUDED
|
|
#define OPM_SPARSE_SYS_HEADER_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);
|
|
|
|
/* ---------------------------------------------------------------------- */
|
|
/* v = zeros([n, 1]) */
|
|
/* ---------------------------------------------------------------------- */
|
|
void
|
|
vector_zero(size_t n, double *v);
|
|
|
|
#ifdef __cplusplus
|
|
}
|
|
#endif
|
|
|
|
#endif /* OPM_SPARSE_SYS_HEADER_INCLUDED */
|