Changeset 3570


Ignore:
Timestamp:
04/19/10 16:26:55 (15 years ago)
Author:
Mathieu Morlighem
Message:

removed exprintf of all ISSMERROR (not needed anymore)
Fixed Discontinuous Galerkin (missing thickness in inflow BC)"

Location:
issm/trunk/src/c
Files:
62 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/Bamgx/objects/Edge.h

    r3566 r3570  
    3939                        int Intersection(const  Edge & e) const {
    4040                                if (!(adj[0]==&e || adj[1]==&e)){
    41                                         ISSMERROR(exprintf("Intersection bug"));
     41                                        ISSMERROR("Intersection bug");
    4242                                }
    4343                                if (adj[0]!=&e && adj[1]!=&e){
    44                                         ISSMERROR(exprintf("adj[0]!=&e && adj[1]!=&e"));
     44                                        ISSMERROR("adj[0]!=&e && adj[1]!=&e");
    4545                                }
    4646                                return adj[0]==&e ? 0 : 1;
  • issm/trunk/src/c/Bamgx/objects/GeometricalEdge.cpp

    r3335 r3570  
    9393                double ca,cb,cta,ctb;
    9494                if ( theta<-1e-12){
    95                         ISSMERROR(exprintf("theta<-1e-12"));
     95                        ISSMERROR("theta<-1e-12");
    9696                }
    9797                if ( theta>1+1e-12){
    98                         ISSMERROR(exprintf("theta>1+1e-12"));
     98                        ISSMERROR("theta>1+1e-12");
    9999                }
    100100                if (TgA())
  • issm/trunk/src/c/Bamgx/objects/Geometry.cpp

    r3407 r3570  
    3333                 subdomains[i].Set(Gh.subdomains[i],Gh,*this);
    3434                if (nbt);  {
    35                         ISSMERROR(exprintf("nbt"));
     35                        ISSMERROR("nbt");
    3636                }
    3737        }
     
    4444
    4545                if (NbRef>0){
    46                         ISSMERROR(exprintf("NbRef>0"));
     46                        ISSMERROR("NbRef>0");
    4747                }
    4848                if(vertices)  delete [] vertices;vertices=0;
     
    8383                }
    8484                if (bamggeom->Edges==NULL){
    85                         ISSMERROR(exprintf("the domain provided does not contain any edge"));
     85                        ISSMERROR("the domain provided does not contain any edge");
    8686                }
    8787
     
    8989                if (bamggeom->Vertices){
    9090                        if(verbose>5) printf("      processing Vertices\n");
    91                         if (bamggeom->VerticesSize[1]!=3) ISSMERROR(exprintf("Vertices should have 3 columns"));
     91                        if (bamggeom->VerticesSize[1]!=3) ISSMERROR("Vertices should have 3 columns");
    9292                        vertices = new GeometricalVertex[nbvx];
    9393                        for (i=0;i<nbv;i++) {
     
    116116                        coefIcoor= (MaxICoor)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
    117117                        if(coefIcoor <=0){
    118                                 ISSMERROR(exprintf("coefIcoor should be positive"));
     118                                ISSMERROR("coefIcoor should be positive");
    119119                        }
    120120                }
    121121                else{
    122                         ISSMERROR(exprintf("No MeshVertex provided"));
     122                        ISSMERROR("No MeshVertex provided");
    123123                }
    124124
     
    129129
    130130                        if(verbose>5) printf("      processing Edges\n");
    131                         if (bamggeom->EdgesSize[1]!=3) ISSMERROR(exprintf("Edges should have 3 columns"));
     131                        if (bamggeom->EdgesSize[1]!=3) ISSMERROR("Edges should have 3 columns");
    132132                        edges = new GeometricalEdge[nbe];
    133133
     
    180180                }
    181181                else{
    182                         ISSMERROR(exprintf("No edges provided"));
     182                        ISSMERROR("No edges provided");
    183183                }
    184184
     
    223223                if (bamggeom->TangentAtEdges){
    224224                        if(verbose>5) printf("      processing TangentAtEdges");
    225                         if (bamggeom->TangentAtEdgesSize[1]!=4) ISSMERROR(exprintf("TangentAtEdges should have 4 columns"));
     225                        if (bamggeom->TangentAtEdgesSize[1]!=4) ISSMERROR("TangentAtEdges should have 4 columns");
    226226                        int n,i,j,k;
    227227                        R2 tg;
     
    234234                                tg.y=bamggeom->TangentAtEdges[k*4+3];
    235235                                if (j!=0 && j!=1){
    236                                         ISSMERROR(exprintf("TangentAtEdges second index should be 1 or 2 only"));
     236                                        ISSMERROR("TangentAtEdges second index should be 1 or 2 only");
    237237                                }
    238238                                edges[i].tg[j] = tg;
     
    243243                if(bamggeom->Corners){
    244244                        if(verbose>5) printf("      processing Corners");
    245                         if (bamggeom->CornersSize[1]!=1) ISSMERROR(exprintf("Corners should have 1 column"));
     245                        if (bamggeom->CornersSize[1]!=1) ISSMERROR("Corners should have 1 column");
    246246                        n=bamggeom->CornersSize[0];
    247247                        for (i=0;i<n;i++) {     
    248248                                j=(int)bamggeom->Corners[i]-1; //for C indexing
    249249                                if (j>nbv-1 || j<0){
    250                                         ISSMERROR(exprintf("Bad corner definition: should in [0 %i]",nbv));
     250                                        ISSMERROR("Bad corner definition: should in [0 %i]",nbv);
    251251                                }
    252252                                vertices[j].SetCorner();
     
    257257                if(bamggeom->RequiredVertices){
    258258                        if(verbose>5) printf("      processing RequiredVertices\n");
    259                         if (bamggeom->RequiredVerticesSize[1]!=1) ISSMERROR(exprintf("RequiredVertices should have 1 column"));
     259                        if (bamggeom->RequiredVerticesSize[1]!=1) ISSMERROR("RequiredVertices should have 1 column");
    260260                        n=bamggeom->RequiredVerticesSize[0];
    261261                        for (i=0;i<n;i++) {     
    262262                                j=(int)bamggeom->RequiredVertices[i]-1; //for C indexing
    263263                                if (j>nbv-1 || j<0){
    264                                         ISSMERROR(exprintf("Bad RequiredVerticess  definition: should in [0 %i]",nbv));
     264                                        ISSMERROR("Bad RequiredVerticess  definition: should in [0 %i]",nbv);
    265265                                }
    266266                                vertices[j].SetRequired();  }
     
    270270                if(bamggeom->RequiredEdges){
    271271                        if(verbose>5) printf("      processing RequiredEdges\n");
    272                         if (bamggeom->RequiredEdgesSize[1]!=1) ISSMERROR(exprintf("RequiredEdges should have 1 column"));
     272                        if (bamggeom->RequiredEdgesSize[1]!=1) ISSMERROR("RequiredEdges should have 1 column");
    273273                        n=bamggeom->RequiredEdgesSize[0];
    274274                        for (i=0;i<n;i++) {     
    275275                                j=(int)bamggeom->RequiredEdges[i]-1; //for C indexing
    276276                                if (j>nbe-1 || j<0){
    277                                         ISSMERROR(exprintf("Bad RequiredEdges definition: should in [0 %i]",nbe));
     277                                        ISSMERROR("Bad RequiredEdges definition: should in [0 %i]",nbe);
    278278                                }
    279279                                edges[j].SetRequired(); 
     
    285285                        long i0,i1,i2,i3;
    286286                        if(verbose>5) printf("      processing SubDomains\n");
    287                         if (bamggeom->SubDomainsSize[1]!=4) ISSMERROR(exprintf("SubDomains should have 4 columns"));
     287                        if (bamggeom->SubDomainsSize[1]!=4) ISSMERROR("SubDomains should have 4 columns");
    288288                        NbSubDomains=bamggeom->SubDomainsSize[0];
    289289                        subdomains = new GeometricalSubDomain[NbSubDomains];
     
    293293                                i2=(int)bamggeom->SubDomains[i*4+2];
    294294                                i3=(int)bamggeom->SubDomains[i*4+3];
    295                                 if (i0!=2) ISSMERROR(exprintf("Bad Subdomain definition: first number should be 2 (for Edges)"));
    296                                 if (i1>nbe || i1<=0) ISSMERROR(exprintf("Bad Subdomain definition: second number should in [1 %i] (edge number)",nbe));
     295                                if (i0!=2) ISSMERROR("Bad Subdomain definition: first number should be 2 (for Edges)");
     296                                if (i1>nbe || i1<=0) ISSMERROR("Bad Subdomain definition: second number should in [1 %i] (edge number)",nbe);
    297297                                subdomains[i].edge=edges + (i1-1);
    298298                                subdomains[i].sens = (int) i2;
     
    466466                                //check that the clostest vertex is not itself...
    467467                                if ( v !=  &(MeshVertex &) vertices[j]){
    468                                         ISSMERROR(exprintf(" v !=  &(MeshVertex &) vertices[j]"));
     468                                        ISSMERROR(" v !=  &(MeshVertex &) vertices[j]");
    469469                                }
    470470                                vertices[i].link = vertices + j;
     
    484484                                if (!vertices[i].IsThe()) printf("  %i and %i\n",i,Number(vertices[i].The()));
    485485                        }
    486                         ISSMERROR(exprintf("See above"));
     486                        ISSMERROR("See above");
    487487                }
    488488
     
    555555                        //check that its length is not 0
    556556                        if(lv10==0) {
    557                                 ISSMERROR(exprintf("Length of edge %i is 0",i));
     557                                ISSMERROR("Length of edge %i is 0",i);
    558558                        }
    559559                        //compute angle in [-Pi Pi]
     
    647647                                long i1 = n1/2 ,j1=n1%2;
    648648                                if( edges[i1].v[j1] != edges[i].v[j]) {
    649                                         ISSMERROR(exprintf("Bug Adj edge"));
     649                                        ISSMERROR("Bug Adj edge");
    650650                                }
    651651                                edges[i1].Adj[j1] = edges + i;
     
    763763                        pon = on;
    764764                        if (k>=100){
    765                                 ISSMERROR(exprintf("k>=100"));
     765                                ISSMERROR("k>=100");
    766766                        }
    767767                        R2 A= (*on)[0];
     
    842842                GeometricalEdge* on=e.onGeometry;
    843843                if (!on){
    844                         ISSMERROR(exprintf("ProjectOnCurve error message: edge provided should be on geometry"));
     844                        ISSMERROR("ProjectOnCurve error message: edge provided should be on geometry");
    845845                }
    846846                if (!e[0].onGeometry ||  !e[1].onGeometry){
    847                         ISSMERROR(exprintf("ProjectOnCurve error message: at least one of the vertex of the edge provided is not on geometry"));
     847                        ISSMERROR("ProjectOnCurve error message: at least one of the vertex of the edge provided is not on geometry");
    848848                }
    849849
     
    886886                                        printf(" 2)  code bug : be sure that we call   Triangles::SetVertexFieldOn() before\n");
    887887                                        printf("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)\n");
    888                                         ISSMERROR(exprintf("see above"));
     888                                        ISSMERROR("see above");
    889889                                  }
    890890                                NbTry++;
     
    906906                                printf(" 2)  code bug : be sure that we call   Triangles::SetVertexFieldOn() before\n");
    907907                                printf("To solve the problem do a coarsening of the geometrical mesh or change the constant value of mxe (dangerous)\n");
    908                                 ISSMERROR(exprintf("see above"));
     908                                ISSMERROR("see above");
    909909                        }
    910910                        GeometricalEdge* tmpge = eg1;
     
    940940                                // search the geometrical edge
    941941                                if (s>1.0){
    942                                         ISSMERROR(exprintf("s>1.0"));
     942                                        ISSMERROR("s>1.0");
    943943                                }
    944944                                double ls= s*ll;
     
    950950                                while (  (l1=lge[i]) < ls ) {
    951951                                        if (i<0 || i>mxe){
    952                                                 ISSMERROR(exprintf("i<0 || i>mxe"));
     952                                                ISSMERROR("i<0 || i>mxe");
    953953                                        }
    954954                                        i++,s0=1-(s1=sensge[i]),l0=l1;
  • issm/trunk/src/c/Bamgx/objects/ListofIntersectionTriangles.cpp

    r3407 r3570  
    6161                                k=(*t)(0) ?  ((  (*t)(1) ? ( (*t)(2) ? -1 : 2) : 1  )) : 0;
    6262                                if (k<0){
    63                                         ISSMERROR(exprintf("k<0"));
     63                                        ISSMERROR("k<0");
    6464                                }
    6565                                ocut = OppositeEdge[k];
     
    105105                                double dij = detj-deti;
    106106                                if (i+j+k != 0 + 1 +2){
    107                                         ISSMERROR(exprintf("i+j+k != 0 + 1 +2"));
     107                                        ISSMERROR("i+j+k != 0 + 1 +2");
    108108                                }
    109109                                ba[j] =  detj/dij;
     
    177177                                        }
    178178                                        else {
    179                                                 ISSMERROR(exprintf("Bug Split Edge"));
     179                                                ISSMERROR("Bug Split Edge");
    180180                                        }
    181181                                }
     
    267267                // check Size
    268268                if (Size<=0){
    269                         ISSMERROR(exprintf("Size<=0"));
     269                        ISSMERROR("Size<=0");
    270270                }
    271271
  • issm/trunk/src/c/Bamgx/objects/ListofIntersectionTriangles.h

    r3407 r3570  
    4343                                        double c01=lEnd-lBegin, c0=(lEnd-s)/c01, c1=(s-lBegin)/c01;
    4444                                        if (lBegin>s || s>lEnd){
    45                                                 ISSMERROR(exprintf("lBegin>s || s>lEnd"));
     45                                                ISSMERROR("lBegin>s || s>lEnd");
    4646                                        }
    4747                                        return e->F(sBegin*c0+sEnd*c1);
     
    9393                                        SegInterpolation * lEn =  new SegInterpolation[MaxNbSeg];
    9494                                        if (!lSegsI || NbSeg>=MaxNbSeg){
    95                                                 ISSMERROR(exprintf("!lSegsI || NbSeg>=MaxNbSeg"));
     95                                                ISSMERROR("!lSegsI || NbSeg>=MaxNbSeg");
    9696                                        }
    9797                                        for (int i=0;i< NbSeg;i++)
  • issm/trunk/src/c/Bamgx/objects/MeshVertex.cpp

    r3407 r3570  
    3737                        j = NextEdge[jc];
    3838                        if (k>=2000){
    39                                 ISSMERROR(exprintf("k>=2000 (Maximum number of iterations reached)"));
     39                                ISSMERROR("k>=2000 (Maximum number of iterations reached)");
    4040                        }
    4141                } while ( tbegin != tria);
     
    104104                                j = NextEdge[jc];
    105105                                if (k>=2000){
    106                                         ISSMERROR(exprintf("k>=2000"));
     106                                        ISSMERROR("k>=2000");
    107107                                }
    108108                        }while ( tbegin != tria);
     
    163163                }
    164164                else{
    165                         ISSMERROR(exprintf("Metrictype %i not supported yet (use 0,1 or 2(default))",Metrictype));
     165                        ISSMERROR("Metrictype %i not supported yet (use 0,1 or 2(default))",Metrictype);
    166166                }
    167167
  • issm/trunk/src/c/Bamgx/objects/Metric.cpp

    r3332 r3570  
    177177                // warning for optimisation S is in [0:0.5] not in [0:1]
    178178                if (i>=512){
    179                         ISSMERROR(exprintf("i>=512"));
     179                        ISSMERROR("i>=512");
    180180                }
    181181                LastMetricInterpole.lab=l;
     
    321321                  }
    322322                if (r>1 || r<0){
    323                         ISSMERROR(exprintf("r>1 || r<0"));
     323                        ISSMERROR("r>1 || r<0");
    324324                }
    325325                return r ;
  • issm/trunk/src/c/Bamgx/objects/QuadTree.cpp

    r3407 r3570  
    8888         root=NewQuadTreeBox();
    8989         if ( MaxISize <= MaxICoor){
    90                  ISSMERROR(exprintf("MaxISize <= MaxICoor"));
     90                 ISSMERROR("MaxISize <= MaxICoor");
    9191         }
    9292         for (int i=0;i<nbv;i++)
     
    156156                //check that l is not 0 (this should not happen as MaxDeep = 30)
    157157                if (level==0){
    158                         ISSMERROR(exprintf("level==0 cannot be true as it has been initialized as MaxISize = %i",MaxISize));
     158                        ISSMERROR("level==0 cannot be true as it has been initialized as MaxISize = %i",MaxISize);
    159159                }
    160160
     
    468468                be = b +ll;
    469469                if (!b){
    470                         ISSMERROR(exprintf("!b"));
     470                        ISSMERROR("!b");
    471471                }
    472472        }
  • issm/trunk/src/c/Bamgx/objects/QuadTree.h

    r3407 r3570  
    6262                        QuadTreeBox* NewQuadTreeBox(){
    6363                                if(! (sb->bc<sb->be)) sb=new StorageQuadTreeBox(lenStorageQuadTreeBox,sb);
    64                                 if (!sb || (sb->bc->n != 0)){ISSMERROR(exprintf("!sb || (sb->bc->n != 0)"));}
     64                                if (!sb || (sb->bc->n != 0)){ISSMERROR("!sb || (sb->bc->n != 0)");}
    6565                                NbQuadTreeBox++;
    6666                                return sb->bc++;
  • issm/trunk/src/c/Bamgx/objects/SetOfE4.cpp

    r3335 r3570  
    7777                //check that nbax <=NbOfEdges
    7878                if (nbax <=NbOfEdges ) {
    79                         ISSMERROR(exprintf("SetOfEdges4::add overflow: NbOfEdges=%i > nbax=%i",NbOfEdges,nbax));
     79                        ISSMERROR("SetOfEdges4::add overflow: NbOfEdges=%i > nbax=%i",NbOfEdges,nbax);
    8080                }
    8181
  • issm/trunk/src/c/Bamgx/objects/SubDomain.cpp

    r3332 r3570  
    1616                *this = Th.subdomains[i];
    1717                if ( head-Th.triangles<0 || head-Th.triangles>=Th.nbt){
    18                         ISSMERROR(exprintf("head-Th.triangles<0 || head-Th.triangles>=Th.nbt"));
     18                        ISSMERROR("head-Th.triangles<0 || head-Th.triangles>=Th.nbt");
    1919                }
    2020                head = ThNew.triangles + Th.Number(head) ;
    2121                if (edge-Th.edges<0 || edge-Th.edges>=Th.nbe);{
    22                         ISSMERROR(exprintf("edge-Th.edges<0 || edge-Th.edges>=Th.nbe"));
     22                        ISSMERROR("edge-Th.edges<0 || edge-Th.edges>=Th.nbe");
    2323                }
    2424                edge = ThNew.edges+ Th.Number(edge);
  • issm/trunk/src/c/Bamgx/objects/Triangle.cpp

    r3407 r3570  
    1515                long nbv = Th->nbv;
    1616                if (i<0 || j<0 || k<0){
    17                         ISSMERROR(exprintf("i<0 || j<0 || k<0"));
     17                        ISSMERROR("i<0 || j<0 || k<0");
    1818                }
    1919                if (i>=nbv || j>=nbv || k>=nbv){
    20                         ISSMERROR(exprintf("i>=nbv || j>=nbv || k>=nbv"));
     20                        ISSMERROR("i>=nbv || j>=nbv || k>=nbv");
    2121                }
    2222                TriaVertices[0]=v+i;
     
    114114                        //check number of iterations
    115115                        if (k>=2000){
    116                                 ISSMERROR(exprintf("too many iteration in Triangle::FindBoundaryEdge (k>=2000)"));
     116                                ISSMERROR("too many iteration in Triangle::FindBoundaryEdge (k>=2000)");
    117117                        }
    118118                } while (this!= t);
     
    137137                do {
    138138                        while (t->swap(j,koption)){
    139                                 if (k>=20000) ISSMERROR(exprintf("k>=20000"));
     139                                if (k>=20000) ISSMERROR("k>=20000");
    140140                                NbSwap++;
    141141                                k++;
  • issm/trunk/src/c/Bamgx/objects/Triangles.cpp

    r3566 r3570  
    122122                                int i2 = Tho.Number(t[2]);
    123123                                if (i0<0 || i1<0 || i2<0){
    124                                         ISSMERROR(exprintf("i0<0 || i1<0 || i2< 0"));
     124                                        ISSMERROR("i0<0 || i1<0 || i2< 0");
    125125                                }
    126126                                if (i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv){
    127                                         ISSMERROR(exprintf("i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv"));
     127                                        ISSMERROR("i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv");
    128128                                }
    129129                                triangles[nbt] = Triangle(this,kk[i0],kk[i1],kk[i2]);
     
    132132                          }
    133133                  if (kt!=nbt){
    134                           ISSMERROR(exprintf("kt!=nbt"));
     134                          ISSMERROR("kt!=nbt");
    135135                  }
    136136                  if (nbt==0 && nbv==0) {
    137                           ISSMERROR(exprintf("All triangles have been removed"));
     137                          ISSMERROR("All triangles have been removed");
    138138                  }
    139139                  delete [] kk;
     
    147147
    148148                  if (!NbSubDomains){
    149                           ISSMERROR(exprintf("NbSubDomains==0"));
     149                          ISSMERROR("NbSubDomains==0");
    150150                  }
    151151                  if (!subdomains[0].head || !subdomains[0].head->link){
    152                           ISSMERROR(exprintf("!subdomains[0].head || !subdomains[0].head->link"));
     152                          ISSMERROR("!subdomains[0].head || !subdomains[0].head->link");
    153153                  }
    154154
     
    334334                }
    335335                else{
    336                         if(verbose>5) ISSMERROR(exprintf("no Vertices found in the initial mesh"));
     336                        if(verbose>5) ISSMERROR("no Vertices found in the initial mesh");
    337337                }
    338338
     
    352352                }
    353353                else{
    354                         if(verbose>5) ISSMERROR(exprintf("no Triangles found in the initial mesh"));
     354                        if(verbose>5) ISSMERROR("no Triangles found in the initial mesh");
    355355                }
    356356
     
    478478                                //Check value
    479479                                if(!(i>=0 && j>=0 && i<nbe && j<Gh.nbe)) {
    480                                         ISSMERROR(exprintf("ReadMesh error: EdgesOnGeometricEdge edge provided (line %i: [%i %i]) is incorrect (must be positive, [0<i<nbe=%i 0<j<Gh.nbe=%i]",i1+1,i+1,j+1,nbe,Gh.nbe));
     480                                        ISSMERROR("ReadMesh error: EdgesOnGeometricEdge edge provided (line %i: [%i %i]) is incorrect (must be positive, [0<i<nbe=%i 0<j<Gh.nbe=%i]",i1+1,i+1,j+1,nbe,Gh.nbe);
    481481                                }
    482482                                edges[i].onGeometry=Gh.edges+j;
     
    504504                                head=(int)bamgmesh->SubDomains[i*3+1]-1;//C indexing
    505505                                sens=(int)bamgmesh->SubDomains[i*3+2];
    506                                 if (i3!=23) ISSMERROR(exprintf("Bad Subdomain definition: first number should be 3"));
    507                                 if (head<0 || head>=nbt) ISSMERROR(exprintf("Bad Subdomain definition: head should in [1 %i] (triangle number)",nbt));
     506                                if (i3!=23) ISSMERROR("Bad Subdomain definition: first number should be 3");
     507                                if (head<0 || head>=nbt) ISSMERROR("Bad Subdomain definition: head should in [1 %i] (triangle number)",nbt);
    508508                                subdomains[i].head = triangles+head;
    509509                        }
     
    552552                                for (j=0;j<3;j++){
    553553                                        int v=Number(triangles[i][j]); //jth vertex of the ith triangle
    554                                         if (k>3*nbt-1 || k<0) ISSMERROR(exprintf("k = %i, nbt = %i",k,nbt));
     554                                        if (k>3*nbt-1 || k<0) ISSMERROR("k = %i, nbt = %i",k,nbt);
    555555                                        next_1[k]=head_1[v];
    556                                         if (v>nbv-1 || v<0)   ISSMERROR(exprintf("v = %i, nbv = %i",v,nbv));
     556                                        if (v>nbv-1 || v<0)   ISSMERROR("v = %i, nbv = %i",v,nbv);
    557557                                        head_1[v]=k++;
    558558                                        connectivitysize_1[v]+=1;
     
    692692                                }
    693693                                if (!stop){
    694                                         ISSMERROR(exprintf("Element holding segment [%i %i] not found...",i1+1,i2+1));
     694                                        ISSMERROR("Element holding segment [%i %i] not found...",i1+1,i2+1);
    695695                                }
    696696                        }
     
    790790                                const VertexOnGeom &v=VerticesOnGeomEdge[i];
    791791                                if (!v.OnGeomEdge()){
    792                                         ISSMERROR(exprintf("A vertices supposed to be OnGeometricEdge is actually not"));
     792                                        ISSMERROR("A vertices supposed to be OnGeometricEdge is actually not");
    793793                                }
    794794                                bamgmesh->VerticesOnGeometricEdge[i*3+0]=Number((MeshVertex*)v)+1; //back to Matlab indexing
     
    873873                        for (j=0;j<2;j++){
    874874                                int v=(int)bamgmesh->ElementEdges[i*i2+j]-1; //back to C indexing
    875                                 if (k>2*i1-1 || k<0) ISSMERROR(exprintf("Index exceed matrix dimensions (k=%i not in [0 %i]",k,2*i1-1));
     875                                if (k>2*i1-1 || k<0) ISSMERROR("Index exceed matrix dimensions (k=%i not in [0 %i]",k,2*i1-1);
    876876                                next_2[k]=head_2[v];
    877                                 if (v>nbv-1 || v<0)   ISSMERROR(exprintf("Index exceed matrix dimensions (v=%i not in [0 %i])",v,nbv-1));
     877                                if (v>nbv-1 || v<0)   ISSMERROR("Index exceed matrix dimensions (v=%i not in [0 %i])",v,nbv-1);
    878878                                head_2[v]=k++;
    879879                                connectivitysize_2[v]+=1;
     
    10081008                //check that hmax is positive
    10091009                if (hmax<=0){
    1010                         ISSMERROR(exprintf("hmax<=0"));
     1010                        ISSMERROR("hmax<=0");
    10111011                }
    10121012
     
    10371037
    10381038                                if (ht<=0 || hn<=0){
    1039                                         ISSMERROR(exprintf("ht<=0 || hn<=0"));
     1039                                        ISSMERROR("ht<=0 || hn<=0");
    10401040                                }
    10411041                                MatVVP2x2 Vp(1/(ht*ht),1/(hn*hn),tg);
     
    10621062                }
    10631063                else{
    1064                         ISSMERROR(exprintf("Hessiantype %i not supported yet (1->use Green formula, 0-> double P2 projection)",Hessiantype));
     1064                        ISSMERROR("Hessiantype %i not supported yet (1->use Green formula, 0-> double P2 projection)",Hessiantype);
    10651065                }
    10661066        }
     
    11041104                //some checks
    11051105                if (( infinitevertexpos <0 ) && (detOld <0) ||  ( infinitevertexpos >=0  ) && (detOld >0) ){
    1106                         ISSMERROR(exprintf("bug in Triangles::Add, bad configuration"));
     1106                        ISSMERROR("bug in Triangles::Add, bad configuration");
    11071107                }
    11081108
     
    11491149                                printf("\nproblem while trying to add:\n");
    11501150                                s.Echo();
    1151                                 ISSMERROR(exprintf("Bug in Triangles::Add points duplicated %i times",nbzerodet));
     1151                                ISSMERROR("Bug in Triangles::Add points duplicated %i times",nbzerodet);
    11521152                        }
    11531153                }
     
    11631163
    11641164                if (nbt>nbtx) {
    1165                         ISSMERROR(exprintf("Not ebough triangles"));
     1165                        ISSMERROR("Not ebough triangles");
    11661166                }
    11671167
     
    12021202
    12031203                        if (!rswap) {
    1204                                 ISSMERROR(exprintf("swap the point s is on a edge"));
     1204                                ISSMERROR("swap the point s is on a edge");
    12051205                        }
    12061206                }
     
    12691269                //check that the mesh is not empty
    12701270                if (nbt<=0 || nbv <=0 ) {
    1271                         ISSMERROR(exprintf("nbt or nbv is negative (Mesh empty?)"));
     1271                        ISSMERROR("nbt or nbv is negative (Mesh empty?)");
    12721272                }
    12731273
     
    12901290                //check that there is no double edge
    12911291                if (nbe !=  edge4->nb()){
    1292                         ISSMERROR(exprintf("Some Double edge in the mesh, the number is %i, nbe4=%i",nbe,edge4->nb()));
     1292                        ISSMERROR("Some Double edge in the mesh, the number is %i, nbe4=%i",nbe,edge4->nb());
    12931293                }
    12941294                //keep nbe in nbeold
     
    13101310                                        //check that it is not an edge on boundary (should not already exist)
    13111311                                        if (triangles[i].TriangleAdj(j) || triangles[st[k]/3].TriangleAdj((int) (st[k]%3))){
    1312                                                 ISSMERROR(exprintf("problem in Geometry reconstruction: an edge on boundary is duplicated (double element?)"));
     1312                                                ISSMERROR("problem in Geometry reconstruction: an edge on boundary is duplicated (double element?)");
    13131313                                        }
    13141314                                        //OK, the element is not on boundary, is belongs to 2 triangles -> build Adjacent triangles list
     
    13281328                                        printf("Edge %i of triangle %i\n",(-st[k]+2)%3,(-st[k]+2)/3);
    13291329                                        printf("Edge %i of triangle %i\n",triangles[(-st[k]+2)/3].NuEdgeTriangleAdj((int)((-st[k]+2)%3)),Number(triangles[(-st[k]+2)/3].TriangleAdj((int)((-st[k]+2)%3))));
    1330                                         ISSMERROR(exprintf("An edge belongs to more than 2 triangles"));
     1330                                        ISSMERROR("An edge belongs to more than 2 triangles");
    13311331                                }       
    13321332                        }
     
    14101410                        //check that we have been through all edges
    14111411                        if (k!=nbe){
    1412                                 ISSMERROR(exprintf("problem in edge construction process: k!=nbe (should not happen)"));
     1412                                ISSMERROR("problem in edge construction process: k!=nbe (should not happen)");
    14131413                        }
    14141414                        //delete edgessave
     
    14551455                                        //check that we have the correct vertex
    14561456                                        if (v!=edges[i0 ].v[j0]){
    1457                                                 ISSMERROR(exprintf("v!=edges[i0 ].v[j0]: this should not happen as the vertex belongs to this edge"));
     1457                                                ISSMERROR("v!=edges[i0 ].v[j0]: this should not happen as the vertex belongs to this edge");
    14581458                                        }
    14591459
     
    14721472                //check that NbSubDomains is empty
    14731473                if (NbSubDomains){
    1474                         ISSMERROR(exprintf("NbSubDomains should be 0"));
     1474                        ISSMERROR("NbSubDomains should be 0");
    14751475                }
    14761476                NbSubDomains=0;
     
    15401540                //check that we have been through all subdomains
    15411541                if (k!= NbSubDomains){
    1542                         ISSMERROR(exprintf("k!= NbSubDomains"));
     1542                        ISSMERROR("k!= NbSubDomains");
    15431543                }
    15441544                //delete colorT and st
     
    16001600                Gh.coefIcoor= (MaxICoor)/(Max(Gh.pmax.x-Gh.pmin.x,Gh.pmax.y-Gh.pmin.y));
    16011601                if (Gh.coefIcoor<=0){
    1602                         ISSMERROR(exprintf("Gh.coefIcoor<=0 in infered Geometry (this should not happen)"));
     1602                        ISSMERROR("Gh.coefIcoor<=0 in infered Geometry (this should not happen)");
    16031603                }
    16041604
     
    16511651                        k = edge4->SortAndAdd(i0,i1);
    16521652                        if (k != i){
    1653                                 ISSMERROR(exprintf("problem in Edge4 construction: k != i"));
     1653                                ISSMERROR("problem in Edge4 construction: k != i");
    16541654                        }
    16551655                }
     
    16801680                        }
    16811681                        else
    1682                          ISSMERROR(exprintf("%i should be >=0"));
     1682                         ISSMERROR("%i should be >=0");
    16831683                  }
    16841684
     
    22372237                                //The vertex cannot be marked more than twice
    22382238                                if (splitvertex[i1]==3 || splitvertex[i2]==3){
    2239                                         ISSMERROR(exprintf("Crossing rifts not supported yet"));
     2239                                        ISSMERROR("Crossing rifts not supported yet");
    22402240                                }
    22412241                        }
     
    23152315                                //printf("%i -> %i %i %i, edge [%i->%i %i->%i]\n",element_renu[Number(ta.t)],Number((*ta.t)[0])+1,Number((*ta.t)[1])+1,Number((*ta.t)[2])+1,i1,j1,i2,j2);
    23162316                                ta = Next(ta).Adj();
    2317                                 if (count++>50) ISSMERROR(exprintf("Maximum number of iteration exceeded"));
     2317                                if (count++>50) ISSMERROR("Maximum number of iteration exceeded");
    23182318                        }while ((tbegin != ta));
    23192319                }
     
    23222322                for(i=0;i<NbCrackedEdges;i++){
    23232323                        if (Edgeflags[i]!=2){
    2324                                 ISSMERROR(exprintf("A problem occured: at least one crack edge (number %i) does not belong to 2 elements",i+1));
     2324                                ISSMERROR("A problem occured: at least one crack edge (number %i) does not belong to 2 elements",i+1);
    23252325                        }
    23262326                }
     
    23482348                        }
    23492349                        if (k!=0) {
    2350                                 ISSMERROR(exprintf("there is %i triangles of mes = 0",k));
     2350                                ISSMERROR("there is %i triangles of mes = 0",k);
    23512351                        }
    23522352
     
    23652365                                                printf("Edge %i: %i %i\n",j,Number(edges[j][0]),Number(edges[j][1]));
    23662366                                        }
    2367                                         ISSMERROR(exprintf("Missing Edge %i, v0=%i,v1=%i",i,Number(edges[i][0]),Number(edges[i][1])));
     2367                                        ISSMERROR("Missing Edge %i, v0=%i,v1=%i",i,Number(edges[i][0]),Number(edges[i][1]));
    23682368                                }
    23692369                        }
    23702370
    23712371                        if (k!=0) {
    2372                                 ISSMERROR(exprintf("There are %i lost edges, the boundary might be crossing",k));
     2372                                ISSMERROR("There are %i lost edges, the boundary might be crossing",k);
    23732373                        }
    23742374                        for (int j=0;j<nbv;j++){
     
    24512451                        it++;} // end while (it<nbt)
    24522452                        if (nbt == NbOutT ||  !NbSubDomTot) {
    2453                                 ISSMERROR(exprintf("The boundary is not close: all triangles are outside"));
     2453                                ISSMERROR("The boundary is not close: all triangles are outside");
    24542454                        }
    24552455
     
    24882488                                                it++;} // end white (it<nbt)
    24892489                                                if (k!=NbSubDomains){
    2490                                                         ISSMERROR(exprintf("k!=NbSubDomains"));
     2490                                                        ISSMERROR("k!=NbSubDomains");
    24912491                                                }
    24922492                                                if(OutSide)
     
    25592559                                        Edge &e = *GeometricalEdgetoEdge[Gh.Number(eg)];
    25602560                                        if (&e==NULL){
    2561                                                 ISSMERROR(exprintf("&e==NULL"));
     2561                                                ISSMERROR("&e==NULL");
    25622562                                        }
    25632563                                        MeshVertex * v0 =  e(0),*v1 = e(1);
     
    25692569                                        subdomains[i].edge = &e;
    25702570                                        if (!t || !sens){
    2571                                                 ISSMERROR(exprintf("!t || !sens"));
     2571                                                ISSMERROR("!t || !sens");
    25722572                                        }
    25732573
     
    25762576                                        while (1) {
    25772577                                                if ( v0 != ta.EdgeVertex(1) ){
    2578                                                         ISSMERROR(exprintf("v0 != ta.EdgeVertex(1)"));
     2578                                                        ISSMERROR("v0 != ta.EdgeVertex(1)");
    25792579                                                }
    25802580                                                if (ta.EdgeVertex(0) == v1) { // ok we find the edge
     
    25842584                                                         subdomains[i].head=t=ta;
    25852585                                                        if(t<triangles || t >= triangles+nbt || t->det < 0 || t->link == 0) {
    2586                                                                 ISSMERROR(exprintf("bad definition of SubSomain %i",i));
     2586                                                                ISSMERROR("bad definition of SubSomain %i",i);
    25872587                                                        }
    25882588                                                        long it = Number(t);
     
    25992599                                                                kkk++;
    26002600                                                                if (mark[Number(tt)]>=0){
    2601                                                                         ISSMERROR(exprintf("mark[Number(tt)]>=0"));
     2601                                                                        ISSMERROR("mark[Number(tt)]>=0");
    26022602                                                                }
    26032603                                                                mark[Number(tt)]=i;
     
    26082608                                                ta = Previous(Adj(ta));         
    26092609                                                if(t == (Triangle *) ta) {
    2610                                                         ISSMERROR(exprintf("bad definition of SubSomain %i",i));
     2610                                                        ISSMERROR("bad definition of SubSomain %i",i);
    26112611                                                }
    26122612                                        }
     
    26392639                else {
    26402640                        if (!quadtree){
    2641                                 ISSMERROR(exprintf("no starting triangle provided and no quadtree available"));
     2641                                ISSMERROR("no starting triangle provided and no quadtree available");
    26422642                        }
    26432643                        MeshVertex *a = quadtree->NearestVertex(B.x,B.y) ;
     
    26472647                                        printf("TriangleContaining vertex number %i, another call to ReMakeTriangleContainingTheVertex was required\n", Number(a));
    26482648                                }
    2649                                 ISSMERROR(exprintf("problem in Triangles::FindTriangleContaining"));
     2649                                ISSMERROR("problem in Triangles::FindTriangleContaining");
    26502650                        }
    26512651                        if (a<vertices || a>=vertices+nbv){
    2652                                 ISSMERROR(exprintf("a<vertices || a>=vertices+nbv"));
     2652                                ISSMERROR("a<vertices || a>=vertices+nbv");
    26532653                        }
    26542654                        t = a->t;
    26552655                        if (t<triangles || t>=triangles+nbt){
    26562656                                a->Echo();
    2657                                 ISSMERROR(exprintf("t<triangles || t>=triangles+nbt"));
     2657                                ISSMERROR("t<triangles || t>=triangles+nbt");
    26582658                        }
    26592659                }
     
    26672667                        int k0=(*t)(0) ?  ((  (*t)(1) ? ( (*t)(2) ? -1 : 2) : 1  )) : 0;
    26682668                        if (k0<0){ // k0 the NULL vertex
    2669                                 ISSMERROR(exprintf("k0<0"));
     2669                                ISSMERROR("k0<0");
    26702670                        }
    26712671                        int k1=NextVertex[k0],k2=PreviousVertex[k0];
     
    26772677                        kkkk++;
    26782678                        if (kkkk>=2){
    2679                                 ISSMERROR(exprintf("kkkk>=2"));
     2679                                ISSMERROR("kkkk>=2");
    26802680                        }
    26812681                }
     
    26872687                        kkkk++;
    26882688                        if (kkkk>=2000){
    2689                                 ISSMERROR(exprintf("kkkk>=2000"));
     2689                                ISSMERROR("kkkk>=2000");
    26902690                        }
    26912691                        j= OppositeVertex[jj];
     
    27082708                        if (k == 2 && BinaryRand()) Exchange(ii[0],ii[1]);
    27092709                        if ( k>=3){
    2710                                 ISSMERROR(exprintf("k>=3"));
     2710                                ISSMERROR("k>=3");
    27112711                        }
    27122712                        TriangleAdjacent t1 = t->Adj(jj=ii[0]);
     
    27642764                VerticesOnGeomVertex = new VertexOnGeom[NbVerticesOnGeomVertex]; 
    27652765                if( NbVerticesOnGeomVertex >= nbvx) {
    2766                         ISSMERROR(exprintf("too many vertices on geometry: %i >= %i",NbVerticesOnGeomVertex,nbvx));
     2766                        ISSMERROR("too many vertices on geometry: %i >= %i",NbVerticesOnGeomVertex,nbvx);
    27672767                }
    27682768
     
    27812781                                }
    27822782                                else{
    2783                                         ISSMERROR(exprintf("Maximum number of vertices (nbvx = %i) too small",nbvx));
     2783                                        ISSMERROR("Maximum number of vertices (nbvx = %i) too small",nbvx);
    27842784                                }
    27852785                               
     
    27972797                //check that edges has been allocated
    27982798                if (edges){
    2799                         ISSMERROR(exprintf("edges is empty"));
     2799                        ISSMERROR("edges is empty");
    28002800                }
    28012801
     
    28432843                                                                        //check that edges has been allocated
    28442844                                                                        if (!edges){
    2845                                                                                 ISSMERROR(exprintf("edges has not been allocated..."));
     2845                                                                                ISSMERROR("edges has not been allocated...");
    28462846                                                                        }
    28472847                                                                        edges[nbe].v[0]=a->to;
     
    28652865                                                                NbNewPoints=0;
    28662866                                                                NbEdgeCurve=0;
    2867                                                                 if (nbvend>=nbvx){
    2868                                                                         ISSMERROR(exprintf("nbvend>=nbvx"));
    2869                                                                 }
     2867                                                                if (nbvend>=nbvx) ISSMERROR("maximum number of vertices too low! Check the domain outline or increase nbvx");
    28702868                                                                lcurve =0;
    28712869                                                                s = lstep;
     
    29402938                                                                                         kk0=kkk,ll0=llk;
    29412939                                                                                }
    2942                                                                                 if (kk1 == kk0){
    2943                                                                                         ISSMERROR(exprintf("kk1 == kk0"));
    2944                                                                                 }
     2940                                                                                ISSMASSERT(kk1!=kk0);
    29452941
    29462942                                                                                double sbb = (ss-ll0  )/(ll1-ll0);
     
    29752971                                                                        k = e->DirAdj[kprev];// next vertices
    29762972                                                                        e = e->Adj[kprev];
    2977                                                                         if (!e){
    2978                                                                                 ISSMERROR(exprintf("!e"));
    2979                                                                         }
     2973                                                                        ISSMASSERT(e);
    29802974                                                                }// for(;;)
    29812975                                                                vb = b->to;
     
    30063000                        if(!step) {
    30073001                                if (edges){
    3008                                         ISSMERROR(exprintf("edges"));
     3002                                        ISSMERROR("edges");
    30093003                                }
    30103004                                if (VerticesOnGeomEdge){
    3011                                         ISSMERROR(exprintf("VerticesOnGeomEdge"));
     3005                                        ISSMERROR("VerticesOnGeomEdge");
    30123006                                }
    30133007                                edges = new Edge[nbex=nbe];
     
    30153009                                 VerticesOnGeomEdge = new VertexOnGeom[NbVerticesOnGeomEdge0];
    30163010                                if (!VerticesOnGeomEdge && NbVerticesOnGeomEdge0!=0){
    3017                                         ISSMERROR(exprintf("!VerticesOnGeomEdge && NbVerticesOnGeomEdge0!=0"));
     3011                                        ISSMERROR("!VerticesOnGeomEdge && NbVerticesOnGeomEdge0!=0");
    30183012                                }
    30193013                                // do the vertex on a geometrical vertex
     
    30213015                        }
    30223016                        else if (NbVerticesOnGeomEdge != NbVerticesOnGeomEdge0){
    3023                                 ISSMERROR(exprintf("NbVerticesOnGeomEdge != NbVerticesOnGeomEdge0"));
     3017                                ISSMERROR("NbVerticesOnGeomEdge != NbVerticesOnGeomEdge0");
    30243018                        }
    30253019                }
     
    30983092                int i;
    30993093                for (i=0;i<Gh.nbv;i++) if (Gh[i].Required()) NbVerticesOnGeomVertex++;
    3100                 if( NbVerticesOnGeomVertex >= nbvx) { ISSMERROR(exprintf("too many vertices on geometry: %i >= %i",NbVerticesOnGeomVertex,nbvx));}
     3094                if( NbVerticesOnGeomVertex >= nbvx) { ISSMERROR("too many vertices on geometry: %i >= %i",NbVerticesOnGeomVertex,nbvx);}
    31013095
    31023096                VerticesOnGeomVertex = new VertexOnGeom[  NbVerticesOnGeomVertex];
     
    31573151                }
    31583152                if (bfind!=Gh.NbOfCurves){
    3159                         ISSMERROR(exprintf("problem generating number of curves (Gh.NbOfCurves=%i bfind=%i)",Gh.NbOfCurves,bfind));
     3153                        ISSMERROR("problem generating number of curves (Gh.NbOfCurves=%i bfind=%i)",Gh.NbOfCurves,bfind);
    31603154                }
    31613155
     
    32493243                                                                                double se= (sNew-L0)/LAB;
    32503244                                                                                if (se<0 || se>=1.000000001){
    3251                                                                                         ISSMERROR(exprintf("Problem creating point on a boundary: se=%g should be in [0 1]",se));
     3245                                                                                        ISSMERROR("Problem creating point on a boundary: se=%g should be in [0 1]",se);
    32523246                                                                                }
    32533247                                                                                se = abscisseInterpole(v0.m,v1.m,AB,se,1);
    32543248                                                                                if (se<0 || se>1){
    3255                                                                                         ISSMERROR(exprintf("Problem creating point on a boundary: se=%g should be in [0 1]",se));
     3249                                                                                        ISSMERROR("Problem creating point on a boundary: se=%g should be in [0 1]",se);
    32563250                                                                                }
    32573251                                                                                se = k1         ? se : 1. - se;
     
    32853279                                                                }
    32863280                                                                if (!ee.adj[k1]) {
    3287                                                                         ISSMERROR(exprintf(" adj edge %i, nbe=%i, Gh.vertices=%i",BTh.Number(ee),nbe,Gh.vertices));
     3281                                                                        ISSMERROR(" adj edge %i, nbe=%i, Gh.vertices=%i",BTh.Number(ee),nbe,Gh.vertices);
    32883282                                                                }
    32893283                                                                pe = ee.adj[k1]; // next edge
     
    33273321                        if(step==0){
    33283322                                if(nbv+NbOfNewPoints > nbvx) {
    3329                                         ISSMERROR(exprintf("too many vertices on geometry: %i >= %i",nbv+NbOfNewPoints,nbvx));
     3323                                        ISSMERROR("too many vertices on geometry: %i >= %i",nbv+NbOfNewPoints,nbvx);
    33303324                                }
    33313325                                edges = new Edge[NbOfNewEdge];
     
    34223416                        //all vertices are aligned!
    34233417                        if  ( ++i >= nbv) {
    3424                                 ISSMERROR(exprintf("all the vertices are aligned"));
     3418                                ISSMERROR("all the vertices are aligned");
    34253419                        }
    34263420                }
     
    35603554                                long  j=vj.ReferenceNumber;
    35613555                                if (&vj!=ordre[j]){
    3562                                         ISSMERROR(exprintf("&vj!= ordre[j]"));
     3556                                        ISSMERROR("&vj!= ordre[j]");
    35633557                                }
    35643558                                if(i!=j){
     
    35703564                                if (tcvj && !tcvj->link){
    35713565                                        tcvj->Echo();
    3572                                         ISSMERROR(exprintf("problem inserting point in InsertNewPoints (tcvj=%p and tcvj->link=%i)",tcvj,tcvj->link));
     3566                                        ISSMERROR("problem inserting point in InsertNewPoints (tcvj=%p and tcvj->link=%i)",tcvj,tcvj->link);
    35733567                                }
    35743568                                quadtree->Add(vj);
     
    35973591
    35983592                if (!Gh.nbe){
    3599                         ISSMERROR(exprintf("!Gh.nbe"));
     3593                        ISSMERROR("!Gh.nbe");
    36003594                }
    36013595                Edge **e= new (Edge* [Gh.nbe]);
     
    36303624                        }
    36313625                }
    3632                 if(kk) ISSMERROR(exprintf("See above"));
     3626                if(kk) ISSMERROR("See above");
    36333627
    36343628                return e;
     
    38193813                                //check i
    38203814                                if (i<0 || i>=nbt ){
    3821                                         ISSMERROR(exprintf("Index problem in NewPoints (i=%i not in [0 %i])",i,nbt-1));
     3815                                        ISSMERROR("Index problem in NewPoints (i=%i not in [0 %i])",i,nbt-1);
    38223816                                }
    38233817                                //change first_np_or_next_t[i]
     
    38703864                                        }
    38713865                                        if (ta.EdgeVertex(0)!=s){
    3872                                                 ISSMERROR(exprintf("ta.EdgeVertex(0)!=s"));
     3866                                                ISSMERROR("ta.EdgeVertex(0)!=s");
    38733867                                        }
    38743868                                        ta = Next(Adj(ta));
     
    39203914                vertices=new MeshVertex[nbvx];
    39213915                if (!vertices){
    3922                         ISSMERROR(exprintf("!vertices"));
     3916                        ISSMERROR("!vertices");
    39233917                }
    39243918                ordre=new (MeshVertex* [nbvx]);
    39253919                if (!ordre){
    3926                         ISSMERROR(exprintf("!ordre"));
     3920                        ISSMERROR("!ordre");
    39273921                }
    39283922                triangles=new Triangle[nbtx];
    39293923                if (!triangles){
    3930                         ISSMERROR(exprintf("!triangles"));
     3924                        ISSMERROR("!triangles");
    39313925                }
    39323926        }
     
    39653959                }
    39663960                else {
    3967                         ISSMERROR(exprintf("ProjectOnCurve On MeshVertex %i forget call to SetVertexFieldOnBTh",BTh.Number(vA)));
     3961                        ISSMERROR("ProjectOnCurve On MeshVertex %i forget call to SetVertexFieldOnBTh",BTh.Number(vA));
    39683962                }
    39693963
     
    39763970                }
    39773971                else {
    3978                         ISSMERROR(exprintf("ProjectOnCurve On MeshVertex %i forget call to SetVertexFieldOnBTh",BTh.Number(vB)));
     3972                        ISSMERROR("ProjectOnCurve On MeshVertex %i forget call to SetVertexFieldOnBTh",BTh.Number(vB));
    39793973                }
    39803974                Edge * e = &BhAB;
    39813975                if (!pA || !pB || !e){
    3982                         ISSMERROR(exprintf("!pA || !pB || !e"));
     3976                        ISSMERROR("!pA || !pB || !e");
    39833977                }
    39843978                // be carefull the back ground edge e is on same geom edge
     
    39863980                //check Is a background Mesh;   
    39873981                if (e<BTh.edges || e>=BTh.edges+BTh.nbe){
    3988                         ISSMERROR(exprintf("e<BTh.edges || e>=BTh.edges+BTh.nbe"));
     3982                        ISSMERROR("e<BTh.edges || e>=BTh.edges+BTh.nbe");
    39893983                }
    39903984                // walk on BTh edge
     
    40094003                  }
    40104004                else{ // do the search by walking
    4011                         ISSMERROR(exprintf("case not supported yet"));
     4005                        ISSMERROR("case not supported yet");
    40124006                  }
    40134007
     
    40384032                                kkk=kkk+1;
    40394033                                if (kkk>=100){
    4040                                         ISSMERROR(exprintf("kkk>=100"));
     4034                                        ISSMERROR("kkk>=100");
    40414035                                }
    40424036                                if (!eee){
    4043                                         ISSMERROR(exprintf("!eee"));
     4037                                        ISSMERROR("!eee");
    40444038                                }
    40454039                                double lg0 = lg;
     
    40504044                                        double thetab = te0*(1-sss)+ sss*iii;
    40514045                                        if (thetab<0 || thetab>1){
    4052                                                 ISSMERROR(exprintf("thetab<0 || thetab>1"));
     4046                                                ISSMERROR("thetab<0 || thetab>1");
    40534047                                        }
    40544048                                        BR = VertexOnEdge(&R,eee,thetab);
     
    40634057                                double lg0 = lg;
    40644058                                if (!eee){
    4065                                         ISSMERROR(exprintf("!eee"));
     4059                                        ISSMERROR("!eee");
    40664060                                }
    40674061                                v1 = pvB;
     
    40744068                                        double thetab = te0*(1-sss)+ sss*tB;
    40754069                                        if (thetab<0 || thetab>1){
    4076                                                 ISSMERROR(exprintf("thetab<0 || thetab>1"));
     4070                                                ISSMERROR("thetab<0 || thetab>1");
    40774071                                        }
    40784072                                        BR = VertexOnEdge(&R,eee,thetab);
     
    40834077
    40844078                  }
    4085                 ISSMERROR(exprintf("Big bug..."));
     4079                ISSMERROR("Big bug...");
    40864080                return 0; // just for the compiler
    40874081        }                 
     
    41244118        }
    41254119        if (kk != nbe){
    4126                 ISSMERROR(exprintf("There are %i double edges in the mesh",kk-nbe));
     4120                ISSMERROR("There are %i double edges in the mesh",kk-nbe);
    41274121        }
    41284122
     
    41564150                        //An edge belongs to 2 triangles
    41574151                        else {
    4158                                 ISSMERROR(exprintf("The edge (%i , %i) belongs to more than 2 triangles",
    4159                                                                 Number(triangles[i][VerticesOfTriangularEdge[j][0]]),Number(triangles[i][VerticesOfTriangularEdge[j][1]])));
     4152                                ISSMERROR("The edge (%i , %i) belongs to more than 2 triangles",Number(triangles[i][VerticesOfTriangularEdge[j][0]]),Number(triangles[i][VerticesOfTriangularEdge[j][1]]));
    41604153                        }
    41614154                }
     
    41954188        }
    41964189        if(k) {
    4197                 ISSMERROR(exprintf("%i boundary edges (from the geometry) are not defined as mesh edges",k));
     4190                ISSMERROR("%i boundary edges (from the geometry) are not defined as mesh edges",k);
    41984191        }
    41994192
     
    42224215        for (i=2;det(ordre[0]->i,ordre[1]->i,ordre[i]->i)==0;)
    42234216         if  (++i>=nbvb) {
    4224                  ISSMERROR(exprintf("ReconstructExistingMesh: All the vertices are aligned"));
     4217                 ISSMERROR("ReconstructExistingMesh: All the vertices are aligned");
    42254218         }
    42264219        //Move this vertex (i) to the 2d position in ordre
     
    42824275        }
    42834276        if(nbloss) {
    4284                 ISSMERROR(exprintf("we lost %i existing edges other %i",nbloss,knbe));
     4277                ISSMERROR("we lost %i existing edges other %i",nbloss,knbe);
    42854278        }
    42864279
     
    43574350        subdomains = savesubdomains;
    43584351        if (k) {
    4359                 ISSMERROR(exprintf("number of triangles edges alone = %i",k));
     4352                ISSMERROR("number of triangles edges alone = %i",k);
    43604353        }
    43614354        FindSubDomain();
     
    43874380
    43884381                                                printf("This edge is on geometry and has no adjacent edge (open curve) and one of the tip is not required\n");
    4389                                                 ISSMERROR(exprintf("See above (might be cryptic...)"));
     4382                                                ISSMERROR("See above (might be cryptic...)");
    43904383                                        }
    43914384                                }
     
    44104403                        t=t0=subdomains[i].head;
    44114404                        if (!t0){ // not empty sub domain
    4412                                 ISSMERROR(exprintf("!t0"));
     4405                                ISSMERROR("!t0");
    44134406                        }
    44144407                        do {
    44154408                                long kt = Number(t);
    44164409                                if (kt<0 || kt >= nbt ){
    4417                                         ISSMERROR(exprintf("kt<0 || kt >= nbt"));
     4410                                        ISSMERROR("kt<0 || kt >= nbt");
    44184411                                }
    44194412                                if (renu[kt]!=-1){
    4420                                         ISSMERROR(exprintf("renu[kt]!=-1"));
     4413                                        ISSMERROR("renu[kt]!=-1");
    44214414                                }
    44224415                                renu[kt]=k++;
     
    44354428                }
    44364429                if (k != nbt){
    4437                         ISSMERROR(exprintf("k != nbt"));
     4430                        ISSMERROR("k != nbt");
    44384431                }
    44394432                // do the change on all the pointeur
     
    45624555        coefIcoor= (MaxICoor)/(Max(pmax.x-pmin.x,pmax.y-pmin.y));
    45634556        if (coefIcoor<=0){
    4564                 ISSMERROR(exprintf("coefIcoor should be positive, a problem in the geometry is likely"));
     4557                ISSMERROR("coefIcoor should be positive, a problem in the geometry is likely");
    45654558        }
    45664559
     
    45834576                                if(Nberr==1)
    45844577                                 if (strfrom){
    4585                                          ISSMERROR(exprintf("Fatal error %s (SetInCoor) : area of Triangle %i < 0",strfrom,i));
     4578                                         ISSMERROR("Fatal error %s (SetInCoor) : area of Triangle %i < 0",strfrom,i);
    45864579                                 }
    45874580                                 else{
    4588                                          ISSMERROR(exprintf("Fatal error (SetInCoor) : area of Triangle %i < 0",i));
     4581                                         ISSMERROR("Fatal error (SetInCoor) : area of Triangle %i < 0",i);
    45894582                                 }
    45904583                        }
     
    47764769                        register Triangle* t= vertices[i].t;
    47774770                        if (!t){
    4778                                 ISSMERROR(exprintf("!t"));
     4771                                ISSMERROR("!t");
    47794772                        }
    47804773                        MeshVertex & vi = vertices[i];
     
    47844777                                ta=Previous(Adj(ta));
    47854778                                if (vertices+i != ta.EdgeVertex(1)){
    4786                                         ISSMERROR(exprintf("vertices+i != ta.EdgeVertex(1)"));
     4779                                        ISSMERROR("vertices+i != ta.EdgeVertex(1)");
    47874780                                }
    47884781                                MeshVertex & vj = *(ta.EdgeVertex(0));
     
    47904783                                        j= &vj-vertices;
    47914784                                        if (j<0 || j >= nbv){
    4792                                                 ISSMERROR(exprintf("j<0 || j >= nbv"));
     4785                                                ISSMERROR("j<0 || j >= nbv");
    47934786                                        }
    47944787                                        R2 Aij = (R2) vj - (R2) vi;
     
    49144907                                        // the first PB is to now a background edge between the 2 vertices
    49154908                                        if (!edgesGtoB){
    4916                                                 ISSMERROR(exprintf("!edgesGtoB"));
     4909                                                ISSMERROR("!edgesGtoB");
    49174910                                        }
    49184911                                        ong= ProjectOnCurve(*edgesGtoB[Gh.Number(edges[i].onGeometry)],
     
    49864979                        Triangle & t = triangles[i];
    49874980                        if (!t.link){
    4988                                 ISSMERROR(exprintf("!t.link"));
     4981                                ISSMERROR("!t.link");
    49894982                        }
    49904983                        for(int j=0;j<3;j++)
     
    50135006                                                                                ||   (cc=Area2( t[0].r , t[1].r , A.r    )) < 0.0)){
    50145007                                                        printf("%i not in triangle %i In= %i %i %i %i %i\n",ke + nbvold,i,!!t.link,aa,bb,cc,dd);
    5015                                                         ISSMERROR(exprintf("Number of triangles with P2 interpolation Problem"));
     5008                                                        ISSMERROR("Number of triangles with P2 interpolation Problem");
    50165009                                                }
    50175010                                        }
     
    50215014                                                                                ||   (cc=Area2( tt[0].r , tt[1].r , A.r     )) < 0)){
    50225015                                                        printf("%i not in triangle %i In= %i %i %i %i %i\n",ke + nbvold,ii,!!tt.link,aa,bb,cc,dd);
    5023                                                         ISSMERROR(exprintf("Number of triangles with P2 interpolation Problem"));
     5016                                                        ISSMERROR("Number of triangles with P2 interpolation Problem");
    50245017                                                }
    50255018                                        }
     
    50695062                                                  } // tt
    50705063                                                else
    5071                                                  ISSMERROR(exprintf("Bug..."));
     5064                                                 ISSMERROR("Bug...");
    50725065                                          } // ke<0           
    50735066                                        else
     
    50825075                          }
    50835076                        if (nbinvisible>=2){
    5084                                 ISSMERROR(exprintf("nbinvisible>=2"));
     5077                                ISSMERROR("nbinvisible>=2");
    50855078                        }
    50865079                        switch (nbsplitedge) {
     
    50945087                        }
    50955088                        if (ksplit[i]<40){
    5096                                 ISSMERROR(exprintf("ksplit[i]<40"));
     5089                                ISSMERROR("ksplit[i]<40");
    50975090                        }
    50985091                  }
     
    51115104                        int  ke=(int) (ksplit[i]%10);
    51125105                        if (kk>=7 || kk<=0){
    5113                                 ISSMERROR(exprintf("kk>=7 || kk<=0"));
     5106                                ISSMERROR("kk>=7 || kk<=0");
    51145107                        }
    51155108
     
    51285121
    51295122                        if (nbmkadj>=10){
    5130                                 ISSMERROR(exprintf("nbmkadj>=10"));
     5123                                ISSMERROR("nbmkadj>=10");
    51315124                        }
    51325125                        // --------------------------
     
    51465139                                                        t1=t0;
    51475140                                                        if (kedge[3*i+i0]<0){
    5148                                                                 ISSMERROR(exprintf("kedge[3*i+i0]<0"));
     5141                                                                ISSMERROR("kedge[3*i+i0]<0");
    51495142                                                        }
    51505143                                                        MeshVertex * v3 = vertices + kedge[3*i+k0];
     
    51625155                                                        t2=t1=t0;
    51635156                                                        if (kedge[3*i+k1]<0){
    5164                                                                 ISSMERROR(exprintf("kedge[3*i+k1]<0"));
     5157                                                                ISSMERROR("kedge[3*i+k1]<0");
    51655158                                                        }
    51665159                                                        if (kedge[3*i+k2]<0){
    5167                                                                 ISSMERROR(exprintf("kedge[3*i+k2]<0"));
     5160                                                                ISSMERROR("kedge[3*i+k2]<0");
    51685161                                                        }
    51695162
     
    51895182                                                        t3=t2=t1=t0;
    51905183                                                        if (kedge[3*i+k0] <0 || kedge[3*i+k1]<0 || kedge[3*i+k2]<0){
    5191                                                                 ISSMERROR(exprintf("kedge[3*i+k0] <0 || kedge[3*i+k1]<0 || kedge[3*i+k2]<0"));
     5184                                                                ISSMERROR("kedge[3*i+k0] <0 || kedge[3*i+k1]<0 || kedge[3*i+k2]<0");
    51925185                                                        }
    51935186                                                        MeshVertex * v12 = vertices + kedge[3*i+k0];
     
    52635256                                }
    52645257                        if (nbmkadj>13){// 13 = 6 + 4 +
    5265                                 ISSMERROR(exprintf("nbmkadj>13"));
     5258                                ISSMERROR("nbmkadj>13");
    52665259                        }
    52675260
     
    53895382                        quadtree->Add(vi);
    53905383                        if (!tcvi || tcvi->det<0){// internal
    5391                                 ISSMERROR(exprintf("!tcvi || tcvi->det < 0"));
     5384                                ISSMERROR("!tcvi || tcvi->det < 0");
    53925385                        }
    53935386                        AddVertex(vi,tcvi,dete);
     
    54265419
    54275420                //check that the subdomain is not empty
    5428                 if (!t0){ ISSMERROR(exprintf("At least one subdomain is empty"));}
     5421                if (!t0){ ISSMERROR("At least one subdomain is empty");}
    54295422
    54305423                //loop
     
    54365429
    54375430                        //check that num is in [0 nbt[
    5438                         if (num<0 || num>=nbt){ ISSMERROR(exprintf("num<0 || num>=nbt"));}
     5431                        if (num<0 || num>=nbt){ ISSMERROR("num<0 || num>=nbt");}
    54395432
    54405433                        //reft of this triangle is the subdomain number
     
    54685461                int dir=0;
    54695462                if (k<0){
    5470                         ISSMERROR(exprintf("k<0"));
     5463                        ISSMERROR("k<0");
    54715464                }
    54725465                int kkk=0; 
     
    54765469                        kkk++;
    54775470                        if (kkk>=1000){
    5478                                 ISSMERROR(exprintf("kkk>=1000"));
     5471                                ISSMERROR("kkk>=1000");
    54795472                        }
    54805473                        MeshVertex  &vI =  *edge.EdgeVertex(0);
     
    54935486                                                        double IJ2 = IJ_IA + IJ_AJ;
    54945487                                                        if (IJ2==0){
    5495                                                                 ISSMERROR(exprintf("IJ2==0"));
     5488                                                                ISSMERROR("IJ2==0");
    54965489                                                        }
    54975490                                                        a= IJ_AJ/IJ2;
     
    55115504                //   int dir=0;
    55125505                if (t->link != 0){
    5513                         ISSMERROR(exprintf("t->link != 0"));
     5506                        ISSMERROR("t->link != 0");
    55145507                }
    55155508                // to have a starting edges
     
    55835576                  }
    55845577                if (cas ==-2){
    5585                         ISSMERROR(exprintf("cas==-2"));
     5578                        ISSMERROR("cas==-2");
    55865579                }
    55875580                // l1 = ||C s1||  , l0 = ||C s0||
     
    56005593                                kkk++;
    56015594                                if (edge.EdgeVertex(0)!=s && kkk>=10000){
    5602                                         ISSMERROR(exprintf("edge.EdgeVertex(0)!=s && kkk>=10000"));
     5595                                        ISSMERROR("edge.EdgeVertex(0)!=s && kkk>=10000");
    56035596                                }
    56045597
     
    56325625
    56335626                        if (!(Triangle *) er){
    5634                                 ISSMERROR(exprintf("!(Triangle *) er"));
     5627                                ISSMERROR("!(Triangle *) er");
    56355628                        }
    56365629                        I2 A((I2)*er.EdgeVertex(0));
     
    56565649        int NbSwap =0;
    56575650        if (!a.t || !b.t){ // the 2 vertex is in a mesh
    5658                 ISSMERROR(exprintf("!a.t || !b.t"));
     5651                ISSMERROR("!a.t || !b.t");
    56595652        }
    56605653        int k=0;
     
    56735666                vbegin =v2;
    56745667                if (!v2){
    5675                         ISSMERROR(exprintf("!v2"));
     5668                        ISSMERROR("!v2");
    56765669                }
    56775670                det2 = det(*v2,a,b);
     
    56905683                        tc = Previous(tc);
    56915684                        if (!v1 || !v2){
    5692                                 ISSMERROR(exprintf("!v1 || !v2"));
     5685                                ISSMERROR("!v1 || !v2");
    56935686                        }
    56945687                        Icoor2 detss = 0,l=0,ks;
    56955688                        while ((ks=SwapForForcingEdge(  va,  vb, tc, detss, det1,det2,NbSwap)))
    56965689                         if(l++ > 10000000) {
    5697                                  ISSMERROR(exprintf("Loop in forcing Egde, nb de swap=%i, nb of try swap (%i) too big",NbSwap,l));
     5690                                 ISSMERROR("Loop in forcing Egde, nb de swap=%i, nb of try swap (%i) too big",NbSwap,l);
    56985691                         }
    56995692                        MeshVertex *aa = tc.EdgeVertex(0), *bb = tc.EdgeVertex(1);
     
    57145707                k++;
    57155708                if (k>=2000){
    5716                         ISSMERROR(exprintf("k>=2000"));
     5709                        ISSMERROR("k>=2000");
    57175710                }
    57185711                if ( vbegin == v2 ) return -1;// error
     
    57865779                short a1=tt1,a2=tt2;// les 2 numero de l arete dans les 2 triangles
    57875780                if ( a1<0 || a1>=3 ){
    5788                         ISSMERROR(exprintf("a1<0 || a1>=3"));
     5781                        ISSMERROR("a1<0 || a1>=3");
    57895782                }
    57905783
     
    57985791                Icoor2 detT = det1+det2;
    57995792                if ((det1<=0 ) || (det2<=0)){
    5800                         ISSMERROR(exprintf("(det1<=0 ) || (det2<=0)"));
     5793                        ISSMERROR("(det1<=0 ) || (det2<=0)");
    58015794                }
    58025795                if ( (detsa>=0) || (detsb<=0) ){ // [a,b] cut infinite line va,bb
    5803                         ISSMERROR(exprintf("(detsa>=0) || (detsb<=0)"));
     5796                        ISSMERROR("(detsa>=0) || (detsb<=0)");
    58045797                }
    58055798                Icoor2 ndet1 = bamg::det(s1,sa,s2);
  • issm/trunk/src/c/BuildNodeSetsx/PartitionSets.cpp

    r3332 r3570  
    7777                                ccount++;
    7878                        }
    79                         if (flags_b_local[i] && flags_c_local[i]) ISSMERROR(exprintf("%s%i%s"," for dof ",i,": breach of exclusive partitioning between sets"));
     79                        if (flags_b_local[i] && flags_c_local[i]) ISSMERROR("%s%i%s"," for dof ",i,": breach of exclusive partitioning between sets");
    8080                       
    8181                        acount++;
  • issm/trunk/src/c/ControlConstrainx/ControlConstrainx.cpp

    r3332 r3570  
    2020                for(i=0;i<numdofnodes;i++){
    2121                        if(isnan(p_g[i])){
    22                                 ISSMERROR(exprintf("NaN found in parameter p_g[%i]",i));
     22                                ISSMERROR("NaN found in parameter p_g[%i]",i);
    2323                        }
    2424                        if(!isnan(cm_min)){
  • issm/trunk/src/c/DataSet/DataSet.cpp

    r3567 r3570  
    291291                }
    292292                else{
    293                         ISSMERROR(exprintf("%s%i"," could not recognize enum type: ",enum_type));
     293                        ISSMERROR("%s%i"," could not recognize enum type: ",enum_type);
    294294                }
    295295
     
    708708        /*Carry out a binary search on the sorted_ids: */
    709709        if(!binary_search(&id_offset,eid, sorted_ids,objects.size())){
    710                 ISSMERROR(exprintf("%s%i"," could not find object with id ",eid));
     710                ISSMERROR("%s%i"," could not find object with id ",eid);
    711711        }
    712712
     
    15011501                        if ((node1 && !node2) || (!node1 && node2)){
    15021502                                /*we are missing one node, not good!*/
    1503                                 ISSMERROR(exprintf("%s%p%s%p"," in Rgb, missing one node. node1: ",node1," node2: ",node2));
     1503                                ISSMERROR("%s%p%s%p"," in Rgb, missing one node. node1: ",node1," node2: ",node2);
    15041504                        }
    15051505
  • issm/trunk/src/c/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp

    r3332 r3570  
    3434        /*Checks*/
    3535        if (data_cols<=0){
    36                 ISSMERROR(exprintf("data provided has a negative number of columns"));
     36                ISSMERROR("data provided has a negative number of columns");
    3737        }
    3838        if (data_rows!=nods_data && data_rows!=nels_data){
    39                 ISSMERROR(exprintf("data provided should have either %i or %i lines (not %i)",nods_data,nels_data,data_rows));
     39                ISSMERROR("data provided should have either %i or %i lines (not %i)",nods_data,nels_data,data_rows);
    4040        }
    4141
     
    101101                        for (j=0;j<data_cols;j++){
    102102                                if (it<0 || it>=nels_data){
    103                                         ISSMERROR(exprintf("Triangle number %i not in [0 %i], because not correctly implemented yet... interpolate on grid first",it,nels_data));
     103                                        ISSMERROR("Triangle number %i not in [0 %i], because not correctly implemented yet... interpolate on grid first",it,nels_data);
    104104                                }
    105105                                data_interp[i*data_cols+j]=data[data_cols*it+j];
  • issm/trunk/src/c/ModelProcessorx/Control/CreateParametersControl.cpp

    r3538 r3570  
    4646                if (strcmp(iomodel->control_type,"drag")==0)   param->SetDouble(0);
    4747                else if (strcmp(iomodel->control_type,"B")==0) param->SetDouble(1);
    48                 else ISSMERROR(exprintf("control_type %s not supported yet!",iomodel->control_type));
     48                else ISSMERROR("control_type %s not supported yet!",iomodel->control_type);
    4949                parameters->AddObject(param);
    5050
  • issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp

    r3567 r3570  
    110110        else{
    111111
    112                 ISSMERROR(exprintf("%s%i%s%i%s"," analysis_type: ",iomodel->analysis_type," sub_analysis_type: ",iomodel->sub_analysis_type," not supported yet!"));
     112                ISSMERROR("%s%i%s%i%s"," analysis_type: ",iomodel->analysis_type," sub_analysis_type: ",iomodel->sub_analysis_type," not supported yet!");
    113113
    114114        }
  • issm/trunk/src/c/ModelProcessorx/Partitioning.cpp

    r3567 r3570  
    207207        }
    208208        else{
    209                 ISSMERROR(exprintf("not implemented yet"));
     209                ISSMERROR("not implemented yet");
    210210        }
    211211
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp

    r3567 r3570  
    6868        }
    6969        else{ //        if (strcmp(meshtype,"2d")==0)
    70                 ISSMERROR(exprintf("not implemented yet"));
     70                ISSMERROR("not implemented yet");
    7171        } //if (strcmp(meshtype,"2d")==0)
    7272
  • issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp

    r3567 r3570  
    2020
    2121        /*Output*/
    22         DataSet*       loads        = NULL;
     22        DataSet* loads=NULL;
    2323
    2424        /*numericalflux intermediary data: */
     
    2727        int  numericalflux_node_ids[MAX_NUMERICALFLUX_NODES];
    2828        int  numericalflux_elem_id;
     29        double numericalflux_h[MAX_NUMERICALFLUX_NODES];
    2930
    3031        /*Create loads: */
     
    3435        IoModelFetchData(&iomodel->edges,&iomodel->numberofedges,NULL,iomodel_handle,"edges");
    3536        IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
     37        IoModelFetchData(&iomodel->thickness,NULL,NULL,iomodel_handle,"thickness");
    3638
    3739        /*First load data:*/
     
    7274                        numericalflux_node_ids[2]=3*(int)e2+pos2;           //ex: 1 2 3
    7375                        numericalflux_node_ids[3]=3*(int)e2+((pos2+1)%3)+1; //ex: 3 1 2
     76
     77                        numericalflux_h[0]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[0]-1] -1];
     78                        numericalflux_h[1]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[1]-1]-1];
     79                        numericalflux_h[2]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[2]-1]-1];
     80                        numericalflux_h[3]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[3]-1]-1];
    7481                }
    7582                else{
     
    8794                        numericalflux_node_ids[0]=3*(int)e1+pos1;
    8895                        numericalflux_node_ids[1]=3*(int)e1+(pos1%3)+1;
     96
     97                        numericalflux_h[0]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[0]-1]-1];
     98                        numericalflux_h[1]=iomodel->thickness[(int)iomodel->elements[numericalflux_node_ids[1]-1]-1];
     99                        numericalflux_h[2]=UNDEF;
     100                        numericalflux_h[3]=UNDEF;
    89101                }
    90102
    91                 loads->AddObject(new Numericalflux(numericalflux_id,numericalflux_type,numericalflux_node_ids,numericalflux_elem_id));
     103                loads->AddObject(new Numericalflux(numericalflux_id,numericalflux_type,numericalflux_node_ids,numericalflux_elem_id,numericalflux_h));
    92104        }
    93105
     
    95107        xfree((void**)&iomodel->edges);
    96108        xfree((void**)&iomodel->elements);
     109        xfree((void**)&iomodel->thickness);
    97110
    98111        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateParametersThermal.cpp

    r3538 r3570  
    6969                param= new Param(count,"t_g",DOUBLEVEC);
    7070                if(iomodel->temperature) param->SetDoubleVec(iomodel->temperature,iomodel->numberofvertices,1);
    71                 else ISSMERROR(exprintf("Missing initial temperature"));
     71                else ISSMERROR("Missing initial temperature");
    7272                parameters->AddObject(param);
    7373
  • issm/trunk/src/c/NodeConnectivityx/NodeConnectivityx.cpp

    r3332 r3570  
    7171         * warn the user to increase the connectivity width: */
    7272        for(i=0;i<nods;i++){
    73                 if (*(connectivity+width*i+maxels)>maxels)ISSMERROR(exprintf("%s%g%s"," max connectivity width reached (",
    74                                         *(connectivity+width*i+maxels),")! increase width of connectivity table"));
     73                if (*(connectivity+width*i+maxels)>maxels)ISSMERROR("%s%g%s"," max connectivity width reached (",*(connectivity+width*i+maxels),")! increase width of connectivity table");
    7574        }
    7675
  • issm/trunk/src/c/Qmux/DakotaResponses.cpp

    r3567 r3570  
    328328                else{
    329329                        if(my_rank==0)printf("%s%s%s"," response descriptor : ",response_descriptor," not supported yet!");
    330                         ISSMERROR(exprintf("%s%s%s"," response descriptor : ",response_descriptor," not supported yet!"));
     330                        ISSMERROR("%s%s%s"," response descriptor : ",response_descriptor," not supported yet!");
    331331                }
    332332        }
  • issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp

    r3567 r3570  
    154154
    155155        }
    156         else ISSMERROR(exprintf("%s%i%s%i%s"," analysis_type ",analysis_type," and sub_analysis_type ",sub_analysis_type," not supported yet!"));
     156        else ISSMERROR("%s%i%s%i%s"," analysis_type ",analysis_type," and sub_analysis_type ",sub_analysis_type," not supported yet!");
    157157       
    158158               
  • issm/trunk/src/c/Solverx/Solverx.cpp

    r3332 r3570  
    102102        KSPGetIterationNumber(ksp,&iteration_number);
    103103        if (iteration_number<0){
    104                 ISSMERROR(exprintf("%s%i"," Solver diverged at iteration number: ",-iteration_number));
     104                ISSMERROR("%s%i"," Solver diverged at iteration number: ",-iteration_number);
    105105        }
    106106        else{
  • issm/trunk/src/c/include/macros.h

    r3499 r3570  
    2020/*Assertion macro*/
    2121#define ISSMASSERT(statement)\
    22   if (!(statement)) ISSMERROR(exprintf("Assertion \"%s\" failed, please report bug to an ISSM developer",#statement))
     22  if (!(statement)) ISSMERROR("Assertion \"%s\" failed, please report bug to an ISSM developer",#statement)
    2323
    2424/*The following macros hide the error exception handling in a matlab module. Just put
  • issm/trunk/src/c/io/FetchParams.cpp

    r3332 r3570  
    6363
    6464                        if (M==0 | N==0){
    65                                 ISSMERROR(exprintf("%s%i (%s) %s%i%s%i%s","array in parameters structure field ",count,name," is of size (",M,",",N,")"));
     65                                ISSMERROR("%s%i (%s) %s%i%s%i%s","array in parameters structure field ",count,name," is of size (",M,",",N,")");
    6666                        }
    6767
     
    147147
    148148                }
    149                 else ISSMERROR(exprintf("%s%i","unknow type in parameters structure field ",i));
     149                else ISSMERROR("%s%i","unknow type in parameters structure field ",i);
    150150        }
    151151
  • issm/trunk/src/c/io/IoModelFetchData.cpp

    r3332 r3570  
    124124        MPI_Bcast(&found,1,MPI_INT,0,MPI_COMM_WORLD);
    125125
    126         if(!found)ISSMERROR(exprintf("%s %s ","could not find data with name",data_name));
     126        if(!found)ISSMERROR("%s %s ","could not find data with name",data_name);
    127127
    128128        return fid;
  • issm/trunk/src/c/io/ParameterInputsInit.cpp

    r3332 r3570  
    8282                                break;
    8383                        default:
    84                                 ISSMERROR(exprintf("%s%i"," unknow parameter input type ",type));
     84                                ISSMERROR("%s%i"," unknow parameter input type ",type);
    8585                }
    8686
  • issm/trunk/src/c/io/WriteParams.cpp

    r3332 r3570  
    142142                                break;
    143143                        default:
    144                                 ISSMERROR(exprintf("%s%i","unknown parameter type: ",param->GetType()));
     144                                ISSMERROR("%s%i","unknown parameter type: ",param->GetType());
    145145                                break;
    146146                }
  • issm/trunk/src/c/io/pfclose.cpp

    r3332 r3570  
    1818        extern int my_rank;
    1919        if(my_rank==0){
    20                 if(fclose(fid)!=0)ISSMERROR(exprintf("%s%s","could not close file ",filename));
     20                if(fclose(fid)!=0)ISSMERROR("%s%s","could not close file ",filename);
    2121        }
    2222}
  • issm/trunk/src/c/io/pfopen.cpp

    r3332 r3570  
    2121        if(my_rank==0){
    2222                fid=fopen(filename,format);
    23                 if(fid==NULL) ISSMERROR(exprintf("%s%s%s","could not open file ",filename," for binary reading or writing"));
     23                if(fid==NULL) ISSMERROR("%s%s%s","could not open file ",filename," for binary reading or writing");
    2424        }
    2525
  • issm/trunk/src/c/objects/BamgOpts.cpp

    r3332 r3570  
    4141        int i;
    4242
    43         if (bamgopts->coef==0) ISSMERROR(exprintf("'coef' should be positive"));
    44         if (bamgopts->maxsubdiv<=1) ISSMERROR(exprintf("'maxsubdiv' should be >1"));
    45         if (bamgopts->Crack!=0  && bamgopts->Crack!=1) ISSMERROR(exprintf("'Crack' supported options are 0 and 1"));
    46         if (bamgopts->Hessiantype!=0  && bamgopts->Hessiantype!=1) ISSMERROR(exprintf("'Hessiantype' supported options are 0 and 1"));
    47         if (bamgopts->Metrictype!=0   && bamgopts->Metrictype!=1 && bamgopts->Metrictype!=2) ISSMERROR(exprintf("'Metrictype' supported options are 0, 1 and 2"));
    48         if (bamgopts->KeepVertices!=0 && bamgopts->KeepVertices!=1) ISSMERROR(exprintf("'KeepVertices' supported options are 0 and 1"));
    49         if (bamgopts->errg<0) ISSMERROR(exprintf("'errg' option should be >0"));
    50         if (bamgopts->nbjacobi<=0) ISSMERROR(exprintf("'nbjacobi' option should be >0"));
    51         if (bamgopts->geometricalmetric!=0  && bamgopts->geometricalmetric!=1) ISSMERROR(exprintf("'geometricalmetric' supported options are 0 and 1"));
    52         if (bamgopts->NbSmooth<=0) ISSMERROR(exprintf("'NbSmooth' option should be >0"));
    53         if (bamgopts->maxnbv<3) ISSMERROR(exprintf("'maxnbv' option should be >3"));
    54         if (bamgopts->hmin<=0) ISSMERROR(exprintf("'hmin' option should be >0"));
    55         if (bamgopts->hmax<=0 || bamgopts->hmax<bamgopts->hmin) ISSMERROR(exprintf("'hmax' option should be between 0 and hmin=%g",bamgopts->hmin));
    56         if (bamgopts->anisomax<1) ISSMERROR(exprintf("'anisomax' option should be >=1"));
    57         if (bamgopts->gradation<1) ISSMERROR(exprintf("'gradation' option should be >=1"));
    58         for (i=0;i<bamgopts->numfields;i++) {if (bamgopts->err[i]<=0) ISSMERROR(exprintf("'err' option should be >0"));};
     43        if (bamgopts->coef==0) ISSMERROR("'coef' should be positive");
     44        if (bamgopts->maxsubdiv<=1) ISSMERROR("'maxsubdiv' should be >1");
     45        if (bamgopts->Crack!=0  && bamgopts->Crack!=1) ISSMERROR("'Crack' supported options are 0 and 1");
     46        if (bamgopts->Hessiantype!=0  && bamgopts->Hessiantype!=1) ISSMERROR("'Hessiantype' supported options are 0 and 1");
     47        if (bamgopts->Metrictype!=0   && bamgopts->Metrictype!=1 && bamgopts->Metrictype!=2) ISSMERROR("'Metrictype' supported options are 0, 1 and 2");
     48        if (bamgopts->KeepVertices!=0 && bamgopts->KeepVertices!=1) ISSMERROR("'KeepVertices' supported options are 0 and 1");
     49        if (bamgopts->errg<0) ISSMERROR("'errg' option should be >0");
     50        if (bamgopts->nbjacobi<=0) ISSMERROR("'nbjacobi' option should be >0");
     51        if (bamgopts->geometricalmetric!=0  && bamgopts->geometricalmetric!=1) ISSMERROR("'geometricalmetric' supported options are 0 and 1");
     52        if (bamgopts->NbSmooth<=0) ISSMERROR("'NbSmooth' option should be >0");
     53        if (bamgopts->maxnbv<3) ISSMERROR("'maxnbv' option should be >3");
     54        if (bamgopts->hmin<=0) ISSMERROR("'hmin' option should be >0");
     55        if (bamgopts->hmax<=0 || bamgopts->hmax<bamgopts->hmin) ISSMERROR("'hmax' option should be between 0 and hmin=%g",bamgopts->hmin);
     56        if (bamgopts->anisomax<1) ISSMERROR("'anisomax' option should be >=1");
     57        if (bamgopts->gradation<1) ISSMERROR("'gradation' option should be >=1");
     58        for (i=0;i<bamgopts->numfields;i++) {if (bamgopts->err[i]<=0) ISSMERROR("'err' option should be >0");};
    5959
    6060}
  • issm/trunk/src/c/objects/Beam.cpp

    r3567 r3570  
    384384                }
    385385                else
    386                         ISSMERROR(exprintf("%s%i%s\n","sub_analysis_type: ",sub_analysis_type," not supported yet"));
     386                        ISSMERROR("%s%i%s\n","sub_analysis_type: ",sub_analysis_type," not supported yet");
    387387        }
    388388        else{
    389                 ISSMERROR(exprintf("%s%i%s\n","analysis: ",analysis_type," not supported yet"));
     389                ISSMERROR("%s%i%s\n","analysis: ",analysis_type," not supported yet");
    390390        }
    391391
     
    435435                }
    436436                else
    437                         ISSMERROR(exprintf("%s%i%s"," analysis ",analysis_type," not supported yet"));
     437                        ISSMERROR("%s%i%s"," analysis ",analysis_type," not supported yet");
    438438        }
    439439        else{
    440                 ISSMERROR(exprintf("%s%i%s"," analysis ",analysis_type," not supported yet"));
     440                ISSMERROR("%s%i%s"," analysis ",analysis_type," not supported yet");
    441441        }
    442442
  • issm/trunk/src/c/objects/DofVec.cpp

    r3567 r3570  
    259259        /*check we are not out of bounds: */
    260260        if (dof>=numdofs){
    261                 ISSMERROR(exprintf("%s%i%s%i\n"," error message: dof ",dof," requested is out of bounds for vector with number of dofs",this->numdofs));
     261                ISSMERROR("%s%i%s%i\n"," error message: dof ",dof," requested is out of bounds for vector with number of dofs",this->numdofs);
    262262        }
    263263
     
    284284        /*check we are not out of bounds: */
    285285        if (dof>=numdofs){
    286                 ISSMERROR(exprintf("%s%i%s%i\n"," error message: dof ",dof," requested is out of bounds for vector with number of dofs",this->numdofs));
     286                ISSMERROR("%s%i%s%i\n"," error message: dof ",dof," requested is out of bounds for vector with number of dofs",this->numdofs);
    287287        }
    288288
     
    306306       
    307307        /*check we are not out of bounds: */
    308         if (requested_numdofs>=this->numdofs) ISSMERROR(exprintf("%s%i%s%i\n"," error message: requested dof list has size  ",requested_numdofs," which is out of bounds for vector with number of dofs",this->numdofs));
     308        if (requested_numdofs>=this->numdofs) ISSMERROR("%s%i%s%i\n"," error message: requested dof list has size  ",requested_numdofs," which is out of bounds for vector with number of dofs",this->numdofs);
    309309
    310310        for(i=0;i<requested_numdofs;i++){
    311                 if (requested_doflist[i]>=this->numdofs)ISSMERROR(exprintf("%s%i%s%i\n"," error message: requested dof   ",requested_doflist[i]," in list, is out of bounds for vector with number of dofs",this->numdofs));
     311                if (requested_doflist[i]>=this->numdofs)ISSMERROR("%s%i%s%i\n"," error message: requested dof   ",requested_doflist[i]," in list, is out of bounds for vector with number of dofs",this->numdofs);
    312312        }
    313313
     
    335335       
    336336        /*check we are not out of bounds: */
    337         if (requested_numdofs>=this->numdofs) ISSMERROR(exprintf("%s%i%s%i\n"," error message: requested dof list has size  ",requested_numdofs," which is out of bounds for vector with number of dofs",this->numdofs));
     337        if (requested_numdofs>=this->numdofs) ISSMERROR("%s%i%s%i\n"," error message: requested dof list has size  ",requested_numdofs," which is out of bounds for vector with number of dofs",this->numdofs);
    338338
    339339        for(i=0;i<requested_numdofs;i++){
    340                 if (requested_doflist[i]>=this->numdofs)ISSMERROR(exprintf("%s%i%s%i\n"," error message: requested dof   ",requested_doflist[i]," in list, is out of bounds for vector with number of dofs",this->numdofs));
     340                if (requested_doflist[i]>=this->numdofs)ISSMERROR("%s%i%s%i\n"," error message: requested dof   ",requested_doflist[i]," in list, is out of bounds for vector with number of dofs",this->numdofs);
    341341        }
    342342
     
    360360
    361361        /*First thing, is the vector of the right size? :*/
    362         if (vector_size!=requested_numdofs*this->numentries)ISSMERROR(exprintf("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries));
     362        if (vector_size!=requested_numdofs*this->numentries)ISSMERROR("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries);
    363363
    364364        /*Create dofset: */
     
    385385
    386386        /*First thing, is the vector of the right size? :*/
    387         if (vector_size!=requested_numdofs*this->numentries)ISSMERROR(exprintf("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries));
     387        if (vector_size!=requested_numdofs*this->numentries)ISSMERROR("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries);
    388388
    389389        /*Create dofset: */
     
    406406
    407407        /*First thing, is the vector of the right size? :*/
    408         if (vector_size!=requested_numdofs*this->numentries)ISSMERROR(exprintf("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries));
     408        if (vector_size!=requested_numdofs*this->numentries)ISSMERROR("%s%i%s%i"," error message: input vector size ",vector_size," is not equal to the expected numdofs*numentries",requested_numdofs*this->numentries);
    409409
    410410        /*Create dofset: */
  • issm/trunk/src/c/objects/Hook.cpp

    r3567 r3570  
    216216                        this->objects[i]=(Object*)dataset->GetObjectById(this->offsets+i,this->ids[i]); //remember the offset for later on.
    217217                        /*check the id is correct!: */
    218                         if (this->objects[i]->GetId()!=this->ids[i]) ISSMERROR(exprintf("%s%i%s%i%s"," wrong id: ",this->objects[i]->GetId()," vs ",this->ids[i],"  in resolved pointer!"));
     218                        if (this->objects[i]->GetId()!=this->ids[i]) ISSMERROR("%s%i%s%i%s"," wrong id: ",this->objects[i]->GetId()," vs ",this->ids[i],"  in resolved pointer!");
    219219                }
    220220        }
     
    281281       
    282282        /*first, check that we only have one T object in our object list: */
    283         if (this->num!=1) ISSMERROR(exprintf("%s%i%s\n"," trying to delivery a single hook object when hook holds ",this->num," objects"));
     283        if (this->num!=1) ISSMERROR("%s%i%s\n"," trying to delivery a single hook object when hook holds ",this->num," objects");
    284284
    285285        /*check NULL: */
  • issm/trunk/src/c/objects/Icefront.cpp

    r3567 r3570  
    1717#include "../include/macros.h"
    1818#include "./Icefront.h"
    19 
    2019               
    2120/*Object constructors and destructor*/
     
    130129                }
    131130                else{
    132                         ISSMERROR(exprintf(" element type %i not supported yet",(int)*(iomodel->elements_type+2*element+0)));
     131                        ISSMERROR(" element type %i not supported yet",(int)*(iomodel->elements_type+2*element+0));
    133132                }
    134133        }
     
    323322
    324323                }
    325                 else ISSMERROR(exprintf("%s%i%s"," sub_analysis_type: ",sub_analysis_type," not supported yet"));
     324                else ISSMERROR("%s%i%s"," sub_analysis_type: ",sub_analysis_type," not supported yet");
    326325
    327326        }
    328327        else{
    329                 ISSMERROR(exprintf("%s%i%s"," analysis ",analysis_type," not supported yet"));
     328                ISSMERROR("%s%i%s"," analysis ",analysis_type," not supported yet");
    330329        }
    331330}
     
    14791478                }
    14801479                else{
    1481                         ISSMERROR(exprintf("fill type %i not supported yet",fill));
     1480                        ISSMERROR("fill type %i not supported yet",fill);
    14821481                }
    14831482
  • issm/trunk/src/c/objects/Input.cpp

    r3567 r3570  
    253253       
    254254        if (!valuesin) ISSMERROR("valuesin has not been allocated properly (pointer is NULL)");
    255         if (type!=DOUBLEVEC) ISSMERROR(exprintf("%s%i%s"," cannot recover values from a ",type," input type"));
     255        if (type!=DOUBLEVEC) ISSMERROR("%s%i%s"," cannot recover values from a ",type," input type");
    256256       
    257257        /*Ok, we are trying to fill values. values is of size ndof*numnodes.  The dofs
  • issm/trunk/src/c/objects/Node.cpp

    r3567 r3570  
    736736                        }
    737737                }
    738                 else ISSMERROR(exprintf("%s%s%s"," field ",field_name," not supported yet!"));
     738                else ISSMERROR("%s%s%s"," field ",field_name," not supported yet!");
    739739
    740740        } //if (extrude)
  • issm/trunk/src/c/objects/Numericalflux.cpp

    r3567 r3570  
    2626/*}}}*/
    2727/*FUNCTION Numericalflux::Numericalflux(char numericalflux_type[NUMERICALFLUXSTRING],int numericalflux_fill...){{{1*/
    28 Numericalflux::Numericalflux(int numericalflux_id,char numericalflux_type[NUMERICALFLUXSTRING], int numericalflux_node_ids[MAX_NUMERICALFLUX_NODES],int numericalflux_element_id){
     28Numericalflux::Numericalflux(int numericalflux_id,char numericalflux_type[NUMERICALFLUXSTRING], int numericalflux_node_ids[MAX_NUMERICALFLUX_NODES],int numericalflux_element_id, double numericalflux_h[MAX_NUMERICALFLUX_NODES]){
    2929
    3030        int i;
     
    4141                node_offsets[i]=UNDEF;
    4242                nodes[i]=NULL;
     43                h[i]=numericalflux_h[i];
    4344        }
    4445
     
    8182        for(i=0;i<MAX_NUMERICALFLUX_NODES;i++)nodes[i]=NULL;
    8283
     84        memcpy(&h,marshalled_dataset,sizeof(h));marshalled_dataset+=sizeof(h);
     85
    8386        /*return: */
    8487        *pmarshalled_dataset=marshalled_dataset;
     
    110113        memcpy(marshalled_dataset,&node_ids,sizeof(node_ids));marshalled_dataset+=sizeof(node_ids);
    111114        memcpy(marshalled_dataset,&node_offsets,sizeof(node_offsets));marshalled_dataset+=sizeof(node_offsets);
     115
     116        memcpy(marshalled_dataset,&h,sizeof(h));marshalled_dataset+=sizeof(h);
    112117
    113118        *pmarshalled_dataset=marshalled_dataset;
     
    124129                sizeof(node_ids)+
    125130                sizeof(node_offsets)+
     131                sizeof(h)+
    126132                sizeof(int); //sizeof(int) for enum type
    127133}
     
    443449        double mean_vy;
    444450        double UdotN;
     451        double thickness;
    445452        double dt;
    446453        int    dofs[1]={0};
     
    489496
    490497                //Get vx, vy and v.n
    491                 GetParameterValue(&vx, &vx_list[0],gauss_coord);
    492                 GetParameterValue(&vy, &vy_list[0],gauss_coord);
     498                GetParameterValue(&vx,       &vx_list[0],gauss_coord);
     499                GetParameterValue(&vy,       &vy_list[0],gauss_coord);
     500                GetParameterValue(&thickness,&h[0],      gauss_coord);
    493501                UdotN=vx*normal[0]+vy*normal[1];
    494502
     
    497505
    498506                /*Compute DL*/
    499                 DL= - gauss_weight*Jdet*dt*UdotN;
     507                DL= - gauss_weight*Jdet*dt*UdotN*thickness;
    500508
    501509                /* Add value into pe_g: */
     
    530538                        if(nodes[i])nodes[i]->Echo();
    531539                }
     540                printf("   h=[%i,%i,%i,%i]\n",h[0],h[1],h[2],h[3]);
    532541        }
    533542        else{
    534                 printf("   node_ids=[%i,%i,%i,%i]\n",node_ids[0],node_ids[1],node_ids[2],node_ids[3]);
    535                 printf("   node_offsets=[%i,%i,%i,%i]\n",node_offsets[0],node_offsets[1],node_offsets[2],node_offsets[3]);
     543                printf("   node_ids=[%i,%i]\n",node_ids[0],node_ids[1]);
     544                printf("   node_offsets=[%i,%i]\n",node_offsets[0],node_offsets[1]);
    536545                for(i=0;i<2;i++){
    537546                        if(nodes[i])nodes[i]->Echo();
    538547                }
     548                printf("   h=[%i,%i]\n",h[0],h[1]);
    539549        }
    540550        return;
     
    559569                printf("   node_ids=[%i,%i,%i,%i]\n",node_ids[0],node_ids[1],node_ids[2],node_ids[3]);
    560570                printf("   node_offsets=[%i,%i,%i,%i]\n",node_offsets[0],node_offsets[1],node_offsets[2],node_offsets[3]);
     571                printf("   h=[%g,%g,%g,%g]\n",h[0],h[1],h[2],h[3]);
    561572        }
    562573        else{
    563574                printf("   node_ids=[%i,%i]\n",node_ids[0],node_ids[1]);
    564575                printf("   node_offsets=[%i,%i]\n",node_offsets[0],node_offsets[1]);
     576                printf("   h=[%g,%g]\n",h[0],h[1]);
    565577        }
    566578
     
    615627                }
    616628        }
    617         else ISSMERROR(exprintf("type %s not supported yet",type));
     629        else ISSMERROR("type %s not supported yet",type);
    618630
    619631        /*Assign output pointers:*/
     
    741753void  Numericalflux::UpdateFromInputs(void* vinputs){
    742754
    743         /*Do nothing...*/
    744 }
    745 /*}}}*/
     755        /*Do nothing FOR NOW...*/
     756}
     757/*}}}*/
  • issm/trunk/src/c/objects/Numericalflux.h

    r3565 r3570  
    3333                int   node_offsets[MAX_NUMERICALFLUX_NODES];
    3434
     35                /*properties: */
     36                double h[MAX_NUMERICALFLUX_NODES]; //thickness
     37
    3538        public:
    3639
    3740                Numericalflux();
    38                 Numericalflux(int numericalflux_id,char numericalflux_type[NUMERICALFLUXSTRING], int numericalflux_node_ids[MAX_NUMERICALFLUX_NODES],int numericalflux_element_id);
     41                Numericalflux(int numericalflux_id,char numericalflux_type[NUMERICALFLUXSTRING], int numericalflux_node_ids[MAX_NUMERICALFLUX_NODES],int numericalflux_element_id,double numericalflux_h[MAX_NUMERICALFLUX_NODES]);
    3942                Object* copy();
    4043                ~Numericalflux();
  • issm/trunk/src/c/objects/Param.cpp

    r3567 r3570  
    2828                (param_type!=DOUBLEVEC) & (param_type!=DOUBLEMAT) & (param_type!=PETSCVEC) & (param_type!=PETSCMAT)
    2929                ){
    30                 ISSMERROR(exprintf("%s%i"," unknow parameter type ",param_type));
     30                ISSMERROR("%s%i"," unknow parameter type ",param_type);
    3131        }
    3232        stringarray=NULL;
     
    9292
    9393                default:
    94                         ISSMERROR(exprintf("%s%i","unknow parameter type ",type));
     94                        ISSMERROR("%s%i","unknow parameter type ",type);
    9595        }
    9696}
     
    210210
    211211                default:
    212                         ISSMERROR(exprintf("%s%i","unknown parameter type",type));
     212                        ISSMERROR("%s%i","unknown parameter type",type);
    213213        }
    214214
     
    311311               
    312312                default:
    313                         ISSMERROR(exprintf("%s%i","unknown parameter type",type));
     313                        ISSMERROR("%s%i","unknown parameter type",type);
    314314                        break;
    315315        }
     
    366366
    367367                default:
    368                         ISSMERROR(exprintf("%s%i","unknown parameter type",type));
     368                        ISSMERROR("%s%i","unknown parameter type",type);
    369369        }
    370370
     
    430430
    431431                default:
    432                         ISSMERROR(exprintf("%s%i","unknow parameter type ",type));
     432                        ISSMERROR("%s%i","unknow parameter type ",type);
    433433        }
    434434}
     
    490490
    491491                default:
    492                         ISSMERROR(exprintf("%s%i","unknow parameter type ",type));
     492                        ISSMERROR("%s%i","unknow parameter type ",type);
    493493        }
    494494}       
     
    518518        Mat  outmat=NULL;
    519519       
    520         if(type!=PETSCMAT)  ISSMERROR(exprintf("%s%i"," error message: trying to recover a Petsc mat from a param with type",type));
     520        if(type!=PETSCMAT)  ISSMERROR("%s%i"," error message: trying to recover a Petsc mat from a param with type",type);
    521521       
    522522        if(mat){
     
    539539void  Param::GetParameterValue(double* pdouble){
    540540
    541         if(type!=DOUBLE)ISSMERROR(exprintf("%s%i"," error message: trying to recover a double from a param with type",type));
     541        if(type!=DOUBLE)ISSMERROR("%s%i"," error message: trying to recover a double from a param with type",type);
    542542        *pdouble=ddouble;
    543543}
     
    546546void  Param::GetParameterValue(int* pinteger){
    547547       
    548         if(type!=DOUBLE)ISSMERROR(exprintf("%s%i"," error message: trying to recover an integer from a param with type",type));
     548        if(type!=DOUBLE)ISSMERROR("%s%i"," error message: trying to recover an integer from a param with type",type);
    549549        *pinteger=(int)ddouble;
    550550}
     
    555555        char*  outstring=NULL;
    556556       
    557         if(type!=STRING)ISSMERROR(exprintf("%s%i"," error message: trying to recover a string from a param with type",type));
     557        if(type!=STRING)ISSMERROR("%s%i"," error message: trying to recover a string from a param with type",type);
    558558        outstring=(char*)xmalloc((strlen(string)+1)*sizeof(char));
    559559        strcpy(outstring,string);
     
    566566
    567567        int i;
    568         if(type!=STRINGARRAY)ISSMERROR(exprintf("%s%i"," error message: trying to recover a string array from a param with type",type));
     568        if(type!=STRINGARRAY)ISSMERROR("%s%i"," error message: trying to recover a string array from a param with type",type);
    569569
    570570        char**  outstringarray=NULL;
     
    585585        double* outdoublearray=NULL;
    586586       
    587         if((type!=DOUBLEVEC) & (type!=DOUBLEMAT)) ISSMERROR(exprintf("%s%i"," error message: trying to recover a double array from a param with type",type));
     587        if((type!=DOUBLEVEC) & (type!=DOUBLEMAT)) ISSMERROR("%s%i"," error message: trying to recover a double array from a param with type",type);
    588588       
    589589        if(type==DOUBLEVEC){
     
    611611        Vec  outvec=NULL;
    612612       
    613         if(type!=PETSCVEC)  ISSMERROR(exprintf("%s%i"," error message: trying to recover a Petsc vec from a param with type",type));
     613        if(type!=PETSCVEC)  ISSMERROR("%s%i"," error message: trying to recover a Petsc vec from a param with type",type);
    614614
    615615        if(vec){
     
    636636/*FUNCTION Param::SetDouble(double value) {{{1*/
    637637void  Param::SetDouble(double value){
    638         if (type!=DOUBLE) ISSMERROR(exprintf("%s%i"," trying to set double for type",type));
     638        if (type!=DOUBLE) ISSMERROR("%s%i"," trying to set double for type",type);
    639639        ddouble=value;
    640640}
     
    642642/*FUNCTION Param::SetDouble(int value) {{{1*/
    643643void  Param::SetDouble(int value){
    644         if (type!=DOUBLE) ISSMERROR(exprintf("%s%i"," trying to set double for type",type));
     644        if (type!=DOUBLE) ISSMERROR("%s%i"," trying to set double for type",type);
    645645        ddouble=(double)value;
    646646}
     
    649649void  Param::SetDoubleMat(double* value,int pM,int pN){
    650650       
    651         if (type!=DOUBLEMAT) ISSMERROR(exprintf("%s%i"," trying to set doublematrix type",type));
     651        if (type!=DOUBLEMAT) ISSMERROR("%s%i"," trying to set doublematrix type",type);
    652652
    653653        this->M=pM;
     
    661661/*FUNCTION Param::SetDoubleVec(double* value,int size) {{{1*/
    662662void  Param::SetDoubleVec(double* value,int size){
    663         if (type!=DOUBLEVEC) ISSMERROR(exprintf("%s%i"," trying to set doublevecfor type",type));
     663        if (type!=DOUBLEVEC) ISSMERROR("%s%i"," trying to set doublevecfor type",type);
    664664       
    665665        M=size;
     
    677677/*FUNCTION Param::SetDoubleVec(double* value,int size, int numberofdofs) {{{1*/
    678678void  Param::SetDoubleVec(double* value,int size, int numberofdofs){
    679         if (type!=DOUBLEVEC) ISSMERROR(exprintf("%s%i"," trying to set doublevecfor type",type));
     679        if (type!=DOUBLEVEC) ISSMERROR("%s%i"," trying to set doublevecfor type",type);
    680680       
    681681        M=size;
     
    689689/*FUNCTION Param::SetString {{{1*/
    690690void  Param::SetString(char* value){
    691         if (type!=STRING) ISSMERROR(exprintf("%s%i"," trying to set string for type",type));
     691        if (type!=STRING) ISSMERROR("%s%i"," trying to set string for type",type);
    692692        strcpy(string,value);
    693693}
     
    697697       
    698698        int i;
    699         if (type!=STRINGARRAY) ISSMERROR(exprintf("%s%i"," trying to set string for type",type));
     699        if (type!=STRINGARRAY) ISSMERROR("%s%i"," trying to set string for type",type);
    700700        M=size;
    701701        stringarray=(char**)xmalloc(M*sizeof(char*));
  • issm/trunk/src/c/objects/ParameterInputs.cpp

    r3389 r3570  
    354354                        /*Now, pick up the parameter corresponding to root: */
    355355                        found=parameters->FindParam(&parameter,NULL,NULL,root);
    356                         if(!found)ISSMERROR(exprintf("%s%s%s"," could not find parameter ",root," for Dakota input update"));
     356                        if(!found)ISSMERROR("%s%s%s"," could not find parameter ",root," for Dakota input update");
    357357                       
    358358                        /*We've got the parameter, we need to update it using qmu_part (a partitioning vector), and the distributed_values: */
  • issm/trunk/src/c/objects/Pengrid.cpp

    r3567 r3570  
    273273        }
    274274        else{
    275                 ISSMERROR(exprintf("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet"));
     275                ISSMERROR("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet");
    276276        }
    277277
     
    377377        }
    378378        else{
    379                 ISSMERROR(exprintf("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet"));
     379                ISSMERROR("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet");
    380380        }
    381381
     
    517517        }
    518518        else{
    519                 ISSMERROR(exprintf("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet"));
     519                ISSMERROR("%s%i%s%i%s","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet");
    520520        }
    521521
  • issm/trunk/src/c/objects/Penta.cpp

    r3567 r3570  
    575575
    576576                }
    577                 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
     577                else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
    578578        }
    579579        else if (analysis_type==SlopecomputeAnalysisEnum){
     
    602602        }
    603603        else{
    604                 ISSMERROR(exprintf("%s%i%s\n","analysis: ",analysis_type," not supported yet"));
     604                ISSMERROR("%s%i%s\n","analysis: ",analysis_type," not supported yet");
    605605        }
    606606
     
    17011701                        CreatePVectorDiagnosticStokes( pg,inputs,analysis_type,sub_analysis_type);
    17021702                }
    1703                 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
     1703                else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
    17041704        }
    17051705        else if (analysis_type==SlopecomputeAnalysisEnum){
     
    17281728        }
    17291729        else{
    1730                 ISSMERROR(exprintf("%s%i%s\n","analysis: ",analysis_type," not supported yet"));
     1730                ISSMERROR("%s%i%s\n","analysis: ",analysis_type," not supported yet");
    17311731        }       
    17321732
     
    27642764
    27652765                }
    2766                 else ISSMERROR(exprintf("%s%s%s"," field ",field_name," not supported yet!"));
     2766                else ISSMERROR("%s%s%s"," field ",field_name," not supported yet!");
    27672767
    27682768        } //if (extrude)
     
    33383338
    33393339        if(*Jdet<0){
    3340                 ISSMERROR(exprintf("%s%i","negative jacobian determinant on element ",id));
     3340                ISSMERROR("%s%i","negative jacobian determinant on element ",id);
    33413341        }
    33423342}
     
    40424042                GradjB( grad_g, inputs,analysis_type,sub_analysis_type);
    40434043        }
    4044         else ISSMERROR(exprintf("%s%s","control type not supported yet: ",control_type));
     4044        else ISSMERROR("%s%s","control type not supported yet: ",control_type);
    40454045}
    40464046/*}}}*/
     
    40754075                return;
    40764076        }
    4077         else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
     4077        else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
    40784078}
    40794079/*}}}*/
  • issm/trunk/src/c/objects/Riftfront.cpp

    r3567 r3570  
    747747                else if(fill==MelangeEnum){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)
    748748                       
    749                         if(!shelf) ISSMERROR(exprintf("%s%i%s","fill type ",fill," not supported on ice sheets yet."));
     749                        if(!shelf) ISSMERROR("%s%i%s","fill type ",fill," not supported on ice sheets yet.");
    750750
    751751                        pressure_litho=rho_ice*gravity*pow(thickness,(double)2)/(double)2;
     
    757757                }
    758758                else{
    759                         ISSMERROR(exprintf("%s%i%s","fill type ",fill," not supported yet."));
     759                        ISSMERROR("%s%i%s","fill type ",fill," not supported yet.");
    760760                }
    761761
  • issm/trunk/src/c/objects/Sing.cpp

    r3567 r3570  
    354354        }
    355355        else{
    356                 ISSMERROR(exprintf("%s%i%s%i%s\n","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet"));
     356                ISSMERROR("%s%i%s%i%s\n","analysis: ",analysis_type," and sub_analysis_type: ",sub_analysis_type," not supported yet");
    357357        }
    358358
     
    391391        }
    392392        else{
    393                 ISSMERROR(exprintf("%s%i%s"," analysis ",analysis_type," not supported yet"));
     393                ISSMERROR("%s%i%s"," analysis ",analysis_type," not supported yet");
    394394        }
    395395
  • issm/trunk/src/c/objects/Tria.cpp

    r3567 r3570  
    550550                }
    551551                else{
    552                         ISSMERROR(exprintf("%s%s","unsupported control type: ",numpar->control_type));
     552                        ISSMERROR("%s%s","unsupported control type: ",numpar->control_type);
    553553                }
    554554
     
    579579                        CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type);
    580580                }
    581                 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
     581                else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
    582582
    583583        }
     
    604604        else{
    605605
    606                 ISSMERROR(exprintf("%s%i%s\n","analysis: ",analysis_type," not supported yet"));
     606                ISSMERROR("%s%i%s\n","analysis: ",analysis_type," not supported yet");
    607607        }
    608608
     
    21872187                        CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type);
    21882188                }
    2189                 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet"));
     2189                else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
    21902190        }
    21912191        else if (analysis_type==SlopecomputeAnalysisEnum){
     
    22102210        }
    22112211        else{
    2212                 ISSMERROR(exprintf("%s%i%s"," analysis ",analysis_type," not supported yet"));
     2212                ISSMERROR("%s%i%s"," analysis ",analysis_type," not supported yet");
    22132213        }
    22142214
     
    35443544        else{
    35453545                /*Not supported yet! : */
    3546                 ISSMERROR(exprintf("%s%g","unsupported type of fit: ",fit));
     3546                ISSMERROR("%s%g","unsupported type of fit: ",fit);
    35473547        }
    35483548
     
    36753675                return 1-((x-x3)*(y2-y3)-(x2-x3)*(y-y3))/area -((x1-x3)*(y-y3)-(x-x3)*(y1-y3))/area;
    36763676        }
    3677         else ISSMERROR(exprintf("%s%i%s\n"," error message: area coordinate ",which_one," done not exist!"));
     3677        else ISSMERROR("%s%i%s\n"," error message: area coordinate ",which_one," done not exist!");
    36783678}
    36793679/*}}}*/
     
    42174217                GradjB( grad_g,inputs,analysis_type,sub_analysis_type);
    42184218        }
    4219         else ISSMERROR(exprintf("%s%s","control type not supported yet: ",control_type));
     4219        else ISSMERROR("%s%s","control type not supported yet: ",control_type);
    42204220}
    42214221/*}}}*/
     
    49164916
    49174917        /*First off, check that this segment belongs to this element: */
    4918         if ((int)*(segment+4)!=this->id)ISSMERROR(exprintf("%s%i%s%i","error message: segment with id ",(int)*(segment+4)," does not belong to element with id:",this->id));
     4918        if ((int)*(segment+4)!=this->id)ISSMERROR("%s%i%s%i","error message: segment with id ",(int)*(segment+4)," does not belong to element with id:",this->id);
    49194919
    49204920        /*Recover segment node locations: */
     
    51325132        else{
    51335133                /*Not supported yet! : */
    5134                 ISSMERROR(exprintf("%s%g","unsupported type of fit: ",fit));
     5134                ISSMERROR("%s%g","unsupported type of fit: ",fit);
    51355135        }
    51365136
  • issm/trunk/src/c/parallel/WriteLockFile.cpp

    r3332 r3570  
    1717        if(my_rank==0){
    1818                fid=fopen(filename,"w");
    19                 if(fid==NULL) ISSMERROR(exprintf("%s%s","error message: could not open lock file ",filename));
     19                if(fid==NULL) ISSMERROR("%s%s","error message: could not open lock file ",filename);
    2020
    2121                /*Close file: */
    22                 if(fclose(fid)!=0) ISSMERROR(exprintf("%s%s","could not close lock file ",filename));
     22                if(fclose(fid)!=0) ISSMERROR("%s%s","could not close lock file ",filename);
    2323        }
    2424
  • issm/trunk/src/c/parallel/transient_core.cpp

    r3332 r3570  
    2626                transient_core_3d(results,model,inputs);
    2727        }
    28         else ISSMERROR(exprintf("%s%i%s"," dimension ",dim," not supported yet!"));
     28        else ISSMERROR("%s%i%s"," dimension ",dim," not supported yet!");
    2929
    3030}
  • issm/trunk/src/c/shared/Dofs/DistributeNumDofs.cpp

    r3567 r3570  
    6262                numdofs=1;
    6363        }
    64         else ISSMERROR(exprintf("%s%i%s"," analysis type: ",analysis_type,"  not implemented yet!"));
     64        else ISSMERROR("%s%i%s"," analysis type: ",analysis_type,"  not implemented yet!");
    6565
    6666        /*Assign output pointers:*/
  • issm/trunk/src/c/shared/Elements/ResolvePointers.cpp

    r3332 r3570  
    5050                        objects[i]=(Object*)dataset->GetObjectById(object_offsets+i,object_ids[i]); //remember the offset for later on.
    5151                        /*check the id is correct!: */
    52                         if (objects[i]->GetId()!=object_ids[i])ISSMERROR(exprintf("%s%i%s%i%s"," wrong id: ",objects[i]->GetId()," vs ",object_ids[i],"  in resolved pointer!"));
     52                        if (objects[i]->GetId()!=object_ids[i])ISSMERROR("%s%i%s%i%s"," wrong id: ",objects[i]->GetId()," vs ",object_ids[i],"  in resolved pointer!");
    5353                }
    5454        }
  • issm/trunk/src/c/shared/Exceptions/exprintf.cpp

    r3332 r3570  
    33 * Instead of returning an int, it will return the char* itself.
    44 * The advantage is to be able to do things like:
    5  * ISSMERROR(exprintf("%s%i\n","test failed for id:",id));
     5 * ErrorException(exprintf("%s%i\n","test failed for id:",id));
    66 */
    77
  • issm/trunk/src/c/shared/Exp/DomainOutlineRead.cpp

    r3332 r3570  
    3737        /*open domain outline file for reading: */
    3838        if ((fid=fopen(domainname,"r"))==NULL){
    39                 ISSMERROR(exprintf("%s%s","could not find domain file ",domainname));
     39                ISSMERROR("%s%s","could not find domain file ",domainname);
    4040                noerr=0; goto cleanupandreturn;
    4141        }
  • issm/trunk/src/c/shared/Matrix/MatrixUtils.cpp

    r3332 r3570  
    270270
    271271        if (!b && nvec) {
    272                 ISSMERROR(exprintf("No right-hand side for nvec=%d.",nvec));
     272                ISSMERROR("No right-hand side for nvec=%d.",nvec);
    273273                noerr=0;
    274274                return noerr;
     
    304304                        xfree((void **)&pivrc);
    305305                        xfree((void **)&pindx);
    306                         ISSMERROR(exprintf("Pivot %f less than machine epsilon",pivot));
     306                        ISSMERROR("Pivot %f less than machine epsilon",pivot);
    307307                        noerr=0;
    308308                        return noerr;
  • issm/trunk/src/c/shared/Numerics/GaussPoints.cpp

    r3466 r3570  
    393393                if (iter >= MAX_GAUS_ITER) {
    394394                        xfree((void **)&work);
    395                         ISSMERROR(exprintf("%s%i"," Max iterations exceeded for l=",MAX_GAUS_ITER));
     395                        ISSMERROR("%s%i"," Max iterations exceeded for l=",MAX_GAUS_ITER);
    396396                }
    397397        }
  • issm/trunk/src/c/toolkits/petsc/patches/NewMat.cpp

    r3332 r3570  
    5858                /*Same as above, except we bypass the default sparsity: */
    5959                if(pconnectivity || pnumberofdofspernode)
    60                         ISSMERROR(exprintf("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ",
    61                                                 psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode));
     60                        ISSMERROR("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ", psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode);
    6261
    6362                sparsity=*psparsity;
     
    7574                 *present per row: */
    7675
    77                 if(psparsity) ISSMERROR(exprintf("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ",
    78                                                 psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode));
     76                if(psparsity) ISSMERROR("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ", psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode);
    7977
    8078                connectivity=*pconnectivity;
     
    8987        }
    9088        else{
    91                 ISSMERROR(exprintf("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ",
    92                                                 psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode));
     89                ISSMERROR("%s%p%s%p%s%p","cannot determine logic with following set of pointers-> sparsity: ", psparsity," connectivity: ",pconnectivity," numberofdofspernode: ",pnumberofdofspernode);
    9390        }
    94 
    9591
    9692        return outmatrix;
  • issm/trunk/src/c/toolkits/petsc/patches/PetscOptionsDetermineSolverType.cpp

    r3332 r3570  
    7676                        if(first[0]!='-'){
    7777                                /*This is not good, the option does not have '-'! Get out*/
    78                                 ISSMERROR(exprintf("%s%s%s","Option ",first," should be preceded by '-'!"));
     78                                ISSMERROR("%s%s%s","Option ",first," should be preceded by '-'!");
    7979                        }
    8080                        break;
  • issm/trunk/src/c/toolkits/petsc/patches/PetscOptionsInsertMultipleString.cpp

    r3332 r3570  
    5555                        if(first[0]!='-'){
    5656                                /*This is not good, the option does not have '-'! Get out*/
    57                                 ISSMERROR(exprintf("%s%s%s","Option ",first," should be preceded by '-'!"));
     57                                ISSMERROR("%s%s%s","Option ",first," should be preceded by '-'!");
    5858                        }
    5959                        /*Reduce first to bare option value*/
  • issm/trunk/src/c/toolkits/petsc/patches/VecMerge.cpp

    r3332 r3570  
    6363                /*Now, modify idxm using the partition vector, and plug values into A*/
    6464                for (i=0;i<range;i++){
    65                         *(idxm+i)=*(row_partition_vector+lower_row+i)-1; //-1 because partition vector comes from Matlab, where indices start at 1.
     65                        *(idxm+i)=(int)*(row_partition_vector+lower_row+i)-1; //-1 because partition vector comes from Matlab, where indices start at 1.
    6666                }
    6767                VecSetValues(A,range,idxm,values,INSERT_VALUES);
Note: See TracChangeset for help on using the changeset viewer.