diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh index 6d8347e250671f8afaadfc6c49d5cbaf236e7063..634ef5c8c79571dffd4b4be791e7fa8627f9ff70 100644 --- a/dune/gfe/hyperbolichalfspacepoint.hh +++ b/dune/gfe/hyperbolichalfspacepoint.hh @@ -44,7 +44,7 @@ class HyperbolicHalfspacePoint return 8.0/15 - 24*(x-1)/35; } else { T d = x*x-1; - return -6*x/(d*d) + (4*x*x+2)*std::acosh(x)/(std::pow(d,5.2)); + return -6*x/(d*d) + (4*x*x+2)*std::acosh(x)/(std::pow(d,2.5)); } }