Changeset 26649 for issm/trunk-jpl/src/jl/solve/gauss.jl

Timestamp:

11/19/21 11:13:48 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: done with GetSolutionFromInputs

File:

• issm/trunk-jpl/src/jl/solve/gauss.jl (modified) (2 diffs)

issm/trunk-jpl/src/jl/solve/gauss.jl

@@ -23 +23 @@
                 coords3 = [0.333333333333333]
         elseif(order==2)
-                weights = [0.577350269189625, 0.577350269189625, 0.577350269189625]
-                coords1 = [0.666666666666667, 0.166666666666667, 0.166666666666667]
-                coords2 = [0.166666666666667, 0.666666666666667, 0.166666666666667]
-                coords3 = [0.166666666666667, 0.166666666666667, 0.666666666666667]
+                weights = [0.577350269189625; 0.577350269189625; 0.577350269189625]
+                coords1 = [0.666666666666667; 0.166666666666667; 0.166666666666667]
+                coords2 = [0.166666666666667; 0.666666666666667; 0.166666666666667]
+                coords3 = [0.166666666666667; 0.166666666666667; 0.666666666666667]
         elseif(order==3)
                 npoints=4
-                weights = [-0.974278579257493, 0.902109795608790, 0.902109795608790, 0.902109795608790]
-                coords1 = [0.333333333333333, 0.600000000000000, 0.200000000000000, 0.200000000000000]
-                coords2 = [0.333333333333333, 0.200000000000000, 0.600000000000000, 0.200000000000000]
-                coords3 = [0.333333333333333, 0.200000000000000, 0.200000000000000, 0.600000000000000]
+                weights = [-0.974278579257493; 0.902109795608790; 0.902109795608790; 0.902109795608790]
+                coords1 = [ 0.333333333333333; 0.600000000000000; 0.200000000000000; 0.200000000000000]
+                coords2 = [ 0.333333333333333; 0.200000000000000; 0.600000000000000; 0.200000000000000]
+                coords3 = [ 0.333333333333333; 0.200000000000000; 0.200000000000000; 0.600000000000000]
         else
                 error("order ",order," not supported yet");
@@ -39 +39 @@
         return GaussTria(npoints,weights,coords1,coords2,coords3)
 end# }}}
+function GaussTria(finiteelement::IssmEnum) #{{{
 
+        if(finiteelement==P0Enum)
+                npoints = 1
+                weights = [1.]
+                coords1 = [0.333333333333333]
+                coords2 = [0.333333333333333]
+                coords3 = [0.333333333333333]
+        elseif(finiteelement==P1Enum)
+                        npoints = 3
+                        weights = 0.333333333333333*ones(3)
+                        coords1 = [1.; 0.; 0.]
+                        coords2 = [0.; 1.; 0.]
+                        coords3 = [0.; 0.; 1.]
+        else
+                error("finite element ", finiteelement," not supported yet");
+        end
+
+        return GaussTria(npoints,weights,coords1,coords2,coords3)
+end# }}}