From 4c8801e20eb1eb1ef2751c436a7b8788940e1125 Mon Sep 17 00:00:00 2001 From: "Michael H. Scott" Date: Sat, 16 Nov 2024 10:51:06 -0800 Subject: [PATCH 1/3] Using 16 digits in locations and weights - generated from gaussq.f --- .../LobattoBeamIntegration.cpp | 262 +++++++++--------- .../forceBeamColumn/RadauBeamIntegration.cpp | 262 +++++++++--------- 2 files changed, 262 insertions(+), 262 deletions(-) diff --git a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp index 3f28105282..365f189e0e 100644 --- a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp @@ -111,87 +111,87 @@ LobattoBeamIntegration::getSectionLocations(int numSections, double *xi) { switch(numSections) { - + case 2: - xi[0] = -1.0; - xi[1] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = 1.00000000000000e+00; break; case 3: - xi[0] = -1.0; - xi[1] = 0.0; - xi[2] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = 0.0; + xi[2] = 1.00000000000000e+00; break; case 4: - xi[0] = -1.0; - xi[1] = -0.44721360; - xi[2] = 0.44721360; - xi[3] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -4.47213595499958e-01; + xi[2] = 4.47213595499958e-01; + xi[3] = 1.00000000000000e+00; break; case 5: - xi[0] = -1.0; - xi[1] = -0.65465367; - xi[2] = 0.0; - xi[3] = 0.65465367; - xi[4] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -6.54653670707977e-01; + xi[2] = 0.0; + xi[3] = 6.54653670707977e-01; + xi[4] = 1.00000000000000e+00; break; case 6: - xi[0] = -1.0; - xi[1] = -0.7650553239; - xi[2] = -0.2852315164; - xi[3] = 0.2852315164; - xi[4] = 0.7650553239; - xi[5] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -7.65055323929465e-01; + xi[2] = -2.85231516480645e-01; + xi[3] = 2.85231516480645e-01; + xi[4] = 7.65055323929465e-01; + xi[5] = 1.00000000000000e+00; break; case 7: - xi[0] = -1.0; - xi[1] = -0.8302238962; - xi[2] = -0.4688487934; - xi[3] = 0.0; - xi[4] = 0.4688487934; - xi[5] = 0.8302238962; - xi[6] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.30223896278567e-01; + xi[2] = -4.68848793470714e-01; + xi[3] = 0.0; + xi[4] = 4.68848793470714e-01; + xi[5] = 8.30223896278567e-01; + xi[6] = 1.00000000000000e+00; break; - + case 8: - xi[0] = -1.0; - xi[1] = -0.8717401485; - xi[2] = -0.5917001814; - xi[3] = -0.2092992179; - xi[4] = 0.2092992179; - xi[5] = 0.5917001814; - xi[6] = 0.8717401485; - xi[7] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.71740148509606e-01; + xi[2] = -5.91700181433143e-01; + xi[3] = -2.09299217902479e-01; + xi[4] = 2.09299217902479e-01; + xi[5] = 5.91700181433143e-01; + xi[6] = 8.71740148509607e-01; + xi[7] = 1.00000000000000e+00; break; case 9: - xi[0] = -1.0; - xi[1] = -0.8997579954; - xi[2] = -0.6771862795; - xi[3] = -0.3631174638; - xi[4] = 0.0; - xi[5] = 0.3631174638; - xi[6] = 0.6771862795; - xi[7] = 0.8997579954; - xi[8] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.99757995411461e-01; + xi[2] = -6.77186279510738e-01; + xi[3] = -3.63117463826178e-01; + xi[4] = 0.0; + xi[5] = 3.63117463826178e-01; + xi[6] = 6.77186279510738e-01; + xi[7] = 8.99757995411461e-01; + xi[8] = 1.00000000000000e+00; break; - + case 10: - xi[0] = -1.0; - xi[1] = -0.9195339082; - xi[2] = -0.7387738651; - xi[3] = -0.4779249498; - xi[4] = -0.1652789577; - xi[5] = 0.1652789577; - xi[6] = 0.4779249498; - xi[7] = 0.7387738651; - xi[8] = 0.9195339082; - xi[9] = 1.0; - break; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.19533908166459e-01; + xi[2] = -7.38773865105505e-01; + xi[3] = -4.77924949810445e-01; + xi[4] = -1.65278957666387e-01; + xi[5] = 1.65278957666387e-01; + xi[6] = 4.77924949810445e-01; + xi[7] = 7.38773865105505e-01; + xi[8] = 9.19533908166459e-01; + xi[9] = 1.00000000000000e+00; + break; default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; @@ -210,86 +210,86 @@ LobattoBeamIntegration::getSectionWeights(int numSections, double L, switch (numSections) { case 2: - wt[0] = 1.0; - wt[1] = 1.0; - break; + wt[0] = 1.00000000000000e+00; + wt[1] = 1.00000000000000e+00; + break; // sum = 2.00000000000000e+00 case 3: - wt[0] = 0.333333333333333; - wt[1] = 1.333333333333333; - wt[2] = 0.333333333333333; - break; + wt[0] = 3.33333333333333e-01; + wt[1] = 1.33333333333333e+00; + wt[2] = 3.33333333333334e-01; + break; // sum = 2.00000000000000e+00 - case 4: - wt[0] = 0.166666666666667; - wt[1] = 0.833333333333333; - wt[2] = 0.833333333333333; - wt[3] = 0.166666666666667; - break; + case 4: + wt[0] = 1.66666666666667e-01; + wt[1] = 8.33333333333333e-01; + wt[2] = 8.33333333333334e-01; + wt[3] = 1.66666666666666e-01; + break; // sum = 2.00000000000000e+00 case 5: - wt[0] = 0.1; - wt[1] = 0.5444444444; - wt[2] = 0.7111111111; - wt[3] = 0.5444444444; - wt[4] = 0.1; - break; + wt[0] = 1.00000000000000e-01; + wt[1] = 5.44444444444445e-01; + wt[2] = 7.11111111111111e-01; + wt[3] = 5.44444444444444e-01; + wt[4] = 1.00000000000000e-01; + break; // sum = 2.00000000000000e+00 - case 6: - wt[0] = 0.06666666667; - wt[1] = 0.3784749562; - wt[2] = 0.5548583770; - wt[3] = 0.5548583770; - wt[4] = 0.3784749562; - wt[5] = 0.06666666667; - break; + case 6: + wt[0] = 6.66666666666666e-02; + wt[1] = 3.78474956297848e-01; + wt[2] = 5.54858377035486e-01; + wt[3] = 5.54858377035487e-01; + wt[4] = 3.78474956297848e-01; + wt[5] = 6.66666666666668e-02; + break; // sum = 2.00000000000000e+00 + + case 7: + wt[0] = 4.76190476190473e-02; + wt[1] = 2.76826047361566e-01; + wt[2] = 4.31745381209863e-01; + wt[3] = 4.87619047619048e-01; + wt[4] = 4.31745381209863e-01; + wt[5] = 2.76826047361567e-01; + wt[6] = 4.76190476190480e-02; + break; // sum = 2.00000000000000e+00 + + case 8: + wt[0] = 3.57142857142856e-02; + wt[1] = 2.10704227143506e-01; + wt[2] = 3.41122692483505e-01; + wt[3] = 4.12458794658704e-01; + wt[4] = 4.12458794658703e-01; + wt[5] = 3.41122692483505e-01; + wt[6] = 2.10704227143506e-01; + wt[7] = 3.57142857142862e-02; + break; // sum = 2.00000000000000e+00 + + case 9: + wt[0] = 2.77777777777778e-02; + wt[1] = 1.65495361560806e-01; + wt[2] = 2.74538712500162e-01; + wt[3] = 3.46428510973046e-01; + wt[4] = 3.71519274376417e-01; + wt[5] = 3.46428510973047e-01; + wt[6] = 2.74538712500161e-01; + wt[7] = 1.65495361560806e-01; + wt[8] = 2.77777777777781e-02; + break; // sum = 2.00000000000000e+00 + + case 10: + wt[0] = 2.22222222222225e-02; + wt[1] = 1.33305990851070e-01; + wt[2] = 2.24889342063126e-01; + wt[3] = 2.92042683679684e-01; + wt[4] = 3.27539761183898e-01; + wt[5] = 3.27539761183897e-01; + wt[6] = 2.92042683679685e-01; + wt[7] = 2.24889342063126e-01; + wt[8] = 1.33305990851070e-01; + wt[9] = 2.22222222222225e-02; + break; // sum = 2.00000000000000e+00 - case 7: - wt[0] = 0.04761904762; - wt[1] = 0.2768260473; - wt[2] = 0.4317453812; - wt[3] = 0.4876190476; - wt[4] = 0.4317453812; - wt[5] = 0.2768260473; - wt[6] = 0.04761904762; - break; - - case 8: - wt[0] = 0.03571428571; - wt[1] = 0.2107042271; - wt[2] = 0.3411226924; - wt[3] = 0.4124587946; - wt[4] = 0.4124587946; - wt[5] = 0.3411226924; - wt[6] = 0.2107042271; - wt[7] = 0.03571428571; - break; - - case 9: - wt[0] = 0.02777777778; - wt[1] = 0.1654953615; - wt[2] = 0.2745387125; - wt[3] = 0.3464285109; - wt[4] = 0.3715192743; - wt[5] = 0.3464285109; - wt[6] = 0.2745387125; - wt[7] = 0.1654953615; - wt[8] = 0.02777777778; - break; - - case 10: - wt[0] = 0.02222222222; - wt[1] = 0.1333059908; - wt[2] = 0.2248893421; - wt[3] = 0.2920426836; - wt[4] = 0.3275397611; - wt[5] = 0.3275397611; - wt[6] = 0.2920426836; - wt[7] = 0.2248893421; - wt[8] = 0.1333059908; - wt[9] = 0.02222222222; - break; - default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; break; diff --git a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp index 793fb7ccbd..d8c81a4caa 100755 --- a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp @@ -131,87 +131,87 @@ RadauBeamIntegration::getSectionLocations(int numSections, double L, case 1: xi[0] = -1.0; break; - + case 2: - xi[0] = -1.0; - xi[1] = 0.3333333333; + xi[0] = -1.00000000000000e+00; + xi[1] = 3.33333333333333e-01; break; case 3: - xi[0] = -1.0; - xi[1] = -0.2898979485; - xi[2] = 0.6898979485; + xi[0] = -1.00000000000000e+00; + xi[1] = -2.89897948556636e-01; + xi[2] = 6.89897948556636e-01; break; case 4: - xi[0] = -1.0; - xi[1] = -0.5753189235; - xi[2] = 0.1810662711; - xi[3] = 0.8228240809; + xi[0] = -1.00000000000000e+00; + xi[1] = -5.75318923521694e-01; + xi[2] = 1.81066271118531e-01; + xi[3] = 8.22824080974592e-01; break; case 5: - xi[0] = -1.0; - xi[1] = -0.7204802713; - xi[2] = -0.1671808647; - xi[3] = 0.4463139727; - xi[4] = 0.8857916077; - break; - + xi[0] = -1.00000000000000e+00; + xi[1] = -7.20480271312439e-01; + xi[2] = -1.67180864737834e-01; + xi[3] = 4.46313972723752e-01; + xi[4] = 8.85791607770964e-01; + break; + case 6: - xi[0] = -1.0; - xi[1] = -0.8029298284; - xi[2] = -0.3909285467; - xi[3] = 0.1240503795; - xi[4] = 0.6039731642; - xi[5] = 0.9203802858; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.02929828402347e-01; + xi[2] = -3.90928546707272e-01; + xi[3] = 1.24050379505228e-01; + xi[4] = 6.03973164252783e-01; + xi[5] = 9.20380285897063e-01; break; - + case 7: - xi[0] = -1.0; - xi[1] = -0.8538913426; - xi[2] = -0.5384677240; - xi[3] = -0.1173430375; - xi[4] = 0.3260306194; - xi[5] = 0.7038428006; - xi[6] = 0.9413671456; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.53891342639482e-01; + xi[2] = -5.38467724060109e-01; + xi[3] = -1.17343037543100e-01; + xi[4] = 3.26030619437692e-01; + xi[5] = 7.03842800663031e-01; + xi[6] = 9.41367145680430e-01; break; - + case 8: - xi[0] = -1.0; - xi[1] = -0.8874748789; - xi[2] = -0.6395186165; - xi[3] = -0.2947505657; - xi[4] = 0.09430725266; - xi[5] = 0.4684203544; - xi[6] = 0.7706418936; - xi[7] = 0.9550412271; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.87474878926156e-01; + xi[2] = -6.39518616526215e-01; + xi[3] = -2.94750565773661e-01; + xi[4] = 9.43072526611107e-02; + xi[5] = 4.68420354430821e-01; + xi[6] = 7.70641893678192e-01; + xi[7] = 9.55041227122575e-01; break; - + case 9: - xi[0] = -1.0; - xi[1] = -0.9107320894; - xi[2] = -0.7112674859; - xi[3] = -0.4263504857; - xi[4] = -0.09037336960; - xi[5] = 0.2561356708; - xi[6] = 0.5713830412; - xi[7] = 0.8173527842; - xi[8] = 0.9644401697; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.10732089420060e-01; + xi[2] = -7.11267485915709e-01; + xi[3] = -4.26350485711139e-01; + xi[4] = -9.03733696068532e-02; + xi[5] = 2.56135670833455e-01; + xi[6] = 5.71383041208739e-01; + xi[7] = 8.17352784200412e-01; + xi[8] = 9.64440169705273e-01; break; - + case 10: - xi[0] = -1.0; - xi[1] = -0.9274843742; - xi[2] = -0.7638420424; - xi[3] = -0.5256460303; - xi[4] = -0.2362344693; - xi[5] = 0.07605919783; - xi[6] = 0.3806648401; - xi[7] = 0.6477666876; - xi[8] = 0.8512252205; - xi[9] = 0.9711751807; - break; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.27484374233581e-01; + xi[2] = -7.63842042420003e-01; + xi[3] = -5.25646030370079e-01; + xi[4] = -2.36234469390588e-01; + xi[5] = 7.60591978379783e-02; + xi[6] = 3.80664840144725e-01; + xi[7] = 6.47766687674009e-01; + xi[8] = 8.51225220581608e-01; + xi[9] = 9.71175180702247e-01; + break; default: break; @@ -232,86 +232,86 @@ RadauBeamIntegration::getSectionWeights(int numSections, double L, break; case 2: - wt[0] = 0.5; - wt[1] = 1.5; - break; + wt[0] = 5.00000000000000e-01; + wt[1] = 1.50000000000000e+00; + break; // sum = 2.00000000000000e+00 case 3: - wt[0] = 0.2222222222; - wt[1] = 1.024971652; - wt[2] = 0.7528061254; - break; + wt[0] = 2.22222222222222e-01; + wt[1] = 1.02497165237684e+00; + wt[2] = 7.52806125400935e-01; + break; // sum = 2.00000000000000e+00 case 4: - wt[0] = 0.125; - wt[1] = 0.6576886399; - wt[2] = 0.7763869376; - wt[3] = 0.4409244223; - break; + wt[0] = 1.25000000000000e-01; + wt[1] = 6.57688639960120e-01; + wt[2] = 7.76386937686343e-01; + wt[3] = 4.40924422353537e-01; + break; // sum = 2.00000000000000e+00 case 5: - wt[0] = 0.08; - wt[1] = 0.4462078021; - wt[2] = 0.6236530459; - wt[3] = 0.5627120302; - wt[4] = 0.2874271215; - break; - + wt[0] = 8.00000000000000e-02; + wt[1] = 4.46207802167141e-01; + wt[2] = 6.23653045951482e-01; + wt[3] = 5.62712030298925e-01; + wt[4] = 2.87427121582452e-01; + break; // sum = 2.00000000000000e+00 + case 6: - wt[0] = 0.05555555555; - wt[1] = 0.3196407532; - wt[2] = 0.4853871884; - wt[3] = 0.5209267831; - wt[4] = 0.4169013343; - wt[5] = 0.2015883852; - break; - + wt[0] = 5.55555555555557e-02; + wt[1] = 3.19640753220511e-01; + wt[2] = 4.85387188468970e-01; + wt[3] = 5.20926783189575e-01; + wt[4] = 4.16901334311908e-01; + wt[5] = 2.01588385253481e-01; + break; // sum = 2.00000000000000e+00 + case 7: - wt[0] = 0.04081632653; - wt[1] = 0.2392274892; - wt[2] = 0.3809498736; - wt[3] = 0.4471098290; - wt[4] = 0.4247037790; - wt[5] = 0.3182042314; - wt[6] = 0.1489884711; - break; - + wt[0] = 4.08163265306122e-02; + wt[1] = 2.39227489225312e-01; + wt[2] = 3.80949873644231e-01; + wt[3] = 4.47109829014566e-01; + wt[4] = 4.24703779005955e-01; + wt[5] = 3.18204231467302e-01; + wt[6] = 1.48988471112020e-01; + break; // sum = 2.00000000000000e+00 + case 8: - wt[0] = 0.03125; - wt[1] = 0.1853581548; - wt[2] = 0.3041306206; - wt[3] = 0.3765175453; - wt[4] = 0.3915721674; - wt[5] = 0.3470147956; - wt[6] = 0.2496479013; - wt[7] = 0.1145088147; - break; - + wt[0] = 3.12500000000000e-02; + wt[1] = 1.85358154802980e-01; + wt[2] = 3.04130620646785e-01; + wt[3] = 3.76517545389118e-01; + wt[4] = 3.91572167452494e-01; + wt[5] = 3.47014795634501e-01; + wt[6] = 2.49647901329865e-01; + wt[7] = 1.14508814744258e-01; + break; // sum = 2.00000000000000e+00 + case 9: - wt[0] = 0.02469135802; - wt[1] = 0.1476540190; - wt[2] = 0.2471893782; - wt[3] = 0.3168437756; - wt[4] = 0.3482730027; - wt[5] = 0.3376939669; - wt[6] = 0.2863866963; - wt[7] = 0.2005532980; - wt[8] = 0.09071450492; - break; - + wt[0] = 2.46913580246915e-02; + wt[1] = 1.47654019046315e-01; + wt[2] = 2.47189378204593e-01; + wt[3] = 3.16843775670438e-01; + wt[4] = 3.48273002772966e-01; + wt[5] = 3.37693966975930e-01; + wt[6] = 2.86386696357232e-01; + wt[7] = 2.00553298024552e-01; + wt[8] = 9.07145049232830e-02; + break; // sum = 2.00000000000000e+00 + case 10: - wt[0] = 0.02; - wt[1] = 0.1202966705; - wt[2] = 0.2042701318; - wt[3] = 0.2681948378; - wt[4] = 0.3058592877; - wt[5] = 0.3135824572; - wt[6] = 0.2906101648; - wt[7] = 0.2391934317; - wt[8] = 0.1643760127; - wt[9] = 0.07361700548; - break; - + wt[0] = 2.00000000000001e-02; + wt[1] = 1.20296670557482e-01; + wt[2] = 2.04270131879001e-01; + wt[3] = 2.68194837841179e-01; + wt[4] = 3.05859287724423e-01; + wt[5] = 3.13582457226938e-01; + wt[6] = 2.90610164832919e-01; + wt[7] = 2.39193431714380e-01; + wt[8] = 1.64376012736922e-01; + wt[9] = 7.36170054867581e-02; + break; // sum = 2.00000000000000e+00 + default: break; } From c7ffd83de916c57bebca9fcee542131b8a64dd59 Mon Sep 17 00:00:00 2001 From: "Michael H. Scott" Date: Sat, 16 Nov 2024 10:51:28 -0800 Subject: [PATCH 2/3] Revert "Using 16 digits in locations and weights - generated from gaussq.f" This reverts commit 4c8801e20eb1eb1ef2751c436a7b8788940e1125. --- .../LobattoBeamIntegration.cpp | 262 +++++++++--------- .../forceBeamColumn/RadauBeamIntegration.cpp | 262 +++++++++--------- 2 files changed, 262 insertions(+), 262 deletions(-) diff --git a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp index 365f189e0e..3f28105282 100644 --- a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp @@ -111,87 +111,87 @@ LobattoBeamIntegration::getSectionLocations(int numSections, double *xi) { switch(numSections) { - + case 2: - xi[0] = -1.00000000000000e+00; - xi[1] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = 1.0; break; case 3: - xi[0] = -1.00000000000000e+00; - xi[1] = 0.0; - xi[2] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = 0.0; + xi[2] = 1.0; break; case 4: - xi[0] = -1.00000000000000e+00; - xi[1] = -4.47213595499958e-01; - xi[2] = 4.47213595499958e-01; - xi[3] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.44721360; + xi[2] = 0.44721360; + xi[3] = 1.0; break; case 5: - xi[0] = -1.00000000000000e+00; - xi[1] = -6.54653670707977e-01; - xi[2] = 0.0; - xi[3] = 6.54653670707977e-01; - xi[4] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.65465367; + xi[2] = 0.0; + xi[3] = 0.65465367; + xi[4] = 1.0; break; case 6: - xi[0] = -1.00000000000000e+00; - xi[1] = -7.65055323929465e-01; - xi[2] = -2.85231516480645e-01; - xi[3] = 2.85231516480645e-01; - xi[4] = 7.65055323929465e-01; - xi[5] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.7650553239; + xi[2] = -0.2852315164; + xi[3] = 0.2852315164; + xi[4] = 0.7650553239; + xi[5] = 1.0; break; case 7: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.30223896278567e-01; - xi[2] = -4.68848793470714e-01; - xi[3] = 0.0; - xi[4] = 4.68848793470714e-01; - xi[5] = 8.30223896278567e-01; - xi[6] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.8302238962; + xi[2] = -0.4688487934; + xi[3] = 0.0; + xi[4] = 0.4688487934; + xi[5] = 0.8302238962; + xi[6] = 1.0; break; - + case 8: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.71740148509606e-01; - xi[2] = -5.91700181433143e-01; - xi[3] = -2.09299217902479e-01; - xi[4] = 2.09299217902479e-01; - xi[5] = 5.91700181433143e-01; - xi[6] = 8.71740148509607e-01; - xi[7] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.8717401485; + xi[2] = -0.5917001814; + xi[3] = -0.2092992179; + xi[4] = 0.2092992179; + xi[5] = 0.5917001814; + xi[6] = 0.8717401485; + xi[7] = 1.0; break; case 9: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.99757995411461e-01; - xi[2] = -6.77186279510738e-01; - xi[3] = -3.63117463826178e-01; - xi[4] = 0.0; - xi[5] = 3.63117463826178e-01; - xi[6] = 6.77186279510738e-01; - xi[7] = 8.99757995411461e-01; - xi[8] = 1.00000000000000e+00; + xi[0] = -1.0; + xi[1] = -0.8997579954; + xi[2] = -0.6771862795; + xi[3] = -0.3631174638; + xi[4] = 0.0; + xi[5] = 0.3631174638; + xi[6] = 0.6771862795; + xi[7] = 0.8997579954; + xi[8] = 1.0; break; - + case 10: - xi[0] = -1.00000000000000e+00; - xi[1] = -9.19533908166459e-01; - xi[2] = -7.38773865105505e-01; - xi[3] = -4.77924949810445e-01; - xi[4] = -1.65278957666387e-01; - xi[5] = 1.65278957666387e-01; - xi[6] = 4.77924949810445e-01; - xi[7] = 7.38773865105505e-01; - xi[8] = 9.19533908166459e-01; - xi[9] = 1.00000000000000e+00; - break; + xi[0] = -1.0; + xi[1] = -0.9195339082; + xi[2] = -0.7387738651; + xi[3] = -0.4779249498; + xi[4] = -0.1652789577; + xi[5] = 0.1652789577; + xi[6] = 0.4779249498; + xi[7] = 0.7387738651; + xi[8] = 0.9195339082; + xi[9] = 1.0; + break; default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; @@ -210,86 +210,86 @@ LobattoBeamIntegration::getSectionWeights(int numSections, double L, switch (numSections) { case 2: - wt[0] = 1.00000000000000e+00; - wt[1] = 1.00000000000000e+00; - break; // sum = 2.00000000000000e+00 + wt[0] = 1.0; + wt[1] = 1.0; + break; case 3: - wt[0] = 3.33333333333333e-01; - wt[1] = 1.33333333333333e+00; - wt[2] = 3.33333333333334e-01; - break; // sum = 2.00000000000000e+00 + wt[0] = 0.333333333333333; + wt[1] = 1.333333333333333; + wt[2] = 0.333333333333333; + break; - case 4: - wt[0] = 1.66666666666667e-01; - wt[1] = 8.33333333333333e-01; - wt[2] = 8.33333333333334e-01; - wt[3] = 1.66666666666666e-01; - break; // sum = 2.00000000000000e+00 + case 4: + wt[0] = 0.166666666666667; + wt[1] = 0.833333333333333; + wt[2] = 0.833333333333333; + wt[3] = 0.166666666666667; + break; case 5: - wt[0] = 1.00000000000000e-01; - wt[1] = 5.44444444444445e-01; - wt[2] = 7.11111111111111e-01; - wt[3] = 5.44444444444444e-01; - wt[4] = 1.00000000000000e-01; - break; // sum = 2.00000000000000e+00 - - case 6: - wt[0] = 6.66666666666666e-02; - wt[1] = 3.78474956297848e-01; - wt[2] = 5.54858377035486e-01; - wt[3] = 5.54858377035487e-01; - wt[4] = 3.78474956297848e-01; - wt[5] = 6.66666666666668e-02; - break; // sum = 2.00000000000000e+00 - - case 7: - wt[0] = 4.76190476190473e-02; - wt[1] = 2.76826047361566e-01; - wt[2] = 4.31745381209863e-01; - wt[3] = 4.87619047619048e-01; - wt[4] = 4.31745381209863e-01; - wt[5] = 2.76826047361567e-01; - wt[6] = 4.76190476190480e-02; - break; // sum = 2.00000000000000e+00 - - case 8: - wt[0] = 3.57142857142856e-02; - wt[1] = 2.10704227143506e-01; - wt[2] = 3.41122692483505e-01; - wt[3] = 4.12458794658704e-01; - wt[4] = 4.12458794658703e-01; - wt[5] = 3.41122692483505e-01; - wt[6] = 2.10704227143506e-01; - wt[7] = 3.57142857142862e-02; - break; // sum = 2.00000000000000e+00 - - case 9: - wt[0] = 2.77777777777778e-02; - wt[1] = 1.65495361560806e-01; - wt[2] = 2.74538712500162e-01; - wt[3] = 3.46428510973046e-01; - wt[4] = 3.71519274376417e-01; - wt[5] = 3.46428510973047e-01; - wt[6] = 2.74538712500161e-01; - wt[7] = 1.65495361560806e-01; - wt[8] = 2.77777777777781e-02; - break; // sum = 2.00000000000000e+00 + wt[0] = 0.1; + wt[1] = 0.5444444444; + wt[2] = 0.7111111111; + wt[3] = 0.5444444444; + wt[4] = 0.1; + break; - case 10: - wt[0] = 2.22222222222225e-02; - wt[1] = 1.33305990851070e-01; - wt[2] = 2.24889342063126e-01; - wt[3] = 2.92042683679684e-01; - wt[4] = 3.27539761183898e-01; - wt[5] = 3.27539761183897e-01; - wt[6] = 2.92042683679685e-01; - wt[7] = 2.24889342063126e-01; - wt[8] = 1.33305990851070e-01; - wt[9] = 2.22222222222225e-02; - break; // sum = 2.00000000000000e+00 + case 6: + wt[0] = 0.06666666667; + wt[1] = 0.3784749562; + wt[2] = 0.5548583770; + wt[3] = 0.5548583770; + wt[4] = 0.3784749562; + wt[5] = 0.06666666667; + break; + case 7: + wt[0] = 0.04761904762; + wt[1] = 0.2768260473; + wt[2] = 0.4317453812; + wt[3] = 0.4876190476; + wt[4] = 0.4317453812; + wt[5] = 0.2768260473; + wt[6] = 0.04761904762; + break; + + case 8: + wt[0] = 0.03571428571; + wt[1] = 0.2107042271; + wt[2] = 0.3411226924; + wt[3] = 0.4124587946; + wt[4] = 0.4124587946; + wt[5] = 0.3411226924; + wt[6] = 0.2107042271; + wt[7] = 0.03571428571; + break; + + case 9: + wt[0] = 0.02777777778; + wt[1] = 0.1654953615; + wt[2] = 0.2745387125; + wt[3] = 0.3464285109; + wt[4] = 0.3715192743; + wt[5] = 0.3464285109; + wt[6] = 0.2745387125; + wt[7] = 0.1654953615; + wt[8] = 0.02777777778; + break; + + case 10: + wt[0] = 0.02222222222; + wt[1] = 0.1333059908; + wt[2] = 0.2248893421; + wt[3] = 0.2920426836; + wt[4] = 0.3275397611; + wt[5] = 0.3275397611; + wt[6] = 0.2920426836; + wt[7] = 0.2248893421; + wt[8] = 0.1333059908; + wt[9] = 0.02222222222; + break; + default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; break; diff --git a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp index d8c81a4caa..793fb7ccbd 100755 --- a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp @@ -131,87 +131,87 @@ RadauBeamIntegration::getSectionLocations(int numSections, double L, case 1: xi[0] = -1.0; break; - + case 2: - xi[0] = -1.00000000000000e+00; - xi[1] = 3.33333333333333e-01; + xi[0] = -1.0; + xi[1] = 0.3333333333; break; case 3: - xi[0] = -1.00000000000000e+00; - xi[1] = -2.89897948556636e-01; - xi[2] = 6.89897948556636e-01; + xi[0] = -1.0; + xi[1] = -0.2898979485; + xi[2] = 0.6898979485; break; case 4: - xi[0] = -1.00000000000000e+00; - xi[1] = -5.75318923521694e-01; - xi[2] = 1.81066271118531e-01; - xi[3] = 8.22824080974592e-01; + xi[0] = -1.0; + xi[1] = -0.5753189235; + xi[2] = 0.1810662711; + xi[3] = 0.8228240809; break; case 5: - xi[0] = -1.00000000000000e+00; - xi[1] = -7.20480271312439e-01; - xi[2] = -1.67180864737834e-01; - xi[3] = 4.46313972723752e-01; - xi[4] = 8.85791607770964e-01; - break; - + xi[0] = -1.0; + xi[1] = -0.7204802713; + xi[2] = -0.1671808647; + xi[3] = 0.4463139727; + xi[4] = 0.8857916077; + break; + case 6: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.02929828402347e-01; - xi[2] = -3.90928546707272e-01; - xi[3] = 1.24050379505228e-01; - xi[4] = 6.03973164252783e-01; - xi[5] = 9.20380285897063e-01; + xi[0] = -1.0; + xi[1] = -0.8029298284; + xi[2] = -0.3909285467; + xi[3] = 0.1240503795; + xi[4] = 0.6039731642; + xi[5] = 0.9203802858; break; - + case 7: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.53891342639482e-01; - xi[2] = -5.38467724060109e-01; - xi[3] = -1.17343037543100e-01; - xi[4] = 3.26030619437692e-01; - xi[5] = 7.03842800663031e-01; - xi[6] = 9.41367145680430e-01; + xi[0] = -1.0; + xi[1] = -0.8538913426; + xi[2] = -0.5384677240; + xi[3] = -0.1173430375; + xi[4] = 0.3260306194; + xi[5] = 0.7038428006; + xi[6] = 0.9413671456; break; - + case 8: - xi[0] = -1.00000000000000e+00; - xi[1] = -8.87474878926156e-01; - xi[2] = -6.39518616526215e-01; - xi[3] = -2.94750565773661e-01; - xi[4] = 9.43072526611107e-02; - xi[5] = 4.68420354430821e-01; - xi[6] = 7.70641893678192e-01; - xi[7] = 9.55041227122575e-01; + xi[0] = -1.0; + xi[1] = -0.8874748789; + xi[2] = -0.6395186165; + xi[3] = -0.2947505657; + xi[4] = 0.09430725266; + xi[5] = 0.4684203544; + xi[6] = 0.7706418936; + xi[7] = 0.9550412271; break; - + case 9: - xi[0] = -1.00000000000000e+00; - xi[1] = -9.10732089420060e-01; - xi[2] = -7.11267485915709e-01; - xi[3] = -4.26350485711139e-01; - xi[4] = -9.03733696068532e-02; - xi[5] = 2.56135670833455e-01; - xi[6] = 5.71383041208739e-01; - xi[7] = 8.17352784200412e-01; - xi[8] = 9.64440169705273e-01; + xi[0] = -1.0; + xi[1] = -0.9107320894; + xi[2] = -0.7112674859; + xi[3] = -0.4263504857; + xi[4] = -0.09037336960; + xi[5] = 0.2561356708; + xi[6] = 0.5713830412; + xi[7] = 0.8173527842; + xi[8] = 0.9644401697; break; - + case 10: - xi[0] = -1.00000000000000e+00; - xi[1] = -9.27484374233581e-01; - xi[2] = -7.63842042420003e-01; - xi[3] = -5.25646030370079e-01; - xi[4] = -2.36234469390588e-01; - xi[5] = 7.60591978379783e-02; - xi[6] = 3.80664840144725e-01; - xi[7] = 6.47766687674009e-01; - xi[8] = 8.51225220581608e-01; - xi[9] = 9.71175180702247e-01; - break; + xi[0] = -1.0; + xi[1] = -0.9274843742; + xi[2] = -0.7638420424; + xi[3] = -0.5256460303; + xi[4] = -0.2362344693; + xi[5] = 0.07605919783; + xi[6] = 0.3806648401; + xi[7] = 0.6477666876; + xi[8] = 0.8512252205; + xi[9] = 0.9711751807; + break; default: break; @@ -232,86 +232,86 @@ RadauBeamIntegration::getSectionWeights(int numSections, double L, break; case 2: - wt[0] = 5.00000000000000e-01; - wt[1] = 1.50000000000000e+00; - break; // sum = 2.00000000000000e+00 + wt[0] = 0.5; + wt[1] = 1.5; + break; case 3: - wt[0] = 2.22222222222222e-01; - wt[1] = 1.02497165237684e+00; - wt[2] = 7.52806125400935e-01; - break; // sum = 2.00000000000000e+00 + wt[0] = 0.2222222222; + wt[1] = 1.024971652; + wt[2] = 0.7528061254; + break; case 4: - wt[0] = 1.25000000000000e-01; - wt[1] = 6.57688639960120e-01; - wt[2] = 7.76386937686343e-01; - wt[3] = 4.40924422353537e-01; - break; // sum = 2.00000000000000e+00 + wt[0] = 0.125; + wt[1] = 0.6576886399; + wt[2] = 0.7763869376; + wt[3] = 0.4409244223; + break; case 5: - wt[0] = 8.00000000000000e-02; - wt[1] = 4.46207802167141e-01; - wt[2] = 6.23653045951482e-01; - wt[3] = 5.62712030298925e-01; - wt[4] = 2.87427121582452e-01; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.08; + wt[1] = 0.4462078021; + wt[2] = 0.6236530459; + wt[3] = 0.5627120302; + wt[4] = 0.2874271215; + break; + case 6: - wt[0] = 5.55555555555557e-02; - wt[1] = 3.19640753220511e-01; - wt[2] = 4.85387188468970e-01; - wt[3] = 5.20926783189575e-01; - wt[4] = 4.16901334311908e-01; - wt[5] = 2.01588385253481e-01; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.05555555555; + wt[1] = 0.3196407532; + wt[2] = 0.4853871884; + wt[3] = 0.5209267831; + wt[4] = 0.4169013343; + wt[5] = 0.2015883852; + break; + case 7: - wt[0] = 4.08163265306122e-02; - wt[1] = 2.39227489225312e-01; - wt[2] = 3.80949873644231e-01; - wt[3] = 4.47109829014566e-01; - wt[4] = 4.24703779005955e-01; - wt[5] = 3.18204231467302e-01; - wt[6] = 1.48988471112020e-01; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.04081632653; + wt[1] = 0.2392274892; + wt[2] = 0.3809498736; + wt[3] = 0.4471098290; + wt[4] = 0.4247037790; + wt[5] = 0.3182042314; + wt[6] = 0.1489884711; + break; + case 8: - wt[0] = 3.12500000000000e-02; - wt[1] = 1.85358154802980e-01; - wt[2] = 3.04130620646785e-01; - wt[3] = 3.76517545389118e-01; - wt[4] = 3.91572167452494e-01; - wt[5] = 3.47014795634501e-01; - wt[6] = 2.49647901329865e-01; - wt[7] = 1.14508814744258e-01; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.03125; + wt[1] = 0.1853581548; + wt[2] = 0.3041306206; + wt[3] = 0.3765175453; + wt[4] = 0.3915721674; + wt[5] = 0.3470147956; + wt[6] = 0.2496479013; + wt[7] = 0.1145088147; + break; + case 9: - wt[0] = 2.46913580246915e-02; - wt[1] = 1.47654019046315e-01; - wt[2] = 2.47189378204593e-01; - wt[3] = 3.16843775670438e-01; - wt[4] = 3.48273002772966e-01; - wt[5] = 3.37693966975930e-01; - wt[6] = 2.86386696357232e-01; - wt[7] = 2.00553298024552e-01; - wt[8] = 9.07145049232830e-02; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.02469135802; + wt[1] = 0.1476540190; + wt[2] = 0.2471893782; + wt[3] = 0.3168437756; + wt[4] = 0.3482730027; + wt[5] = 0.3376939669; + wt[6] = 0.2863866963; + wt[7] = 0.2005532980; + wt[8] = 0.09071450492; + break; + case 10: - wt[0] = 2.00000000000001e-02; - wt[1] = 1.20296670557482e-01; - wt[2] = 2.04270131879001e-01; - wt[3] = 2.68194837841179e-01; - wt[4] = 3.05859287724423e-01; - wt[5] = 3.13582457226938e-01; - wt[6] = 2.90610164832919e-01; - wt[7] = 2.39193431714380e-01; - wt[8] = 1.64376012736922e-01; - wt[9] = 7.36170054867581e-02; - break; // sum = 2.00000000000000e+00 - + wt[0] = 0.02; + wt[1] = 0.1202966705; + wt[2] = 0.2042701318; + wt[3] = 0.2681948378; + wt[4] = 0.3058592877; + wt[5] = 0.3135824572; + wt[6] = 0.2906101648; + wt[7] = 0.2391934317; + wt[8] = 0.1643760127; + wt[9] = 0.07361700548; + break; + default: break; } From 8f92d32e102ea6c77d2e3b134ce31da0c56128c9 Mon Sep 17 00:00:00 2001 From: "Michael H. Scott" Date: Sun, 17 Nov 2024 15:08:03 -0800 Subject: [PATCH 3/3] Revert "Revert "Using 16 digits in locations and weights - generated from gaussq.f"" This reverts commit c7ffd83de916c57bebca9fcee542131b8a64dd59. --- .../LobattoBeamIntegration.cpp | 262 +++++++++--------- .../forceBeamColumn/RadauBeamIntegration.cpp | 262 +++++++++--------- 2 files changed, 262 insertions(+), 262 deletions(-) diff --git a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp index 3f28105282..365f189e0e 100644 --- a/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/LobattoBeamIntegration.cpp @@ -111,87 +111,87 @@ LobattoBeamIntegration::getSectionLocations(int numSections, double *xi) { switch(numSections) { - + case 2: - xi[0] = -1.0; - xi[1] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = 1.00000000000000e+00; break; case 3: - xi[0] = -1.0; - xi[1] = 0.0; - xi[2] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = 0.0; + xi[2] = 1.00000000000000e+00; break; case 4: - xi[0] = -1.0; - xi[1] = -0.44721360; - xi[2] = 0.44721360; - xi[3] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -4.47213595499958e-01; + xi[2] = 4.47213595499958e-01; + xi[3] = 1.00000000000000e+00; break; case 5: - xi[0] = -1.0; - xi[1] = -0.65465367; - xi[2] = 0.0; - xi[3] = 0.65465367; - xi[4] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -6.54653670707977e-01; + xi[2] = 0.0; + xi[3] = 6.54653670707977e-01; + xi[4] = 1.00000000000000e+00; break; case 6: - xi[0] = -1.0; - xi[1] = -0.7650553239; - xi[2] = -0.2852315164; - xi[3] = 0.2852315164; - xi[4] = 0.7650553239; - xi[5] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -7.65055323929465e-01; + xi[2] = -2.85231516480645e-01; + xi[3] = 2.85231516480645e-01; + xi[4] = 7.65055323929465e-01; + xi[5] = 1.00000000000000e+00; break; case 7: - xi[0] = -1.0; - xi[1] = -0.8302238962; - xi[2] = -0.4688487934; - xi[3] = 0.0; - xi[4] = 0.4688487934; - xi[5] = 0.8302238962; - xi[6] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.30223896278567e-01; + xi[2] = -4.68848793470714e-01; + xi[3] = 0.0; + xi[4] = 4.68848793470714e-01; + xi[5] = 8.30223896278567e-01; + xi[6] = 1.00000000000000e+00; break; - + case 8: - xi[0] = -1.0; - xi[1] = -0.8717401485; - xi[2] = -0.5917001814; - xi[3] = -0.2092992179; - xi[4] = 0.2092992179; - xi[5] = 0.5917001814; - xi[6] = 0.8717401485; - xi[7] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.71740148509606e-01; + xi[2] = -5.91700181433143e-01; + xi[3] = -2.09299217902479e-01; + xi[4] = 2.09299217902479e-01; + xi[5] = 5.91700181433143e-01; + xi[6] = 8.71740148509607e-01; + xi[7] = 1.00000000000000e+00; break; case 9: - xi[0] = -1.0; - xi[1] = -0.8997579954; - xi[2] = -0.6771862795; - xi[3] = -0.3631174638; - xi[4] = 0.0; - xi[5] = 0.3631174638; - xi[6] = 0.6771862795; - xi[7] = 0.8997579954; - xi[8] = 1.0; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.99757995411461e-01; + xi[2] = -6.77186279510738e-01; + xi[3] = -3.63117463826178e-01; + xi[4] = 0.0; + xi[5] = 3.63117463826178e-01; + xi[6] = 6.77186279510738e-01; + xi[7] = 8.99757995411461e-01; + xi[8] = 1.00000000000000e+00; break; - + case 10: - xi[0] = -1.0; - xi[1] = -0.9195339082; - xi[2] = -0.7387738651; - xi[3] = -0.4779249498; - xi[4] = -0.1652789577; - xi[5] = 0.1652789577; - xi[6] = 0.4779249498; - xi[7] = 0.7387738651; - xi[8] = 0.9195339082; - xi[9] = 1.0; - break; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.19533908166459e-01; + xi[2] = -7.38773865105505e-01; + xi[3] = -4.77924949810445e-01; + xi[4] = -1.65278957666387e-01; + xi[5] = 1.65278957666387e-01; + xi[6] = 4.77924949810445e-01; + xi[7] = 7.38773865105505e-01; + xi[8] = 9.19533908166459e-01; + xi[9] = 1.00000000000000e+00; + break; default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; @@ -210,86 +210,86 @@ LobattoBeamIntegration::getSectionWeights(int numSections, double L, switch (numSections) { case 2: - wt[0] = 1.0; - wt[1] = 1.0; - break; + wt[0] = 1.00000000000000e+00; + wt[1] = 1.00000000000000e+00; + break; // sum = 2.00000000000000e+00 case 3: - wt[0] = 0.333333333333333; - wt[1] = 1.333333333333333; - wt[2] = 0.333333333333333; - break; + wt[0] = 3.33333333333333e-01; + wt[1] = 1.33333333333333e+00; + wt[2] = 3.33333333333334e-01; + break; // sum = 2.00000000000000e+00 - case 4: - wt[0] = 0.166666666666667; - wt[1] = 0.833333333333333; - wt[2] = 0.833333333333333; - wt[3] = 0.166666666666667; - break; + case 4: + wt[0] = 1.66666666666667e-01; + wt[1] = 8.33333333333333e-01; + wt[2] = 8.33333333333334e-01; + wt[3] = 1.66666666666666e-01; + break; // sum = 2.00000000000000e+00 case 5: - wt[0] = 0.1; - wt[1] = 0.5444444444; - wt[2] = 0.7111111111; - wt[3] = 0.5444444444; - wt[4] = 0.1; - break; + wt[0] = 1.00000000000000e-01; + wt[1] = 5.44444444444445e-01; + wt[2] = 7.11111111111111e-01; + wt[3] = 5.44444444444444e-01; + wt[4] = 1.00000000000000e-01; + break; // sum = 2.00000000000000e+00 - case 6: - wt[0] = 0.06666666667; - wt[1] = 0.3784749562; - wt[2] = 0.5548583770; - wt[3] = 0.5548583770; - wt[4] = 0.3784749562; - wt[5] = 0.06666666667; - break; + case 6: + wt[0] = 6.66666666666666e-02; + wt[1] = 3.78474956297848e-01; + wt[2] = 5.54858377035486e-01; + wt[3] = 5.54858377035487e-01; + wt[4] = 3.78474956297848e-01; + wt[5] = 6.66666666666668e-02; + break; // sum = 2.00000000000000e+00 + + case 7: + wt[0] = 4.76190476190473e-02; + wt[1] = 2.76826047361566e-01; + wt[2] = 4.31745381209863e-01; + wt[3] = 4.87619047619048e-01; + wt[4] = 4.31745381209863e-01; + wt[5] = 2.76826047361567e-01; + wt[6] = 4.76190476190480e-02; + break; // sum = 2.00000000000000e+00 + + case 8: + wt[0] = 3.57142857142856e-02; + wt[1] = 2.10704227143506e-01; + wt[2] = 3.41122692483505e-01; + wt[3] = 4.12458794658704e-01; + wt[4] = 4.12458794658703e-01; + wt[5] = 3.41122692483505e-01; + wt[6] = 2.10704227143506e-01; + wt[7] = 3.57142857142862e-02; + break; // sum = 2.00000000000000e+00 + + case 9: + wt[0] = 2.77777777777778e-02; + wt[1] = 1.65495361560806e-01; + wt[2] = 2.74538712500162e-01; + wt[3] = 3.46428510973046e-01; + wt[4] = 3.71519274376417e-01; + wt[5] = 3.46428510973047e-01; + wt[6] = 2.74538712500161e-01; + wt[7] = 1.65495361560806e-01; + wt[8] = 2.77777777777781e-02; + break; // sum = 2.00000000000000e+00 + + case 10: + wt[0] = 2.22222222222225e-02; + wt[1] = 1.33305990851070e-01; + wt[2] = 2.24889342063126e-01; + wt[3] = 2.92042683679684e-01; + wt[4] = 3.27539761183898e-01; + wt[5] = 3.27539761183897e-01; + wt[6] = 2.92042683679685e-01; + wt[7] = 2.24889342063126e-01; + wt[8] = 1.33305990851070e-01; + wt[9] = 2.22222222222225e-02; + break; // sum = 2.00000000000000e+00 - case 7: - wt[0] = 0.04761904762; - wt[1] = 0.2768260473; - wt[2] = 0.4317453812; - wt[3] = 0.4876190476; - wt[4] = 0.4317453812; - wt[5] = 0.2768260473; - wt[6] = 0.04761904762; - break; - - case 8: - wt[0] = 0.03571428571; - wt[1] = 0.2107042271; - wt[2] = 0.3411226924; - wt[3] = 0.4124587946; - wt[4] = 0.4124587946; - wt[5] = 0.3411226924; - wt[6] = 0.2107042271; - wt[7] = 0.03571428571; - break; - - case 9: - wt[0] = 0.02777777778; - wt[1] = 0.1654953615; - wt[2] = 0.2745387125; - wt[3] = 0.3464285109; - wt[4] = 0.3715192743; - wt[5] = 0.3464285109; - wt[6] = 0.2745387125; - wt[7] = 0.1654953615; - wt[8] = 0.02777777778; - break; - - case 10: - wt[0] = 0.02222222222; - wt[1] = 0.1333059908; - wt[2] = 0.2248893421; - wt[3] = 0.2920426836; - wt[4] = 0.3275397611; - wt[5] = 0.3275397611; - wt[6] = 0.2920426836; - wt[7] = 0.2248893421; - wt[8] = 0.1333059908; - wt[9] = 0.02222222222; - break; - default: opserr << "LobattoBeamIntegration -- max # integration points is 10\n"; break; diff --git a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp index 793fb7ccbd..d8c81a4caa 100755 --- a/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp +++ b/SRC/element/forceBeamColumn/RadauBeamIntegration.cpp @@ -131,87 +131,87 @@ RadauBeamIntegration::getSectionLocations(int numSections, double L, case 1: xi[0] = -1.0; break; - + case 2: - xi[0] = -1.0; - xi[1] = 0.3333333333; + xi[0] = -1.00000000000000e+00; + xi[1] = 3.33333333333333e-01; break; case 3: - xi[0] = -1.0; - xi[1] = -0.2898979485; - xi[2] = 0.6898979485; + xi[0] = -1.00000000000000e+00; + xi[1] = -2.89897948556636e-01; + xi[2] = 6.89897948556636e-01; break; case 4: - xi[0] = -1.0; - xi[1] = -0.5753189235; - xi[2] = 0.1810662711; - xi[3] = 0.8228240809; + xi[0] = -1.00000000000000e+00; + xi[1] = -5.75318923521694e-01; + xi[2] = 1.81066271118531e-01; + xi[3] = 8.22824080974592e-01; break; case 5: - xi[0] = -1.0; - xi[1] = -0.7204802713; - xi[2] = -0.1671808647; - xi[3] = 0.4463139727; - xi[4] = 0.8857916077; - break; - + xi[0] = -1.00000000000000e+00; + xi[1] = -7.20480271312439e-01; + xi[2] = -1.67180864737834e-01; + xi[3] = 4.46313972723752e-01; + xi[4] = 8.85791607770964e-01; + break; + case 6: - xi[0] = -1.0; - xi[1] = -0.8029298284; - xi[2] = -0.3909285467; - xi[3] = 0.1240503795; - xi[4] = 0.6039731642; - xi[5] = 0.9203802858; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.02929828402347e-01; + xi[2] = -3.90928546707272e-01; + xi[3] = 1.24050379505228e-01; + xi[4] = 6.03973164252783e-01; + xi[5] = 9.20380285897063e-01; break; - + case 7: - xi[0] = -1.0; - xi[1] = -0.8538913426; - xi[2] = -0.5384677240; - xi[3] = -0.1173430375; - xi[4] = 0.3260306194; - xi[5] = 0.7038428006; - xi[6] = 0.9413671456; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.53891342639482e-01; + xi[2] = -5.38467724060109e-01; + xi[3] = -1.17343037543100e-01; + xi[4] = 3.26030619437692e-01; + xi[5] = 7.03842800663031e-01; + xi[6] = 9.41367145680430e-01; break; - + case 8: - xi[0] = -1.0; - xi[1] = -0.8874748789; - xi[2] = -0.6395186165; - xi[3] = -0.2947505657; - xi[4] = 0.09430725266; - xi[5] = 0.4684203544; - xi[6] = 0.7706418936; - xi[7] = 0.9550412271; + xi[0] = -1.00000000000000e+00; + xi[1] = -8.87474878926156e-01; + xi[2] = -6.39518616526215e-01; + xi[3] = -2.94750565773661e-01; + xi[4] = 9.43072526611107e-02; + xi[5] = 4.68420354430821e-01; + xi[6] = 7.70641893678192e-01; + xi[7] = 9.55041227122575e-01; break; - + case 9: - xi[0] = -1.0; - xi[1] = -0.9107320894; - xi[2] = -0.7112674859; - xi[3] = -0.4263504857; - xi[4] = -0.09037336960; - xi[5] = 0.2561356708; - xi[6] = 0.5713830412; - xi[7] = 0.8173527842; - xi[8] = 0.9644401697; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.10732089420060e-01; + xi[2] = -7.11267485915709e-01; + xi[3] = -4.26350485711139e-01; + xi[4] = -9.03733696068532e-02; + xi[5] = 2.56135670833455e-01; + xi[6] = 5.71383041208739e-01; + xi[7] = 8.17352784200412e-01; + xi[8] = 9.64440169705273e-01; break; - + case 10: - xi[0] = -1.0; - xi[1] = -0.9274843742; - xi[2] = -0.7638420424; - xi[3] = -0.5256460303; - xi[4] = -0.2362344693; - xi[5] = 0.07605919783; - xi[6] = 0.3806648401; - xi[7] = 0.6477666876; - xi[8] = 0.8512252205; - xi[9] = 0.9711751807; - break; + xi[0] = -1.00000000000000e+00; + xi[1] = -9.27484374233581e-01; + xi[2] = -7.63842042420003e-01; + xi[3] = -5.25646030370079e-01; + xi[4] = -2.36234469390588e-01; + xi[5] = 7.60591978379783e-02; + xi[6] = 3.80664840144725e-01; + xi[7] = 6.47766687674009e-01; + xi[8] = 8.51225220581608e-01; + xi[9] = 9.71175180702247e-01; + break; default: break; @@ -232,86 +232,86 @@ RadauBeamIntegration::getSectionWeights(int numSections, double L, break; case 2: - wt[0] = 0.5; - wt[1] = 1.5; - break; + wt[0] = 5.00000000000000e-01; + wt[1] = 1.50000000000000e+00; + break; // sum = 2.00000000000000e+00 case 3: - wt[0] = 0.2222222222; - wt[1] = 1.024971652; - wt[2] = 0.7528061254; - break; + wt[0] = 2.22222222222222e-01; + wt[1] = 1.02497165237684e+00; + wt[2] = 7.52806125400935e-01; + break; // sum = 2.00000000000000e+00 case 4: - wt[0] = 0.125; - wt[1] = 0.6576886399; - wt[2] = 0.7763869376; - wt[3] = 0.4409244223; - break; + wt[0] = 1.25000000000000e-01; + wt[1] = 6.57688639960120e-01; + wt[2] = 7.76386937686343e-01; + wt[3] = 4.40924422353537e-01; + break; // sum = 2.00000000000000e+00 case 5: - wt[0] = 0.08; - wt[1] = 0.4462078021; - wt[2] = 0.6236530459; - wt[3] = 0.5627120302; - wt[4] = 0.2874271215; - break; - + wt[0] = 8.00000000000000e-02; + wt[1] = 4.46207802167141e-01; + wt[2] = 6.23653045951482e-01; + wt[3] = 5.62712030298925e-01; + wt[4] = 2.87427121582452e-01; + break; // sum = 2.00000000000000e+00 + case 6: - wt[0] = 0.05555555555; - wt[1] = 0.3196407532; - wt[2] = 0.4853871884; - wt[3] = 0.5209267831; - wt[4] = 0.4169013343; - wt[5] = 0.2015883852; - break; - + wt[0] = 5.55555555555557e-02; + wt[1] = 3.19640753220511e-01; + wt[2] = 4.85387188468970e-01; + wt[3] = 5.20926783189575e-01; + wt[4] = 4.16901334311908e-01; + wt[5] = 2.01588385253481e-01; + break; // sum = 2.00000000000000e+00 + case 7: - wt[0] = 0.04081632653; - wt[1] = 0.2392274892; - wt[2] = 0.3809498736; - wt[3] = 0.4471098290; - wt[4] = 0.4247037790; - wt[5] = 0.3182042314; - wt[6] = 0.1489884711; - break; - + wt[0] = 4.08163265306122e-02; + wt[1] = 2.39227489225312e-01; + wt[2] = 3.80949873644231e-01; + wt[3] = 4.47109829014566e-01; + wt[4] = 4.24703779005955e-01; + wt[5] = 3.18204231467302e-01; + wt[6] = 1.48988471112020e-01; + break; // sum = 2.00000000000000e+00 + case 8: - wt[0] = 0.03125; - wt[1] = 0.1853581548; - wt[2] = 0.3041306206; - wt[3] = 0.3765175453; - wt[4] = 0.3915721674; - wt[5] = 0.3470147956; - wt[6] = 0.2496479013; - wt[7] = 0.1145088147; - break; - + wt[0] = 3.12500000000000e-02; + wt[1] = 1.85358154802980e-01; + wt[2] = 3.04130620646785e-01; + wt[3] = 3.76517545389118e-01; + wt[4] = 3.91572167452494e-01; + wt[5] = 3.47014795634501e-01; + wt[6] = 2.49647901329865e-01; + wt[7] = 1.14508814744258e-01; + break; // sum = 2.00000000000000e+00 + case 9: - wt[0] = 0.02469135802; - wt[1] = 0.1476540190; - wt[2] = 0.2471893782; - wt[3] = 0.3168437756; - wt[4] = 0.3482730027; - wt[5] = 0.3376939669; - wt[6] = 0.2863866963; - wt[7] = 0.2005532980; - wt[8] = 0.09071450492; - break; - + wt[0] = 2.46913580246915e-02; + wt[1] = 1.47654019046315e-01; + wt[2] = 2.47189378204593e-01; + wt[3] = 3.16843775670438e-01; + wt[4] = 3.48273002772966e-01; + wt[5] = 3.37693966975930e-01; + wt[6] = 2.86386696357232e-01; + wt[7] = 2.00553298024552e-01; + wt[8] = 9.07145049232830e-02; + break; // sum = 2.00000000000000e+00 + case 10: - wt[0] = 0.02; - wt[1] = 0.1202966705; - wt[2] = 0.2042701318; - wt[3] = 0.2681948378; - wt[4] = 0.3058592877; - wt[5] = 0.3135824572; - wt[6] = 0.2906101648; - wt[7] = 0.2391934317; - wt[8] = 0.1643760127; - wt[9] = 0.07361700548; - break; - + wt[0] = 2.00000000000001e-02; + wt[1] = 1.20296670557482e-01; + wt[2] = 2.04270131879001e-01; + wt[3] = 2.68194837841179e-01; + wt[4] = 3.05859287724423e-01; + wt[5] = 3.13582457226938e-01; + wt[6] = 2.90610164832919e-01; + wt[7] = 2.39193431714380e-01; + wt[8] = 1.64376012736922e-01; + wt[9] = 7.36170054867581e-02; + break; // sum = 2.00000000000000e+00 + default: break; }