Stability Heatmap

lbmweights/mrt_moments.py

@@ -126,7 +126,7 @@ def get_v37_moments():
     print(calc_weights)
     print(actual_weights)
     assert (np.array(calc_weights) - actual_weights <= 1e-15).all()
-    selected_moments = gram_schmidt(selected_moments, lattice.velocities, lattice.weights)
+    #selected_moments = gram_schmidt(selected_moments, lattice.velocities, lattice.weights)
     M = sp.Matrix(get_moments_of_continuous_maxwellian_equilibrium(moments=selected_moments, dim=2, order=4, c_s_sq=lattice.c_s_sq))
     pprint(m.inv() * M)