compute_cell_contrib(): Advance derivative pointer for each connection.

Failing to do this operation resulted in incorrect matrices in cases
with anisotropic tensors and/or non-cube cells.  The error has been
present since the inception of this implementation.
This commit is contained in:
Bård Skaflestad 2012-05-15 14:34:56 +02:00
parent 8ed2e18e15
commit 83274c4027

View File

@ -606,6 +606,8 @@ compute_cell_contrib(struct UnstructuredGrid *G ,
pimpl->ratio->mat_row[ 0 ] += s * dt * dF1;
pimpl->ratio->mat_row[ off ] += s * dt * dF2;
dv += 2 * np; /* '2' == number of one-sided derivatives. */
}
}
}