Changeset 26649 for issm/trunk-jpl/src/jl/solve/inputs.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/inputs.jl (modified) (5 diffs)

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

@@ -13 +13 @@
         interp::IssmEnum
         values::Vector{Float64}
+        element_values::Vector{Float64}
 end# }}}
 
@@ -23 +24 @@
 
 #Inputs functions
+function GetInput(inputs::Inputs,enum::IssmEnum) #{{{
+
+        #Does this input exist
+        if !haskey(inputs.lookup,enum)
+                error("Input ",enum," not found")
+        end
+
+        #return input
+        return inputs.lookup[enum]
+
+end#}}}
 function SetInput(inputs::Inputs,enum::IssmEnum,index::Int64,value::Bool) #{{{
 
@@ -48 +60 @@
                 @assert inputs.numberofelements > 0
                 if interp==P0Enum
-                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofelements))
+                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofelements),Vector{Float64}(undef,1))
                 elseif interp==P1Enum
-                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofvertices))
+                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofvertices),Vector{Float64}(undef,3))
                 else
                         error("not supported yet")
@@ -71 +83 @@
                 @assert inputs.numberofvertices>0
                 if interp==P1Enum
-                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofvertices))
+                        inputs.lookup[enum] = ElementInput(enum,interp,zeros(inputs.numberofvertices),Vector{Float64}(undef,3))
                 else
                         error("not supported yet")
@@ -85 +97 @@
         input.values[indices] = values
 end#}}}
+function GetInputValue(input::ElementInput,gauss::GaussTria,i::Int64)
+
+        if input.interp==P0Enum
+                return input.element_value
+        elseif input.interp==P1Enum
+                value = input.element_values[1]*gauss.coords1[i] +  input.element_values[2]*gauss.coords2[i] +  input.element_values[3]*gauss.coords3[i]
+        else
+                error("not implemented yet")
+        end
+
+        return value
+
+end#}}}