Changeset 1861
- Timestamp:
- 08/25/09 10:00:32 (16 years ago)
- Location:
- issm/trunk/src/c/parallel
- Files:
-
- 2 added
- 2 deleted
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/objectivefunctionC.cpp
r1854 r1861 46 46 int numberofnodes; 47 47 48 /* thermalstatic: */49 int thermalstatic=0;48 /*steadystate: */ 49 int steadystate=0; 50 50 int isstokes=0; 51 DataSet* results_ thermalstatic=NULL;51 DataSet* results_steadystate=NULL; 52 52 int dofs01[2]={0,1}; 53 53 … … 72 72 femmodel->parameters->FindParam((void*)&numberofnodes,"numberofnodes"); 73 73 femmodel->parameters->FindParam((void*)&numberofdofspernode,"numberofdofspernode"); 74 femmodel->parameters->FindParam((void*)& thermalstatic,"thermalstatic");74 femmodel->parameters->FindParam((void*)&steadystate,"steadystate"); 75 75 femmodel->parameters->FindParam((void*)&isstokes,"isstokes"); 76 76 … … 88 88 inputs->Add(control_type,param_g_copy,1,numberofnodes); 89 89 90 if(! thermalstatic){90 if(!steadystate){ 91 91 //Run diagnostic with updated parameters. 92 92 diagnostic_core_nonlinear(&u_g,NULL,NULL,NULL,femmodel,inputs,DiagnosticAnalysisEnum(),sub_analysis_type); … … 94 94 } 95 95 else{ 96 //Run full thermalstaticsolution with updated parameters.97 results_ thermalstatic=new DataSet(ResultsEnum());98 thermalstatic_core(results_thermalstatic,model,inputs);96 //Run full steadystate solution with updated parameters. 97 results_steadystate=new DataSet(ResultsEnum()); 98 steadystate_core(results_steadystate,model,inputs); 99 99 100 100 //get u_g 101 results_ thermalstatic->FindResult(&u_g,"u_g");102 delete results_ thermalstatic;101 results_steadystate->FindResult(&u_g,"u_g"); 102 delete results_steadystate; 103 103 104 104 //extract the correct number of dofs (3 or 4) -
issm/trunk/src/c/parallel/parallel.h
r1854 r1861 18 18 void thermal_core_nonlinear(Vec* ptg,double* pmelting_offset,FemModel* fem,ParameterInputs* inputs,int analysis_type,int sub_analysis_type); 19 19 20 void thermalstatic_core(DataSet* results,Model* model, ParameterInputs* inputs);20 void steadystate_core(DataSet* results,Model* model, ParameterInputs* inputs); 21 21 22 22 void diagnostic_core_nonlinear(Vec* pug,Mat* pK_ff0,Mat* pK_fs0, DataSet* loads, FemModel* fem,ParameterInputs* inputs,int analysis_type,int sub_analysis_type);
Note:
See TracChangeset
for help on using the changeset viewer.