Changeset 16137 for issm/trunk/src/c/analyses/objectivefunction.cpp
- Timestamp:
- 09/16/13 09:43:55 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:ignore
-
old new 1 nightlylog 2 configure.sh 1 3 par 2 4 ad
-
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 15397-15401,15403-15487,15489-15701,15704-15735,15737-16076,16082-16133
- Property svn:ignore
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/c
- Property svn:ignore
-
old new 14 14 probe.results 15 15 stXXXX* 16 16 .deps 17 .dirstamp
-
- Property svn:ignore
-
issm/trunk/src/c/analyses
-
Property svn:ignore
set to
.deps
.dirstamp
-
Property svn:ignore
set to
-
issm/trunk/src/c/analyses/objectivefunction.cpp
r15339 r16137 24 24 /*parameters: */ 25 25 int solution_type,analysis_type; 26 bool is stokes= false;26 bool isFS = false; 27 27 bool conserve_loads = true; 28 28 FemModel *femmodel = NULL; … … 32 32 33 33 /*Recover parameters: */ 34 femmodel->parameters->FindParam(&is stokes,FlowequationIsstokesEnum);34 femmodel->parameters->FindParam(&isFS,FlowequationIsFSEnum); 35 35 femmodel->parameters->FindParam(&analysis_type,AnalysisTypeEnum); 36 36 femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum); 37 37 38 38 /*set analysis type to compute velocity: */ 39 if (solution_type==SteadystateSolutionEnum || solution_type== DiagnosticSolutionEnum){40 femmodel->SetCurrentConfiguration( DiagnosticHorizAnalysisEnum);39 if (solution_type==SteadystateSolutionEnum || solution_type==StressbalanceSolutionEnum){ 40 femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum); 41 41 } 42 42 else if (solution_type==BalancethicknessSolutionEnum){ … … 53 53 InputControlUpdatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,search_scalar,false); 54 54 55 /*Run diagnosticwith updated inputs: */55 /*Run stressbalance with updated inputs: */ 56 56 if (solution_type==SteadystateSolutionEnum){ 57 diagnostic_core(femmodel); //We need a 3D velocity!! (vz is required for the next thermal run)57 stressbalance_core(femmodel); //We need a 3D velocity!! (vz is required for the next thermal run) 58 58 } 59 else if (solution_type== DiagnosticSolutionEnum){59 else if (solution_type==StressbalanceSolutionEnum){ 60 60 solutionsequence_nonlinear(femmodel,conserve_loads); 61 61 }
Note:
See TracChangeset
for help on using the changeset viewer.