Mechanics of Defect Evolution Library
Loading...
Searching...
No Matches
GaussLegendre_2_37.h
Go to the documentation of this file.
1 /* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
2 *
3 *
4 * MoDELib is distributed without any warranty under the
5 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
6 */
7
8#ifndef model_GAUSSLEGENDRE_2_37_H_
9#define model_GAUSSLEGENDRE_2_37_H_
10
11namespace model{
12
14 template <>
15 struct GaussLegendre<2,37>{
16 // from http://people.sc.fsu.edu/~jburkardt/datasets/quadrature_rules_tri/quadrature_rules_tri.html
17 static Eigen::Matrix<double,3,37> abcsissasAndWeights(){
18 Eigen::Matrix<double,37,3> U;
19 U<<3.333333333333330e-01,3.333333333333330e-01,2.586988303287200e-02,
20 9.502756629241060e-01,2.486216853794700e-02,4.003899777782500e-03,
21 2.486216853794700e-02,9.502756629241060e-01,4.003899777782500e-03,
22 2.486216853794700e-02,2.486216853794700e-02,4.003899777782500e-03,
23 1.716149149238350e-01,4.141925425380820e-01,2.343444949091100e-02,
24 4.141925425380820e-01,1.716149149238350e-01,2.343444949091100e-02,
25 4.141925425380820e-01,4.141925425380820e-01,2.343444949091100e-02,
26 5.394122436771900e-01,2.302938781614050e-01,2.329547009198800e-02,
27 2.302938781614050e-01,5.394122436771900e-01,2.329547009198800e-02,
28 2.302938781614050e-01,2.302938781614050e-01,2.329547009198800e-02,
29 7.721600366765330e-01,1.139199816617340e-01,1.550847165689800e-02,
30 1.139199816617340e-01,7.721600366765330e-01,1.550847165689800e-02,
31 1.139199816617340e-01,1.139199816617340e-01,1.550847165689800e-02,
32 9.085399949835000e-03,4.954573000250820e-01,5.395806368315500e-03,
33 4.954573000250820e-01,9.085399949835000e-03,5.395806368315500e-03,
34 4.954573000250820e-01,4.954573000250820e-01,5.395806368315500e-03,
35 6.227729030588700e-02,4.688613548470560e-01,1.609776712121600e-02,
36 4.688613548470560e-01,6.227729030588700e-02,1.609776712121600e-02,
37 4.688613548470560e-01,4.688613548470560e-01,1.609776712121600e-02,
38 2.207628965362400e-02,8.513065041743481e-01,7.722917105351000e-03,
39 2.207628965362400e-02,1.266172061720270e-01,7.722917105351000e-03,
40 8.513065041743481e-01,2.207628965362400e-02,7.722917105351000e-03,
41 8.513065041743481e-01,1.266172061720270e-01,7.722917105351000e-03,
42 1.266172061720270e-01,2.207628965362400e-02,7.722917105351000e-03,
43 1.266172061720270e-01,8.513065041743481e-01,7.722917105351000e-03,
44 1.862052280252100e-02,6.894419707285910e-01,8.911494961589500e-03,
45 1.862052280252100e-02,2.919375064688880e-01,8.911494961589500e-03,
46 6.894419707285910e-01,1.862052280252100e-02,8.911494961589500e-03,
47 6.894419707285910e-01,2.919375064688880e-01,8.911494961589500e-03,
48 2.919375064688880e-01,1.862052280252100e-02,8.911494961589500e-03,
49 2.919375064688880e-01,6.894419707285910e-01,8.911494961589500e-03,
50 9.650648129215900e-02,6.358678594338730e-01,1.851934184069250e-02,
51 9.650648129215900e-02,2.676256592739680e-01,1.851934184069250e-02,
52 6.358678594338730e-01,9.650648129215900e-02,1.851934184069250e-02,
53 6.358678594338730e-01,2.676256592739680e-01,1.851934184069250e-02,
54 2.676256592739680e-01,9.650648129215900e-02,1.851934184069250e-02,
55 2.676256592739680e-01,6.358678594338730e-01,1.851934184069250e-02;
56// 0.333333333333333, 0.333333333333333, 0.051739766065744,
57// 0.950275662924106, 0.024862168537947, 0.008007799555565,
58// 0.024862168537947, 0.950275662924106, 0.008007799555565,
59// 0.024862168537947, 0.024862168537947, 0.008007799555565,
60// 0.171614914923835, 0.414192542538082, 0.046868898981822,
61// 0.414192542538082, 0.171614914923835, 0.046868898981822,
62// 0.414192542538082, 0.414192542538082, 0.046868898981822,
63// 0.539412243677190, 0.230293878161405, 0.046590940183976,
64// 0.230293878161405, 0.539412243677190, 0.046590940183976,
65// 0.230293878161405, 0.230293878161405, 0.046590940183976,
66// 0.772160036676533, 0.113919981661734, 0.031016943313796,
67// 0.113919981661734, 0.772160036676533, 0.031016943313796,
68// 0.113919981661734, 0.113919981661734, 0.031016943313796,
69// 0.009085399949835, 0.495457300025082, 0.010791612736631,
70// 0.495457300025082, 0.009085399949835, 0.010791612736631,
71// 0.495457300025082, 0.495457300025082, 0.010791612736631,
72// 0.062277290305887, 0.468861354847056, 0.032195534242432,
73// 0.468861354847056, 0.062277290305887, 0.032195534242432,
74// 0.468861354847056, 0.468861354847056, 0.032195534242432,
75// 0.022076289653624, 0.851306504174348, 0.015445834210702,
76// 0.022076289653624, 0.126617206172027, 0.015445834210702,
77// 0.851306504174348, 0.022076289653624, 0.015445834210702,
78// 0.851306504174348, 0.126617206172027, 0.015445834210702,
79// 0.126617206172027, 0.022076289653624, 0.015445834210702,
80// 0.126617206172027, 0.851306504174348, 0.015445834210702,
81// 0.018620522802521, 0.689441970728591, 0.017822989923179,
82// 0.018620522802521, 0.291937506468888, 0.017822989923179,
83// 0.689441970728591, 0.018620522802521, 0.017822989923179,
84// 0.689441970728591, 0.291937506468888, 0.017822989923179,
85// 0.291937506468888, 0.018620522802521, 0.017822989923179,
86// 0.291937506468888, 0.689441970728591, 0.017822989923179,
87// 0.096506481292159, 0.635867859433873, 0.037038683681385,
88// 0.096506481292159, 0.267625659273968, 0.037038683681385,
89// 0.635867859433873, 0.096506481292159, 0.037038683681385,
90// 0.635867859433873, 0.267625659273968, 0.037038683681385,
91// 0.267625659273968, 0.096506481292159, 0.037038683681385,
92// 0.267625659273968, 0.635867859433873, 0.037038683681385;
93 return U.transpose();
94 }
95 };
96/*************************************************/
97}
98#endif
99
static Eigen::Matrix< double, 3, 37 > abcsissasAndWeights()
Class template defining the GaussLegendre rules for determination of quadrature abscissas and weights...