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 2acc750d93
commit 20cf0eda36

View File

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