Fix for broadcasting.

trixi-framework · Oct 4, 2023 · ea115cf · ea115cf
1 parent 2e0f348
commit ea115cf
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/solvers/dgmulti/dg.jl b/src/solvers/dgmulti/dg.jl
@@ -302,7 +302,9 @@ function calc_volume_integral!(du, u, mesh::DGMultiMesh,
     @threaded for e in eachelement(mesh, dg, cache)
         flux_values = local_values_threaded[Threads.threadid()]
         for i in eachdim(mesh)
-            flux_values .= flux.(view(u_values, :, e), i, equations)
+            for j in eachindex(flux_values)
+                flux_values[j] = flux(u_values[j, e], i, equations)
+            end
             for j in eachdim(mesh)
                 apply_to_each_field(mul_by_accum!(weak_differentiation_matrices[j],
                                                   dxidxhatj[i, j][1, e]),