OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Follow file conventions in mult_csr_matrix().

Browse Source
This commit is contained in:
Bård Skaflestad 2012-05-09 19:07:59 +02:00
parent fbb8ae77b1
commit 3c567c409c
1 changed files with 14 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/core/pressure/tpfa/ifs_tpfa.c
View File

@ -11,20 +11,26 @@
#include <opm/core/pressure/tpfa/ifs_tpfa.h>
/* ---------------------------------------------------------------------- */
static void
mult_csr_matrix(const struct CSRMatrix* A,
const double* u,
double* v)
mult_csr_matrix(const struct CSRMatrix *A,
const double *u,
double *v)
/* ---------------------------------------------------------------------- */
{
size_t i, j;
for (j = 0; j < A->m; ++j) {
v[j] = 0;
for (i = (size_t) (A->ia[j]); i < (size_t) (A->ia[j+1]); ++i) {
v[j] += A->sa[i]*u[A->ja[i]];
int j;
size_t i;
for (i = 0, j = 0; i < A->m; i++) {
v[i] = 0.0;
for (; j < A->ia[i + 1]; j++) {
v[i] += A->sa[j] * u[ A->ja[j] ];
}
}
}
struct ifs_tpfa_impl {
double *fgrav; /* Accumulated grav contrib/face */