Ignore:
Timestamp:
12/24/13 01:40:32 (11 years ago)
Author:
Mathieu Morlighem
Message:

NEW: working on P2xP4 finite elements

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/gauss/GaussPenta.cpp

    r17041 r17057  
    657657                        }
    658658                        break;
     659                case P2xP4Enum:
     660                        switch(iv){
     661                                case 0: coord1=1.; coord2=0.; coord3=0.; coord4=-1.; break;
     662                                case 1: coord1=0.; coord2=1.; coord3=0.; coord4=-1.; break;
     663                                case 2: coord1=0.; coord2=0.; coord3=1.; coord4=-1.; break;
     664                                case 3: coord1=1.; coord2=0.; coord3=0.; coord4=+1.; break;
     665                                case 4: coord1=0.; coord2=1.; coord3=0.; coord4=+1.; break;
     666                                case 5: coord1=0.; coord2=0.; coord3=1.; coord4=+1.; break;
     667
     668                                case 6: coord1=1.; coord2=0.; coord3=0.; coord4=0.; break;
     669                                case 7: coord1=0.; coord2=1.; coord3=0.; coord4=0.; break;
     670                                case 8: coord1=0.; coord2=0.; coord3=1.; coord4=0.; break;
     671
     672                                case  9: coord1=0.; coord2=.5; coord3=.5; coord4=-1.;break;
     673                                case 10: coord1=.5; coord2=0.; coord3=.5; coord4=-1.;break;
     674                                case 11: coord1=.5; coord2=.5; coord3=0.; coord4=-1.;break;
     675                                case 12: coord1=0.; coord2=.5; coord3=.5; coord4=+1.;break;
     676                                case 13: coord1=.5; coord2=0.; coord3=.5; coord4=+1.;break;
     677                                case 14: coord1=.5; coord2=.5; coord3=0.; coord4=+1.;break;
     678
     679                                case 15: coord1=1.; coord2=0.; coord3=0.; coord4=-.5; break;
     680                                case 16: coord1=0.; coord2=1.; coord3=0.; coord4=-.5; break;
     681                                case 17: coord1=0.; coord2=0.; coord3=1.; coord4=-.5; break;
     682                                case 18: coord1=1.; coord2=0.; coord3=0.; coord4=+.5; break;
     683                                case 19: coord1=0.; coord2=1.; coord3=0.; coord4=+.5; break;
     684                                case 20: coord1=0.; coord2=0.; coord3=1.; coord4=+.5; break;
     685
     686                                case 21: coord1=0.; coord2=.5; coord3=.5; coord4=-.5;break;
     687                                case 22: coord1=.5; coord2=0.; coord3=.5; coord4=-.5;break;
     688                                case 23: coord1=.5; coord2=.5; coord3=0.; coord4=-.5;break;
     689                                case 24: coord1=0.; coord2=.5; coord3=.5; coord4=0.;break;
     690                                case 25: coord1=.5; coord2=0.; coord3=.5; coord4=0.;break;
     691                                case 26: coord1=.5; coord2=.5; coord3=0.; coord4=0.;break;
     692                                case 27: coord1=0.; coord2=.5; coord3=.5; coord4=+.5;break;
     693                                case 28: coord1=.5; coord2=0.; coord3=.5; coord4=+.5;break;
     694                                case 29: coord1=.5; coord2=.5; coord3=0.; coord4=+.5;break;
     695                                default: _error_("node index should be in [0 14]");
     696                        }
     697                        break;
    659698                default: _error_("Finite element "<<EnumToStringx(finiteelement)<<" not supported");
    660699        }
Note: See TracChangeset for help on using the changeset viewer.