Changeset 11419
- Timestamp:
- 02/13/12 16:14:21 (13 years ago)
- Location:
- issm/trunk-jpl-damage/src/c
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl-damage/src/c/objects/Elements/Penta.cpp ¶
r11417 r11419 769 769 } 770 770 /*}}}*/ 771 <<<<<<< .working772 =======773 771 /*FUNCTION Penta::CreateJacobianMatrix{{{1*/ 774 772 void Penta::CreateJacobianMatrix(Mat Jff){ … … 804 802 } 805 803 /*}}}*/ 806 >>>>>>> .merge-right.r11410807 804 /*FUNCTION Penta::DeepEcho{{{1*/ 808 805 void Penta::DeepEcho(void){ … … 4339 4336 if(converged){ 4340 4337 /*Convert enthalpy into temperature and water fraction*/ 4341 <<<<<<< .working4342 4338 for(i=0;i<numdof;i++) matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]); 4343 =======4344 4339 for(i=0;i<numdof;i++){ 4345 4340 matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]); … … 4347 4342 //if(waterfraction[i]>1) _error_("Water fraction >1 found in solution vector"); 4348 4343 } 4349 >>>>>>> .merge-right.r114104350 4344 4351 4345 this->inputs->AddInput(new PentaP1Input(EnthalpyEnum,values)); -
TabularUnified issm/trunk-jpl-damage/src/c/objects/Elements/Penta.h ¶
r11417 r11419 231 231 ElementMatrix* CreateKMatrixDiagnosticVertVolume(void); 232 232 ElementMatrix* CreateKMatrixDiagnosticVertSurface(void); 233 <<<<<<< .working234 =======235 233 ElementMatrix* CreateJacobianDiagnosticHoriz(void); 236 234 ElementMatrix* CreateJacobianDiagnosticMacayeal2d(void); 237 235 ElementMatrix* CreateJacobianDiagnosticPattyn(void); 238 >>>>>>> .merge-right.r11410239 236 ElementMatrix* CreateJacobianDiagnosticStokes(void); 240 237 void InputUpdateFromSolutionDiagnosticHoriz( double* solutiong); -
TabularUnified issm/trunk-jpl-damage/src/c/objects/Elements/Tria.h ¶
r11417 r11419 84 84 void CreateKMatrix(Mat Kff, Mat Kfs,Vec df); 85 85 void CreatePVector(Vec pf); 86 <<<<<<< .working87 =======88 86 void CreateJacobianMatrix(Mat Jff); 89 >>>>>>> .merge-right.r1141090 87 int GetNodeIndex(Node* node); 91 88 int Sid(); -
TabularUnified issm/trunk-jpl-damage/src/c/objects/KML/KML_File.cpp ¶
r11417 r11419 66 66 KML_Object::DeepEcho(indent); 67 67 68 <<<<<<< .working69 68 /* loop over the kml objects for the file */ 70 69 … … 82 81 _printf_(flag,"%s kmlobj: [empty]\n" ,indent); 83 82 84 =======85 >>>>>>> .merge-right.r1141086 83 return; 87 84 } -
TabularUnified issm/trunk-jpl-damage/src/c/objects/KML/KML_Object.cpp ¶
r11417 r11419 26 26 attrib =new DataSet; 27 27 commnt =new DataSet; 28 <<<<<<< .working29 =======30 28 kmlobj =new DataSet; 31 >>>>>>> .merge-right.r1141032 29 33 30 } … … 44 41 commnt =NULL; 45 42 } 46 <<<<<<< .working47 =======48 43 if (kmlobj) { 49 44 delete kmlobj; 50 45 kmlobj =NULL; 51 46 } 52 >>>>>>> .merge-right.r1141053 47 54 48 } … … 63 57 _printf_(flag," attrib: (size=%d)\n" ,attrib->Size()); 64 58 _printf_(flag," commnt: (size=%d)\n" ,commnt->Size()); 65 <<<<<<< .working66 =======67 59 _printf_(flag," kmlobj: (size=%d)\n" ,kmlobj->Size()); 68 >>>>>>> .merge-right.r1141069 60 70 61 return; … … 106 97 _printf_(flag,"%s commnt: [empty]\n" ,indent); 107 98 108 <<<<<<< .working109 =======110 99 /* loop over the unknown objects for the object */ 111 100 … … 122 111 _printf_(flag,"%s kmlobj: [empty]\n" ,indent); 123 112 124 >>>>>>> .merge-right.r11410125 113 return; 126 114 } … … 134 122 ; 135 123 136 <<<<<<< .working137 =======138 124 memcpy(indent2,indent,(strlen(indent)+1)*sizeof(char)); 139 125 strcat(indent2," "); … … 144 130 } 145 131 146 >>>>>>> .merge-right.r11410147 132 return; 148 133 } … … 273 258 kobj->Read(fid,kstr); 274 259 kmlobj ->AddObject((Object*)kobj); 275 >>>>>>> .merge-right.r11410276 260 } 277 261 -
TabularUnified issm/trunk-jpl-damage/src/c/objects/KML/KML_Object.h ¶
r11417 r11419 21 21 DataSet* attrib; 22 22 DataSet* commnt; 23 <<<<<<< .working24 =======25 23 DataSet* kmlobj; 26 >>>>>>> .merge-right.r1141027 24 28 25 /*KML_Object constructors, destructors {{{1*/ -
TabularUnified issm/trunk-jpl-damage/src/c/objects/KML/KML_Unknown.cpp ¶
r11417 r11419 98 98 void KML_Unknown::Write(FILE* filout,const char* indent){ 99 99 100 <<<<<<< .working101 =======102 100 char* valuei; 103 101 char* vtoken; 104 102 char nl[]={'\n','\0'}; 105 103 106 >>>>>>> .merge-right.r11410107 104 fprintf(filout,"%s<%s",indent,name); 108 105 WriteAttrib(filout," "); … … 110 107 WriteCommnt(filout,indent); 111 108 112 <<<<<<< .working113 109 KML_Object::Write(filout,indent); 114 110 115 111 if (value ) 116 112 fprintf(filout,"%s %s\n",indent,value); 117 =======118 113 if (value ) { 119 114 valuei=(char *) xmalloc((strlen(value)+1)*sizeof(char)); … … 125 120 while (vtoken=strtok(NULL,nl)) 126 121 fprintf(filout,"%s %s\n",indent,vtoken); 127 >>>>>>> .merge-right.r11410128 122 129 <<<<<<< .working130 =======131 123 xfree((void**)&valuei); 132 124 } … … 134 126 KML_Object::Write(filout,indent); 135 127 136 >>>>>>> .merge-right.r11410137 128 fprintf(filout,"%s</%s>\n",indent,name); 138 129 … … 166 157 _error_("KML_Unknown::Read -- Unexpected closing tag %s.\n",kstri); 167 158 168 <<<<<<< .working169 159 else if (strncmp(kstri,"<",1)) 170 160 KMLFileTokenParse( value ,NULL,0, 171 161 NULL, 172 162 fid); 173 =======174 163 else if (strncmp(kstri,"<",1)) { 175 164 if (value) { … … 183 172 } 184 173 } 185 >>>>>>> .merge-right.r11410186 174 187 175 else if (!strncmp(kstri,"<",1)) -
TabularUnified issm/trunk-jpl-damage/src/c/objects/Loads/Icefront.h ¶
r11417 r11419 76 76 void CreateKMatrix(Mat Kff, Mat Kfs); 77 77 void CreatePVector(Vec pf); 78 <<<<<<< .working79 =======80 78 void CreateJacobianMatrix(Mat Jff); 81 >>>>>>> .merge-right.r1141082 79 void PenaltyCreateKMatrix(Mat Kff, Mat kfs, double kmax); 83 80 void PenaltyCreatePVector(Vec pf, double kmax); -
TabularUnified issm/trunk-jpl-damage/src/c/objects/Loads/Penpair.cpp ¶
r11417 r11419 214 214 } 215 215 /*}}}1*/ 216 <<<<<<< .working217 =======218 216 /*FUNCTION Penpair::CreateJacobianMatrix{{{1*/ 219 217 void Penpair::CreateJacobianMatrix(Mat Jff){ … … 221 219 } 222 220 /*}}}1*/ 223 >>>>>>> .merge-right.r11410224 221 /*FUNCTION Penpair::PenaltyCreateKMatrix {{{1*/ 225 222 void Penpair::PenaltyCreateKMatrix(Mat Kff, Mat Kfs,double kmax){ -
TabularUnified issm/trunk-jpl-damage/src/c/solutions/controltao_core.cpp ¶
r11417 r11419 25 25 26 26 /*TAO*/ 27 <<<<<<< .working28 27 int ierr,numberofvertices; 29 28 AppCtx user; … … 32 31 Vec XL = NULL; 33 32 Vec XU = NULL; 34 =======35 33 int ierr; 36 34 int num_controls,solution_type; … … 43 41 Vec XL = NULL; 44 42 Vec XU = NULL; 45 >>>>>>> .merge-right.r1141046 43 47 44 /*Initialize TAO*/ … … 51 48 if(ierr) _error_("Could not initialize Tao"); 52 49 53 <<<<<<< .working54 =======55 50 /*Recover some parameters*/ 56 51 femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum); … … 61 56 maxiter=nsteps*(int)dummy[0]; xfree((void**)&dummy); 62 57 63 >>>>>>> .merge-right.r1141064 58 /*Initialize TAO*/ 65 59 TaoCreate(PETSC_COMM_WORLD,&tao); … … 121 115 VecFree(&gradient); 122 116 CostFunctionx(fcn,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters); 123 <<<<<<< .working124 =======125 117 126 118 /*Compute gradient*/ … … 132 124 xfree((void**)&cost_functions); 133 125 xfree((void**)&cost_functionsd); 134 >>>>>>> .merge-right.r11410135 126 return 0; 136 127 }
Note:
See TracChangeset
for help on using the changeset viewer.