Only sum inside pressure contribution

WaterLily-jl · Dec 16, 2023 · bd88748 · bd88748
1 parent 7cd9fcc
commit bd88748
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/Metrics.jl b/src/Metrics.jl
@@ -82,7 +82,7 @@ Surface normal integral of field `p` over the `body`.
 """
 function ∮nds(p::AbstractArray{T,N},df::AbstractArray{T},body::AbstractBody,t=0) where {T,N}
     @loop df[I,:] = p[I]*nds(body,loc(0,I,T),t) over I ∈ inside(p)
-    reshape(sum(df,dims=1:N),N) |> Array
+    [sum(@inbounds(df[inside(p),i])) for i ∈ 1:N] |> Array
 end
 @inline function nds(body::AbstractBody,x,t)
     d,n,_ = measure(body,x,t)

diff --git a/test/runtests.jl b/test/runtests.jl
@@ -203,7 +203,7 @@ using StaticArrays
 function get_flow(N,f)
     a = Flow((N,N),(1.,0.);f,T=Float32)
     @inside a.p[I] = loc(0, I)[2]
-    sdf(x,t) = √sum(abs2,x.-(N÷2+0.5))-N÷4
+    sdf(x,t) = √sum(abs2,x.-(N/2))-N÷4
     map(x,t) = x.-SVector(t,0)
     body = AutoBody(sdf,map)
     WaterLily.measure!(a,body)
@@ -264,7 +264,7 @@ import WaterLily: ×
         N = 32
         a,body = get_flow(N,f)
         force = WaterLily.∮nds(a.p,a.V,body)
-        @test sum(abs,force/(π*(N÷4)^2) - [0,1]) < 1e-2
+        @test sum(abs,force/(π*(N/4)^2) - [0,1]) < 2e-3
     end
 end