bug fix

Author: carlodev

src/Commons/Equations/StabilizationOperations.jl

@@ -62,12 +62,8 @@ Bazilevs, Y., Calo, V. M., Cottrell, J. A., Hughes, T. J. R., Reali, A., & Scova
 """
 function continuity_stabilization(uu, stab_coeff::TensorStabilization,simcase::SimulationCase)
      @unpack   gg = stab_coeff
-     @sunpack τm_comp = simcase
-     if τm_comp > 0
-        return (uu ⋅ uu) * momentum_stabilization(uu, stab_coeff, simcase) #More SUPG-style
-    elseif τm_comp <0 
-        return 1 / (momentum_stabilization(uu,stab_coeff,simcase) ⋅ gg) #Standard VMS
-    end
+
+     return 1 / (momentum_stabilization(uu,stab_coeff,simcase) ⋅ gg) #Standard VMS
 
 end
 

src/Commons/ParametersDef/StabilizationStruct.jl

@@ -27,7 +27,8 @@ end
 
 @with_kw struct TensorFormulation <: StabilizationFormulation
     r::Int64 = 2 
-    Ci::Vector{Int64} = [4,36]
+    Ci::Vector{Real} = [4,36]
+    @assert (Ci[1]>=0 && Ci[2]>=0) "Ci values must be non-negative"
 end