Changeset 2046
- Timestamp:
- 08/28/09 11:00:21 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Tria.cpp
r1904 r2046 1647 1647 GetNodalFunctionsDerivativesBasic(&dh1dh2dh3_basic[0][0],xyz_list, gauss_l1l2l3); 1648 1648 1649 #ifdef _ISSM_ ISSM_DEBUG_1649 #ifdef _ISSM_DEBUG_ 1650 1650 for (i=0;i<3;i++){ 1651 1651 printf("Node %i dh/dx=%lf dh/dy=%lf \n",i,dh1dh2dh3_basic[0][i],dh1dh2dh3_basic[1][i]); … … 1777 1777 GetNodalFunctions(&l1l2l3[0],gauss_l1l2l3); 1778 1778 1779 #ifdef _ISSM_ ISSM_DEBUG_1779 #ifdef _ISSM_DEBUG_ 1780 1780 for (i=0;i<3;i++){ 1781 1781 printf("Node %i h=%lf \n",i,l1l2l3[i]); … … 2128 2128 /* Get Jacobian determinant: */ 2129 2129 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3); 2130 #ifdef _ISSM_ ISSM_DEBUG_2130 #ifdef _ISSM_DEBUG_ 2131 2131 printf("Element id %i Jacobian determinant: %g\n",GetId(),Jdet); 2132 2132 #endif … … 2365 2365 GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3); 2366 2366 GetParameterValue(&drag, &k[0],gauss_l1l2l3); 2367 #ifdef _ISSM_ ISSM_DEBUG_2367 #ifdef _ISSM_DEBUG_ 2368 2368 printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag); 2369 2369 #endif … … 2372 2372 GetParameterValue(&lambda, &adjx_list[0],gauss_l1l2l3); 2373 2373 GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3); 2374 #ifdef _ISSM_ ISSM_DEBUG_2374 #ifdef _ISSM_DEBUG_ 2375 2375 printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu); 2376 2376 #endif … … 2379 2379 GetParameterValue(&vx, &vx_list[0],gauss_l1l2l3); 2380 2380 GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3); 2381 #ifdef _ISSM_ ISSM_DEBUG_2381 #ifdef _ISSM_DEBUG_ 2382 2382 printf("Velocity vector %20.20lf %20.20lf\n",vx,vy); 2383 2383 #endif … … 2385 2385 /* Get Jacobian determinant: */ 2386 2386 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3); 2387 #ifdef _ISSM_ ISSM_DEBUG_2387 #ifdef _ISSM_DEBUG_ 2388 2388 printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet); 2389 2389 #endif … … 2563 2563 GetParameterValue(&alpha_complement, &alpha_complement_list[0],gauss_l1l2l3); 2564 2564 GetParameterValue(&drag, &k[0],gauss_l1l2l3); 2565 #ifdef _ISSM_ ISSM_DEBUG_2565 #ifdef _ISSM_DEBUG_ 2566 2566 printf("Drag complement: %20.20lf Drag: %20.20lf\n",alpha_complement,drag); 2567 2567 #endif … … 2571 2571 GetParameterValue(&mu, &adjy_list[0],gauss_l1l2l3); 2572 2572 GetParameterValue(&xi, &adjz_list[0],gauss_l1l2l3); 2573 #ifdef _ISSM_ ISSM_DEBUG_2573 #ifdef _ISSM_DEBUG_ 2574 2574 printf("Adjoint vector %20.20lf %20.20lf\n",lambda,mu); 2575 2575 #endif … … 2579 2579 GetParameterValue(&vy, &vy_list[0],gauss_l1l2l3); 2580 2580 GetParameterValue(&vz, &vz_list[0],gauss_l1l2l3); 2581 #ifdef _ISSM_ ISSM_DEBUG_2581 #ifdef _ISSM_DEBUG_ 2582 2582 printf("Velocity vector %20.20lf %20.20lf\n",vx,vy); 2583 2583 … … 2592 2592 /* Get Jacobian determinant: */ 2593 2593 GetJacobianDeterminant3d(&Jdet, &xyz_list[0][0],gauss_l1l2l3); 2594 #ifdef _ISSM_ ISSM_DEBUG_2594 #ifdef _ISSM_DEBUG_ 2595 2595 printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet); 2596 2596 #endif … … 2743 2743 /* Get gaussian points and weights (make this a statically initialized list of points? fstd): */ 2744 2744 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_ 2746 2746 for (i=0;i<num_gauss;i++){ 2747 2747 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)); … … 2777 2777 /* Get nodal functions value at gaussian point:*/ 2778 2778 GetNodalFunctions(l1l2l3, gauss_l1l2l3); 2779 #ifdef _ISSM_ ISSM_DEBUG_2779 #ifdef _ISSM_DEBUG_ 2780 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]); 2781 2781 #endif … … 2913 2913 GaussTria( &num_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &gauss_weights, 2); 2914 2914 2915 #ifdef _ISSM_ ISSM_DEBUG_2915 #ifdef _ISSM_DEBUG_ 2916 2916 for (i=0;i<num_gauss;i++){ 2917 2917 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)); … … 2929 2929 /* Get Jacobian determinant: */ 2930 2930 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3); 2931 #ifdef _ISSM_ ISSM_DEBUG_2931 #ifdef _ISSM_DEBUG_ 2932 2932 printf("Element id %i Jacobian determinant: %lf\n",GetId(),Jdet); 2933 2933 #endif
Note:
See TracChangeset
for help on using the changeset viewer.