- Timestamp:
- 12/04/14 10:34:36 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp
r18057 r18930 6 6 7 7 /*Model processing*/ 8 void L2ProjectionEPLAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ 9 10 /*No constraints*/ 11 }/*}}}*/ 12 void L2ProjectionEPLAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/ 13 14 /*No loads*/ 15 }/*}}}*/ 16 void L2ProjectionEPLAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/ 17 /*Now, do we really want DC?*/ 18 int hydrology_model; 19 iomodel->Constant(&hydrology_model,HydrologyModelEnum); 20 if(hydrology_model!=HydrologydcEnum) return; 21 22 /*Do we want an efficient layer*/ 23 bool isefficientlayer; 24 iomodel->Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 25 if(!isefficientlayer) return; 26 27 if(iomodel->domaintype==Domain3DEnum){ 28 iomodel->FetchData(1,MeshVertexonbaseEnum); 29 } 30 else if(iomodel->domaintype==Domain2DverticalEnum){ 31 iomodel->FetchData(1,MeshVertexonbaseEnum); 32 } 33 ::CreateNodes(nodes,iomodel,L2ProjectionEPLAnalysisEnum,P1Enum); 34 iomodel->DeleteData(1,MeshVertexonbaseEnum); 35 }/*}}}*/ 8 36 int L2ProjectionEPLAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/ 9 37 return 1; 10 }/*}}}*/11 void L2ProjectionEPLAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/12 38 }/*}}}*/ 13 39 void L2ProjectionEPLAnalysis::UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/ … … 41 67 } 42 68 }/*}}}*/ 43 void L2ProjectionEPLAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/ 44 /*Now, do we really want DC?*/ 45 int hydrology_model; 46 iomodel->Constant(&hydrology_model,HydrologyModelEnum); 47 if(hydrology_model!=HydrologydcEnum) return; 48 49 /*Do we want an efficient layer*/ 50 bool isefficientlayer; 51 iomodel->Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 52 if(!isefficientlayer) return; 53 54 if(iomodel->domaintype==Domain3DEnum){ 55 iomodel->FetchData(1,MeshVertexonbaseEnum); 56 } 57 else if(iomodel->domaintype==Domain2DverticalEnum){ 58 iomodel->FetchData(1,MeshVertexonbaseEnum); 59 } 60 ::CreateNodes(nodes,iomodel,L2ProjectionEPLAnalysisEnum,P1Enum); 61 iomodel->DeleteData(1,MeshVertexonbaseEnum); 62 }/*}}}*/ 63 void L2ProjectionEPLAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ 64 65 /*No constraints*/ 66 }/*}}}*/ 67 void L2ProjectionEPLAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/ 68 69 /*No loads*/ 69 void L2ProjectionEPLAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/ 70 70 }/*}}}*/ 71 71 … … 226 226 return pe; 227 227 }/*}}}*/ 228 void L2ProjectionEPLAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/228 void L2ProjectionEPLAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/ 229 229 _error_("not implemented yet"); 230 230 }/*}}}*/ 231 void L2ProjectionEPLAnalysis::GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index){/*{{{*/231 void L2ProjectionEPLAnalysis::GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index){/*{{{*/ 232 232 _error_("Not implemented yet"); 233 233 }/*}}}*/ 234 void L2ProjectionEPLAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/234 void L2ProjectionEPLAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/ 235 235 int inputenum,domaintype; 236 236 … … 250 250 } 251 251 }/*}}}*/ 252 void L2ProjectionEPLAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/252 void L2ProjectionEPLAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/ 253 253 /*Default, do nothing*/ 254 254 return;
Note:
See TracChangeset
for help on using the changeset viewer.