Changeset 16152
- Timestamp:
- 09/17/13 14:35:17 (12 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/bamg/Mesh.cpp
r15219 r16152 122 122 } 123 123 for (i=0;i<Tho.nbt;i++) 124 if( reft[i] >=0 && flag[i]) 125 { 124 if(reft[i] >=0 && flag[i]){ 126 125 const Triangle & t = Tho.triangles[i]; 127 126 int i0 = Tho.GetId(t[0]); 128 127 int i1 = Tho.GetId(t[1]); 129 128 int i2 = Tho.GetId(t[2]); 130 if 129 if(i0<0 || i1<0 || i2<0){ 131 130 delete [] refv; 132 131 _error_("i0<0 || i1<0 || i2< 0"); 133 132 } 134 if (i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv){ 133 if(i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv){ 134 delete [] refv; 135 135 _error_("i0>=Tho.nbv || i1>=Tho.nbv || i2>=Tho.nbv"); 136 136 } … … 4138 4138 first_np_or_next_t1[i]=-1; 4139 4139 kk=0; 4140 while (Head0>=0&& kk++<100){4140 while(Head0>=0&& kk++<100){ 4141 4141 kch=0; 4142 for 4142 for(i=Head0;i>=0;i=first_np_or_next_t0[ip=i],first_np_or_next_t0[ip]=-1) { 4143 4143 // pour tous les triangles autour du sommet s 4144 4144 register Triangle* t= vertices[i].t; -
issm/trunk-jpl/src/c/shared/Bamg/HeapSort.h
r15072 r16152 85 85 } 86 86 } 87 88 /*Make cppcheck happy*/ 89 *porder=order; 87 90 }/*}}}*/ 88 91 -
issm/trunk-jpl/src/c/shared/Exceptions/Exceptions.cpp
r15839 r16152 15 15 ErrorException::ErrorException(const string &what_arg){/*{{{*/ 16 16 17 what_str=what_arg; 18 file_name=""; 19 function_name=""; 20 file_line=0; 17 what_str = what_arg; 18 file_name = ""; 19 function_name = ""; 20 file_line = 0; 21 21 22 }/*}}}*/ 22 23 ErrorException::ErrorException(const string& what_file, const string& what_function,int what_line, const string& what_arg){/*{{{*/ 23 24 24 what_str=what_arg; 25 file_name=what_file; 26 function_name=what_function; 27 file_line=what_line; 25 what_str = what_arg; 26 file_name = what_file; 27 function_name = what_function; 28 file_line = what_line; 29 28 30 }/*}}}*/ 29 31 ErrorException::~ErrorException() throw(){/*{{{*/ -
issm/trunk-jpl/src/c/shared/Exceptions/exceptions.h
r15104 r16152 83 83 public: 84 84 ErrorException(const string &what_arg); //for windows 85 ErrorException(const string &what_file,const string& what_function,int what_line,const string& what_arg);//UNIX85 ErrorException(const string &what_file,const string& what_function,int what_line,const string& what_arg);//UNIX 86 86 ~ErrorException() throw(); 87 87 virtual const char *what() const throw(); -
issm/trunk-jpl/src/c/shared/LatLong/Ll2xyx.cpp
r15741 r16152 46 46 int i,iret=0; 47 47 double delta,slat; 48 double cde,re,ex2,ex;48 double re,ex2,ex; 49 49 double latitude,longitude; 50 50 double T,rho,sl,tc,mc; … … 55 55 slat = standard_parallel; 56 56 57 /* Conversion constant from degrees to radians */58 cde = 57.29577951;59 57 /* Radius of the earth in meters */ 60 58 re = 6378.273*1.e3; -
issm/trunk-jpl/src/c/shared/LatLong/Xy2llx.cpp
r15741 r16152 48 48 int i,iret=0; 49 49 double delta,slat; 50 double cde,re,ex2,ex;50 double re,ex2,ex; 51 51 double sl,rho,cm,T,chi; 52 52 … … 56 56 slat = standard_parallel; 57 57 58 /* Conversion constant from degrees to radians */59 cde = 57.29577951;60 58 /* Radius of the earth in meters */ 61 59 re = 6378.273e+3; -
issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp
r15741 r16152 45 45 static IssmPDouble* xip [MAX_LINE_GAUS_PTS]={xi1 ,xi2 ,xi3 ,xi4 }; 46 46 47 static int np[MAX_LINE_GAUS_PTS]={sizeof(wgt1 )/sizeof(IssmPDouble),48 sizeof(wgt2 )/sizeof(IssmPDouble),49 sizeof(wgt3 )/sizeof(IssmPDouble),50 sizeof(wgt4 )/sizeof(IssmPDouble)};51 52 // _printf0_("Gauss-Legendre recurrence coefficients ngaus=" << ngaus << "\n");53 47 *pxgaus =xNew<IssmPDouble>(ngaus); 54 48 *pxwgt =xNew<IssmPDouble>(ngaus); 55 49 56 50 /* check to see if Gauss points need to be calculated */ 57 if 51 if(ngaus <= MAX_LINE_GAUS_PTS) { 58 52 59 53 /* copy the points from the static arrays (noting that the pointers … … 61 55 have to know to not free them) */ 62 56 63 for 57 for(i=0; i<ngaus; i++) { 64 58 (*pxgaus)[i]=xip [ngaus-1][i]; 65 59 (*pxwgt )[i]=wgtp[ngaus-1][i]; 66 60 } 67 61 } 68 else 62 else{ 69 63 70 64 /* calculate the Gauss points using recurrence relations */ … … 1259 1253 0.250000000000000, 0.500000000000000, 0.166666666666667, 1260 1254 0.166666666666667, 0.166666666666667}; 1261 static IssmPDouble l23[]={1262 0.250000000000000, 0.166666666666667, 0.500000000000000,1263 0.166666666666667, 0.166666666666667};1264 1255 static IssmPDouble l33[]={ 1265 1256 0.250000000000000, 0.166666666666667, 0.166666666666667, … … 1511 1502 static IssmPDouble* xip [MAX_LINE_GLOB_PTS]={xi1 ,xi2 ,xi3 ,xi4 ,xi5 }; 1512 1503 1513 static int np[MAX_LINE_GLOB_PTS]={sizeof(wgt1 )/sizeof(IssmPDouble),1514 sizeof(wgt2 )/sizeof(IssmPDouble),1515 sizeof(wgt3 )/sizeof(IssmPDouble),1516 sizeof(wgt4 )/sizeof(IssmPDouble),1517 sizeof(wgt5 )/sizeof(IssmPDouble)};1518 1519 // _printf0_("Gauss-Lobatto recurrence coefficients ngaus=" << ngaus << "\n");1520 1504 *pxgaus =xNew<IssmPDouble>(ngaus); 1521 1505 *pxwgt =xNew<IssmPDouble>(ngaus); -
issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h
r16144 r16152 504 504 IssmAbsVec<IssmDouble>* Solve(IssmAbsVec<IssmDouble>* pfin, Parameters* parameters){ 505 505 506 int i;507 508 506 /*output: */ 509 507 IssmMpiVec<IssmDouble>* uf=NULL; … … 523 521 #endif 524 522 return (IssmAbsVec<IssmDouble>*)uf; 525 break; 526 } 523 } 527 524 case GslEnum: { 528 525 … … 535 532 536 533 return (IssmAbsVec<IssmDouble>*)uf; 537 break;538 534 #else 539 535 _error_("GSL support not compiled in!"); 540 536 #endif 541 }537 } 542 538 default: 543 539 _error_("solver type not supported yet!"); -
issm/trunk-jpl/src/c/toolkits/mumps/MpiDenseMumpsSolve.cpp
r16149 r16152 110 110 /*Variables*/ 111 111 ISSM_MPI_Comm comm; 112 int my_rank;113 112 int num_procs; 114 113 int i,j; … … 125 124 126 125 /*Communicator info */ 127 my_rank = IssmComm::GetRank();128 126 num_procs = IssmComm::GetSize(); 129 127 comm = IssmComm::GetComm();
Note:
See TracChangeset
for help on using the changeset viewer.