Added error checks to stdderiv.h
This commit is contained in:
@@ -228,12 +228,28 @@ value(in VectorField uu)
|
|||||||
Preprocessed Vector
|
Preprocessed Vector
|
||||||
gradient(in ScalarField vertex)
|
gradient(in ScalarField vertex)
|
||||||
{
|
{
|
||||||
|
assert(AC_dsx > 0);
|
||||||
|
assert(AC_dsy > 0);
|
||||||
|
assert(AC_dsz > 0);
|
||||||
|
|
||||||
|
assert(AC_inv_dsx > 0);
|
||||||
|
assert(AC_inv_dsy > 0);
|
||||||
|
assert(AC_inv_dsz > 0);
|
||||||
|
|
||||||
return (Vector){derx(vertexIdx, vertex), dery(vertexIdx, vertex), derz(vertexIdx, vertex)};
|
return (Vector){derx(vertexIdx, vertex), dery(vertexIdx, vertex), derz(vertexIdx, vertex)};
|
||||||
}
|
}
|
||||||
|
|
||||||
Preprocessed Matrix
|
Preprocessed Matrix
|
||||||
hessian(in ScalarField vertex)
|
hessian(in ScalarField vertex)
|
||||||
{
|
{
|
||||||
|
assert(AC_dsx > 0);
|
||||||
|
assert(AC_dsy > 0);
|
||||||
|
assert(AC_dsz > 0);
|
||||||
|
|
||||||
|
assert(AC_inv_dsx > 0);
|
||||||
|
assert(AC_inv_dsy > 0);
|
||||||
|
assert(AC_inv_dsz > 0);
|
||||||
|
|
||||||
Matrix mat;
|
Matrix mat;
|
||||||
|
|
||||||
mat.row[0] = (Vector){derxx(vertexIdx, vertex), derxy(vertexIdx, vertex),
|
mat.row[0] = (Vector){derxx(vertexIdx, vertex), derxy(vertexIdx, vertex),
|
||||||
|
Reference in New Issue
Block a user