Changeset 9206

Timestamp:

08/08/11 15:21:04 (14 years ago)

Author:

Eric.Larour

Message:

New NodalValue module

Location:

issm/trunk/src/c

Files:

• Makefile.am (modified) (2 diffs)
• modules/NodalValuex (added)
• modules/NodalValuex/NodalValuex.cpp (added)
• modules/NodalValuex/NodalValuex.h (added)
• modules/Responsex/Responsex.cpp (modified) (1 diff)
• modules/modules.h (modified) (1 diff)
• objects/Elements/Element.h (modified) (1 diff)
• objects/Elements/Penta.cpp (modified) (1 diff)
• objects/Elements/Penta.h (modified) (1 diff)
• objects/Elements/Tria.cpp (modified) (1 diff)
• objects/Elements/Tria.h (modified) (1 diff)

issm/trunk/src/c/Makefile.am

@@ -494 +494 @@
                                         ./modules/NodesDofx/NodesDofx.h\
                                         ./modules/NodesDofx/NodesDofx.cpp\
+                                        ./modules/NodalValuex/NodalValuex.h\
+                                        ./modules/NodalValuex/NodalValuex.cpp\
                                         ./modules/TriaSearchx/TriaSearchx.h\
                                         ./modules/TriaSearchx/TriaSearchx.cpp\
@@ -1183 +1185 @@
                                         ./modules/NodesDofx/NodesDofx.h\
                                         ./modules/NodesDofx/NodesDofx.cpp\
+                                        ./modules/NodalValuex/NodalValuex.h\
+                                        ./modules/NodalValuex/NodalValuex.cpp\
                                         ./modules/TriaSearchx/TriaSearchx.h\
                                         ./modules/TriaSearchx/TriaSearchx.cpp\

issm/trunk/src/c/modules/Responsex/Responsex.cpp

@@ -42 +42 @@
                 case DragCoefficientAbsGradientEnum:DragCoefficientAbsGradientx(responses, elements,nodes, vertices, loads, materials, parameters,process_units,weight_index); break;
                 case RheologyBbarEnum:RheologyBbarx(responses, elements,nodes, vertices, loads, materials, parameters,process_units); break;
+                case DragCoefficientEnum:NodalValuex(responses, DragCoefficientEnum,elements,nodes, vertices, loads, materials, parameters,process_units); break;
                 default: _error_(" response descriptor \"%s\" not supported yet!",response_descriptor); break;
         }

issm/trunk/src/c/modules/modules.h

@@ -79 +79 @@
 #include "./ModelProcessorx/ModelProcessorx.h"
 #include "./ParsePetscOptionsx/ParsePetscOptionsx.h"
+#include "./NodalValuex/NodalValuex.h"
 #include "./NodeConnectivityx/NodeConnectivityx.h"
 #include "./NodesDofx/NodesDofx.h"

issm/trunk/src/c/objects/Elements/Element.h

@@ -96 +96 @@
                 virtual int    UpdatePotentialSheetUngrounding(double* potential_sheet_ungrounding,Vec vec_nodes_on_iceshelf,double* nodes_on_iceshelf)=0;
                 virtual double RheologyBbarx(void)=0;
+                virtual int    NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units)=0;
 
                 /*Implementation: */

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -6940 +6940 @@
         extern int my_rank;
         return my_rank; 
+}
+/*}}}*/
+/*FUNCTION Penta::NodalValue {{{1*/
+int    Penta::NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units){
+
+        int i;
+        int found=0;
+        double value;
+        Input* data=NULL;
+        GaussPenta* gauss=NULL;
+
+        /*First, serarch the input: */
+        data=inputs->GetInput(natureofdataenum); 
+
+        /*figure out if we have the vertex id: */
+        found=0;
+        for(i=0;i<NUMVERTICES;i++){
+                if(index==nodes[i]->GetVertexId()){
+                        /*Do we have natureofdataenum in our inputs? :*/
+                        if(data){
+                                /*ok, we are good. retrieve value of input at vertex :*/
+                                gauss=new GaussPenta(); gauss->GaussVertex(i);
+                                data->GetParameterValue(&value,gauss);
+                                found=1;
+                                break;
+                        }
+                }
+        }
+
+        if(found)*pvalue=value;
+        return found;
 }
 /*}}}*/

issm/trunk/src/c/objects/Elements/Penta.h

@@ -146 +146 @@
                 double RheologyBbarx(void);
                 void   ViscousHeatingCreateInput(void);
+                int    NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units);
                 /*}}}*/
                 /*Penta specific routines:{{{1*/

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -4525 +4525 @@
 }
 /*}}}*/
+/*FUNCTION Tria::NodalValue {{{1*/
+int    Tria::NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units){
+
+        int i;
+        int found=0;
+        double value;
+        Input* data=NULL;
+        GaussTria *gauss                            = NULL;
+
+        /*First, serarch the input: */
+        data=inputs->GetInput(natureofdataenum); 
+
+        /*figure out if we have the vertex id: */
+        found=0;
+        for(i=0;i<NUMVERTICES;i++){
+                if(index==nodes[i]->GetVertexId()){
+                        /*Do we have natureofdataenum in our inputs? :*/
+                        if(data){
+                                /*ok, we are good. retrieve value of input at vertex :*/
+                                gauss=new GaussTria(); gauss->GaussVertex(i);
+                                data->GetParameterValue(&value,gauss);
+                                found=1;
+                                break;
+                        }
+                }
+        }
+
+        if(found)*pvalue=value;
+        return found;
+}
+/*}}}*/
 /*FUNCTION Tria::PatchFill{{{1*/
 void  Tria::PatchFill(int* prow, Patch* patch){

issm/trunk/src/c/objects/Elements/Tria.h

@@ -148 +148 @@
                 int*   GetHorizontalNeighboorSids(void);
                 double RheologyBbarx(void);
+                int    NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units);
                 /*}}}*/
                 /*Tria specific routines:{{{1*/