Ignore:
Timestamp:
09/14/10 16:08:56 (15 years ago)
Author:
Mathieu Morlighem
Message:

New Icefront MUCH cleaner. Enjoy

File:
1 edited

Legend:

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

    r5797 r5808  
    179179
    180180        /*Quads: get the gauss points using the product of two line rules  */
    181         if(index1==0 && index2==1 && index3==3 && index4==4){
     181        if(index1==0 && index2==1 && index3==4 && index4==3){
    182182                for(i=0;i<order_horiz;i++){
    183183                        for(j=0;j<order_vert;j++){
    184                                 coords1[i*order_horiz+j]=  0.5*(1-seg_horiz_coords[i]);
    185                                 coords2[i*order_horiz+j]=1-0.5*(1-seg_horiz_coords[i]);
    186                                 coords3[i*order_horiz+j]=0.0;
    187                                 coords4[i*order_horiz+j]=seg_vert_coords[j];
    188                                 weights[i*order_horiz+j]=seg_horiz_weights[i]*seg_vert_weights[j];
     184                                coords1[i*order_vert+j]=  0.5*(1-seg_horiz_coords[i]);
     185                                coords2[i*order_vert+j]=1-0.5*(1-seg_horiz_coords[i]);
     186                                coords3[i*order_vert+j]=0.0;
     187                                coords4[i*order_vert+j]=seg_vert_coords[j];
     188                                weights[i*order_vert+j]=seg_horiz_weights[i]*seg_vert_weights[j];
    189189                        }
    190190                }
    191191        }
    192         else if(index1==1 && index2==2 && index3==4 && index4==5){
     192        else if(index1==1 && index2==2 && index3==5 && index4==4){
    193193                for(i=0;i<order_horiz;i++){
    194194                        for(j=0;j<order_vert;j++){
    195                                 coords1[i*order_horiz+j]=0.0;
    196                                 coords2[i*order_horiz+j]=  0.5*(1-seg_horiz_coords[i]);
    197                                 coords3[i*order_horiz+j]=1-0.5*(1-seg_horiz_coords[i]);
    198                                 coords4[i*order_horiz+j]=seg_vert_coords[j];
    199                                 weights[i*order_horiz+j]=seg_horiz_weights[i]*seg_vert_weights[j];
     195                                coords1[i*order_vert+j]=0.0;
     196                                coords2[i*order_vert+j]=  0.5*(1-seg_horiz_coords[i]);
     197                                coords3[i*order_vert+j]=1-0.5*(1-seg_horiz_coords[i]);
     198                                coords4[i*order_vert+j]=seg_vert_coords[j];
     199                                weights[i*order_vert+j]=seg_horiz_weights[i]*seg_vert_weights[j];
    200200                        }
    201201                }
    202202        }
    203         else if(index1==2 && index2==0 && index3==5 && index4==3){
     203        else if(index1==2 && index2==0 && index3==3 && index4==5){
    204204                for(i=0;i<order_horiz;i++){
    205205                        for(j=0;j<order_vert;j++){
    206                                 coords1[i*order_horiz+j]=1-0.5*(1-seg_horiz_coords[i]);
    207                                 coords2[i*order_horiz+j]=0.0;
    208                                 coords3[i*order_horiz+j]=  0.5*(1-seg_horiz_coords[i]);
    209                                 coords4[i*order_horiz+j]=seg_vert_coords[j];
    210                                 weights[i*order_horiz+j]=seg_horiz_weights[i]*seg_vert_weights[j];
     206                                coords1[i*order_vert+j]=1-0.5*(1-seg_horiz_coords[i]);
     207                                coords2[i*order_vert+j]=0.0;
     208                                coords3[i*order_vert+j]=  0.5*(1-seg_horiz_coords[i]);
     209                                coords4[i*order_vert+j]=seg_vert_coords[j];
     210                                weights[i*order_vert+j]=seg_horiz_weights[i]*seg_vert_weights[j];
    211211                        }
    212212                }
     
    216216        }
    217217
     218        /*clean-up*/
     219        xfree((void**)&seg_horiz_coords);
     220        xfree((void**)&seg_horiz_weights);
     221        xfree((void**)&seg_vert_coords);
     222        xfree((void**)&seg_vert_weights);
    218223}
    219224/*}}}*/
Note: See TracChangeset for help on using the changeset viewer.