Bugfix - use of wrong default algorithms for surface and volume

JuliaGeometry · Sep 8, 2024 · f45d7bb · f45d7bb
1 parent 248b32a
commit f45d7bb
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/integral_aliases.jl b/src/integral_aliases.jl
@@ -84,7 +84,7 @@ function surfaceintegral(
     Dim = Meshes.paramdim(geometry)
 
     if Dim == 2
-        return integral(f, geometry, GaussKronrod())
+        return integral(f, geometry, HAdaptiveCubature())
     else
         error("Performing a surface integral on a geometry with $Dim parametric dimensions not supported.")
     end
@@ -145,7 +145,7 @@ function volumeintegral(
     Dim = Meshes.paramdim(geometry)
 
     if Dim == 3
-        return integral(f, geometry, GaussKronrod())
+        return integral(f, geometry, HAdaptiveCubature())
     else
         error("Performing a volume integral on a geometry with $Dim parametric dimensions not supported.")
     end