Changeset 2046


Ignore:
Timestamp:
08/28/09 11:00:21 (15 years ago)
Author:
Eric.Larour
Message:

wrong debug flag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Tria.cpp

    r1904 r2046  
    16471647        GetNodalFunctionsDerivativesBasic(&dh1dh2dh3_basic[0][0],xyz_list, gauss_l1l2l3);
    16481648
    1649         #ifdef _ISSM_ISSM_DEBUG_
     1649        #ifdef _ISSM_DEBUG_
    16501650        for (i=0;i<3;i++){
    16511651                printf("Node %i  dh/dx=%lf dh/dy=%lf \n",i,dh1dh2dh3_basic[0][i],dh1dh2dh3_basic[1][i]);
     
    17771777        GetNodalFunctions(&l1l2l3[0],gauss_l1l2l3);
    17781778
    1779         #ifdef _ISSM_ISSM_DEBUG_
     1779        #ifdef _ISSM_DEBUG_
    17801780        for (i=0;i<3;i++){
    17811781                printf("Node %i  h=%lf \n",i,l1l2l3[i]);
     
    21282128                /* Get Jacobian determinant: */
    21292129                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2130                 #ifdef _ISSM_ISSM_DEBUG_
     2130                #ifdef _ISSM_DEBUG_
    21312131                printf("Element id %i Jacobian determinant: %g\n",GetId(),Jdet);
    21322132                #endif
     
    23652365                GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3);
    23662366                GetParameterValue(&drag, &k[0],gauss_l1l2l3);
    2367                 #ifdef _ISSM_ISSM_DEBUG_
     2367                #ifdef _ISSM_DEBUG_
    23682368                        printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag);
    23692369                #endif
     
    23722372                GetParameterValue(&lambda, &adjx_list[0],gauss_l1l2l3);
    23732373                GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3);
    2374                 #ifdef _ISSM_ISSM_DEBUG_
     2374                #ifdef _ISSM_DEBUG_
    23752375                        printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu);
    23762376                #endif
     
    23792379                GetParameterValue(&vx, &vx_list[0],gauss_l1l2l3);
    23802380                GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3);
    2381                 #ifdef _ISSM_ISSM_DEBUG_
     2381                #ifdef _ISSM_DEBUG_
    23822382                        printf("Velocity vector %20.20lf %20.20lf\n",vx,vy);
    23832383                #endif
     
    23852385                /* Get Jacobian determinant: */
    23862386                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2387                 #ifdef _ISSM_ISSM_DEBUG_
     2387                #ifdef _ISSM_DEBUG_
    23882388                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    23892389                #endif
     
    25632563                GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3);
    25642564                GetParameterValue(&drag, &k[0],gauss_l1l2l3);
    2565 #ifdef _ISSM_ISSM_DEBUG_
     2565#ifdef _ISSM_DEBUG_
    25662566                printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag);
    25672567#endif
     
    25712571                GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3);
    25722572                GetParameterValue(&xi, &adjz_list[0],gauss_l1l2l3);
    2573 #ifdef _ISSM_ISSM_DEBUG_
     2573#ifdef _ISSM_DEBUG_
    25742574                printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu);
    25752575#endif
     
    25792579                GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3);
    25802580                GetParameterValue(&vz, &vz_list[0],gauss_l1l2l3);
    2581 #ifdef _ISSM_ISSM_DEBUG_
     2581#ifdef _ISSM_DEBUG_
    25822582                printf("Velocity vector %20.20lf %20.20lf\n",vx,vy);
    25832583
     
    25922592                /* Get Jacobian determinant: */
    25932593                GetJacobianDeterminant3d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2594 #ifdef _ISSM_ISSM_DEBUG_
     2594#ifdef _ISSM_DEBUG_
    25952595                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    25962596#endif
     
    27432743        /* Get gaussian points and weights (make this a statically initialized list of points? fstd): */
    27442744        GaussTria( &num_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &gauss_weights, 4);
    2745         #ifdef _ISSM_ISSM_DEBUG_
     2745        #ifdef _ISSM_DEBUG_
    27462746        for (i=0;i<num_gauss;i++){
    27472747                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));
     
    27772777                /* Get nodal functions value at gaussian point:*/
    27782778                GetNodalFunctions(l1l2l3, gauss_l1l2l3);
    2779                 #ifdef _ISSM_ISSM_DEBUG_
     2779                #ifdef _ISSM_DEBUG_
    27802780                        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]);
    27812781                #endif
     
    29132913        GaussTria( &num_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &gauss_weights, 2);
    29142914       
    2915         #ifdef _ISSM_ISSM_DEBUG_
     2915        #ifdef _ISSM_DEBUG_
    29162916        for (i=0;i<num_gauss;i++){
    29172917                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));
     
    29292929                /* Get Jacobian determinant: */
    29302930                GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
    2931                 #ifdef _ISSM_ISSM_DEBUG_
     2931                #ifdef _ISSM_DEBUG_
    29322932                printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet);
    29332933                #endif
Note: See TracChangeset for help on using the changeset viewer.