fix scale_by_mass_matrix!

JoshuaLampert · Dec 11, 2024 · f017619 · f017619
1 parent 40f52cb
commit f017619
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/equations/equations.jl b/src/equations/equations.jl
@@ -567,7 +567,7 @@ function solve_system_matrix!(dv, system_matrix, rhs,
         if issuccess(factorization)
             scale_by_mass_matrix!(rhs, D1)
             # see https://github.com/JoshuaLampert/DispersiveShallowWater.jl/issues/122
-            dv .= factorization \ rhs
+            dv .= system_matrix \ rhs[2:(end - 1)]
         else
             # The factorization may fail if the time step is too large
             # and h becomes negative.

diff --git a/src/equations/svaerd_kalisch_1d.jl b/src/equations/svaerd_kalisch_1d.jl
@@ -502,7 +502,7 @@ function rhs!(dq, q, t, mesh, equations::SvaerdKalischEquations1D,
     @trixi_timeit timer() "solving elliptic system" begin
         tmp1 .= dv
         solve_system_matrix!((@view dv[2:(end - 1)]), system_matrix,
-                             (@view tmp1[2:(end - 1)]), equations, D1,
+                             tmp1, equations, D1,
                              cache)
         dv[1] = dv[end] = zero(eltype(dv))
     end