diff --git a/Plugin/ShapeFunctions.h b/Plugin/ShapeFunctions.h index c4ebe19e8e9d556566c5bf624e852afbea0fac4e..68b1d1a3dcc906769793bfd260138a2088a6e2c0 100644 --- a/Plugin/ShapeFunctions.h +++ b/Plugin/ShapeFunctions.h @@ -172,9 +172,9 @@ public: inline int getNumGaussPoints(){ return 3; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u3[3] = {0.16666666666666,0.66666666666666,0.16666666666666}; - static double v3[3] = {0.16666666666666,0.16666666666666,0.66666666666666}; - static double p3[3] = {0.16666666666666,0.16666666666666,0.16666666666666}; + double u3[3] = {0.16666666666666,0.66666666666666,0.16666666666666}; + double v3[3] = {0.16666666666666,0.16666666666666,0.66666666666666}; + double p3[3] = {0.16666666666666,0.16666666666666,0.16666666666666}; if(num < 0 || num > 2) return; u = u3[num]; v = v3[num]; @@ -228,9 +228,9 @@ public: inline int getNumGaussPoints(){ return 4; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u4[4] = {0.577350269189,-0.577350269189,0.577350269189,-0.577350269189}; - static double v4[4] = {0.577350269189,0.577350269189,-0.577350269189,-0.577350269189}; - static double p4[4] = {1.,1.,1.,1.}; + double u4[4] = {0.577350269189,-0.577350269189,0.577350269189,-0.577350269189}; + double v4[4] = {0.577350269189,0.577350269189,-0.577350269189,-0.577350269189}; + double p4[4] = {1.,1.,1.,1.}; if(num < 0 || num > 3) return; u = u4[num]; v = v4[num]; @@ -287,10 +287,10 @@ public: inline int getNumGaussPoints(){ return 4; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u4[4] = {0.138196601125,0.138196601125,0.138196601125,0.585410196625}; - static double v4[4] = {0.138196601125,0.138196601125,0.585410196625,0.138196601125}; - static double w4[4] = {0.138196601125,0.585410196625,0.138196601125,0.138196601125}; - static double p4[4] = {0.0416666666667,0.0416666666667,0.0416666666667,0.0416666666667}; + double u4[4] = {0.138196601125,0.138196601125,0.138196601125,0.585410196625}; + double v4[4] = {0.138196601125,0.138196601125,0.585410196625,0.138196601125}; + double w4[4] = {0.138196601125,0.585410196625,0.138196601125,0.138196601125}; + double p4[4] = {0.0416666666667,0.0416666666667,0.0416666666667,0.0416666666667}; if(num < 0 || num > 3) return; u = u4[num]; v = v4[num]; @@ -327,14 +327,14 @@ public: inline int getNumGaussPoints(){ return 6; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u6[6] = { 0.40824826, 0.40824826, -0.40824826, - -0.40824826, -0.81649658, 0.81649658}; - static double v6[6] = { 0.70710678, -0.70710678, 0.70710678, - -0.70710678, 0., 0.}; - static double w6[6] = {-0.57735027, -0.57735027, 0.57735027, - 0.57735027, -0.57735027, 0.57735027}; - static double p6[6] = { 1.3333333333, 1.3333333333, 1.3333333333, - 1.3333333333, 1.3333333333, 1.3333333333}; + double u6[6] = { 0.40824826, 0.40824826, -0.40824826, + -0.40824826, -0.81649658, 0.81649658}; + double v6[6] = { 0.70710678, -0.70710678, 0.70710678, + -0.70710678, 0., 0.}; + double w6[6] = {-0.57735027, -0.57735027, 0.57735027, + 0.57735027, -0.57735027, 0.57735027}; + double p6[6] = { 1.3333333333, 1.3333333333, 1.3333333333, + 1.3333333333, 1.3333333333, 1.3333333333}; if(num < 0 || num > 5) return; u = u6[num]; v = v6[num]; @@ -395,14 +395,14 @@ public: inline int getNumGaussPoints(){ return 6; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u6[6] = {0.166666666666666, 0.333333333333333, 0.166666666666666, - 0.166666666666666, 0.333333333333333, 0.166666666666666}; - static double v6[6] = {0.166666666666666, 0.166666666666666, 0.333333333333333, - 0.166666666666666, 0.166666666666666, 0.333333333333333}; - static double w6[6] = {-0.577350269189, -0.577350269189, -0.577350269189, - 0.577350269189, 0.577350269189, 0.577350269189}; - static double p6[6] = {0.166666666666666,0.166666666666666,0.166666666666666, - 0.166666666666666,0.166666666666666,0.166666666666666,}; + double u6[6] = {0.166666666666666, 0.333333333333333, 0.166666666666666, + 0.166666666666666, 0.333333333333333, 0.166666666666666}; + double v6[6] = {0.166666666666666, 0.166666666666666, 0.333333333333333, + 0.166666666666666, 0.166666666666666, 0.333333333333333}; + double w6[6] = {-0.577350269189, -0.577350269189, -0.577350269189, + 0.577350269189, 0.577350269189, 0.577350269189}; + double p6[6] = {0.166666666666666,0.166666666666666,0.166666666666666, + 0.166666666666666,0.166666666666666,0.166666666666666,}; if(num < 0 || num > 5) return; u = u6[num]; v = v6[num]; @@ -455,22 +455,22 @@ public: inline int getNumGaussPoints(){ return 8; } void getGaussPoint(int num, double &u, double &v, double &w, double &weight) { - static double u8[8] = {0.3595161057791018,0.09633205020967324, - 0.3595161057791018,0.09633205020967324, - 0.6920507403468987,0.1854344369976602, - 0.6920507403468987,0.1854344369976602}; - static double v8[8] = {0.3595161057791018,0.3595161057791018, - 0.09633205020967324,0.09633205020967324, - 0.6920507403468987,0.6920507403468987, - 0.1854344369976602,0.1854344369976602}; - static double w8[8] = {0.544151844011225,0.544151844011225, - 0.544151844011225,0.544151844011225, - 0.122514822655441,0.122514822655441, - 0.122514822655441,0.122514822655441}; - static double p8[8] = {0.02519647051995625,0.02519647051995625, - 0.02519647051995625,0.02519647051995625, - 0.058136862813377,0.058136862813377, - 0.058136862813377,0.058136862813377}; + double u8[8] = {0.3595161057791018,0.09633205020967324, + 0.3595161057791018,0.09633205020967324, + 0.6920507403468987,0.1854344369976602, + 0.6920507403468987,0.1854344369976602}; + double v8[8] = {0.3595161057791018,0.3595161057791018, + 0.09633205020967324,0.09633205020967324, + 0.6920507403468987,0.6920507403468987, + 0.1854344369976602,0.1854344369976602}; + double w8[8] = {0.544151844011225,0.544151844011225, + 0.544151844011225,0.544151844011225, + 0.122514822655441,0.122514822655441, + 0.122514822655441,0.122514822655441}; + double p8[8] = {0.02519647051995625,0.02519647051995625, + 0.02519647051995625,0.02519647051995625, + 0.058136862813377,0.058136862813377, + 0.058136862813377,0.058136862813377}; if(num < 0 || num > 7) return; u = u8[num]; v = v8[num];