Skip to content
Snippets Groups Projects
Select Git revision
  • 9ca500f78387ba9ff289b58d0c2683810ef2dae6
  • master default protected
  • dof-renumbering
  • test-dof-hash
  • gdemesy-master-patch-30528
  • eval-space-time
  • oscillating_multiharm
  • MH_movement
  • axisqu
  • write_vtu_and_ensight_formats
  • movingband
  • CP_1972_add_vtu_file_writing
  • mortar
  • fast_freq_sweep_Resolution
  • applyresolvent_again
  • marteaua-master-patch-54323
  • patch-1
  • binde-master-patch-08072
  • binde-master-patch-52461
  • BCGSL
  • resolvent
  • getdp_3_5_0
  • getdp_3_4_0
  • getdp_3_3_0
  • getdp_3_2_0
  • getdp_3_1_0
  • getdp_3_0_4
  • getdp_3_0_3
  • getdp_3_0_2
  • getdp_3_0_1
  • getdp_3_0_0
  • onelab_mobile_2.1.0
  • getdp_2_11_3 protected
  • getdp_2_11_2 protected
  • getdp_2_11_1 protected
  • getdp_2_11_0 protected
  • getdp_2_10_0 protected
  • getdp_2_9_2 protected
  • getdp_2_9_1 protected
  • getdp_2_9_0 protected
  • getdp_2_8_0 protected
41 results

Gauss_Tetrahedron.h

Blame
  • Gauss_Tetrahedron.h 5.03 KiB
    // GetDP - Copyright (C) 1997-2018 P. Dular and C. Geuzaine, University of Liege
    //
    // See the LICENSE.txt file for license information. Please report all
    // bugs and problems to the public mailing list <getdp@onelab.info>.
    
    /* 1 integration point */
    
    static double xtet1[1] = {0.25};
    static double ytet1[1] = {0.25};
    static double ztet1[1] = {0.25};
    static double ptet1[1] = {0.166666666667};
    
    /* 4 integration points */
    
    static double xtet4[4] = {0.138196601125,0.138196601125,0.138196601125,0.585410196625};
    static double ytet4[4] = {0.138196601125,0.138196601125,0.585410196625,0.138196601125};
    static double ztet4[4] = {0.138196601125,0.585410196625,0.138196601125,0.138196601125};
    static double ptet4[4] = {0.0416666666667,0.0416666666667,0.0416666666667,0.0416666666667};
    
    /* 5 integration points */
    
    static double xtet5[5] = {0.25,0.166666666667,0.166666666667,0.166666666667,0.5};
    static double ytet5[5] = {0.25,0.166666666667,0.166666666667,0.5,0.166666666667};
    static double ztet5[5] = {0.25,0.166666666667,0.5,0.166666666667,0.166666666667};
    static double ptet5[5] = {-0.133333333333,0.075,0.075,0.075,0.075};
    
    /* 15 integration points */
    
    static double xtet15[15] = {0.25,
    			    0.0919710780526,0.0919710780526,0.0919710780526,0.724086765841,
    			    0.319793627829,0.319793627829,0.319793627829,0.0406191165118,
    			    0.0563508326895,0.0563508326895,0.44364916731,0.0563508326895,
    			    0.44364916731,0.44364916731};
    static double ytet15[15] = {0.25,
    			    0.0919710780526,0.0919710780526,0.724086765841,0.0919710780526,
    			    0.319793627829,0.319793627829,0.0406191165118,0.319793627829,
    			    0.0563508326895,0.44364916731,0.0563508326895,0.44364916731,
    			    0.0563508326895,0.44364916731};
    static double ztet15[15] = {0.25,
    			    0.0919710780526,0.724086765841,0.0919710780526,0.0919710780526,
    			    0.319793627829,0.0406191165118,0.319793627829,0.319793627829,
    			    0.44364916731,0.0563508326895,0.0563508326895,0.44364916731,
    			    0.44364916731,0.0563508326895};
    static double ptet15[15] = {0.0197530864198,
    			    0.0119895139632,0.0119895139632,0.0119895139632,0.0119895139632,  
    			    0.011511367871,0.011511367871,0.011511367871,0.011511367871,
    			    0.00881834215168,0.00881834215168,0.00881834215168,0.00881834215168,
    			    0.00881834215168,0.00881834215168};
    
    /* 16 integration points */
    
    #define a16  0.0503737941001228 / 6.0  
    #define b16  0.0665420686332923 / 6.0
    #define c16  0.7716429020672371  
    #define d16  0.0761190326442543  
    #define e16  0.1197005277978019  
    #define f16  0.0718316452676693  
    #define g16  0.4042339134672644  
    
    static double xtet16[16] = {c16,d16,d16,d16,e16,f16,e16,f16,
    			    g16,g16,g16,g16,e16,f16,g16,g16};
    static double ytet16[16] = {d16,c16,d16,d16,f16,e16,g16,g16,
    			    g16,g16,e16,f16,g16,g16,e16,f16};
    static double ztet16[16] = {d16,d16,c16,d16,g16,g16,g16,g16,
    			    e16,f16,f16,e16,f16,e16,g16,g16};
    static double ptet16[16] = {a16,a16,a16,a16,b16,b16,b16,b16,
    			    b16,b16,b16,b16,b16,b16,b16,b16};
    #undef a16  
    #undef b16  
    #undef c16