1 // GetDP - Copyright (C) 1997-2021 P. Dular and C. Geuzaine, University of Liege 2 // 3 // See the LICENSE.txt file for license information. Please report all 4 // issues on https://gitlab.onelab.info/getdp/getdp/issues. 5 6 /* 1 integration point */ 7 8 static double xt1[1] = {0.333333333333333}; 9 static double yt1[1] = {0.333333333333333}; 10 static double pt1[1] = {0.5}; 11 12 /* 3 integration points */ 13 14 static double xt3[3] = {0.16666666666666,0.66666666666666,0.16666666666666}; 15 static double yt3[3] = {0.16666666666666,0.16666666666666,0.66666666666666}; 16 static double pt3[3] = {0.16666666666666,0.16666666666666,0.16666666666666}; 17 18 /* 4 integration points */ 19 20 static double xt4[4] = {0.333333333333333,0.6,0.2,0.2}; 21 static double yt4[4] = {0.3333333333333333,0.2,0.6,0.2}; 22 static double pt4[4] = {-0.28125,.260416666666,.260416666666,.260416666666}; 23 24 /* 6 integration points */ 25 26 static double xt6[6] = {0.816847572980459,0.091576213509771,0.091576213509771, 27 0.108103018168070,0.445948490915965,0.445948490915965}; 28 static double yt6[6] = {0.091576213509771,0.816847572980459,0.091576213509771, 29 0.445948490915965,0.108103018168070,0.445948490915965}; 30 static double pt6[6] = {0.054975871827661,0.054975871827661,0.054975871827661, 31 0.111690794839,0.111690794839,0.111690794839}; 32 33 /* 7 integration points */ 34 35 static double xt7[7] = {0.333333333333333,0.797426985353087,0.101286507323456, 36 0.101286507323456,0.470142064105115,0.059715871789770, 37 0.470142064105115}; 38 static double yt7[7] = {0.333333333333333,0.101286507323456,0.797426985353087, 39 0.101286507323456,0.059715871789770,0.470142064105115, 40 0.470142064105115}; 41 static double pt7[7] = {0.112500000000000,0.062969590272414,0.062969590272414, 42 0.062969590272414,0.066197076394253,0.066197076394253, 43 0.066197076394253}; 44 /* 12 integration points */ 45 46 static double xt12[12] = {0.873821971016996,0.063089014491502,0.063089014491502, 47 0.501426509658179,0.249286745170910,0.249286745170910, 48 0.636502499121399,0.310352451033785,0.636502499121399, 49 0.310352451033785,0.053145049844816,0.053145049844816}; 50 static double yt12[12] = {0.063089014491502,0.873821971016996,0.063089014491502, 51 0.249286745170910,0.501426509658179,0.249286745170910, 52 0.310352451033785,0.636502499121399,0.053145049844816, 53 0.053145049844816,0.310352451033785,0.636502499121399}; 54 static double pt12[12] = {0.025422453185104,0.025422453185104,0.025422453185104, 55 0.058393137863189,0.058393137863189,0.058393137863189, 56 0.041425537809187,0.041425537809187,0.041425537809187, 57 0.041425537809187,0.041425537809187,0.041425537809187}; 58 59 /* 13 integration points */ 60 61 static double xt13[13] = {0.333333333333333,0.479308067841920,0.260345966079040, 62 0.260345966079040,0.869739794195568,0.065130102902216, 63 0.065130102902216,0.048690315425316,0.312865496004874, 64 0.638444188569810,0.048690315425316,0.312865496004874, 65 0.638444188569810}; 66 static double yt13[13] = {0.333333333333333,0.260345966079040,0.479308067841920, 67 0.260345966079040,0.065130102902216,0.869739794195568, 68 0.065130102902216,0.312865496004874,0.048690315425316, 69 0.048690315425316,0.638444188569810,0.638444188569810, 70 0.312865496004874}; 71 static double pt13[13] = {-0.074785022233841,0.087807628716604,0.087807628716604, 72 0.087807628716604,0.026673617804419,0.026673617804419, 73 0.026673617804419,0.038556880445128,0.038556880445128, 74 0.038556880445128,0.038556880445128,0.038556880445128, 75 0.038556880445128}; 76 77 /* 16 integration points */ 78 79 static double xt16[16] = {0.333333333333333,0.081414823414554,0.459292588292723, 80 0.459292588292723,0.658861384496480,0.170569307751760, 81 0.170569307751760,0.898905543365938,0.050547228317031, 82 0.050547228317031,0.008394777409958,0.728492392955404, 83 0.263112829634638,0.008394777409958,0.263112829634638, 84 0.728492392955404}; 85 static double yt16[16] = {0.333333333333333,0.459292588292723,0.081414823414554, 86 0.459292588292723,0.170569307751760,0.658861384496480, 87 0.170569307751760,0.050547228317031,0.898905543365938, 88 0.050547228317031,0.728492392955404,0.008394777409958, 89 0.008394777409958,0.263112829634638,0.728492392955404, 90 0.263112829634638}; 91 static double pt16[16] = {0.072157803838894,0.047545817133643,0.047545817133643, 92 0.047545817133643,0.051608685267359,0.051608685267359, 93 0.051608685267359,0.016229248811599,0.016229248811599, 94 0.016229248811599,0.013615157087217,0.013615157087217, 95 0.013615157087217,0.013615157087217,0.013615157087217, 96 0.013615157087217}; 97 98 /* 99 GAUSS TRIANGLE WITH 1/R SINGULARITY OVER NODE (0,0,0) 100 101 ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia, 102 Some numerical integration formulae over triangles and 103 squares with a 1/R singularity, 104 Appl. Math. Modelling, Vol 5, June 1981, pp 209--211 105 */ 106 107 /* 1 integration point */ 108 109 static double xts1[1] = {0.25}; 110 static double yts1[1] = {0.25}; 111 static double pts1[1] = {1.24645048}; 112 113 /* 3 integration points */ 114 115 static double xts3[3] = {0.16666667,0.81742619,0.18257381}; 116 static double yts3[3] = {0.16666667,0.18257381,0.81742619}; 117 static double pts3[3] = {0.93483790,0.15580629,0.15580629}; 118 119 /* 4 integration points */ 120 121 static double xts4[4] = {0.16385495,0.61114353,0.04756957,0.17753138}; 122 static double yts4[4] = {0.04756957,0.17753138,0.16385495,0.61114353}; 123 static double pts4[4] = {0.31161231,0.31161293,0.31161231,0.31161293}; 124