Changeset 15849 for issm/trunk-jpl/src/c/modules/SurfaceAverageVelMisfitx/SurfaceAverageVelMisfitx.cpp
- Timestamp:
- 08/20/13 14:50:44 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/SurfaceAverageVelMisfitx/SurfaceAverageVelMisfitx.cpp
r15838 r15849 9 9 #include "../SurfaceAreax/SurfaceAreax.h" 10 10 11 void SurfaceAverageVelMisfitx( IssmDouble* pJ, Elements* elements,Nodes* nodes, Vertices* vertices, Loads* loads, Materials* materials,Parameters* parameters,int weight_index){11 void SurfaceAverageVelMisfitx(IssmDouble* pJ,FemModel* femmodel,int weight_index){ 12 12 13 13 /*Intermediary*/ 14 int i;15 14 Element* element=NULL; 16 15 17 16 /*output: */ 18 IssmDouble J =0;17 IssmDouble J = 0.; 19 18 IssmDouble J_sum; 20 19 21 20 /*Compute surface area and add to elements inputs */ 22 SurfaceAreax(NULL, elements,nodes,vertices, loads,materials,parameters);21 SurfaceAreax(NULL,femmodel); 23 22 24 23 /*Compute Misfit: */ 25 for (i=0;i<elements->Size();i++){26 element=dynamic_cast<Element*>( elements->GetObjectByOffset(i));24 for(int i=0;i<femmodel->elements->Size();i++){ 25 element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i)); 27 26 J+=element->SurfaceAverageVelMisfit(weight_index); 28 27 }
Note:
See TracChangeset
for help on using the changeset viewer.