Changeset 1904


Ignore:
Timestamp:
08/25/09 17:56:24 (15 years ago)
Author:
Eric.Larour
Message:

Valgrind memory checker now working.

Location:
issm/trunk/src
Files:
44 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/ContourToMeshx/ContourToMeshx.cpp

    r291 r1904  
    3030        /*Loop through all contours: */
    3131        for (i=0;i<numcontours;i++){
    32                 #ifdef _DEBUG_
     32                #ifdef _ISSM_DEBUG_
    3333                        printf("\nHandling contour %i/%i\n",i,numcontours);
    3434                #endif

  • issm/trunk/src/c/ContourToNodesx/ContourToNodesx.cpp

    r1105 r1904  
    2626        /*Loop through all contours: */
    2727        for (i=0;i<numcontours;i++){
    28                 #ifdef _DEBUG_
     28                #ifdef _ISSM_DEBUG_
    2929                        printf("\nHandling contour %i/%i\n",i,numcontours);
    3030                #endif

  • issm/trunk/src/c/DataSet/DataSet.cpp

    r1898 r1904  
    214214        }
    215215
    216         #ifdef _DEBUG_
     216        #ifdef _ISSM_DEBUG_
    217217        _printf_("Number of objects in dataset being demarshalled: %i\n",numobjects);
    218218        #endif
    … …  
    641641        NodeRank(ranks);
    642642
    643         #ifdef _DEBUG_
     643        #ifdef _ISSM_DEBUG_
    644644        for(i=0;i<numberofnodes;i++){
    645645                _printf_("%i\n",ranks[i]);
    … …  
    653653        MPI_Allreduce ( (void*)ranks,(void*)minranks,numberofnodes,MPI_INT,MPI_MIN,MPI_COMM_WORLD);
    654654
    655         #ifdef _DEBUG_
     655        #ifdef _ISSM_DEBUG_
    656656        for(i=0;i<numberofnodes;i++){
    657657                _printf_("%i\n",minranks[i]);

  • issm/trunk/src/c/GriddataMeshToGridx/GriddataMeshToGridx.cpp

    r819 r1904  
    137137        }
    138138        for (i=0;i<numcontours;i++){
    139                 #ifdef _DEBUG_
     139                #ifdef _ISSM_DEBUG_
    140140                        _printf_("Handling contour %i/%i\n",i,numcontours);
    141141                #endif

  • issm/trunk/src/c/HoleFillerx/HoleFillerx.cpp

    r1454 r1904  
    4444        int             imageoutsize;
    4545
    46         #ifdef _DEBUG_
     46        #ifdef _ISSM_DEBUG_
    4747                if ( smooth == 1 ){
    4848                        printf("Data patches will be SMOOTHED. \n");
    … …  
    9393        #endif
    9494
    95         #ifdef _DEBUG_
     95        #ifdef _ISSM_DEBUG_
    9696                printf( "\n" );
    9797                printf( "Iterations %5ld   and %5ld.", count-1, count );
    … …  
    119119        #endif
    120120
    121         #ifdef _DEBUG_
     121        #ifdef _ISSM_DEBUG_
    122122                printf( "\b\b\b\b\b\b\b\b\b\b\b\b\b\b" );
    123123                printf( "\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b" );
    … …  
    367367                }
    368368
    369                 #ifdef _DEBUG_
     369                #ifdef _ISSM_DEBUG_
    370370                        if ( !(i%100) ){
    371371                        fprintf(stdout, "\b\b\b\b\b\b%5ld ", i );   /* Count lines on output screen */
    … …  
    375375        }
    376376
    377         #ifdef _DEBUG_
     377        #ifdef _ISSM_DEBUG_
    378378                fprintf(stdout, "\b\b\b\b\b\b%5ld ", lines );
    379379                fflush(stdout);
    … …  
    389389
    390390        time(&t2);     
    391         #ifdef _DEBUG_
     391        #ifdef _ISSM_DEBUG_
    392392                printf ( "\n\nEnd                  "); printf( ctime(&t2) );
    393393                printf( "\n                           Total time:        ");
    … …  
    403403        end:
    404404
    405         #ifdef _DEBUG_
     405        #ifdef _ISSM_DEBUG_
    406406                printf( "Done.\a\n" );
    407407        #endif

  • issm/trunk/src/c/Makefile.am

    r1901 r1904  
    302302
    303303
    304 libISSM_a_CXXFLAGS = -fPIC -DMATLAB -D_SERIAL_ -ansi -D_GNU_SOURCE -fno-omit-frame-pointer -pthread -D_CPP_
     304libISSM_a_CXXFLAGS = -fPIC -DMATLAB -D_SERIAL_ -ansi -D_GNU_SOURCE -fno-omit-frame-pointer -pthread -D_CPP_  -g -O0
    305305if LARGEARRAYS
    306306libISSM_a_CXXFLAGS += -D__GCC4BUILD__ 

  • issm/trunk/src/c/MeshPartitionx/MeshPartitionx.cpp

    r582 r1904  
    8888        }
    8989
    90         #ifdef _DEBUG_
     90        #ifdef _ISSM_DEBUG_
    9191        if(my_rank==0){
    9292                for (i=0;i<numberofelements;i++){

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

    r1832 r1904  
    505505                VecAssemblyEnd(gridborder);
    506506
    507                 #ifdef _DEBUG_
     507                #ifdef _ISSM_DEBUG_
    508508                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    509509                #endif
    … …  
    511511                VecToMPISerial(&my_bordergrids,gridborder);
    512512
    513                 #ifdef _DEBUG_
     513                #ifdef _ISSM_DEBUG_
    514514                if(my_rank==0){
    515515                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp

    r1835 r1904  
    184184        VecAssemblyEnd(gridborder);
    185185
    186         #ifdef _DEBUG_
     186        #ifdef _ISSM_DEBUG_
    187187        VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    188188        #endif
    … …  
    190190        VecToMPISerial(&my_bordergrids,gridborder);
    191191
    192         #ifdef _DEBUG_
     192        #ifdef _ISSM_DEBUG_
    193193        if(my_rank==0){
    194194                for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateElementsNodesAndMaterialsDiagnosticStokes.cpp

    r1891 r1904  
    334334                VecAssemblyEnd(gridborder);
    335335
    336                 #ifdef _DEBUG_
     336                #ifdef _ISSM_DEBUG_
    337337                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    338338                #endif
    … …  
    340340                VecToMPISerial(&my_bordergrids,gridborder);
    341341
    342                 #ifdef _DEBUG_
     342                #ifdef _ISSM_DEBUG_
    343343                if(my_rank==0){
    344344                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateElementsNodesAndMaterialsDiagnosticVert.cpp

    r1834 r1904  
    284284                VecAssemblyEnd(gridborder);
    285285
    286                 #ifdef _DEBUG_
     286                #ifdef _ISSM_DEBUG_
    287287                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    288288                #endif
    … …  
    290290                VecToMPISerial(&my_bordergrids,gridborder);
    291291
    292                 #ifdef _DEBUG_
     292                #ifdef _ISSM_DEBUG_
    293293                if(my_rank==0){
    294294                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

    r1834 r1904  
    301301                VecAssemblyEnd(gridborder);
    302302
    303                 #ifdef _DEBUG_
     303                #ifdef _ISSM_DEBUG_
    304304                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    305305                #endif
    … …  
    307307                VecToMPISerial(&my_bordergrids,gridborder);
    308308
    309                 #ifdef _DEBUG_
     309                #ifdef _ISSM_DEBUG_
    310310                if(my_rank==0){
    311311                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp

    r1834 r1904  
    350350                VecAssemblyEnd(gridborder);
    351351
    352                 #ifdef _DEBUG_
     352                #ifdef _ISSM_DEBUG_
    353353                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    354354                #endif
    … …  
    356356                VecToMPISerial(&my_bordergrids,gridborder);
    357357
    358                 #ifdef _DEBUG_
     358                #ifdef _ISSM_DEBUG_
    359359                if(my_rank==0){
    360360                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/Qmu/CreateParametersQmu.cpp

    r1836 r1904  
    146146                parameters->AddObject(param);
    147147
    148                 #ifdef _DEBUG_
     148                #ifdef _ISSM_DEBUG_
    149149                        for(i=0;i<iomodel->numberofvariables;i++){
    150150                                _printf_("variable descriptor %s\n",variabledescriptors[i]);

  • issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateElementsNodesAndMaterialsSlopeCompute.cpp

    r1833 r1904  
    312312                VecAssemblyEnd(gridborder);
    313313
    314                 #ifdef _DEBUG_
     314                #ifdef _ISSM_DEBUG_
    315315                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    316316                #endif
    … …  
    318318                VecToMPISerial(&my_bordergrids,gridborder);
    319319
    320                 #ifdef _DEBUG_
     320                #ifdef _ISSM_DEBUG_
    321321                if(my_rank==0){
    322322                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

    r1834 r1904  
    304304                VecAssemblyEnd(gridborder);
    305305
    306                 #ifdef _DEBUG_
     306                #ifdef _ISSM_DEBUG_
    307307                VecView(gridborder,PETSC_VIEWER_STDOUT_WORLD);
    308308                #endif
    … …  
    310310                VecToMPISerial(&my_bordergrids,gridborder);
    311311
    312                 #ifdef _DEBUG_
     312                #ifdef _ISSM_DEBUG_
    313313                if(my_rank==0){
    314314                        for (i=0;i<iomodel->numberofnodes;i++){

  • issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp

    r1805 r1904  
    263263        #endif
    264264
    265         #ifdef _DEBUG_
     265        #ifdef _ISSM_DEBUG_
    266266                printf("Maximum rift penetration2: %g\n",max_penetration);
    267267        #endif

  • issm/trunk/src/c/PenaltySystemMatricesx/PenaltySystemMatricesx.cpp

    r586 r1904  
    2828        /*Now, figure out maximum value of K_gg, so that we can penalize it correctly: */
    2929        MatNorm(Kgg,NORM_INFINITY,&kmax);
    30         #ifdef _DEBUG_
     30        #ifdef _ISSM_DEBUG_
    3131         _printf_("   K_gg infinity norm: %g\n",kmax);
    3232        #endif
    … …  
    4747        }
    4848
    49         #ifdef _DEBUG_
     49        #ifdef _ISSM_DEBUG_
    5050        MatNorm(Kgg,NORM_INFINITY,&kmax2);
    5151         _printf_("   K_gg infinity norm after penalties: %g\n",kmax2);

  • issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp

    r1881 r1904  
    6969        model->FindParam(&qmu_npart,"qmu_npart");
    7070        model->FindParam(&qmu_part,"qmu_part");
    71         #ifdef _DEBUG_
     71        #ifdef _ISSM_DEBUG_
    7272        for(i=0;i<numresponses;i++){
    7373                PetscSynchronizedPrintf(MPI_COMM_WORLD,"response descriptor %i: %s\n",i,responses_descriptors[i]);
    … …  
    8181        MPI_Bcast(variables,numvariables,MPI_DOUBLE,0,MPI_COMM_WORLD);
    8282
    83         #ifdef _DEBUG_
     83        #ifdef _ISSM_DEBUG_
    8484        for(i=0;i<numvariables;i++){
    8585                PetscSynchronizedPrintf(MPI_COMM_WORLD,"variable %i: %g\n",i,variables[i]);
    … …  
    103103        }
    104104
    105         #ifdef _DEBUG_
     105        #ifdef _ISSM_DEBUG_
    106106        for(i=0;i<numvariables;i++){
    107107                PetscSynchronizedPrintf(MPI_COMM_WORLD,"variable descriptor %i: %s value: %g\n",i,variables_descriptors[i],variables[i]);

  • issm/trunk/src/c/Solverx/Solverx.cpp

    r1 r1904  
    7575               
    7676
    77         #ifdef _DEBUG_
     77        #ifdef _ISSM_DEBUG_
    7878                KSPView(ksp,PETSC_VIEWER_STDOUT_WORLD);
    7979                KSPGetInitialGuessNonzero(ksp,&flag);
    … …  
    8989        }
    9090        else{
    91                 #ifdef _DEBUG_
     91                #ifdef _ISSM_DEBUG_
    9292                _printf_("%s%i%s\n","Solver converged after ",iteration_number," iterations");
    9393                #endif

  • issm/trunk/src/c/objects/Matice.cpp

    r1439 r1904  
    181181                }
    182182        }
    183         #ifdef _DEBUG_
    184         _printf_("Viscosity %lf\n",viscosity);
     183        #ifdef _ISSM_DEBUG_
     184        printf("Viscosity %lf\n",viscosity);
    185185        #endif
    186186
    … …  
    241241        }
    242242               
    243         #ifdef _DEBUG_
    244         _printf_("viscosity_complement %lf\n",mu2);
     243        #ifdef _ISSM_DEBUG_
     244        printf("viscosity_complement %lf\n",viscosity_complement);
    245245        #endif
    246246
    … …  
    306306                }
    307307        }
    308         #ifdef _DEBUG_
    309         _printf_("Viscosity %lf\n",viscosity3d);
     308        #ifdef _ISSM_DEBUG_
     309        printf("Viscosity %lf\n",viscosity3d);
    310310        #endif
    311311
    … …  
    375375        }
    376376
    377         #ifdef _DEBUG_
    378         _printf_("Viscosity %lf\n",viscosity3d);
     377        #ifdef _ISSM_DEBUG_
     378        printf("Viscosity %lf\n",viscosity3d);
    379379        #endif
    380380

  • issm/trunk/src/c/objects/ParameterInputs.cpp

    r962 r1904  
    329329                        }
    330330
    331                         #ifdef _DEBUG_
     331                        #ifdef _ISSM_DEBUG_
    332332                                PetscSynchronizedPrintf(MPI_COMM_WORLD,"Parameter vetor:");
    333333                                PetscSynchronizedFlush(MPI_COMM_WORLD);

  • issm/trunk/src/c/objects/Penta.cpp

    r1901 r1904  
    710710
    711711        GaussPenta( &num_area_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &area_gauss_weights, &fourth_gauss_vert_coord,&vert_gauss_weights,order_area_gauss,num_vert_gauss);
    712         #ifdef _DEBUG_
     712        #ifdef _ISSM_DEBUG_
    713713        for (i=0;i<num_area_gauss;i++){
    714                 _printf_("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
     714                printf("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
    715715        }
    716716        for (i=0;i<num_vert_gauss;i++){
    717                 _printf_("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
     717                printf("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
    718718        }
    719719        #endif
    … …  
    14341434        GetB(&B[0][0], xyz_list, gauss_l1l2l3l4);
    14351435
    1436         #ifdef _DEBUG_
    1437         _printf_("B for grid1 : [ %lf   %lf  \n",B[0][0],B[0][1]);
    1438         _printf_("              [ %lf   %lf  \n",B[1][0],B[1][1]);
    1439         _printf_("              [ %lf   %lf  ]\n",B[2][0],B[2][1]);
    1440         _printf_("              [ %lf   %lf  ]\n",B[3][0],B[3][1]);
    1441         _printf_("              [ %lf   %lf  ]\n",B[4][0],B[4][1]);
    1442        
    1443         _printf_("B for grid2 : [ %lf   %lf  \n",B[0][2],B[0][3]);
    1444         _printf_("              [ %lf   %lf  \n",B[1][2],B[1][3]);
    1445         _printf_("              [ %lf   %lf  ]\n",B[2][2],B[2][3]);
    1446         _printf_("              [ %lf   %lf  ]\n",B[3][2],B[3][3]);
    1447         _printf_("              [ %lf   %lf  ]\n",B[4][2],B[4][3]);
    1448        
    1449         _printf_("B for grid3 : [ %lf   %lf  \n", B[0][4],B[0][5]);
    1450         _printf_("              [ %lf   %lf  \n", B[1][4],B[1][5]);
    1451         _printf_("              [ %lf   %lf  ]\n",B[2][4],B[2][5]);
    1452         _printf_("              [ %lf   %lf  ]\n",B[3][4],B[3][5]);
    1453         _printf_("              [ %lf   %lf  ]\n",B[4][4],B[4][5]);
    1454        
    1455         _printf_("B for grid4 : [ %lf   %lf  \n", B[0][6],B[0][7]);
    1456         _printf_("              [ %lf   %lf  \n", B[1][6],B[1][7]);
    1457         _printf_("              [ %lf   %lf  ]\n",B[2][6],B[2][7]);
    1458         _printf_("              [ %lf   %lf  ]\n",B[3][6],B[3][7]);
    1459         _printf_("              [ %lf   %lf  ]\n",B[4][6],B[4][7]);
     1436        #ifdef _ISSM_DEBUG_
     1437        printf("B for grid1 : [ %lf   %lf  \n",B[0][0],B[0][1]);
     1438        printf("              [ %lf   %lf  \n",B[1][0],B[1][1]);
     1439        printf("              [ %lf   %lf  ]\n",B[2][0],B[2][1]);
     1440        printf("              [ %lf   %lf  ]\n",B[3][0],B[3][1]);
     1441        printf("              [ %lf   %lf  ]\n",B[4][0],B[4][1]);
     1442       
     1443        printf("B for grid2 : [ %lf   %lf  \n",B[0][2],B[0][3]);
     1444        printf("              [ %lf   %lf  \n",B[1][2],B[1][3]);
     1445        printf("              [ %lf   %lf  ]\n",B[2][2],B[2][3]);
     1446        printf("              [ %lf   %lf  ]\n",B[3][2],B[3][3]);
     1447        printf("              [ %lf   %lf  ]\n",B[4][2],B[4][3]);
     1448       
     1449        printf("B for grid3 : [ %lf   %lf  \n", B[0][4],B[0][5]);
     1450        printf("              [ %lf   %lf  \n", B[1][4],B[1][5]);
     1451        printf("              [ %lf   %lf  ]\n",B[2][4],B[2][5]);
     1452        printf("              [ %lf   %lf  ]\n",B[3][4],B[3][5]);
     1453        printf("              [ %lf   %lf  ]\n",B[4][4],B[4][5]);
     1454       
     1455        printf("B for grid4 : [ %lf   %lf  \n", B[0][6],B[0][7]);
     1456        printf("              [ %lf   %lf  \n", B[1][6],B[1][7]);
     1457        printf("              [ %lf   %lf  ]\n",B[2][6],B[2][7]);
     1458        printf("              [ %lf   %lf  ]\n",B[3][6],B[3][7]);
     1459        printf("              [ %lf   %lf  ]\n",B[4][6],B[4][7]);
    14601460                               
    1461         _printf_("B for grid5 : [ %lf   %lf  \n", B[0][8],B[0][9]);
    1462         _printf_("              [ %lf   %lf  \n", B[1][8],B[1][9]);
    1463         _printf_("              [ %lf   %lf  ]\n",B[2][8],B[2][9]);
    1464         _printf_("              [ %lf   %lf  ]\n",B[3][8],B[3][9]);
    1465         _printf_("              [ %lf   %lf  ]\n",B[4][8],B[4][9]);
    1466 
    1467         _printf_("B for grid6 : [ %lf   %lf  \n", B[0][10],B[0][11]);
    1468         _printf_("              [ %lf   %lf  \n", B[1][10],B[1][11]);
    1469         _printf_("              [ %lf   %lf  ]\n",B[2][10],B[2][11]);
    1470         _printf_("              [ %lf   %lf  ]\n",B[3][10],B[3][11]);
    1471         _printf_("              [ %lf   %lf  ]\n",B[4][10],B[4][11]);
    1472 
    1473         for (i=0;i<numgrids;i++){
    1474                 _printf_("Velocity for grid %i %lf %lf\n",i,*(vxvy_list+2*i+0),*(vxvy_list+2*i+1));
    1475         }
     1461        printf("B for grid5 : [ %lf   %lf  \n", B[0][8],B[0][9]);
     1462        printf("              [ %lf   %lf  \n", B[1][8],B[1][9]);
     1463        printf("              [ %lf   %lf  ]\n",B[2][8],B[2][9]);
     1464        printf("              [ %lf   %lf  ]\n",B[3][8],B[3][9]);
     1465        printf("              [ %lf   %lf  ]\n",B[4][8],B[4][9]);
     1466
     1467        printf("B for grid6 : [ %lf   %lf  \n", B[0][10],B[0][11]);
     1468        printf("              [ %lf   %lf  \n", B[1][10],B[1][11]);
     1469        printf("              [ %lf   %lf  ]\n",B[2][10],B[2][11]);
     1470        printf("              [ %lf   %lf  ]\n",B[3][10],B[3][11]);
     1471        printf("              [ %lf   %lf  ]\n",B[4][10],B[4][11]);
     1472
    14761473        #endif
    14771474
    … …  
    15101507        GetNodalFunctionsDerivativesBasic(&dh1dh2dh3dh4dh5dh6_basic[0][0],xyz_list, gauss_l1l2l3l4);
    15111508
    1512         #ifdef _DEBUG_
     1509        #ifdef _ISSM_DEBUG_
    15131510        for (i=0;i<numgrids;i++){
    1514                 _printf_("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
     1511                printf("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
    15151512        }
    15161513        #endif
    … …  
    15641561        GetNodalFunctionsDerivativesBasic(&dh1dh2dh3dh4dh5dh6_basic[0][0],xyz_list, gauss_l1l2l3l4);
    15651562
    1566         #ifdef _DEBUG_
     1563        #ifdef _ISSM_DEBUG_
    15671564        for (i=0;i<numgrids;i++){
    1568                 _printf_("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
     1565                printf("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
    15691566        }
    15701567        #endif
    … …  
    17041701        *(J+NDOF3*2+2)=sqrt3/12.0*(z1+z2-2*z3-z4-z5+2*z6)*eta+1.0/4.0*(z1-z2-z4+z5)*xi+1.0/4.0*(-z1+z5-z2+z4);
    17051702
    1706         #ifdef _DEBUG_
     1703        #ifdef _ISSM_DEBUG_
    17071704        for(i=0;i<3;i++){
    17081705                for (j=0;j<3;j++){
    … …  
    18691866
    18701867                GaussPenta( &num_area_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &area_gauss_weights, &fourth_gauss_vert_coord,&vert_gauss_weights,order_area_gauss,num_vert_gauss);
    1871                 #ifdef _DEBUG_
     1868                #ifdef _ISSM_DEBUG_
    18721869                for (i=0;i<num_area_gauss;i++){
    1873                         _printf_("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
     1870                        printf("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
    18741871                }
    18751872                for (i=0;i<num_vert_gauss;i++){
    1876                         _printf_("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
     1873                        printf("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
    18771874                }
    18781875                #endif
    … …  
    21512148        GetNodalFunctionsDerivativesBasic(&dh1dh2dh3dh4dh5dh6_basic[0][0],xyz_list, gauss_l1l2l3l4);
    21522149
    2153         #ifdef _DEBUG_
     2150        #ifdef _ISSM_DEBUG_
    21542151        for (i=0;i<numgrids;i++){
    2155                 _printf_("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
     2152                printf("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf\n",i,dh1dh2dh3dh4dh5dh6_basic[0][i],dh1dh2dh3dh4dh5dh6_basic[1][i],dh1dh2dh3dh4dh5dh6_basic[2][i]);
    21562153        }
    21572154        #endif
    … …  
    22612258
    22622259        GaussPenta( &num_area_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &area_gauss_weights, &fourth_gauss_vert_coord,&vert_gauss_weights,order_area_gauss,num_vert_gauss);
    2263         #ifdef _DEBUG_
     2260        #ifdef _ISSM_DEBUG_
    22642261        for (i=0;i<num_area_gauss;i++){
    2265                 _printf_("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
     2262                printf("Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
    22662263        }
    22672264        for (i=0;i<num_vert_gauss;i++){
    2268                 _printf_("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
     2265                printf("Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
    22692266        }
    22702267        #endif
    … …  
    22932290                        /* Get Jacobian determinant: */
    22942291                        GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss_l1l2l3l4);
    2295                         #ifdef _DEBUG_
    2296                         _printf_("Element id %i Jacobian determinant: %lf\n",PentaElementGetID(this),Jdet);
     2292                        #ifdef _ISSM_DEBUG_
     2293                        printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    22972294                        #endif
    22982295               
    … …  
    26292626        GetNodalFunctions(l1l6, gauss_coord);
    26302627       
    2631         #ifdef _DEBUG_
     2628        #ifdef _ISSM_DEBUG_
    26322629        for (i=0;i<7;i++){
    26332630                printf("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf \n",i,dh1dh7_basic[0][i],dh1dh7_basic[1][i],dh1dh7_basic[2][i]);
    … …  
    27112708        GetNodalFunctions(l1l6, gauss_coord);
    27122709       
    2713         #ifdef _DEBUG_
     2710        #ifdef _ISSM_DEBUG_
    27142711        for (i=0;i<6;i++){
    27152712                printf("Node %i  dh/dx=%lf dh/dy=%lf dh/dz=%lf \n",i,dh1dh7_basic[0][i],dh1dh7_basic[1][i]);

  • issm/trunk/src/c/objects/Riftfront.cpp

    r1881 r1904  
    341341                 *contact slip friction. */
    342342                 
    343                 #ifdef _DEBUG_
     343                #ifdef _ISSM_DEBUG_
    344344                printf("Dealing with grid pair (%i,%i)\n",nodes[0]->GetId(),nodes[1]->GetId());
    345345                #endif
    … …  
    350350                thickness=h[0];
    351351
    352                 #ifdef _DEBUG_
     352                #ifdef _ISSM_DEBUG_
    353353                        printf("Thickness at grid (%i,%i): %lg\n",nodes[0]->GetId(),nodes[1]->GetID(),thickness);
    354354                #endif
    … …  
    446446                 * and we want to avoid zigzagging of the loads, we want lump the loads onto grids, not onto surfaces between grids.:*/
    447447       
    448                 #ifdef _DEBUG_
     448                #ifdef _ISSM_DEBUG_
    449449                _printf_("Grids  (%i,%i) are free of constraints\n",nodes[0]->GetId(),nodes[1]->GetID());
    450450                #endif

  • issm/trunk/src/c/objects/Tria.cpp

    r1901 r1904  
    284284
    285285}
    286 int    Tria::GetId(void){ return id; }
     286int    Tria::GetId(){ return id; }
    287287
    288288int    Tria::MyRank(void){
    … …  
    15591559        GetB(&B[0][0], xyz_list, gauss_l1l2l3);
    15601560
    1561         #ifdef _DEBUG_
    1562         printf("B for grid1 : [ %lf   %lf  \n",B[0][0],B[0][1]);
    1563         printf("              [ %lf   %lf  \n",B[1][0],B[1][1]);
    1564         printf("              [ %lf   %lf  ]\n",B[2][0],B[2][1]);
    1565 
    1566         printf("B for grid2 : [ %lf   %lf  \n",B[0][2],B[0][3]);
    1567         printf("              [ %lf   %lf  \n",B[1][2],B[1][3]);
    1568         printf("              [ %lf   %lf  ]\n",B[2][2],B[2][3]);
    1569 
    1570         printf("B for grid3 : [ %lf   %lf  \n",B[0][4],B[0][5]);
    1571         printf("              [ %lf   %lf  \n",B[1][4],B[1][5]);
    1572         printf("              [ %lf   %lf  ]\n",B[2][4],B[2][5]);
    1573                
    1574         for (i=0;i<numgrids;i++){
    1575                 printf("Velocity for grid %i %lf %lf\n",i,*(vxvy_list+2*i+0),*(vxvy_list+2*i+1));
    1576         }
    1577         #endif
    1578 
    15791561        /*Multiply B by velocity, to get strain rate: */
    15801562        MatrixMultiply( &B[0][0],3,NDOF2*numgrids,0,
    … …  
    16651647        GetNodalFunctionsDerivativesBasic(&dh1dh2dh3_basic[0][0],xyz_list, gauss_l1l2l3);
    16661648
    1667         #ifdef _DEBUG_
     1649        #ifdef _ISSM_ISSM_DEBUG_
    16681650        for (i=0;i<3;i++){
    16691651                printf("Node %i  dh/dx=%lf dh/dy=%lf \n",i,dh1dh2dh3_basic[0][i],dh1dh2dh3_basic[1][i]);
    … …  
    17951777        GetNodalFunctions(&l1l2l3[0],gauss_l1l2l3);
    17961778
    1797         #ifdef _DEBUG_
     1779        #ifdef _ISSM_ISSM_DEBUG_
    17981780        for (i=0;i<3;i++){
    17991781                printf("Node %i  h=%lf \n",i,l1l2l3[i]);
    … …  
    21462128                /* Get Jacobian determinant: */
    21472129                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2148                 #ifdef _DEBUG_
    2149                 printf("Element id %i Jacobian determinant: %lf\n",TriaElementGetID(this),Jdet);
     2130                #ifdef _ISSM_ISSM_DEBUG_
     2131                printf("Element id %i Jacobian determinant: %g\n",GetId(),Jdet);
    21502132                #endif
    21512133               
    … …  
    23832365                GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3);
    23842366                GetParameterValue(&drag, &k[0],gauss_l1l2l3);
    2385                 #ifdef _DEBUG_
     2367                #ifdef _ISSM_ISSM_DEBUG_
    23862368                        printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag);
    23872369                #endif
    … …  
    23902372                GetParameterValue(&lambda, &adjx_list[0],gauss_l1l2l3);
    23912373                GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3);
    2392                 #ifdef _DEBUG_
     2374                #ifdef _ISSM_ISSM_DEBUG_
    23932375                        printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu);
    23942376                #endif
    … …  
    23972379                GetParameterValue(&vx, &vx_list[0],gauss_l1l2l3);
    23982380                GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3);
    2399                 #ifdef _DEBUG_
     2381                #ifdef _ISSM_ISSM_DEBUG_
    24002382                        printf("Velocity vector %20.20lf %20.20lf\n",vx,vy);
    24012383                #endif
    … …  
    24032385                /* Get Jacobian determinant: */
    24042386                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2405                 #ifdef _DEBUG_
    2406                 printf("Element id %i Jacobian determinant: %lf\n",TriaElementGetID(this),Jdet);
     2387                #ifdef _ISSM_ISSM_DEBUG_
     2388                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    24072389                #endif
    24082390               
    … …  
    25812563                GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3);
    25822564                GetParameterValue(&drag, &k[0],gauss_l1l2l3);
    2583 #ifdef _DEBUG_
     2565#ifdef _ISSM_ISSM_DEBUG_
    25842566                printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag);
    25852567#endif
    … …  
    25892571                GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3);
    25902572                GetParameterValue(&xi, &adjz_list[0],gauss_l1l2l3);
    2591 #ifdef _DEBUG_
     2573#ifdef _ISSM_ISSM_DEBUG_
    25922574                printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu);
    25932575#endif
    … …  
    25972579                GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3);
    25982580                GetParameterValue(&vz, &vz_list[0],gauss_l1l2l3);
    2599 #ifdef _DEBUG_
     2581#ifdef _ISSM_ISSM_DEBUG_
    26002582                printf("Velocity vector %20.20lf %20.20lf\n",vx,vy);
    26012583
    … …  
    26102592                /* Get Jacobian determinant: */
    26112593                GetJacobianDeterminant3d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2612 #ifdef _DEBUG_
    2613                 printf("Element id %i Jacobian determinant: %lf\n",TriaElementGetID(this),Jdet);
     2594#ifdef _ISSM_ISSM_DEBUG_
     2595                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    26142596#endif
    26152597
    … …  
    27612743        /* Get gaussian points and weights (make this a statically initialized list of points? fstd): */
    27622744        GaussTria( &num_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &gauss_weights, 4);
    2763         #ifdef _DEBUG_
     2745        #ifdef _ISSM_ISSM_DEBUG_
    27642746        for (i=0;i<num_gauss;i++){
    27652747                printf("Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(gauss_weights+i));
    … …  
    27952777                /* Get nodal functions value at gaussian point:*/
    27962778                GetNodalFunctions(l1l2l3, gauss_l1l2l3);
    2797                 #ifdef _DEBUG_
    2798                         _printf_("viscositycomp %g thickness %g dvx [%g %g] dvy [%g %g]  dadjx [%g %g] dadjy[%g %g]\n",viscosity_complement,thickness,dvx[0],dvx[1],dvy[0],dvy[1],dadjx[0],dadjx[1],dadjy[0],dadjy[1]);
     2779                #ifdef _ISSM_ISSM_DEBUG_
     2780                        printf("viscositycomp %g thickness %g dvx [%g %g] dvy [%g %g]  dadjx [%g %g] dadjy[%g %g]\n",viscosity_complement,thickness,dvx[0],dvx[1],dvy[0],dvy[1],dadjx[0],dadjx[1],dadjy[0],dadjy[1]);
    27992781                #endif
    28002782
    … …  
    29312913        GaussTria( &num_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &gauss_weights, 2);
    29322914       
    2933         #ifdef _DEBUG_
     2915        #ifdef _ISSM_ISSM_DEBUG_
    29342916        for (i=0;i<num_gauss;i++){
    29352917                printf("Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(gauss_weights+i));
    … …  
    29472929                /* Get Jacobian determinant: */
    29482930                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2949                 #ifdef _DEBUG_
    2950                 printf("Element id %i Jacobian determinant: %lf\n",TriaElementGetID(this),Jdet);
     2931                #ifdef _ISSM_ISSM_DEBUG_
     2932                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    29512933                #endif
    29522934               

  • issm/trunk/src/c/shared/Exp/DomainOutlineRead.cpp

    r1 r1904  
    6060                nprof++;
    6161        }
    62         #ifdef _DEBUG_
     62        #ifdef _ISSM_DEBUG_
    6363                printf("Number of profiles in domain outline file: %i\n",nprof);
    6464        #endif
    … …  
    115115
    116116       
    117         #ifdef _DEBUG_
     117        #ifdef _ISSM_DEBUG_
    118118        for (i=0;i<nprof;i++){
    119119                printf("Profile #%i\n",i);

  • issm/trunk/src/c/shared/Sorting/binary_search.cpp

    r1 r1904  
    4040        else{
    4141                while((beg <= end) && (*mid != target)){
    42                         #ifdef _DEBUG_
     42                        #ifdef _ISSM_DEBUG_
    4343                                printf("1: %i %i %i\n",*beg,*mid,*(end-1));
    4444                        #endif
    … …  
    5252                                mid = beg + (end-beg)/2;  //new middle
    5353                        }
    54                         #ifdef _DEBUG_
     54                        #ifdef _ISSM_DEBUG_
    5555                                printf("2: %i %i %i\n",*beg,*mid,*(end-1));
    5656                        #endif

  • issm/trunk/src/c/shared/TriMesh/TriMeshUtils.cpp

    r1 r1904  
    159159                if (el2!=-1){
    160160                        /*el and el2 are on a segment rift, facing one another, plug them into riftsegments_uncompressed: */
    161                         #ifdef _DEBUG_
     161                        #ifdef _ISSM_DEBUG_
    162162                                printf("Elements %i and %i are on a rift\n",el+1,el2+1);
    163163                        #endif
    … …  
    184184        }
    185185
    186         #ifdef _DEBUG_
     186        #ifdef _ISSM_DEBUG_
    187187        for (i=0;i<nriftsegs;i++){
    188188                printf("Elements %i and %i are on a rift, joined by grids %i and %i\n",*(riftsegments+4*i+0)+1,*(riftsegments+4*i+1)+1,*(riftsegments+4*i+2),*(riftsegments+4*i+3));
    … …  
    215215        /*Build a list of all the elements connected to this grid: */
    216216        GridElementsList(&GridElements,&NumGridElements,grid,index,nel);
    217         #ifdef _DEBUG_
     217        #ifdef _ISSM_DEBUG_
    218218                printf("Connected elements for grid %i\n",grid);
    219219                for (k=0;k<NumGridElements;k++){
    … …  
    230230                                                                                                                           done rotating*/
    231231        GridElementListOnOneSideOfRift[1]=*(riftsegments+4*segmentnumber+1);
    232         #ifdef _DEBUG_
     232        #ifdef _ISSM_DEBUG_
    233233        printf("Starting with elements %i and %i for grid %i\n",GridElementListOnOneSideOfRift[0]+1,GridElementListOnOneSideOfRift[1]+1,grid);
    234234        #endif
    … …  
    238238                 * equal to GridElementListOnOneSideOfRift[counter-1]*/
    239239                for (k=0;k<NumGridElements;k++){
    240                         #ifdef _DEBUG_
     240                        #ifdef _ISSM_DEBUG_
    241241                        printf("k: %i GridElements[k]: %i GridElementListOnOneSideOfRift[counter-1]: %i GridElementListOnOneSideOfRift[counter] %i Neighboor %i\n",k,GridElements[k]+1,
    242242                                GridElementListOnOneSideOfRift[counter-1]+1,GridElementListOnOneSideOfRift[counter]+1,IsNeighbor(GridElements[k],GridElementListOnOneSideOfRift[counter],index));
    … …  
    266266                        GridElementListOnOneSideOfRift[l]=GridElementListOnOneSideOfRift[l+1];
    267267                }
    268                 #ifdef _DEBUG_
     268                #ifdef _ISSM_DEBUG_
    269269                printf("Grid %i is owned by  the following elements  on the same side of the rift: \n",grid);
    270270                for (l=0;l<NumGridElementListOnOneSideOfRift;l++){
    … …  
    313313                for (j=0;j<nsegs;j++){
    314314                        if (*(segments+3*j+2)==(el1+1)){
    315                                 #ifdef _DEBUG_
     315                                #ifdef _ISSM_DEBUG_
    316316                                        printf("Segment %i is the same as rift segment %i\n",j,i);
    317317                                #endif
    … …  
    348348                        }
    349349                        if (*(segments+3*j+2)==(el2+1)){
    350                                 #ifdef _DEBUG_
     350                                #ifdef _ISSM_DEBUG_
    351351                                        printf("Segment %i is the same as rift segment %i\n",j,i);
    352352                                #endif
    … …  
    441441        int n;
    442442        int el=-1;
    443         #ifdef _DEBUG_
     443        #ifdef _ISSM_DEBUG_
    444444                printf("Looking for %lf %lf\n",A,B);
    445445        #endif
    … …  
    447447                if (((*(index+3*n+0)==A)  || (*(index+3*n+1)==A) || (*(index+3*n+2)==A) ) && ((*(index+3*n+0)==B)  || (*(index+3*n+1)==B) || (*(index+3*n+2)==B))){
    448448                        el=n;
    449                         #ifdef _DEBUG_
     449                        #ifdef _ISSM_DEBUG_
    450450                                printf("Found them: %lf %lf %lf\n",*(index+3*n+0),*(index+3*n+1),*(index+3*n+2));
    451451                        #endif
    … …  
    938938                }
    939939
    940                 #ifdef _DEBUG_
     940                #ifdef _ISSM_DEBUG_
    941941                        printf("Tips for rift#%i   (%i-%i)\n",i,tip1,tip2);
    942942                #endif
    … …  
    977977                }
    978978
    979                 #ifdef _DEBUG_
     979                #ifdef _ISSM_DEBUG_
    980980                for (j=0;j<numsegs;j++){
    981981                        printf("%i\n",order[j]);

  • issm/trunk/src/c/toolkits/petsc/patches/MatPartition.cpp

    r1 r1904  
    6666                upper_row--;
    6767                range=upper_row-lower_row+1;
    68                 #ifdef _DEBUG_
     68                #ifdef _ISSM_DEBUG_
    6969                        PetscSynchronizedPrintf(MPI_COMM_WORLD,"My rank %i Range %i",my_rank,range);
    7070                        PetscSynchronizedFlush(MPI_COMM_WORLD);
    … …  
    8888                /*Now each node has a node_rows vectors holding which rows they should extract from matrixA. Create an Index Set from node_rows.*/
    8989                ISCreateGeneral(MPI_COMM_WORLD,count,node_rows,&row_index);
    90                 #ifdef _DEBUG_
     90                #ifdef _ISSM_DEBUG_
    9191                        ISView(row_index,PETSC_VIEWER_STDOUT_WORLD);
    9292                #endif
    … …  
    9898                }
    9999                ISCreateGeneral(MPI_COMM_WORLD,col_partition_vector_size,node_cols,&col_index);
    100                 #ifdef _DEBUG_
     100                #ifdef _ISSM_DEBUG_
    101101                        ISView(col_index,PETSC_VIEWER_STDOUT_WORLD);
    102102                #endif

  • issm/trunk/src/c/toolkits/petsc/patches/MatlabMatrixToDoubleMatrix.cpp

    r1 r1904  
    8686                }
    8787               
    88                 #ifdef _DEBUG_
     88                #ifdef _ISSM_DEBUG_
    8989                for(i=0;i<rows;i++){
    9090                        for(j=0;j<cols;j++){

  • issm/trunk/src/c/toolkits/petsc/patches/MatlabMatrixToPetscMatrix.cpp

    r1 r1904  
    8080                cols=mxGetN(mxmatrix);
    8181
    82                 #ifdef _DEBUG_
     82                #ifdef _ISSM_DEBUG_
    8383                for(i=0;i<rows;i++){
    8484                        for(j=0;j<cols;j++){

  • issm/trunk/src/c/toolkits/petsc/patches/PetscMatrixToMatlabMatrix.cpp

    r376 r1904  
    8585                        MatGetRow(matrix,i,&ncols,&columns,&column_values);
    8686
    87                         #ifdef _DEBUG_
     87                        #ifdef _ISSM_DEBUG_
    8888                        for(j=0;j<ncols;j++){
    8989                                printf("%i %i: %g\n",i,columns[j],column_values[j]);

  • issm/trunk/src/c/toolkits/petsc/patches/VecMerge.cpp

    r1 r1904  
    3434        VecGetSize(B,&MB);
    3535
    36         #ifdef _DEBUG_
     36        #ifdef _ISSM_DEBUG_
    3737                OutputServerMessages(0,0,"\n%s%s\n",STROFFSET,"Row partition vector:");
    3838                if (my_rank==0){
    … …  
    7373        VecAssemblyEnd(A);
    7474
    75         #ifdef _DEBUG_
     75        #ifdef _ISSM_DEBUG_
    7676                _printf_("Vector B:\n");
    7777                fflush(stdout);

  • issm/trunk/src/c/toolkits/petsc/patches/VecPartition.cpp

    r1 r1904  
    5858                upper_row--;
    5959                range=upper_row-lower_row+1;
    60                 #ifdef _DEBUG_
     60                #ifdef _ISSM_DEBUG_
    6161                        PetscSynchronizedPrintf(MPI_COMM_WORLD,"My rank %i Range %i\n",my_rank,range);
    6262                        PetscSynchronizedFlush(MPI_COMM_WORLD);
    … …  
    8686                }
    8787
    88                 #ifdef _DEBUG_
     88                #ifdef _ISSM_DEBUG_
    8989                        PetscSynchronizedPrintf(MPI_COMM_WORLD,"My rank: %i My count: %i node_rows: %p values: %p\n",my_rank,count,node_rows,values);
    9090                        PetscSynchronizedFlush(MPI_COMM_WORLD);
    … …  
    123123        }
    124124       
    125         #ifdef _DEBUG_
     125        #ifdef _ISSM_DEBUG_
    126126                _printf_("Vector A:\n");
    127127                fflush(stdout);

  • issm/trunk/src/c/toolkits/petsc/patches/VecToMPISerial.cpp

    r1 r1904  
    6464                if (my_rank==0){
    6565                        MPI_Recv(buffer,3,MPI_INT,i,1,MPI_COMM_WORLD,&status);
    66                         #ifdef _DEBUG_
     66                        #ifdef _ISSM_DEBUG_
    6767                                _printf_("Received from node %i: %i-%i\n",buffer[0],buffer[1],buffer[2]);
    6868                        #endif

  • issm/trunk/src/c/toolkits/plapack/patches/CyclicalFactorization.cpp

    r1439 r1904  
    3838        *pnprows=nprows;
    3939        *pnpcols=npcols;
    40         #ifdef _DEBUG_
     40        #ifdef _ISSM_DEBUG_
    4141                _printf_("Decomposition: %i-%i\n",nprows,npcols);
    4242        #endif
    … …  
    5656        for (i=0;i<input;i++){
    5757                SmallestPrimeFactor(&prime_factor,*(decomp+i));
    58                 #ifdef _DEBUG_
     58                #ifdef _ISSM_DEBUG_
    5959                        _printf_("Smallest prime factor for term %i : %i\n",i,prime_factor);
    6060                #endif
    … …  
    6868                }
    6969        }
    70         #ifdef _DEBUG_
     70        #ifdef _ISSM_DEBUG_
    7171                _printf_("Prime factor decomposition for integer %i: \n",input);
    7272                for(i=0;i<*pdecomp_size;i++){

  • issm/trunk/src/c/toolkits/plapack/patches/PlapackInvertMatrix.cpp

    r831 r1904  
    119119        PLA_API_end();
    120120
    121         #ifdef _DEBUG_
     121        #ifdef _ISSM_DEBUG_
    122122                PLA_Global_show("Matrix A",a," %lf","Done with A");
    123123                MatView(*A,PETSC_VIEWER_STDOUT_WORLD);
    … …  
    131131        PlapackToPetsc(inv_A,local_mA,local_nA,mA,nA,type,a,templ,nprows,npcols,nb);
    132132
    133         #ifdef _DEBUG_
     133        #ifdef _ISSM_DEBUG_
    134134                PLA_Global_show("Inverse of A",a," %lf","Done...");
    135135                MatView(*inv_A,PETSC_VIEWER_STDOUT_WORLD);

  • issm/trunk/src/m/classes/@model/model.m

    r1899 r1904  
    202202        md.element_debug=0;
    203203        md.element_debugid=NaN;
     204        md.mem_debug=0;
    204205
    205206        %Results fields

  • issm/trunk/src/m/classes/public/BuildQueueingScriptGeneric.m

    r1790 r1904  
    44%   Usage:
    55%      BuildQueueingScriptGeneric(md,executionpath,codepath)
     6
     7global ISSM_DIR
    68
    79%Open queuing script file
    … …  
    1517fprintf(fid,'#!/bin/sh\n');
    1618fprintf(fid,'rm -rf %s/%s.lock\n',executionpath,md.name);
     19
     20if md.mem_debug==0,
    1721fprintf(fid,'mpirun -np %i %s/%s.exe %s %s.bin %s.outbin %s.lock  2> %s.errlog >%s.outlog & ',md.np,codepath,AnalysisTypeFromEnum(md.analysis_type),executionpath,md.name,md.name,md.name,md.name,md.name);
     22else
     23fprintf(fid,'LD_PRELOAD=%s mpirun -np %i %s --leak-check=full %s/%s.exe %s %s.bin %s.outbin %s.lock  2> %s.errlog >%s.outlog & ',[ISSM_DIR '/externalpackages/valgrind/install/lib/libmpidebug.so'],md.np,[ISSM_DIR '/externalpackages/valgrind/install/bin/valgrind'],codepath,AnalysisTypeFromEnum(md.analysis_type),executionpath,md.name,md.name,md.name,md.name,md.name);
     24end
     25
    1826fclose(fid);

  • issm/trunk/src/mex/ContourToMesh/ContourToMesh.cpp

    r292 r1904  
    8282        }
    8383
    84         #ifdef _DEBUG_
     84        #ifdef _ISSM_DEBUG_
    8585        for(i=0;i<numcontours;i++){
    8686                printf("\nContour echo: contour number  %i / %i\n",i+1,numcontours);

  • issm/trunk/src/mex/ContourToNodes/ContourToNodes.cpp

    r1103 r1904  
    7575        }
    7676
    77         #ifdef _DEBUG_
     77        #ifdef _ISSM_DEBUG_
    7878        for(i=0;i<numcontours;i++){
    7979                printf("\nContour echo: contour number  %i / %i\n",i+1,numcontours);

  • issm/trunk/src/mex/TriMesh/TriMesh.cpp

    r1130 r1904  
    9898        }
    9999       
    100         #ifdef _DEBUG_
     100        #ifdef _ISSM_DEBUG_
    101101        printf("Domain name: %s\n",domainname);
    102102        printf("Rift name: %s\n",riftname);
    … …  
    122122        }
    123123
    124         #ifdef _DEBUG_
     124        #ifdef _ISSM_DEBUG_
    125125        printf("# points in domain outline: %i\n",nprof);
    126126        #endif
    … …  
    137137
    138138
    139         #ifdef _DEBUG_
     139        #ifdef _ISSM_DEBUG_
    140140        printf("# points in rift outline: %i\n",numrifts);
    141141        #endif

  • issm/trunk/src/mex/TriMeshProcessRifts/TriMeshProcessRifts.cpp

    r1805 r1904  
    142142        }
    143143
    144         #ifdef _DEBUG_
     144        #ifdef _ISSM_DEBUG_
    145145        printf("Index: \n");
    146146        for (i=0;i<nel;i++){

  • issm/trunk/src/mex/TriMeshRefine/TriMeshRefine.cpp

    r1 r1904  
    160160
    161161
    162         #ifdef _DEBUG_
     162        #ifdef _ISSM_DEBUG_
    163163        /*printf("Index: \n");
    164164        for (i=0;i<nel;i++){
Note: See TracChangeset for help on using the changeset viewer.