Changeset 14384


Ignore:
Timestamp:
02/23/13 16:55:59 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: getting rid of GetJacobian3d, which did not make sense (this is done by PentaRef), renamed GetJacobian2d accordingly

Location:
issm/trunk-jpl/src/c/classes/objects/Elements
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp

    r14378 r14384  
    301301
    302302                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    303                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0], gauss);
     303                GetJacobianDeterminant(&Jdet, &xyz_list[0][0], gauss);
    304304
    305305                D_scalar=latentheat/heatcapacity*gauss->weight*Jdet;
    … …  
    378378                gauss->GaussPoint(ig);
    379379
    380                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     380                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    381381                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    382382
    … …  
    494494                vyaverage_input->GetInputValue(&vy,gauss);
    495495
    496                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     496                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    497497                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    498498
    … …  
    547547                gauss->GaussPoint(ig);
    548548
    549                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     549                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    550550                D=gauss->weight*Jdet;
    551551
    … …  
    668668                gauss->GaussPoint(ig);
    669669
    670                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     670                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    671671                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    672672
    … …  
    716716                gauss->GaussPoint(ig);
    717717
    718                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     718                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    719719                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    720720
    … …  
    766766                gauss->GaussPoint(ig);
    767767
    768                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     768                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    769769                GetNodalFunctions(basis, gauss);
    770770
    … …  
    12071207
    12081208                /*Compute fraction of grounded element*/
    1209                 GetJacobianDeterminant2d(&area_init, xyz_list,gauss);
    1210                 GetJacobianDeterminant2d(&area_grounded, &xyz_bis[0][0],gauss);
     1209                GetJacobianDeterminant(&area_init, xyz_list,gauss);
     1210                GetJacobianDeterminant(&area_grounded, &xyz_bis[0][0],gauss);
    12111211                if(mainlyfloating==true) area_grounded=area_init-area_grounded;
    12121212                phi=area_grounded/area_init;
    … …  
    30803080                gauss->GaussPoint(ig);
    30813081
    3082                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3082                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    30833083                GetBMacAyeal(&B[0][0], &xyz_list[0][0], gauss);
    30843084                GetBprimeMacAyeal(&Bprime[0][0], &xyz_list[0][0], gauss);
    … …  
    31643164
    31653165                GetL(&L[0][0], &xyz_list[0][0], gauss,NDOF2);
    3166                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3166                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    31673167                DL_scalar=alpha2*gauss->weight*Jdet;
    31683168                for (i=0;i<2;i++) DL[i][i]=DL_scalar;
    … …  
    32353235                gauss->GaussPoint(ig);
    32363236
    3237                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3237                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    32383238                GetNodalFunctions(basis, gauss);
    32393239
    … …  
    33423342                gauss->GaussPoint(ig);
    33433343
    3344                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3344                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    33453345                GetNodalFunctionsDerivatives(&dphi[0][0],&xyz_list[0][0],gauss);
    33463346
    … …  
    36113611
    36123612                /* Get Jacobian determinant: */
    3613                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3613                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    36143614
    36153615                /*Get all parameters at gaussian point*/
    … …  
    38213821                gauss->GaussPoint(ig);
    38223822
    3823                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3823                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    38243824                GetNodalFunctionsDerivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
    38253825                weights_input->GetInputValue(&weight,gauss,weight_index);
    … …  
    38613861                gauss->GaussPoint(ig);
    38623862
    3863                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3863                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    38643864                GetNodalFunctionsDerivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
    38653865                weights_input->GetInputValue(&weight,gauss,weight_index);
    … …  
    39193919                material->GetViscosityComplement(&viscosity_complement,&epsilon[0]);
    39203920
    3921                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3921                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    39223922                GetNodalFunctions(basis,gauss);
    39233923
    … …  
    39763976                material->GetViscosityZComplement(&viscosity_complement,&epsilon[0]);
    39773977
    3978                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     3978                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    39793979                GetNodalFunctions(basis,gauss);
    39803980
    … …  
    40334033                gauss->GaussPoint(ig);
    40344034
    4035                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4035                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    40364036                GetNodalFunctions(basis, gauss);
    40374037
    … …  
    41044104                gauss->GaussPoint(ig);
    41054105
    4106                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4106                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    41074107                GetNodalFunctionsDerivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
    41084108                weights_input->GetInputValue(&weight,gauss,weight_index);
    … …  
    41664166                gauss->GaussPoint(ig);
    41674167
    4168                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4168                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    41694169                GetNodalFunctions(basis, gauss);
    41704170
    … …  
    42094209                gauss->GaussPoint(ig);
    42104210
    4211                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4211                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    42124212                GetNodalFunctions(basis, gauss);
    42134213
    … …  
    42674267                gauss->GaussPoint(ig);
    42684268
    4269                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4269                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    42704270                GetNodalFunctions(basis, gauss);
    42714271                GetNodalFunctionsDerivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
    … …  
    43764376
    43774377                /* Get Jacobian determinant: */
    4378                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4378                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    43794379
    43804380                /*Get all parameters at gaussian point*/
    … …  
    44214421
    44224422                /* Get Jacobian determinant: */
    4423                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4423                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    44244424
    44254425                /*Get all parameters at gaussian point*/
    … …  
    44814481
    44824482                /* Get Jacobian determinant: */
    4483                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4483                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    44844484
    44854485                /*Get all parameters at gaussian point*/
    … …  
    45434543
    45444544                /* Get Jacobian determinant: */
    4545                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4545                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    45464546
    45474547                /*Get all parameters at gaussian point*/
    … …  
    46034603
    46044604                /* Get Jacobian determinant: */
    4605                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4605                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    46064606
    46074607                /*Get all parameters at gaussian point*/
    … …  
    46644664
    46654665                /* Get Jacobian determinant: */
    4666                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4666                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    46674667
    46684668                /*Get all parameters at gaussian point*/
    … …  
    47224722
    47234723                /* Get Jacobian determinant: */
    4724                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4724                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    47254725
    47264726                /*Get all parameters at gaussian point*/
    … …  
    47684768
    47694769                /* Get Jacobian determinant: */
    4770                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4770                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    47714771
    47724772                /*Get all parameters at gaussian point*/
    … …  
    48194819
    48204820                /* Get Jacobian determinant: */
    4821                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4821                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    48224822
    48234823                /*Get all parameters at gaussian point*/
    … …  
    48664866
    48674867                /* Get Jacobian determinant: */
    4868                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4868                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    48694869
    48704870                /*Get parameters at gauss point*/
    … …  
    49214921                gauss->GaussPoint(ig);
    49224922
    4923                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     4923                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    49244924                GetNodalFunctions(basis, gauss);
    49254925                GetNodalFunctionsDerivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
    … …  
    50185018
    50195019                /* Get Jacobian determinant: */
    5020                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     5020                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    50215021
    50225022                /*Get all parameters at gaussian point*/
    … …  
    52005200
    52015201                /* Get Jacobian determinant: */
    5202                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     5202                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    52035203
    52045204                /*Get all parameters at gaussian point*/
    … …  
    53665366
    53675367                /* Get Jacobian determinant: */
    5368                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     5368                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    53695369
    53705370                /*Get all parameters at gaussian point*/
    … …  
    54395439                gauss->GaussPoint(ig);
    54405440
    5441                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss);
     5441                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    54425442                GetNodalFunctionsDerivatives(&dphi[0][0],&xyz_list[0][0],gauss);
    54435443
    … …  
    57005700                gauss->GaussPoint(ig);
    57015701
    5702                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     5702                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    57035703                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    57045704
    … …  
    57935793                gauss->GaussPoint(ig);
    57945794
    5795                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     5795                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    57965796                GetNodalFunctions(basis, gauss);
    57975797
    … …  
    60796079                gauss->GaussPoint(ig);
    60806080
    6081                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     6081                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    60826082                GetBPrognostic(&B[0][0], &xyz_list[0][0], gauss);
    60836083                GetBprimePrognostic(&Bprime[0][0], &xyz_list[0][0], gauss);
    … …  
    61736173                gauss->GaussPoint(ig);
    61746174
    6175                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     6175                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    61766176                /*WARNING: B and Bprime are inverted compared to usual prognostic!!!!*/
    61776177                GetBPrognostic(&Bprime[0][0], &xyz_list[0][0], gauss);
    … …  
    62426242                dhdt_input->GetInputValue(&dhdt_g,gauss);
    62436243
    6244                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     6244                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    62456245                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    62466246
    … …  
    62856285                dhdt_input->GetInputValue(&dhdt_g,gauss);
    62866286
    6287                 GetJacobianDeterminant2d(&Jdettria, &xyz_list[0][0],gauss);
     6287                GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
    62886288                GetL(&L[0], &xyz_list[0][0], gauss,NDOF1);
    62896289

  • issm/trunk-jpl/src/c/classes/objects/Elements/TriaRef.cpp

    r13622 r14384  
    341341}
    342342/*}}}*/
    343 /*FUNCTION TriaRef::GetJacobianDeterminant2d{{{*/
    344 void TriaRef::GetJacobianDeterminant2d(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss){
     343/*FUNCTION TriaRef::GetJacobianDeterminant{{{*/
     344void TriaRef::GetJacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss){
    345345        /*The Jacobian determinant is constant over the element, discard the gaussian points.
    346346         * J is assumed to have been allocated of size NDOF2xNDOF2.*/
    … …  
    352352        /*Get Determinant*/
    353353        Matrix2x2Determinant(Jdet,&J[0][0]);
    354         if(*Jdet<0) _error_("negative jacobian determinant!");
    355 
    356 }
    357 /*}}}*/
    358 /*FUNCTION TriaRef::GetJacobianDeterminant3d {{{*/
    359 void TriaRef::GetJacobianDeterminant3d(IssmDouble*  Jdet, IssmDouble* xyz_list,GaussTria* gauss){
    360         /*The Jacobian determinant is constant over the element, discard the gaussian points.
    361          * J is assumed to have been allocated of size NDOF2xNDOF2.*/
    362 
    363         IssmDouble x1,x2,x3,y1,y2,y3,z1,z2,z3;
    364 
    365         x1=*(xyz_list+3*0+0);
    366         y1=*(xyz_list+3*0+1);
    367         z1=*(xyz_list+3*0+2);
    368         x2=*(xyz_list+3*1+0);
    369         y2=*(xyz_list+3*1+1);
    370         z2=*(xyz_list+3*1+2);
    371         x3=*(xyz_list+3*2+0);
    372         y3=*(xyz_list+3*2+1);
    373         z3=*(xyz_list+3*2+2);
    374 
    375         *Jdet=SQRT3/6.0*pow(pow(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1)),2.0)+pow(((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1)),2.0)+pow(((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1)),2.0),0.5);
    376354        if(*Jdet<0) _error_("negative jacobian determinant!");
    377355

  • issm/trunk-jpl/src/c/classes/objects/Elements/TriaRef.h

    r13623 r14384  
    3232                void GetJacobian(IssmDouble* J, IssmDouble* xyz_list,GaussTria* gauss);
    3333                void GetSegmentJacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss);
    34                 void GetJacobianDeterminant2d(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss);
    35                 void GetJacobianDeterminant3d(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss);
     34                void GetJacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,GaussTria* gauss);
    3635                void GetJacobianInvert(IssmDouble*  Jinv, IssmDouble* xyz_list,GaussTria* gauss);
    3736                void GetNodalFunctions(IssmDouble* l1l2l3,GaussTria* gauss);
Note: See TracChangeset for help on using the changeset viewer.