/*
* Copyright (C) 1997-2006 P. Dular, C. Geuzaine
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
* USA.
*
* Please report all bugs and problems to <getdp@geuz.org>.
*/
/* 6 integration points (DHATT ET TOUZOT) */
double upri6[6] = {0.166666666666666, 0.333333333333333, 0.166666666666666,
0.166666666666666, 0.333333333333333, 0.166666666666666};
double vpri6[6] = {0.166666666666666, 0.166666666666666, 0.333333333333333,
0.166666666666666, 0.166666666666666, 0.333333333333333};
double wpri6[6] = {-0.577350269189, -0.577350269189, -0.577350269189,
0.577350269189, 0.577350269189, 0.577350269189};
double ppri6[6] = {0.166666666666666,0.166666666666666,0.166666666666666,
0.166666666666666,0.166666666666666,0.166666666666666,};
/* 9 integration points (DHATT ET TOUZOT) */
double upri9[9] = {0.166666666666666, 0.333333333333333, 0.166666666666666,
0.166666666666666, 0.333333333333333, 0.166666666666666,
0.166666666666666, 0.333333333333333, 0.166666666666666};
double vpri9[9] = {0.166666666666666, 0.166666666666666, 0.333333333333333,
0.166666666666666, 0.166666666666666, 0.333333333333333,
0.166666666666666, 0.166666666666666, 0.333333333333333};
double wpri9[9] = {-0.774596669241483, -0.774596669241483, -0.774596669241483,
0. , 0. , 0. ,
0.774596669241483, 0.774596669241483, 0.774596669241483};
double ppri9[9] = {0.0925925925929,0.0925925925929,0.0925925925929,
0.1481481481481,0.1481481481481,0.1481481481481,
0.0925925925929,0.0925925925929,0.0925925925929};
/* 21 integration points (DHATT ET TOUZOT) */
double upri21[21] = {0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456};
double vpri21[21] = {0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535};
double wpri21[21] = {-0.774596669241483,-0.774596669241483,
-0.774596669241483,-0.774596669241483,
-0.774596669241483,-0.774596669241483,
-0.774596669241483,
0. , 0. ,
0. , 0. ,
0. , 0. ,
0. ,
0.774596669241483, 0.774596669241483,
0.774596669241483, 0.774596669241483,
0.774596669241483, 0.774596669241483,
0.774596669241483};
double ppri21[21] = {0.625000000000000E-01, 0.367761535523628E-01,
0.367761535523628E-01, 0.367761535523628E-01,
0.349831057068964E-01, 0.349831057068964E-01,
0.349831057068964E-01,
0.100000000000000E+00, 0.588418456837804E-01,
0.588418456837804E-01, 0.588418456837804E-01,
0.559729691310342E-01, 0.559729691310342E-01,
0.559729691310342E-01,
0.625000000000000E-01, 0.367761535523628E-01,
0.367761535523628E-01, 0.367761535523628E-01,
0.349831057068964E-01, 0.349831057068964E-01,
0.349831057068964E-01};
/* 42 integration points (DHATT ET TOUZOT) */
double upri42[42] = {0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456,
0.333333333333333, 0.470142064105115,
0.05971587179 , 0.470142064105115,
0.101286507323456, 0.79742698535 ,
0.101286507323456};
double vpri42[42] = {0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535 ,
0.333333333333333, 0.470142064105115,
0.470142064105115, 0.05971587179 ,
0.101286507323456, 0.101286507323456,
0.79742698535};
double wpri42[42] = {-0.238619186083197,
-0.238619186083197,-0.238619186083197,
-0.238619186083197,-0.238619186083197,
-0.238619186083197,-0.238619186083197,
0.238619186083197,
0.238619186083197, 0.238619186083197,
0.238619186083197, 0.238619186083197,
0.238619186083197, 0.238619186083197,
-0.661209386466265,
-0.661209386466265,-0.661209386466265,
-0.661209386466265,-0.661209386466265,
-0.661209386466265,-0.661209386466265,
0.661209386466265,
0.661209386466265, 0.661209386466265,
0.661209386466265, 0.661209386466265,
0.661209386466265, 0.661209386466265,
-0.932469514203152,
-0.932469514203152,-0.932469514203152,
-0.932469514203152,-0.932469514203152,
-0.932469514203152,-0.932469514203152,
0.932469514203152,
0.932469514203152, 0.932469514203152,
0.932469514203152, 0.932469514203152,
0.932469514203152, 0.932469514203152};
double ppri42[42] = {0.5264031763942775E-01, 0.3097453447284393E-01,
0.3097453447284393E-01, 0.3097453447284393E-01,
0.2946434874279526E-01, 0.2946434874279526E-01,
0.2946434874279526E-01,
0.5264031763942775E-01, 0.3097453447284393E-01,
0.3097453447284393E-01, 0.3097453447284393E-01,
0.2946434874279526E-01, 0.2946434874279526E-01,
0.2946434874279526E-01,
0.4058567696791564E-01, 0.2388136141117854E-01,
0.2388136141117854E-01, 0.2388136141117854E-01,
0.2271700844087269E-01, 0.2271700844087269E-01,
0.2271700844087269E-01,
0.4058567696791564E-01, 0.2388136141117854E-01,
0.2388136141117854E-01, 0.2388136141117854E-01,
0.2271700844087269E-01, 0.2271700844087269E-01,
0.2271700844087269E-01,
0.1927400539265662E-01, 0.1134118051023053E-01,
0.1134118051023053E-01, 0.1134118051023053E-01,
0.1078823308874557E-01, 0.1078823308874557E-01,
0.1078823308874557E-01,
0.1927400539265662E-01, 0.1134118051023053E-01,
0.1134118051023053E-01, 0.1134118051023053E-01,
0.1078823308874557E-01, 0.1078823308874557E-01,
0.1078823308874557E-01};
syntax highlighted by Code2HTML, v. 0.9.1