default integration rules for lineintegral, surfaceintegral, and …

…`volumeintegral`
JuliaGeometry · Sep 26, 2024 · a76adb7 · a76adb7
1 parent 1032805
commit a76adb7
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/integral_aliases.jl b/src/integral_aliases.jl
@@ -17,7 +17,7 @@ Rule types available:
 function lineintegral(
         f::F,
         geometry::G,
-        rule::I;
+        rule::I = GaussKronrod();
         kwargs...
 ) where {F <: Function, G <: Meshes.Geometry, I <: IntegrationRule}
     N = Meshes.paramdim(geometry)
@@ -48,7 +48,7 @@ Algorithm types available:
 function surfaceintegral(
         f::F,
         geometry::G,
-        rule::I;
+        rule::I = HAdaptiveCubature();
         kwargs...
 ) where {F <: Function, G <: Meshes.Geometry, I <: IntegrationRule}
     N = Meshes.paramdim(geometry)
@@ -79,7 +79,7 @@ Algorithm types available:
 function volumeintegral(
         f::F,
         geometry::G,
-        rule::I;
+        rule::I = HAdaptiveCubature();
         kwargs...
 ) where {F <: Function, G <: Meshes.Geometry, I <: IntegrationRule}
     N = Meshes.paramdim(geometry)