fix gpu again

Author: LasNikas

src/schemes/boundary/open_boundary/mirroring.jl

@@ -10,7 +10,7 @@ to retrieve first order kernel and particle consistency.
 """
 struct FirstOrderMirroring{ELTYPE}
     firstorder_tolerance::ELTYPE
-    function FirstOrderMirroring(; firstorder_tolerance::Real=1/1000)
+    function FirstOrderMirroring(; firstorder_tolerance::Real=1e-3)
         return new{typeof(firstorder_tolerance)}(firstorder_tolerance)
     end
 end
@@ -27,7 +27,7 @@ the corrected gradient as proposed by [Negi et al. (2022)](@cite Negi2022).
 """
 struct SimpleMirroring{ELTYPE}
     firstorder_tolerance::ELTYPE
-    function SimpleMirroring(; firstorder_tolerance::Real=1/1000)
+    function SimpleMirroring(; firstorder_tolerance::Real=1e-3)
         return new{typeof(firstorder_tolerance)}(firstorder_tolerance)
     end
 end
@@ -63,7 +63,7 @@ end
 
 function BoundaryModelTafuni(;
                              mirror_method=FirstOrderMirroring(;
-                                                               firstorder_tolerance=1/1000))
+                                                               firstorder_tolerance=1e-3))
     return BoundaryModelTafuni(mirror_method)
 end
 

test/examples/gpu.jl

@@ -319,7 +319,9 @@ end
                                                              joinpath(examples_dir(),
                                                                       "fluid",
                                                                       "pipe_flow_2d.jl"),
-                                                             open_boundary_model=BoundaryModelTafuni(),
+                                                             open_boundary_model=BoundaryModelTafuni(;
+                                                                                                     mirror_method=FirstOrderMirroring(;
+                                                                                                                                       firstorder_tolerance=0.001f0)),
                                                              boundary_type_in=BidirectionalFlow(),
                                                              boundary_type_out=BidirectionalFlow(),
                                                              reference_density_in=nothing,
@@ -340,7 +342,9 @@ end
                                                                       "pipe_flow_2d.jl"),
                                                              wcsph=true, sound_speed=20.0f0,
                                                              pressure=0.0f0,
-                                                             open_boundary_model=BoundaryModelTafuni(),
+                                                             open_boundary_model=BoundaryModelTafuni(;
+                                                                                                     mirror_method=FirstOrderMirroring(;
+                                                                                                                                       firstorder_tolerance=0.001f0)),
                                                              boundary_type_in=BidirectionalFlow(),
                                                              boundary_type_out=BidirectionalFlow(),
                                                              reference_density_in=nothing,