Changeset 18075


Ignore:
Timestamp:
05/31/14 04:19:56 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: minor cosmetics

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

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

    r18068 r18075  
    58145814
    58155815        int num_dof=2;
    5816         IssmDouble L1L2l3[3];
     5816        IssmDouble basis[3];
    58175817
    58185818        /*Cast gauss to GaussPenta*/
    … …  
    58205820        GaussPenta* gauss = dynamic_cast<GaussPenta*>(gauss_in);
    58215821
    5822         /*Get L1L2l3 in actual coordinate system: */
    5823         L1L2l3[0]=gauss->coord1*(1-gauss->coord4)/2.0;
    5824         L1L2l3[1]=gauss->coord2*(1-gauss->coord4)/2.0;
    5825         L1L2l3[2]=gauss->coord3*(1-gauss->coord4)/2.0;
     5822        /*Get basis in actual coordinate system: */
     5823        basis[0]=gauss->coord1*(1-gauss->coord4)/2.0;
     5824        basis[1]=gauss->coord2*(1-gauss->coord4)/2.0;
     5825        basis[2]=gauss->coord3*(1-gauss->coord4)/2.0;
    58265826
    58275827        /*Build LFS: */
    58285828        for(int i=0;i<3;i++){
    5829                 LFS[num_dof*3*0+num_dof*i+0] = L1L2l3[i];
     5829                LFS[num_dof*3*0+num_dof*i+0] = basis[i];
    58305830                LFS[num_dof*3*0+num_dof*i+1] = 0;
    58315831                LFS[num_dof*3*1+num_dof*i+0] = 0;
    5832                 LFS[num_dof*3*1+num_dof*i+1] = L1L2l3[i];
    5833                 LFS[num_dof*3*2+num_dof*i+0] = L1L2l3[i];
     5832                LFS[num_dof*3*1+num_dof*i+1] = basis[i];
     5833                LFS[num_dof*3*2+num_dof*i+0] = basis[i];
    58345834                LFS[num_dof*3*2+num_dof*i+1] = 0;
    58355835                LFS[num_dof*3*3+num_dof*i+0] = 0;
    5836                 LFS[num_dof*3*3+num_dof*i+1] = L1L2l3[i];
    5837                 LFS[num_dof*3*4+num_dof*i+0] = L1L2l3[i];
     5836                LFS[num_dof*3*3+num_dof*i+1] = basis[i];
     5837                LFS[num_dof*3*4+num_dof*i+0] = basis[i];
    58385838                LFS[num_dof*3*4+num_dof*i+1] = 0;
    58395839                LFS[num_dof*3*5+num_dof*i+0] = 0;
    5840                 LFS[num_dof*3*5+num_dof*i+1] = L1L2l3[i];
    5841                 LFS[num_dof*3*6+num_dof*i+0] = L1L2l3[i];
     5840                LFS[num_dof*3*5+num_dof*i+1] = basis[i];
     5841                LFS[num_dof*3*6+num_dof*i+0] = basis[i];
    58425842                LFS[num_dof*3*6+num_dof*i+1] = 0;
    58435843                LFS[num_dof*3*7+num_dof*i+0] = 0;
    5844                 LFS[num_dof*3*7+num_dof*i+1] = L1L2l3[i];
     5844                LFS[num_dof*3*7+num_dof*i+1] = basis[i];
    58455845        }
    58465846}/*}}}*/
    … …  
    58625862        int num_dof_vel=3*7;
    58635863        int num_dof_total=3*7+1*6;
    5864         IssmDouble L1L2l3[3];
     5864        IssmDouble basis[3];
    58655865        IssmDouble dbasis[3][6];
    58665866
    … …  
    58695869        GaussPenta* gauss = dynamic_cast<GaussPenta*>(gauss_in);
    58705870
    5871         /*Get L1L2l3 in actual coordinate system: */
    5872         L1L2l3[0]=gauss->coord1*(1-gauss->coord4)/2.0;
    5873         L1L2l3[1]=gauss->coord2*(1-gauss->coord4)/2.0;
    5874         L1L2l3[2]=gauss->coord3*(1-gauss->coord4)/2.0;
     5871        /*Get basis in actual coordinate system: */
     5872        basis[0]=gauss->coord1*(1-gauss->coord4)/2.0;
     5873        basis[1]=gauss->coord2*(1-gauss->coord4)/2.0;
     5874        basis[2]=gauss->coord3*(1-gauss->coord4)/2.0;
    58755875
    58765876        element->NodalFunctionsP1Derivatives(&dbasis[0][0],xyz_list,gauss);
    … …  
    58785878        /*Build LprimeFS: */
    58795879        for(int i=0;i<3;i++){
    5880                 LprimeFS[num_dof_total*0+num_dof*i+0] = L1L2l3[i];
     5880                LprimeFS[num_dof_total*0+num_dof*i+0] = basis[i];
    58815881                LprimeFS[num_dof_total*0+num_dof*i+1] = 0.;
    58825882                LprimeFS[num_dof_total*0+num_dof*i+2] = 0.;
    58835883                LprimeFS[num_dof_total*1+num_dof*i+0] = 0.;
    5884                 LprimeFS[num_dof_total*1+num_dof*i+1] = L1L2l3[i];
     5884                LprimeFS[num_dof_total*1+num_dof*i+1] = basis[i];
    58855885                LprimeFS[num_dof_total*1+num_dof*i+2] = 0.;
    58865886                LprimeFS[num_dof_total*2+num_dof*i+0] = 0.;
    58875887                LprimeFS[num_dof_total*2+num_dof*i+1] = 0.;
    5888                 LprimeFS[num_dof_total*2+num_dof*i+2] = L1L2l3[i];
     5888                LprimeFS[num_dof_total*2+num_dof*i+2] = basis[i];
    58895889                LprimeFS[num_dof_total*3+num_dof*i+0] = 0.;
    58905890                LprimeFS[num_dof_total*3+num_dof*i+1] = 0.;
    5891                 LprimeFS[num_dof_total*3+num_dof*i+2] = L1L2l3[i];
     5891                LprimeFS[num_dof_total*3+num_dof*i+2] = basis[i];
    58925892                LprimeFS[num_dof_total*4+num_dof*i+0] = 0.;
    58935893                LprimeFS[num_dof_total*4+num_dof*i+1] = 0.;
    … …  
    59365936                LprimeFS[num_dof_total*4+num_dof_vel+i] = 0.;
    59375937                LprimeFS[num_dof_total*5+num_dof_vel+i] = 0.;
    5938                 LprimeFS[num_dof_total*6+num_dof_vel+i] = L1L2l3[i];
    5939                 LprimeFS[num_dof_total*7+num_dof_vel+i] = L1L2l3[i];
     5938                LprimeFS[num_dof_total*6+num_dof_vel+i] = basis[i];
     5939                LprimeFS[num_dof_total*7+num_dof_vel+i] = basis[i];
    59405940        }
    59415941        for(int i=3;i<6;i++){
    … …  
    59625962
    59635963        int num_dof=3;
    5964         IssmDouble L1L2l3[3];
     5964        IssmDouble basis[3];
    59655965
    59665966        /*Cast gauss to GaussPenta*/
    … …  
    59685968        GaussPenta* gauss = dynamic_cast<GaussPenta*>(gauss_in);
    59695969
    5970         /*Get L1L2l3 in actual coordinate system: */
    5971         L1L2l3[0]=gauss->coord1*(1-gauss->coord4)/2.0;
    5972         L1L2l3[1]=gauss->coord2*(1-gauss->coord4)/2.0;
    5973         L1L2l3[2]=gauss->coord3*(1-gauss->coord4)/2.0;
     5970        /*Get basis in actual coordinate system: */
     5971        basis[0]=gauss->coord1*(1-gauss->coord4)/2.0;
     5972        basis[1]=gauss->coord2*(1-gauss->coord4)/2.0;
     5973        basis[2]=gauss->coord3*(1-gauss->coord4)/2.0;
    59745974
    59755975        /*Build LFS: */
    59765976        for(int i=0;i<3;i++){
    5977                 LFS[num_dof*3*0+num_dof*i+0] = L1L2l3[i];
     5977                LFS[num_dof*3*0+num_dof*i+0] = basis[i];
    59785978                LFS[num_dof*3*0+num_dof*i+1] = 0.;
    59795979                LFS[num_dof*3*0+num_dof*i+2] = 0.;
    59805980                LFS[num_dof*3*1+num_dof*i+0] = 0.;
    5981                 LFS[num_dof*3*1+num_dof*i+1] = L1L2l3[i];
     5981                LFS[num_dof*3*1+num_dof*i+1] = basis[i];
    59825982                LFS[num_dof*3*1+num_dof*i+2] = 0.;
    59835983                LFS[num_dof*3*2+num_dof*i+0] = 0.;
    59845984                LFS[num_dof*3*2+num_dof*i+1] = 0.;
    5985                 LFS[num_dof*3*2+num_dof*i+2] = L1L2l3[i];
     5985                LFS[num_dof*3*2+num_dof*i+2] = basis[i];
    59865986                LFS[num_dof*3*3+num_dof*i+0] = 0.;
    59875987                LFS[num_dof*3*3+num_dof*i+1] = 0.;
    5988                 LFS[num_dof*3*3+num_dof*i+2] = L1L2l3[i];
     5988                LFS[num_dof*3*3+num_dof*i+2] = basis[i];
    59895989        }
    59905990}/*}}}*/
    … …  
    60006000         */
    60016001        int num_dof=2;
    6002         IssmDouble L1L2l3[3];
     6002        IssmDouble basis[3];
    60036003
    60046004        /*Cast gauss to GaussPenta*/
    … …  
    60066006        GaussPenta* gauss = dynamic_cast<GaussPenta*>(gauss_in);
    60076007
    6008         /*Get L1L2l3 in actual coordinate system: */
    6009         L1L2l3[0]=gauss->coord1*(1-gauss->coord4)/2.0;
    6010         L1L2l3[1]=gauss->coord2*(1-gauss->coord4)/2.0;
    6011         L1L2l3[2]=gauss->coord3*(1-gauss->coord4)/2.0;
     6008        /*Get basis in actual coordinate system: */
     6009        basis[0]=gauss->coord1*(1-gauss->coord4)/2.0;
     6010        basis[1]=gauss->coord2*(1-gauss->coord4)/2.0;
     6011        basis[2]=gauss->coord3*(1-gauss->coord4)/2.0;
    60126012
    60136013        /*Build LprimeFS: */
    60146014        for(int i=0;i<3;i++){
    6015                 LprimeFS[num_dof*3*0+num_dof*i+0] = L1L2l3[i];
     6015                LprimeFS[num_dof*3*0+num_dof*i+0] = basis[i];
    60166016                LprimeFS[num_dof*3*0+num_dof*i+1] = 0.;
    60176017                LprimeFS[num_dof*3*1+num_dof*i+0] = 0.;
    6018                 LprimeFS[num_dof*3*1+num_dof*i+1] = L1L2l3[i];
    6019                 LprimeFS[num_dof*3*2+num_dof*i+0] = L1L2l3[i];
     6018                LprimeFS[num_dof*3*1+num_dof*i+1] = basis[i];
     6019                LprimeFS[num_dof*3*2+num_dof*i+0] = basis[i];
    60206020                LprimeFS[num_dof*3*2+num_dof*i+1] = 0.;
    60216021                LprimeFS[num_dof*3*3+num_dof*i+0] = 0.;
    6022                 LprimeFS[num_dof*3*3+num_dof*i+1] = L1L2l3[i];
     6022                LprimeFS[num_dof*3*3+num_dof*i+1] = basis[i];
    60236023        }
    60246024}/*}}}*/

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

    r17962 r18075  
    883883        GetJacobianInvert(&Jinv[0][0], xyz_list, gauss);
    884884
    885         /*Build dh1dh3:
     885        /*Build basis function derivatives:
    886886         *
    887887         * [dhi/dx]= Jinv*[dhi/dr]

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

    r17926 r18075  
    2828                void GetNodalFunctionsDerivativesReference(IssmDouble* dbasis,Gauss* gauss);
    2929                void GetNodalFunctionsDerivativesReference(IssmDouble* dbasis,Gauss* gauss,int finiteelement);
    30                 void GetNodalFunctionsP1(IssmDouble* l1l6, Gauss* gauss);
    31                 void GetNodalFunctionsP1Derivatives(IssmDouble* dh1dh6,IssmDouble* xyz_list, Gauss* gauss);
    32                 void GetNodalFunctionsMINIDerivatives(IssmDouble* dh1dh7,IssmDouble* xyz_list, Gauss* gauss);
     30                void GetNodalFunctionsP1(IssmDouble* basis, Gauss* gauss);
     31                void GetNodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list, Gauss* gauss);
     32                void GetNodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list, Gauss* gauss);
    3333                void GetNodalFunctionsP1DerivativesReference(IssmDouble* dl1dl6,Gauss* gauss);
    3434                void GetNodalFunctionsMINIDerivativesReference(IssmDouble* dl1dl7,Gauss* gauss);
Note: See TracChangeset for help on using the changeset viewer.