Changeset 16158
- Timestamp:
- 09/17/13 16:00:44 (12 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/ResetBoundaryConditions.cpp
r16149 r16158 10 10 11 11 /*variables: */ 12 Vector<IssmDouble>* yg = NULL; 13 Nodes *nodes = NULL; 12 Vector<IssmDouble>* yg = NULL; 14 13 15 14 if(VerboseSolution()) _printf0_(" updating boundary conditions...\n"); -
issm/trunk-jpl/src/c/analyses/controltao_core.cpp
r16047 r16158 29 29 AppCtx user; 30 30 TaoSolver tao = 0; 31 IssmDouble *dummy = NULL;32 31 int *control_list = NULL; 33 32 Vector<IssmDouble> *X = NULL; … … 163 162 AppCtx *user = (AppCtx *)userCtx; 164 163 FemModel *femmodel = user->femmodel; 165 Element *element = NULL;166 164 int *responses = NULL; 167 165 -
issm/trunk-jpl/src/c/analyses/dakota_core.cpp
r16127 r16158 161 161 char *dakota_output_file = NULL; 162 162 char *dakota_error_file = NULL; 163 int status = 0;164 163 Dakota::ModelLIter ml_iter; 165 Parameters *parameters = NULL;166 164 167 165 /*Recover dakota_input_file, dakota_output_file and dakota_error_file, in the parameters dataset in parallel */ -
issm/trunk-jpl/src/c/analyses/gia_core.cpp
r15849 r16158 11 11 void gia_core(FemModel* femmodel){ 12 12 13 int i; 14 Vector<IssmDouble>* wg = NULL; 15 Vector<IssmDouble>* dwdtg = NULL; 16 IssmDouble* x = NULL; 17 IssmDouble* y = NULL; 13 Vector<IssmDouble> *wg = NULL; 14 Vector<IssmDouble> *dwdtg = NULL; 15 IssmDouble *x = NULL; 16 IssmDouble *y = NULL; 18 17 19 18 /*parameters: */ -
issm/trunk-jpl/src/c/analyses/stressbalance_core.cpp
r15849 r16158 78 78 } 79 79 80 if (dim==3 & (isSIA || isSSA || isL1L2 || isHO)){80 if (dim==3 && (isSIA || isSSA || isL1L2 || isHO)){ 81 81 if(VerboseSolution()) _printf0_(" computing vertical velocities\n"); 82 82 femmodel->SetCurrentConfiguration(StressbalanceVerticalAnalysisEnum); -
issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp
r15067 r16158 535 535 register int l=0; // level 536 536 register BamgQuadtreeBox * b; 537 Icoor1 h=MaxISize;538 537 Icoor1 hb = MaxISize; 539 538 Icoor1 i0=0,j0=0; … … 549 548 ii[0]=i0; 550 549 jj[0]=j0; 551 h=hb;552 550 do { 553 551 b= pb[l]; … … 560 558 { 561 559 R2 XY(X,b->v[k]->r); 562 double dd; 563 if( (dd= LengthInterpole(Mx(XY), b->v[k]->m(XY))) < seuil ){ 560 if(LengthInterpole(Mx(XY), b->v[k]->m(XY)) < seuil){ 564 561 return b->v[k]; 565 562 } -
issm/trunk-jpl/src/c/bamg/Direction.cpp
r12821 r16158 27 27 int r =1; 28 28 if (dir!= MaxICoor) { 29 Icoor2 x(dir/2),y1(MaxICoor/2-Abs(x)),y( dir%2?-y1:y1);29 Icoor2 x(dir/2),y1(MaxICoor/2-Abs(x)),y((dir%2)?-y1:y1); 30 30 r = (x*i + y*j) >=0; 31 31 } -
issm/trunk-jpl/src/c/bamg/Geometry.cpp
r15104 r16158 62 62 63 63 int verbose; 64 nbv=0;65 nbe=0;66 64 nbcurves=0; 67 65 … … 286 284 int nbreqv=0; 287 285 int nbtan=0; 288 int nbcracked=0;289 286 int i,count; 290 287 … … 378 375 if (nbtan){ 379 376 bamggeom->TangentAtEdges=xNew<double>(4*nbtan); 380 count=0;381 377 for (i=0;i<nbe;i++){ 382 378 if (edges[i].TgA() && edges[i][0].Corner()){ … … 392 388 bamggeom->TangentAtEdges[4*i+3]=edges[i].tg[1].y; 393 389 } 394 count=count+1;395 390 } 396 391 } -
issm/trunk-jpl/src/c/bamg/ListofIntersectionTriangles.cpp
r15100 r16158 48 48 49 49 SegInterpolation* SegI=lSegsI; 50 SegI=lSegsI;51 50 lSegsI[NbSeg].last=Size; 52 51 int EndSeg=Size; … … 225 224 Triangle *tbegin, *t; 226 225 227 long int verbose=2;228 226 Icoor2 deta[3], deti,detj; 229 227 double ba[3]; 230 int nbt =0,ifirst=-1,ilast;228 int ifirst=-1,ilast; 231 229 int i0,i1,i2; 232 230 int ocut,i,j,k=-1; … … 288 286 if (t == tbegin) { // 289 287 double ba,bb; 290 long int verbose=2;291 288 AdjacentTriangle edge=CloseBoundaryEdge(a,t,ba,bb); 292 289 BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); -
issm/trunk-jpl/src/c/bamg/Mesh.cpp
r16156 r16158 269 269 void Mesh::ReadMesh(int* index,double* x,double* y,int nods,int nels){ 270 270 271 double Hmin = HUGE_VAL;// the infinie value272 271 long i1,i2,i3; 273 272 long i; -
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r16151 r16158 43 43 Penta::Penta(int penta_id, int penta_sid, int index, IoModel* iomodel,int nummodels) 44 44 :PentaRef(nummodels) 45 ,ElementHook(nummodels,index+1,6,iomodel) //index+1: material id, iomodel->numberofelements+1: matpar id 46 { //i is the element index 47 48 int i; 45 ,ElementHook(nummodels,index+1,6,iomodel){ 46 49 47 int penta_elements_ids[2]; 50 48 51 49 /*Checks in debugging mode*/ 52 /*{{{*/53 50 _assert_(iomodel->Data(MeshUpperelementsEnum)); 54 51 _assert_(iomodel->Data(MeshLowerelementsEnum)); 55 /*}}}*/56 52 57 53 /*id: */ … … 4544 4540 4545 4541 /*Intermediaries*/ 4546 int i, j,found=0;4542 int i,found=0; 4547 4543 int stabilization; 4548 4544 IssmDouble Jdet,phi,dt; … … 6526 6522 int* pdoflist=NULL; 6527 6523 IssmDouble FSreconditioning; 6528 GaussPenta *gauss;6529 6524 6530 6525 /*Fetch number of nodes and dof for this finite element*/ … … 8586 8581 8587 8582 /*Intermediaries */ 8588 int i,j;8589 8583 int approximation; 8590 8584 IssmDouble viscosity,Jdet; … … 8632 8626 material->GetViscosity3dFS(&viscosity,&epsilon[0]); 8633 8627 8634 for(i =0;i<NUMVERTICES;i++){8628 for(int i=0;i<NUMVERTICES;i++){ 8635 8629 pe->values[i*NDOF3+0]+=-Jdet*gauss->weight*viscosity*dw[0]*dbasis[2][i]; 8636 8630 pe->values[i*NDOF3+1]+=-Jdet*gauss->weight*viscosity*dw[1]*dbasis[2][i]; -
issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp
r15872 r16158 846 846 * where h is the interpolation function for node i. 847 847 */ 848 int i; 849 int num_dof=4; 850 848 849 int num_dof=4; 851 850 IssmDouble L1L2l3[NUMNODESP1_2d]; 852 851 IssmDouble dbasis[3][NUMNODESP1]; -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r16153 r16158 24 24 /*FUNCTION Tria::Tria(){{{*/ 25 25 Tria::Tria(){ 26 27 int i;28 29 26 this->nodes = NULL; 30 27 this->vertices = NULL; -
issm/trunk-jpl/src/c/toolkits/issm/Bucket.h
r15838 r16158 131 131 void SpawnBucketsPerCpu(DataSet* bucketsofcpu_i,int rank_i,int* rowranks){ /*{{{*/ 132 132 133 int i,j;134 135 133 /*go through our idxm index of rows this bucket owns, and spawn buckets 136 134 *if these rows belong to cpu rank_i. Use rowranks to determine this.*/ 137 for(i =0;i<m;i++){135 for(int i=0;i<m;i++){ 138 136 if (rowranks[idxm[i]]==rank_i){ 139 137 /*This row belongs to cpu rank_i, so spawn a bucket with this row, and add it to the bucketsofcpu_i dataset: */ … … 147 145 } 148 146 149 }; 150 /*}}}*/ 147 }; /*}}}*/ 151 148 int BucketType(void){ /*{{{*/ 152 149 -
issm/trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h
r15838 r16158 149 149 150 150 int *row_indices_forcpu = NULL; 151 int *col_indices_forcpu = NULL;152 151 int *modes_forcpu = NULL; 153 152 doubletype *values_forcpu = NULL;
Note:
See TracChangeset
for help on using the changeset viewer.