stability False for P2

python/fluid.py

@@ -157,6 +157,8 @@ class FluidProblem :
         self.order_fields = order
         self.order_porosity = [1]
         self.order_concentration = self.order_porosity
+        if np.max(order) > 1:
+            stability = False
         if dim == 2 :
             self._lib = lib2
         elif dim == 3 :