Scaled wedge 11 (#68)

Added back a scaling factor in the wedge 11. This is reflected by scaling down the precomputed coeffs by 2.
AlgebraicJulia · Feb 21, 2024 · a5553f4 · a5553f4
1 parent a7c3639
commit a5553f4
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/FastDEC.jl b/src/FastDEC.jl
@@ -133,7 +133,7 @@ function dec_p_wedge_product(::Type{Tuple{1,1}}, sd::EmbeddedDeltaDualComplex2D{
 
     shift = ntriangles(sd)
     @inbounds for i in 1:ntriangles(sd)
-        area = areas[i]
+        area = 2 * areas[i]
         coeffs[1, i] = (d_areas[i] + d_areas[i+shift]) / area
         coeffs[2, i] = (d_areas[i+2*shift] + d_areas[i+3*shift]) / area
         coeffs[3, i] = (d_areas[i+4*shift] + d_areas[i+5*shift]) / area

diff --git a/test/Operators.jl b/test/Operators.jl
@@ -188,7 +188,7 @@ end
 
         @test all(dec_wedge_product(Tuple{0, 2}, sd)(V_1, T_2) / 2 .== ∧(Tuple{0, 2}, sd, V_1, T_2))
 
-        @test all(dec_wedge_product(Tuple{1, 1}, sd)(E_1, E_2) / 2 .== ∧(Tuple{1, 1}, sd, E_1, E_2))
+        @test all(dec_wedge_product(Tuple{1, 1}, sd)(E_1, E_2) .== ∧(Tuple{1, 1}, sd, E_1, E_2))
     end
 end
-Original file line number
+Diff line change
@@ Expand Up / @@ -188,7 +188,7 @@ end @@
             @test all(dec_wedge_product(Tuple{0, 2}, sd)(V_1, T_2) / 2 .== ∧(Tuple{0, 2}, sd, V_1, T_2))
-            @test all(dec_wedge_product(Tuple{1, 1}, sd)(E_1, E_2) / 2 .== ∧(Tuple{1, 1}, sd, E_1, E_2))
+            @test all(dec_wedge_product(Tuple{1, 1}, sd)(E_1, E_2) .== ∧(Tuple{1, 1}, sd, E_1, E_2))
         end
     end
@@ Expand Down @@