diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh index aa77c3ce2a5226a000f75b9bd78af8871960bbfa..307ad7f20a932ef7f49672e42df35e274864b3bb 100644 --- a/dune/gfe/hyperbolichalfspacepoint.hh +++ b/dune/gfe/hyperbolichalfspacepoint.hh @@ -387,7 +387,7 @@ public: } else if (i==N-1 and j==N-1 and k==N-1) { - dFdqdqdq[i][j][k] = (2.0*p[N-1] -1 -4*p[N-1] -3*diffNormSquared)/(p[N-1]*Dune::Power<3>::eval(q[N-1])); + dFdqdqdq[i][j][k] = (-1 -6*p[N-1] + 4*q[N-1] -3*diffNormSquared)/(p[N-1]*Dune::Power<3>::eval(q[N-1])); }