1C> \ingroup nwxc 2C> @{ 3C> 4C> \file nwxcm_c_perdew86.F 5C> The nwxcm_c_perdew86 functional 6C> 7C> @} 8C> 9C> \ingroup nwxc_priv 10C> @{ 11C> 12C> \brief Evaluate the nwxcm_c_perdew86 functional [1] 13C> 14C> \f{eqnarray*}{ 15C> {\it t_1} &=& \rho_\beta+\rho_\alpha\\\\ 16C> {\it t_2} &=& {{1}\over{{\it t_1}}}\\\\ 17C> {\it t_3} &=& \rho_\alpha-\rho_\beta\\\\ 18C> {\it t_4} &=& {{1}\over{{\it t_1}^{{{2}\over{3}}}}}\\\\ 19C> {\it t_5} &=& \left| {\it t_1}\right| \\\\ 20C> {\it t_6} &=& {{1}\over{{\it t_5}^{{{1}\over{3}}}}}\\\\ 21C> {\it t_7} &=& {{0.01443307452126544\,{\it t_6} 22C> +2.843543831490385 \times 10^{-6}\,{\it t_4}+0.002568} 23C> \over{0.01763993811759022\,{\it t_4}\,{\it t_6} 24C> +5.411317332115466\,{\it t_6}+0.1816419932959077\,{\it t_4} 25C> +1.0}}+0.001667\\\\ 26C> {\it t_8} &=& \sigma_{\beta\beta}+2.0\,\sigma_{\alpha\beta} 27C> +\sigma_{\alpha\alpha}\\\\ 28C> {\it t_9} &=& {{1}\over{\rho_s^{{{2}\over{3}}}}}\\\\ 29C> {\it t_{10}} &=& \left| \rho_s\right| \\\\ 30C> {\it t_{11}} &=& {{1}\over{{\it t_{10}}^{{{1}\over{3}}}}}\\\\ 31C> {\it t_{12}} &=& {{0.01443307452126544\,{\it t_{11}} 32C> +2.843543831490385 \times 10^{-6}\,{\it t_9}+0.002568} 33C> \over{0.01763993811759022\,{\it t_9}\,{\it t_{11}} 34C> +5.411317332115466\,{\it t_{11}}+0.1816419932959077\,{\it t_9} 35C> +1.0}}+0.001667\\\\ 36C> f &=& {{0.7937005259840999\,{\it t_2}\,{\it t_7}\,{\it t_6} 37C> \,{\it t_8}}\over{\sqrt{0.3149802624737183\,\left(1.0\,{ 38C> \it t_3}\,{\it t_2}+1.0\right)^{{{5}\over{3}}} 39C> +0.3149802624737183\,\left(1.0-1.0\,{\it t_3}\,{ 40C> \it t_2}\right)^{{{5}\over{3}}}} 41C> \,e^{{{8.129082500000001 \times 10^{-4}\,\sqrt{{ 42C> \it t_8}}}\over{{\it t_7}\,{\it t_5}^{1.166666666666667}}}}}}\\\\ 43C> g &=& 0\\\\ 44C> G &=& {{0.7937005259840999\,{\it t_{12}}\,{\it t_{11}} 45C> \,\sigma_{ss}}\over{\rho_s\,e^{{{8.129082500000001 \times 10^{ 46C> -4}\,\sqrt{\sigma_{ss}}}\over{{\it t_{12}}\,{ 47C> \it t_{10}}^{1.166666666666667}}}}}}\\\\ 48C> \f} 49C> 50C> Code generated with Maxima 5.34.0 [2] 51C> driven by autoxc [3]. 52C> 53C> ### References ### 54C> 55C> [1] JP Perdew, Phys.Rev. B 33, 8822 (1986) , DOI: 56C> <a href="https://doi.org/10.1103/PhysRevB.33.8822 "> 57C> 10.1103/PhysRevB.33.8822 </a> 58C> 59C> [2] Maxima, a computer algebra system, 60C> <a href="http://maxima.sourceforge.net/"> 61C> http://maxima.sourceforge.net/</a> 62C> 63C> [3] autoxc, revision 27097 2015-05-08 64C> 65 subroutine nwxcm_c_perdew86(param,tol_rho,ipol,nq,wght, 66 +rho,rgamma,fnc,Amat,Cmat) 67c $Id: $ 68#ifdef NWXC_QUAD_PREC 69 implicit real(kind=selected_real_kind(30))(a-h,o-z),integer(i-n) 70 integer, parameter :: rk=selected_real_kind(30) 71#else 72 implicit real(kind=selected_real_kind(15))(a-h,o-z),integer(i-n) 73 integer, parameter :: rk=selected_real_kind(15) 74#endif 75#include "nwxc_param.fh" 76 double precision param(*) !< [Input] Parameters of functional 77 double precision tol_rho !< [Input] The lower limit on the density 78 integer ipol !< [Input] The number of spin channels 79 integer nq !< [Input] The number of points 80 double precision wght !< [Input] The weight of the functional 81 double precision rho(nq,*) !< [Input] The density 82 double precision rgamma(nq,*) !< [Input] The norm of the density 83 !< gradients 84 double precision fnc(nq) !< [Output] The value of the functional 85c 86c Sampling Matrices for the XC Kernel 87c 88 double precision Amat(nq,*) !< [Output] The derivative wrt rho 89 double precision Cmat(nq,*) !< [Output] The derivative wrt rgamma 90 integer iq 91 double precision tmp 92 double precision rhoa,rhob 93 double precision gammaaa,gammaab,gammabb 94 double precision taua,taub 95 double precision nwxcm_heaviside 96 external nwxcm_heaviside 97CDIR$ NOVECTOR 98 do iq = 1, nq 99 if (ipol.eq.1) then 100 rhoa = 0.5d0*rho(iq,R_T) 101 gammaaa = 0.25d0*rgamma(iq,G_TT) 102 if (rhoa.gt.tol_rho) then 103 t1 = 3.968502629920499d-1 104 t2 = 1/rhoa 105 t3 = gammaaa**5.0d-1 106 t4 = 6.299605249474366d-1 107 t5 = 1/rhoa**6.666666666666666d-1 108 t6 = 7.937005259840998d-1 109 t7 = abs(rhoa) 110 t8 = 1/t7**3.333333333333333d-1 111 t9 = 1.443307452126544d-2*t6*t8+2.8435438314903855d-6*t4*t5+ 112 1 2.568d-3 113 t10 = 5.411317332115466d+0*t6*t8+8.819969058795109d-3*t5*t8+ 114 1 1.816419932959077d-1*t4*t5+1.0d+0 115 t11 = 1/t10 116 t12 = t11*t9+1.6670000000000001d-3 117 t13 = 1/t12 118 t14 = 1/t7**1.166666666666667d+0 119 t15 = exp(-7.242189178907065d-4*t13*t14*t3) 120 t16 = 1/t7 121 t17 = 3.149802624737183d-1 122 t18 = 1/rhoa**1.6666666666666669d+0 123 t19 = 1/rhoa**3.333333333333333d-1 124 t20 = t11*(-4.8110248404218126d-3*t1*t16*t19-1.8956958876602 125 1 57d-6*t17*t18)-(-2.939989686265036d-3*t18*t8-1.4699948431 126 2 32518d-3*t16*t2-1.8037724440384886d+0*t1*t16*t19-1.210946 127 3 6219727179d-1*t17*t18)*t9/t10**2 128 t21 = 1/t7**1.5000000000000002d+0 129 fnc(iq) = 4.000000000000001d+0*gammaaa*t1*t12*t2*exp(-7.2421 130 1 89178907064d-4*t13*t14*t3)*t8*wght+fnc(iq) 131 Amat(iq,D1_RA) = 6.299605249474365d-1*t2*t6*t8*(3.1748021039 132 1 363994d+0*gammaaa*t12*t15*(4.224610354362455d-4*rhoa*t13* 133 2 t3/t7**3.166666666666667d+0+7.242189178907065d-4*t14*t20* 134 3 t3/t12**2)*wght+3.1748021039363994d+0*gammaaa*t15*t20*wgh 135 4 t)-1.0d+0*gammaaa*t12*t15*t6*t8*wght/rhoa**2-2.6456684199 136 5 46998d-1*gammaaa*t12*t15*t16*wght/rhoa**1.333333333333333 137 6 3d+0+Amat(iq,D1_RA) 138 Cmat(iq,D1_GAA) = 5.0d-1*t12*t15*t2*t6*t8*wght-1.43703234014 139 1 37232d-4*t15*t2*t21*t3*wght+Cmat(iq,D1_GAA) 140 Cmat(iq,D1_GAB) = 1.0d+0*t12*t15*t2*t6*t8*wght-2.87406468028 141 1 74464d-4*t15*t2*t21*t3*wght+Cmat(iq,D1_GAB) 142 endif ! rhoa.gt.tol_rho 143 else ! ipol.eq.1 144 rhoa = rho(iq,R_A) 145 rhob = rho(iq,R_B) 146 gammaaa = rgamma(iq,G_AA) 147 gammaab = rgamma(iq,G_AB) 148 gammabb = rgamma(iq,G_BB) 149 if (rhoa.gt.tol_rho.and.rhob.gt.tol_rho) then 150 t1 = gammabb+2.0d+0*gammaab+gammaaa 151 t2 = rhob+rhoa 152 t3 = 1/t2 153 t4 = rhoa-rhob 154 t5 = 1.0d+0-1.0d+0*t3*t4 155 t6 = 1.0d+0*t3*t4+1.0d+0 156 t7 = (3.149802624737183d-1*t6**1.6666666666666669d+0+3.14980 157 1 2624737183d-1*t5**1.6666666666666669d+0)**5.0d-1 158 t8 = 1/t7 159 t9 = t1**5.0d-1 160 t10 = 1/t2**6.666666666666666d-1 161 t11 = abs(t2) 162 t12 = 1/t11**3.333333333333333d-1 163 t13 = 1.443307452126544d-2*t12+2.8435438314903855d-6*t10+2.5 164 1 68d-3 165 t14 = 1.7639938117590215d-2*t10*t12+5.411317332115466d+0*t12 166 1 +1.816419932959077d-1*t10+1.0d+0 167 t15 = 1/t14 168 t16 = t13*t15+1.6670000000000001d-3 169 t17 = 1/t16 170 t18 = 1/t11**1.166666666666667d+0 171 t19 = exp(-8.129082500000001d-4*t17*t18*t9) 172 t20 = 1/t11 173 t21 = -2.6456684199469993d-1*t1*t16*t19*t20*t8*wght/t2**1.33 174 1 33333333333333d+0 175 t22 = 1/t2**2 176 t23 = 1.0d+0*t22*t4 177 t24 = -1.0d+0*t3 178 t25 = t5**6.666666666666666d-1 179 t26 = -1.0d+0*t22*t4 180 t27 = 1.0d+0*t3 181 t28 = t6**6.666666666666666d-1 182 t29 = 1/t7**3 183 t30 = -7.937005259840999d-1*t1*t12*t16*t19*t22*t8*wght 184 t31 = 1/t2**1.6666666666666669d+0 185 t32 = 1/t2**3.333333333333333d-1 186 t33 = t15*(-4.8110248404218126d-3*t20*t32-1.895695887660257d 187 1 -6*t31)-t13*(-1.8037724440384886d+0*t20*t32-1.17599587450 188 2 60142d-2*t12*t31-1.2109466219727179d-1*t31-5.879979372530 189 3 072d-3*t20*t3)/t14**2 190 t34 = t12*t3*t8*(7.937005259840999d-1*t1*t16*t19*(8.12908250 191 1 0000001d-4*t18*t33*t9/t16**2+9.483929583333336d-4*t17*t2* 192 2 t9/t11**3.166666666666667d+0)*wght+7.937005259840999d-1*t 193 3 1*t19*t33*wght) 194 t35 = 1/t11**1.5000000000000002d+0 195 t36 = 7.937005259840999d-1*t12*t16*t19*t3*t8*wght-3.22602852 196 1 8009071d-4*t19*t3*t35*t8*t9*wght 197 fnc(iq) = 7.937005259840999d-1*t1*t12*t16*t19*t3*t8*wght+fnc 198 1 (iq) 199 Amat(iq,D1_RA) = -3.9685026299204995d-1*t1*t12*t16*t19*(5.24 200 1 9671041228638d-1*(t27+t26)*t28+5.249671041228638d-1*(t24+ 201 2 t23)*t25)*t29*t3*wght+t34+t30+t21+Amat(iq,D1_RA) 202 Amat(iq,D1_RB) = -3.9685026299204995d-1*t1*t12*t16*t19*(5.24 203 1 9671041228638d-1*(t26+t24)*t28+5.249671041228638d-1*t25*( 204 2 t27+t23))*t29*t3*wght+t34+t30+t21+Amat(iq,D1_RB) 205 Cmat(iq,D1_GAA) = t36+Cmat(iq,D1_GAA) 206 Cmat(iq,D1_GAB) = -6.452057056018142d-4*t19*t3*t35*t8*t9*wgh 207 1 t+1.5874010519681997d+0*t12*t16*t19*t3*t8*wght+Cmat(iq,D1 208 2 _GAB) 209 Cmat(iq,D1_GBB) = t36+Cmat(iq,D1_GBB) 210 elseif (rhoa.gt.tol_rho.and.rhob.le.tol_rho) then 211 t1 = 1/rhoa 212 t2 = gammaaa**5.0d-1 213 t3 = 1/rhoa**6.666666666666666d-1 214 t4 = abs(rhoa) 215 t5 = 1/t4**3.333333333333333d-1 216 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 217 1 -3 218 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 219 1 16419932959077d-1*t3+1.0d+0 220 t8 = 1/t7 221 t9 = t6*t8+1.6670000000000001d-3 222 t10 = 1/t9 223 t11 = 1/t4**1.166666666666667d+0 224 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 225 t13 = 1/t4 226 t14 = 1/rhoa**1.6666666666666669d+0 227 t15 = 1/rhoa**3.333333333333333d-1 228 t16 = (-4.8110248404218126d-3*t13*t15-1.895695887660257d-6*t 229 1 14)*t8-(-1.1759958745060142d-2*t14*t5-1.8037724440384886d 230 2 +0*t13*t15-1.2109466219727179d-1*t14-5.879979372530072d-3 231 3 *t1*t13)*t6/t7**2 232 fnc(iq) = 7.937005259840999d-1*gammaaa*t1*t12*t5*t9*wght+fnc 233 1 (iq) 234 Amat(iq,D1_RA) = t1*t5*(7.937005259840999d-1*gammaaa*t12*(8. 235 1 129082500000001d-4*t11*t16*t2/t9**2+9.483929583333336d-4* 236 2 rhoa*t10*t2/t4**3.166666666666667d+0)*t9*wght+7.937005259 237 3 840999d-1*gammaaa*t12*t16*wght)-7.937005259840999d-1*gamm 238 4 aaa*t12*t5*t9*wght/rhoa**2-2.6456684199469993d-1*gammaaa* 239 5 t12*t13*t9*wght/rhoa**1.3333333333333333d+0+Amat(iq,D1_RA 240 6 ) 241 Cmat(iq,D1_GAA) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 242 1 26028528009071d-4*t1*t12*t2*wght/t4**1.5000000000000002d+ 243 2 0+Cmat(iq,D1_GAA) 244 elseif (rhoa.le.tol_rho.and.rhob.gt.tol_rho) then 245 t1 = 1/rhob 246 t2 = gammabb**5.0d-1 247 t3 = 1/rhob**6.666666666666666d-1 248 t4 = abs(rhob) 249 t5 = 1/t4**3.333333333333333d-1 250 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 251 1 -3 252 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 253 1 16419932959077d-1*t3+1.0d+0 254 t8 = 1/t7 255 t9 = t6*t8+1.6670000000000001d-3 256 t10 = 1/t9 257 t11 = 1/t4**1.166666666666667d+0 258 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 259 t13 = 1/t4 260 t14 = 1/rhob**1.6666666666666669d+0 261 t15 = 1/rhob**3.333333333333333d-1 262 t16 = (-4.8110248404218126d-3*t13*t15-1.895695887660257d-6*t 263 1 14)*t8-(-1.1759958745060142d-2*t14*t5-1.8037724440384886d 264 2 +0*t13*t15-1.2109466219727179d-1*t14-5.879979372530072d-3 265 3 *t1*t13)*t6/t7**2 266 fnc(iq) = 7.937005259840999d-1*gammabb*t1*t12*t5*t9*wght+fnc 267 1 (iq) 268 Amat(iq,D1_RB) = t1*t5*(7.937005259840999d-1*gammabb*t12*(8. 269 1 129082500000001d-4*t11*t16*t2/t9**2+9.483929583333336d-4* 270 2 rhob*t10*t2/t4**3.166666666666667d+0)*t9*wght+7.937005259 271 3 840999d-1*gammabb*t12*t16*wght)-7.937005259840999d-1*gamm 272 4 abb*t12*t5*t9*wght/rhob**2-2.6456684199469993d-1*gammabb* 273 5 t12*t13*t9*wght/rhob**1.3333333333333333d+0+Amat(iq,D1_RB 274 6 ) 275 Cmat(iq,D1_GBB) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 276 1 26028528009071d-4*t1*t12*t2*wght/t4**1.5000000000000002d+ 277 2 0+Cmat(iq,D1_GBB) 278 endif ! rhoa.gt.tol_rho.and.rhob.gt.tol_rho 279 endif ! ipol.eq.1 280 enddo ! iq 281 end 282C> 283C> \brief Evaluate the nwxcm_c_perdew86 functional [1] 284C> 285C> \f{eqnarray*}{ 286C> {\it t_1} &=& \rho_\beta+\rho_\alpha\\\\ 287C> {\it t_2} &=& {{1}\over{{\it t_1}}}\\\\ 288C> {\it t_3} &=& \rho_\alpha-\rho_\beta\\\\ 289C> {\it t_4} &=& {{1}\over{{\it t_1}^{{{2}\over{3}}}}}\\\\ 290C> {\it t_5} &=& \left| {\it t_1}\right| \\\\ 291C> {\it t_6} &=& {{1}\over{{\it t_5}^{{{1}\over{3}}}}}\\\\ 292C> {\it t_7} &=& {{0.01443307452126544\,{\it t_6} 293C> +2.843543831490385 \times 10^{-6}\,{\it t_4}+0.002568} 294C> \over{0.01763993811759022\,{\it t_4}\,{\it t_6} 295C> +5.411317332115466\,{\it t_6}+0.1816419932959077\,{\it t_4} 296C> +1.0}}+0.001667\\\\ 297C> {\it t_8} &=& \sigma_{\beta\beta}+2.0\,\sigma_{\alpha\beta} 298C> +\sigma_{\alpha\alpha}\\\\ 299C> {\it t_9} &=& {{1}\over{\rho_s^{{{2}\over{3}}}}}\\\\ 300C> {\it t_{10}} &=& \left| \rho_s\right| \\\\ 301C> {\it t_{11}} &=& {{1}\over{{\it t_{10}}^{{{1}\over{3}}}}}\\\\ 302C> {\it t_{12}} &=& {{0.01443307452126544\,{\it t_{11}} 303C> +2.843543831490385 \times 10^{-6}\,{\it t_9}+0.002568} 304C> \over{0.01763993811759022\,{\it t_9}\,{\it t_{11}} 305C> +5.411317332115466\,{\it t_{11}}+0.1816419932959077\,{\it t_9} 306C> +1.0}}+0.001667\\\\ 307C> f &=& {{0.7937005259840999\,{\it t_2}\,{\it t_7}\,{\it t_6} 308C> \,{\it t_8}}\over{\sqrt{0.3149802624737183\,\left(1.0\,{ 309C> \it t_3}\,{\it t_2}+1.0\right)^{{{5}\over{3}}} 310C> +0.3149802624737183\,\left(1.0-1.0\,{\it t_3}\,{ 311C> \it t_2}\right)^{{{5}\over{3}}}} 312C> \,e^{{{8.129082500000001 \times 10^{-4}\,\sqrt{{ 313C> \it t_8}}}\over{{\it t_7}\,{\it t_5}^{1.166666666666667}}}}}}\\\\ 314C> g &=& 0\\\\ 315C> G &=& {{0.7937005259840999\,{\it t_{12}}\,{\it t_{11}} 316C> \,\sigma_{ss}}\over{\rho_s\,e^{{{8.129082500000001 \times 10^{ 317C> -4}\,\sqrt{\sigma_{ss}}}\over{{\it t_{12}}\,{ 318C> \it t_{10}}^{1.166666666666667}}}}}}\\\\ 319C> \f} 320C> 321C> Code generated with Maxima 5.34.0 [2] 322C> driven by autoxc [3]. 323C> 324C> ### References ### 325C> 326C> [1] JP Perdew, Phys.Rev. B 33, 8822 (1986) , DOI: 327C> <a href="https://doi.org/10.1103/PhysRevB.33.8822 "> 328C> 10.1103/PhysRevB.33.8822 </a> 329C> 330C> [2] Maxima, a computer algebra system, 331C> <a href="http://maxima.sourceforge.net/"> 332C> http://maxima.sourceforge.net/</a> 333C> 334C> [3] autoxc, revision 27097 2015-05-08 335C> 336 subroutine nwxcm_c_perdew86_d2(param,tol_rho,ipol,nq,wght, 337 +rho,rgamma,fnc,Amat,Amat2,Cmat,Cmat2) 338c $Id: $ 339#ifdef NWXC_QUAD_PREC 340 implicit real(kind=selected_real_kind(30))(a-h,o-z),integer(i-n) 341 integer, parameter :: rk=selected_real_kind(30) 342#else 343 implicit real(kind=selected_real_kind(15))(a-h,o-z),integer(i-n) 344 integer, parameter :: rk=selected_real_kind(15) 345#endif 346#include "nwxc_param.fh" 347 double precision param(*) !< [Input] Parameters of functional 348 double precision tol_rho !< [Input] The lower limit on the density 349 integer ipol !< [Input] The number of spin channels 350 integer nq !< [Input] The number of points 351 double precision wght !< [Input] The weight of the functional 352 double precision rho(nq,*) !< [Input] The density 353 double precision rgamma(nq,*) !< [Input] The norm of the density 354 !< gradients 355 double precision fnc(nq) !< [Output] The value of the functional 356c 357c Sampling Matrices for the XC Kernel 358c 359 double precision Amat(nq,*) !< [Output] The derivative wrt rho 360 double precision Cmat(nq,*) !< [Output] The derivative wrt rgamma 361c 362c Sampling Matrices for the XC Kernel 363c 364 double precision Amat2(nq,*) !< [Output] The 2nd derivative wrt rho 365 double precision Cmat2(nq,*) !< [Output] The 2nd derivative wrt rgamma 366 !< and possibly rho 367 integer iq 368 double precision tmp 369 double precision rhoa,rhob 370 double precision gammaaa,gammaab,gammabb 371 double precision taua,taub 372 double precision nwxcm_heaviside 373 external nwxcm_heaviside 374CDIR$ NOVECTOR 375 do iq = 1, nq 376 if (ipol.eq.1) then 377 rhoa = 0.5d0*rho(iq,R_T) 378 gammaaa = 0.25d0*rgamma(iq,G_TT) 379 if (rhoa.gt.tol_rho) then 380 t1 = 3.968502629920499d-1 381 t2 = 1/rhoa 382 t3 = gammaaa**5.0d-1 383 t4 = 6.299605249474366d-1 384 t5 = 1/rhoa**6.666666666666666d-1 385 t6 = 7.937005259840998d-1 386 t7 = abs(rhoa) 387 t8 = 1/t7**3.333333333333333d-1 388 t9 = 1.443307452126544d-2*t6*t8+2.8435438314903855d-6*t4*t5+ 389 1 2.568d-3 390 t10 = 5.411317332115466d+0*t6*t8+8.819969058795109d-3*t5*t8+ 391 1 1.816419932959077d-1*t4*t5+1.0d+0 392 t11 = 1/t10 393 t12 = t11*t9+1.6670000000000001d-3 394 t13 = 1/t12 395 t14 = 1/t7**1.166666666666667d+0 396 t15 = 1.9842513149602492d-1 397 t16 = 1/rhoa**1.3333333333333333d+0 398 t17 = exp(-7.242189178907065d-4*t13*t14*t3) 399 t18 = 1/t7 400 t19 = rhoa**2 401 t20 = 1/t19 402 t21 = 3.149802624737183d-1 403 t22 = 1/rhoa**1.6666666666666669d+0 404 t23 = 1/rhoa**3.333333333333333d-1 405 t24 = -2.939989686265036d-3*t22*t8-1.8037724440384886d+0*t1* 406 1 t18*t23-1.2109466219727179d-1*t21*t22-1.469994843132518d- 407 2 3*t18*t2 408 t25 = 1/t10**2 409 t26 = -4.8110248404218126d-3*t1*t18*t23-1.895695887660257d-6 410 1 *t21*t22 411 t27 = t11*t26-t24*t25*t9 412 t28 = 1/t7**3.166666666666667d+0 413 t29 = 1/t12**2 414 t30 = 7.242189178907065d-4*t14*t27*t29*t3+4.224610354362455d 415 1 -4*rhoa*t13*t28*t3 416 t31 = 1/t7**1.5000000000000002d+0 417 t32 = 4.409447366578332d-1*gammaaa*t12*t17*t18*wght/rhoa**2. 418 1 3333333333333334d+0 419 t33 = 1/rhoa**3 420 t34 = 1.5749013123685918d-1 421 t35 = 1/rhoa**2.6666666666666666d+0 422 t36 = -t25*(2.4499914052208632d-3*t35*t8+2.01824436995453d-1 423 1 *t34*t35+1.9599931241766905d-3*t18*t20+2.4050299253846513 424 2 d+0*t15*t16*t18)*t9+2*t24**2*t9/t10**3+t11*(3.15949314610 425 3 04283d-6*t34*t35+6.414699787229083d-3*t15*t16*t18)-2*t24* 426 4 t25*t26 427 t37 = 6.299605249474365d-1*t2*t6*t8*(3.1748021039363994d+0*g 428 1 ammaaa*t12*t17*(-6.688966394407221d-4*t13*t19*t3/t7**5.16 429 2 6666666666667d+0+7.242189178907065d-4*t14*t29*t3*t36-8.44 430 3 922070872491d-4*rhoa*t27*t28*t29*t3+2.1123051771812276d-4 431 4 *t13*t28*t3-1.448437835781413d-3*t14*t27**2*t3/t12**3)*wg 432 5 ht+3.1748021039363994d+0*gammaaa*t17*t36*wght+3.174802103 433 6 9363994d+0*gammaaa*t12*t17*t30**2*wght+6.349604207872799d 434 7 +0*gammaaa*t17*t27*t30*wght) 435 t38 = 3.1498026247371824d-1*t20*t6*t8*(-6.349604207872799d+0 436 1 *gammaaa*t12*t17*t30*wght-6.349604207872799d+0*gammaaa*t1 437 2 7*t27*wght) 438 t39 = 4.19973683298291d-1*t15*t16*t18*(-3.1748021039363994d+ 439 1 0*gammaaa*t12*t17*t30*wght-3.1748021039363994d+0*gammaaa* 440 2 t17*t27*wght) 441 t40 = 1.259921049894873d+0*t15*t16*t18*(-1.0582673679787998d 442 1 +0*gammaaa*t12*t17*t30*wght-1.0582673679787998d+0*gammaaa 443 2 *t17*t27*wght) 444 t41 = 1/t7**2.166666666666667d+0 445 t42 = 1/t3 446 t43 = 6.299605249474365d-1*t2*t6*t8*(3.1748021039363994d+0*g 447 1 ammaaa*t12*t17*(9.052736473633831d-5*t14*t27*t29*t42+5.28 448 2 0762942953069d-5*rhoa*t13*t28*t42)*wght-2.874064680287446 449 3 4d-4*t14*t17*t3*t30*wght+7.937005259840999d-1*t12*t17*t30 450 4 *wght-2.8740646802874464d-4*t13*t14*t17*t27*t3*wght+7.937 451 5 005259840999d-1*t17*t27*wght)-2.5d-1*t12*t17*t20*t6*t8*wg 452 6 ht+6.035157649089219d-5*t1*t16*t17*t3*t41*wght+7.18516170 453 7 0718616d-5*t17*t20*t3*t31*wght-3.333333333333333d-1*t12*t 454 8 15*t16*t17*t18*wght 455 t44 = 1/t7**2.666666666666667d+0 456 t45 = 1.300907507941046d-8*t13*t17*t2*t44*wght-5.38887127553 457 1 8962d-5*t17*t2*t31*t42*wght 458 fnc(iq) = 4.000000000000001d+0*gammaaa*t1*t12*t2*exp(-7.2421 459 1 89178907064d-4*t13*t14*t3)*t8*wght+fnc(iq) 460 Amat(iq,D1_RA) = 6.299605249474365d-1*t2*t6*t8*(3.1748021039 461 1 363994d+0*gammaaa*t12*t17*t30*wght+3.1748021039363994d+0* 462 2 gammaaa*t17*t27*wght)-1.0d+0*gammaaa*t12*t17*t20*t6*t8*wg 463 3 ht-1.3333333333333333d+0*gammaaa*t12*t15*t16*t17*t18*wght 464 4 +Amat(iq,D1_RA) 465 Cmat(iq,D1_GAA) = 5.0d-1*t12*t17*t2*t6*t8*wght-1.43703234014 466 1 37232d-4*t17*t2*t3*t31*wght+Cmat(iq,D1_GAA) 467 Cmat(iq,D1_GAB) = 1.0d+0*t12*t17*t2*t6*t8*wght-2.87406468028 468 1 74464d-4*t17*t2*t3*t31*wght+Cmat(iq,D1_GAB) 469 Amat2(iq,D2_RA_RA) = 7.222222222222222d-1*gammaaa*t12*t17*t3 470 1 3*t6*t8*wght+t40+t39+t38+t37+t32+Amat2(iq,D2_RA_RA) 471 Amat2(iq,D2_RA_RB) = 1.2777777777777777d+0*gammaaa*t12*t17*t 472 1 33*t6*t8*wght+t40+t39+t38+t37+t32+Amat2(iq,D2_RA_RB) 473 Cmat2(iq,D2_RA_GAA) = t43+Cmat2(iq,D2_RA_GAA) 474 Cmat2(iq,D2_RA_GAB) = 6.299605249474365d-1*t2*t6*t8*(3.17480 475 1 21039363994d+0*gammaaa*t12*t17*(1.8105472947267662d-4*t14 476 2 *t27*t29*t42+1.0561525885906138d-4*rhoa*t13*t28*t42)*wght 477 3 -5.748129360574893d-4*t14*t17*t3*t30*wght+1.5874010519681 478 4 997d+0*t12*t17*t30*wght-5.748129360574893d-4*t13*t14*t17* 479 5 t27*t3*wght+1.5874010519681997d+0*t17*t27*wght)-5.0d-1*t1 480 6 2*t17*t20*t6*t8*wght+1.2070315298178438d-4*t1*t16*t17*t3* 481 7 t41*wght+1.4370323401437232d-4*t17*t20*t3*t31*wght-6.6666 482 8 66666666666d-1*t12*t15*t16*t17*t18*wght+Cmat2(iq,D2_RA_GA 483 9 B) 484 Cmat2(iq,D2_RA_GBB) = t43+Cmat2(iq,D2_RA_GBB) 485 Cmat2(iq,D2_GAA_GAA) = t45+Cmat2(iq,D2_GAA_GAA) 486 Cmat2(iq,D2_GAA_GAB) = 2.601815015882092d-8*t13*t17*t2*t44*w 487 1 ght-1.0777742551077923d-4*t17*t2*t31*t42*wght+Cmat2(iq,D2 488 2 _GAA_GAB) 489 Cmat2(iq,D2_GAA_GBB) = t45+Cmat2(iq,D2_GAA_GBB) 490 Cmat2(iq,D2_GAB_GAB) = 5.203630031764184d-8*t13*t17*t2*t44*w 491 1 ght-2.1555485102155847d-4*t17*t2*t31*t42*wght+Cmat2(iq,D2 492 2 _GAB_GAB) 493 endif ! rhoa.gt.tol_rho 494 else ! ipol.eq.1 495 rhoa = rho(iq,R_A) 496 rhob = rho(iq,R_B) 497 gammaaa = rgamma(iq,G_AA) 498 gammaab = rgamma(iq,G_AB) 499 gammabb = rgamma(iq,G_BB) 500 if (rhoa.gt.tol_rho.and.rhob.gt.tol_rho) then 501 t1 = gammabb+2.0d+0*gammaab+gammaaa 502 t2 = rhob+rhoa 503 t3 = 1/t2 504 t4 = rhoa-rhob 505 t5 = 1.0d+0-1.0d+0*t3*t4 506 t6 = 1.0d+0*t3*t4+1.0d+0 507 t7 = (3.149802624737183d-1*t6**1.6666666666666669d+0+3.14980 508 1 2624737183d-1*t5**1.6666666666666669d+0)**5.0d-1 509 t8 = 1/t7 510 t9 = t1**5.0d-1 511 t10 = 1/t2**6.666666666666666d-1 512 t11 = abs(t2) 513 t12 = 1/t11**3.333333333333333d-1 514 t13 = 1.443307452126544d-2*t12+2.8435438314903855d-6*t10+2.5 515 1 68d-3 516 t14 = 1.7639938117590215d-2*t10*t12+5.411317332115466d+0*t12 517 1 +1.816419932959077d-1*t10+1.0d+0 518 t15 = 1/t14 519 t16 = t13*t15+1.6670000000000001d-3 520 t17 = 1/t16 521 t18 = 1/t11**1.166666666666667d+0 522 t19 = exp(-8.129082500000001d-4*t17*t18*t9) 523 t20 = 1/t2**1.3333333333333333d+0 524 t21 = 1/t11 525 t22 = -2.6456684199469993d-1*t1*t16*t19*t20*t21*t8*wght 526 t23 = t2**2 527 t24 = 1/t23 528 t25 = 1.0d+0*t24*t4 529 t26 = -1.0d+0*t3 530 t27 = t26+t25 531 t28 = t5**6.666666666666666d-1 532 t29 = -1.0d+0*t24*t4 533 t30 = 1.0d+0*t3 534 t31 = t30+t29 535 t32 = t6**6.666666666666666d-1 536 t33 = 5.249671041228638d-1*t31*t32+5.249671041228638d-1*t27* 537 1 t28 538 t34 = 1/t7**3 539 t35 = -7.937005259840999d-1*t1*t12*t16*t19*t24*t8*wght 540 t36 = 1/t2**1.6666666666666669d+0 541 t37 = 1/t2**3.333333333333333d-1 542 t38 = -1.8037724440384886d+0*t21*t37-1.1759958745060142d-2*t 543 1 12*t36-1.2109466219727179d-1*t36-5.879979372530072d-3*t21 544 2 *t3 545 t39 = 1/t14**2 546 t40 = -4.8110248404218126d-3*t21*t37-1.895695887660257d-6*t3 547 1 6 548 t41 = t15*t40-t13*t38*t39 549 t42 = 1/t11**3.166666666666667d+0 550 t43 = 1/t16**2 551 t44 = 8.129082500000001d-4*t18*t41*t43*t9+9.483929583333336d 552 1 -4*t17*t2*t42*t9 553 t45 = 7.937005259840999d-1*t1*t16*t19*t44*wght+7.93700525984 554 1 0999d-1*t1*t19*t41*wght 555 t46 = t12*t3*t45*t8 556 t47 = t30+t25 557 t48 = t29+t26 558 t49 = 5.249671041228638d-1*t32*t48+5.249671041228638d-1*t28* 559 1 t47 560 t50 = 1/t11**1.5000000000000002d+0 561 t51 = 7.937005259840999d-1*t12*t16*t19*t3*t8*wght-3.22602852 562 1 8009071d-4*t19*t3*t50*t8*t9*wght 563 t52 = 8.818894733156665d-1*t1*t16*t19*t21*t8*wght/t2**2.3333 564 1 333333333334d+0 565 t53 = 1/t7**5 566 t54 = 1/t2**3 567 t55 = 1.5874010519681997d+0*t1*t12*t16*t19*t54*t8*wght 568 t56 = 1/t2**2.6666666666666666d+0 569 t57 = -t13*t39*(1.9599931241766905d-2*t12*t56+2.018244369954 570 1 53d-1*t56+1.5679944993413522d-2*t21*t24+2.405029925384651 571 2 3d+0*t20*t21)+t15*(3.1594931461004283d-6*t56+6.4146997872 572 3 29083d-3*t20*t21)-2*t38*t39*t40+2*t13*t38**2/t14**3 573 t58 = t12*t3*t8*(7.937005259840999d-1*t1*t16*t19*(8.12908250 574 1 0000001d-4*t18*t43*t57*t9-1.8967859166666673d-3*t2*t41*t4 575 2 2*t43*t9+9.483929583333336d-4*t17*t42*t9-1.62581650000000 576 3 02d-3*t18*t41**2*t9/t16**3-3.003244368055557d-3*t17*t23*t 577 4 9/t11**5.166666666666667d+0)*wght+7.937005259840999d-1*t1 578 5 *t19*t57*wght+7.937005259840999d-1*t1*t16*t19*t44**2*wght 579 6 +1.5874010519681997d+0*t1*t19*t41*t44*wght) 580 t59 = t12*t24*t8*(-1.5874010519681997d+0*t1*t16*t19*t44*wght 581 1 -1.5874010519681997d+0*t1*t19*t41*wght) 582 t60 = 3.333333333333333d-1*t20*t21*t8*(-7.937005259840999d-1 583 1 *t1*t16*t19*t44*wght-7.937005259840999d-1*t1*t19*t41*wght 584 2 ) 585 t61 = t20*t21*t8*(-2.6456684199469993d-1*t1*t16*t19*t44*wght 586 1 -2.6456684199469993d-1*t1*t19*t41*wght) 587 t62 = -3.9685026299204995d-1*t1*t19*t33*t41*wght 588 t63 = 1/t5**3.333333333333333d-1 589 t64 = -2.0d+0*t4*t54 590 t65 = 2.0d+0*t24 591 t66 = 1/t6**3.333333333333333d-1 592 t67 = 2.0d+0*t4*t54 593 t68 = -2.0d+0*t24 594 t69 = -3.9685026299204995d-1*t1*t16*t19*t33*t44*wght 595 t70 = -5.0d-1*t12*t3*t34*t45*t49 596 t71 = 1/t11**2.166666666666667d+0 597 t72 = 1.0753428426696904d-4*t19*t20*t71*t8*t9*wght 598 t73 = 3.226028528009071d-4*t19*t24*t50*t8*t9*wght 599 t74 = -2.6456684199469993d-1*t16*t19*t20*t21*t8*wght 600 t75 = -7.937005259840999d-1*t12*t16*t19*t24*t8*wght 601 t76 = 1/t9 602 t77 = t12*t3*t8*(-3.226028528009071d-4*t18*t19*t44*t9*wght-3 603 1 .226028528009071d-4*t17*t18*t19*t41*t9*wght+7.93700525984 604 2 0999d-1*t1*t16*t19*(4.0645412500000005d-4*t18*t41*t43*t76 605 3 +4.741964791666668d-4*t17*t2*t42*t76)*wght+7.937005259840 606 4 999d-1*t16*t19*t44*wght+7.937005259840999d-1*t19*t41*wght 607 5 ) 608 t78 = 1.6130142640045356d-4*t19*t3*t33*t34*t50*t9*wght-3.968 609 1 5026299204995d-1*t12*t16*t19*t3*t33*t34*wght+t77+t75+t74+ 610 2 t73+t72 611 t79 = 2.1506856853393808d-4*t19*t20*t71*t8*t9*wght 612 t80 = 6.452057056018142d-4*t19*t24*t50*t8*t9*wght 613 t81 = -5.291336839893999d-1*t16*t19*t20*t21*t8*wght 614 t82 = -1.5874010519681997d+0*t12*t16*t19*t24*t8*wght 615 t83 = t12*t3*t8*(-6.452057056018142d-4*t18*t19*t44*t9*wght-6 616 1 .452057056018142d-4*t17*t18*t19*t41*t9*wght+7.93700525984 617 2 0999d-1*t1*t16*t19*(8.129082500000001d-4*t18*t41*t43*t76+ 618 3 9.483929583333336d-4*t17*t2*t42*t76)*wght+1.5874010519681 619 4 997d+0*t16*t19*t44*wght+1.5874010519681997d+0*t19*t41*wgh 620 5 t) 621 t84 = 1.6130142640045356d-4*t19*t3*t34*t49*t50*t9*wght-3.968 622 1 5026299204995d-1*t12*t16*t19*t3*t34*t49*wght+t77+t75+t74+ 623 2 t73+t72 624 t85 = 1/t11**2.666666666666667d+0 625 t86 = 1.311232602576965d-7*t17*t19*t3*t8*t85*wght-4.83904279 626 1 20136064d-4*t19*t3*t50*t76*t8*wght 627 t87 = 2.62246520515393d-7*t17*t19*t3*t8*t85*wght-9.678085584 628 1 027213d-4*t19*t3*t50*t76*t8*wght 629 fnc(iq) = 7.937005259840999d-1*t1*t12*t16*t19*t3*t8*wght+fnc 630 1 (iq) 631 Amat(iq,D1_RA) = -3.9685026299204995d-1*t1*t12*t16*t19*t3*t3 632 1 3*t34*wght+t46+t35+t22+Amat(iq,D1_RA) 633 Amat(iq,D1_RB) = -3.9685026299204995d-1*t1*t12*t16*t19*t3*t3 634 1 4*t49*wght+t46+t35+t22+Amat(iq,D1_RB) 635 Cmat(iq,D1_GAA) = t51+Cmat(iq,D1_GAA) 636 Cmat(iq,D1_GAB) = -6.452057056018142d-4*t19*t3*t50*t8*t9*wgh 637 1 t+1.5874010519681997d+0*t12*t16*t19*t3*t8*wght+Cmat(iq,D1 638 2 _GAB) 639 Cmat(iq,D1_GBB) = t51+Cmat(iq,D1_GBB) 640 Amat2(iq,D2_RA_RA) = t12*t3*t34*(-3.9685026299204995d-1*t1*t 641 1 16*t19*(5.249671041228638d-1*t32*(t68+t67)+3.499780694152 642 2 4257d-1*t31**2*t66+5.249671041228638d-1*t28*(t65+t64)+3.4 643 3 997806941524257d-1*t27**2*t63)*wght+t69+t62)+5.9527539448 644 4 80749d-1*t1*t12*t16*t19*t3*t33**2*t53*wght+7.937005259840 645 5 999d-1*t1*t12*t16*t19*t24*t33*t34*wght+2.6456684199469993 646 6 d-1*t1*t16*t19*t20*t21*t33*t34*wght+t61+t60+t59+t58+t55+t 647 7 52-5.0d-1*t12*t3*t33*t34*t45+Amat2(iq,D2_RA_RA) 648 Amat2(iq,D2_RA_RB) = t12*t3*t34*(-3.9685026299204995d-1*t1*t 649 1 16*t19*(3.4997806941524257d-1*t31*t48*t66+3.4997806941524 650 2 257d-1*t27*t47*t63+1.0499342082457278d+0*t32*t4*t54-1.049 651 3 9342082457278d+0*t28*t4*t54)*wght+t69+t62)+t12*t24*t34*(3 652 4 .9685026299204995d-1*t1*t16*t19*t49*wght+3.96850262992049 653 5 95d-1*t1*t16*t19*t33*wght)+t20*t21*t34*(1.322834209973499 654 6 7d-1*t1*t16*t19*t49*wght+1.3228342099734997d-1*t1*t16*t19 655 7 *t33*wght)+5.952753944880749d-1*t1*t12*t16*t19*t3*t33*t49 656 8 *t53*wght+t70+t61+t60+t59+t58+t55+t52+Amat2(iq,D2_RA_RB) 657 Amat2(iq,D2_RB_RB) = t12*t3*t34*(-3.9685026299204995d-1*t1*t 658 1 16*t19*(5.249671041228638d-1*t28*(t68+t64)+5.249671041228 659 2 638d-1*t32*(t67+t65)+3.4997806941524257d-1*t48**2*t66+3.4 660 3 997806941524257d-1*t47**2*t63)*wght-3.9685026299204995d-1 661 4 *t1*t16*t19*t44*t49*wght-3.9685026299204995d-1*t1*t19*t41 662 5 *t49*wght)+5.952753944880749d-1*t1*t12*t16*t19*t3*t49**2* 663 6 t53*wght+7.937005259840999d-1*t1*t12*t16*t19*t24*t34*t49* 664 7 wght+2.6456684199469993d-1*t1*t16*t19*t20*t21*t34*t49*wgh 665 8 t+t70+t61+t60+t59+t58+t55+t52+Amat2(iq,D2_RB_RB) 666 Cmat2(iq,D2_RA_GAA) = t78+Cmat2(iq,D2_RA_GAA) 667 Cmat2(iq,D2_RA_GAB) = 3.226028528009071d-4*t19*t3*t33*t34*t5 668 1 0*t9*wght-7.937005259840999d-1*t12*t16*t19*t3*t33*t34*wgh 669 2 t+t83+t82+t81+t80+t79+Cmat2(iq,D2_RA_GAB) 670 Cmat2(iq,D2_RA_GBB) = t78+Cmat2(iq,D2_RA_GBB) 671 Cmat2(iq,D2_RB_GAA) = t84+Cmat2(iq,D2_RB_GAA) 672 Cmat2(iq,D2_RB_GAB) = 3.226028528009071d-4*t19*t3*t34*t49*t5 673 1 0*t9*wght-7.937005259840999d-1*t12*t16*t19*t3*t34*t49*wgh 674 2 t+t83+t82+t81+t80+t79+Cmat2(iq,D2_RB_GAB) 675 Cmat2(iq,D2_RB_GBB) = t84+Cmat2(iq,D2_RB_GBB) 676 Cmat2(iq,D2_GAA_GAA) = t86+Cmat2(iq,D2_GAA_GAA) 677 Cmat2(iq,D2_GAA_GAB) = t87+Cmat2(iq,D2_GAA_GAB) 678 Cmat2(iq,D2_GAA_GBB) = t86+Cmat2(iq,D2_GAA_GBB) 679 Cmat2(iq,D2_GAB_GAB) = 5.24493041030786d-7*t17*t19*t3*t8*t85 680 1 *wght-1.9356171168054426d-3*t19*t3*t50*t76*t8*wght+Cmat2( 681 2 iq,D2_GAB_GAB) 682 Cmat2(iq,D2_GAB_GBB) = t87+Cmat2(iq,D2_GAB_GBB) 683 Cmat2(iq,D2_GBB_GBB) = t86+Cmat2(iq,D2_GBB_GBB) 684 elseif (rhoa.gt.tol_rho.and.rhob.le.tol_rho) then 685 t1 = 1/rhoa 686 t2 = gammaaa**5.0d-1 687 t3 = 1/rhoa**6.666666666666666d-1 688 t4 = abs(rhoa) 689 t5 = 1/t4**3.333333333333333d-1 690 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 691 1 -3 692 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 693 1 16419932959077d-1*t3+1.0d+0 694 t8 = 1/t7 695 t9 = t6*t8+1.6670000000000001d-3 696 t10 = 1/t9 697 t11 = 1/t4**1.166666666666667d+0 698 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 699 t13 = 1/rhoa**1.3333333333333333d+0 700 t14 = 1/t4 701 t15 = rhoa**2 702 t16 = 1/t15 703 t17 = 1/rhoa**1.6666666666666669d+0 704 t18 = 1/rhoa**3.333333333333333d-1 705 t19 = -1.1759958745060142d-2*t17*t5-1.8037724440384886d+0*t1 706 1 4*t18-1.2109466219727179d-1*t17-5.879979372530072d-3*t1*t 707 2 14 708 t20 = 1/t7**2 709 t21 = -4.8110248404218126d-3*t14*t18-1.895695887660257d-6*t1 710 1 7 711 t22 = t21*t8-t19*t20*t6 712 t23 = 1/t4**3.166666666666667d+0 713 t24 = 1/t9**2 714 t25 = 8.129082500000001d-4*t11*t2*t22*t24+9.483929583333336d 715 1 -4*rhoa*t10*t2*t23 716 t26 = 1/t4**1.5000000000000002d+0 717 t27 = 1/rhoa**2.6666666666666666d+0 718 t28 = (3.1594931461004283d-6*t27+6.414699787229083d-3*t13*t1 719 1 4)*t8+2*t19**2*t6/t7**3-t20*(1.9599931241766905d-2*t27*t5 720 2 +2.01824436995453d-1*t27+1.5679944993413522d-2*t14*t16+2. 721 3 4050299253846513d+0*t13*t14)*t6-2*t19*t20*t21 722 t29 = 1/t2 723 fnc(iq) = 7.937005259840999d-1*gammaaa*t1*t12*t5*t9*wght+fnc 724 1 (iq) 725 Amat(iq,D1_RA) = t1*t5*(7.937005259840999d-1*gammaaa*t12*t25 726 1 *t9*wght+7.937005259840999d-1*gammaaa*t12*t22*wght)-7.937 727 2 005259840999d-1*gammaaa*t12*t16*t5*t9*wght-2.645668419946 728 3 9993d-1*gammaaa*t12*t13*t14*t9*wght+Amat(iq,D1_RA) 729 Cmat(iq,D1_GAA) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 730 1 26028528009071d-4*t1*t12*t2*t26*wght+Cmat(iq,D1_GAA) 731 Amat2(iq,D2_RA_RA) = t1*t5*(7.937005259840999d-1*gammaaa*t12 732 1 *(-1.6258165000000002d-3*t11*t2*t22**2/t9**3-3.0032443680 733 2 55557d-3*t10*t15*t2/t4**5.166666666666667d+0+8.1290825000 734 3 00001d-4*t11*t2*t24*t28-1.8967859166666673d-3*rhoa*t2*t22 735 4 *t23*t24+9.483929583333336d-4*t10*t2*t23)*t9*wght+7.93700 736 5 5259840999d-1*gammaaa*t12*t25**2*t9*wght+7.93700525984099 737 6 9d-1*gammaaa*t12*t28*wght+1.5874010519681997d+0*gammaaa*t 738 7 12*t22*t25*wght)+t13*t14*(-2.6456684199469993d-1*gammaaa* 739 8 t12*t25*t9*wght-2.6456684199469993d-1*gammaaa*t12*t22*wgh 740 9 t)+3.333333333333333d-1*t13*t14*(-7.937005259840999d-1*ga 741 : mmaaa*t12*t25*t9*wght-7.937005259840999d-1*gammaaa*t12*t2 742 ; 2*wght)+t16*t5*(-1.5874010519681997d+0*gammaaa*t12*t25*t9 743 < *wght-1.5874010519681997d+0*gammaaa*t12*t22*wght)+1.58740 744 = 10519681997d+0*gammaaa*t12*t5*t9*wght/rhoa**3+8.818894733 745 > 156665d-1*gammaaa*t12*t14*t9*wght/rhoa**2.333333333333333 746 ? 4d+0+Amat2(iq,D2_RA_RA) 747 Cmat2(iq,D2_RA_GAA) = t1*t5*(7.937005259840999d-1*gammaaa*t1 748 1 2*(4.0645412500000005d-4*t11*t22*t24*t29+4.74196479166666 749 2 8d-4*rhoa*t10*t23*t29)*t9*wght+7.937005259840999d-1*t12*t 750 3 25*t9*wght-3.226028528009071d-4*t11*t12*t2*t25*wght-3.226 751 4 028528009071d-4*t10*t11*t12*t2*t22*wght+7.937005259840999 752 5 d-1*t12*t22*wght)-7.937005259840999d-1*t12*t16*t5*t9*wght 753 6 -2.6456684199469993d-1*t12*t13*t14*t9*wght+1.075342842669 754 7 6904d-4*t12*t13*t2*wght/t4**2.166666666666667d+0+3.226028 755 8 528009071d-4*t12*t16*t2*t26*wght+Cmat2(iq,D2_RA_GAA) 756 Cmat2(iq,D2_GAA_GAA) = 1.311232602576965d-7*t1*t10*t12*wght/ 757 1 t4**2.666666666666667d+0-4.8390427920136064d-4*t1*t12*t26 758 2 *t29*wght+Cmat2(iq,D2_GAA_GAA) 759 elseif (rhoa.le.tol_rho.and.rhob.gt.tol_rho) then 760 t1 = 1/rhob 761 t2 = gammabb**5.0d-1 762 t3 = 1/rhob**6.666666666666666d-1 763 t4 = abs(rhob) 764 t5 = 1/t4**3.333333333333333d-1 765 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 766 1 -3 767 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 768 1 16419932959077d-1*t3+1.0d+0 769 t8 = 1/t7 770 t9 = t6*t8+1.6670000000000001d-3 771 t10 = 1/t9 772 t11 = 1/t4**1.166666666666667d+0 773 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 774 t13 = 1/rhob**1.3333333333333333d+0 775 t14 = 1/t4 776 t15 = rhob**2 777 t16 = 1/t15 778 t17 = 1/rhob**1.6666666666666669d+0 779 t18 = 1/rhob**3.333333333333333d-1 780 t19 = -1.1759958745060142d-2*t17*t5-1.8037724440384886d+0*t1 781 1 4*t18-1.2109466219727179d-1*t17-5.879979372530072d-3*t1*t 782 2 14 783 t20 = 1/t7**2 784 t21 = -4.8110248404218126d-3*t14*t18-1.895695887660257d-6*t1 785 1 7 786 t22 = t21*t8-t19*t20*t6 787 t23 = 1/t4**3.166666666666667d+0 788 t24 = 1/t9**2 789 t25 = 8.129082500000001d-4*t11*t2*t22*t24+9.483929583333336d 790 1 -4*rhob*t10*t2*t23 791 t26 = 1/t4**1.5000000000000002d+0 792 t27 = 1/rhob**2.6666666666666666d+0 793 t28 = (3.1594931461004283d-6*t27+6.414699787229083d-3*t13*t1 794 1 4)*t8+2*t19**2*t6/t7**3-t20*(1.9599931241766905d-2*t27*t5 795 2 +2.01824436995453d-1*t27+1.5679944993413522d-2*t14*t16+2. 796 3 4050299253846513d+0*t13*t14)*t6-2*t19*t20*t21 797 t29 = 1/t2 798 fnc(iq) = 7.937005259840999d-1*gammabb*t1*t12*t5*t9*wght+fnc 799 1 (iq) 800 Amat(iq,D1_RB) = t1*t5*(7.937005259840999d-1*gammabb*t12*t25 801 1 *t9*wght+7.937005259840999d-1*gammabb*t12*t22*wght)-7.937 802 2 005259840999d-1*gammabb*t12*t16*t5*t9*wght-2.645668419946 803 3 9993d-1*gammabb*t12*t13*t14*t9*wght+Amat(iq,D1_RB) 804 Cmat(iq,D1_GBB) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 805 1 26028528009071d-4*t1*t12*t2*t26*wght+Cmat(iq,D1_GBB) 806 Amat2(iq,D2_RB_RB) = t1*t5*(7.937005259840999d-1*gammabb*t12 807 1 *(-1.6258165000000002d-3*t11*t2*t22**2/t9**3-3.0032443680 808 2 55557d-3*t10*t15*t2/t4**5.166666666666667d+0+8.1290825000 809 3 00001d-4*t11*t2*t24*t28-1.8967859166666673d-3*rhob*t2*t22 810 4 *t23*t24+9.483929583333336d-4*t10*t2*t23)*t9*wght+7.93700 811 5 5259840999d-1*gammabb*t12*t25**2*t9*wght+7.93700525984099 812 6 9d-1*gammabb*t12*t28*wght+1.5874010519681997d+0*gammabb*t 813 7 12*t22*t25*wght)+t13*t14*(-2.6456684199469993d-1*gammabb* 814 8 t12*t25*t9*wght-2.6456684199469993d-1*gammabb*t12*t22*wgh 815 9 t)+3.333333333333333d-1*t13*t14*(-7.937005259840999d-1*ga 816 : mmabb*t12*t25*t9*wght-7.937005259840999d-1*gammabb*t12*t2 817 ; 2*wght)+t16*t5*(-1.5874010519681997d+0*gammabb*t12*t25*t9 818 < *wght-1.5874010519681997d+0*gammabb*t12*t22*wght)+1.58740 819 = 10519681997d+0*gammabb*t12*t5*t9*wght/rhob**3+8.818894733 820 > 156665d-1*gammabb*t12*t14*t9*wght/rhob**2.333333333333333 821 ? 4d+0+Amat2(iq,D2_RB_RB) 822 Cmat2(iq,D2_RB_GBB) = t1*t5*(7.937005259840999d-1*gammabb*t1 823 1 2*(4.0645412500000005d-4*t11*t22*t24*t29+4.74196479166666 824 2 8d-4*rhob*t10*t23*t29)*t9*wght+7.937005259840999d-1*t12*t 825 3 25*t9*wght-3.226028528009071d-4*t11*t12*t2*t25*wght-3.226 826 4 028528009071d-4*t10*t11*t12*t2*t22*wght+7.937005259840999 827 5 d-1*t12*t22*wght)-7.937005259840999d-1*t12*t16*t5*t9*wght 828 6 -2.6456684199469993d-1*t12*t13*t14*t9*wght+1.075342842669 829 7 6904d-4*t12*t13*t2*wght/t4**2.166666666666667d+0+3.226028 830 8 528009071d-4*t12*t16*t2*t26*wght+Cmat2(iq,D2_RB_GBB) 831 Cmat2(iq,D2_GBB_GBB) = 1.311232602576965d-7*t1*t10*t12*wght/ 832 1 t4**2.666666666666667d+0-4.8390427920136064d-4*t1*t12*t26 833 2 *t29*wght+Cmat2(iq,D2_GBB_GBB) 834 endif ! rhoa.gt.tol_rho.and.rhob.gt.tol_rho 835 endif ! ipol.eq.1 836 enddo ! iq 837 end 838C> 839C> \brief Evaluate the nwxcm_c_perdew86 functional [1] 840C> 841C> \f{eqnarray*}{ 842C> {\it t_1} &=& \rho_\beta+\rho_\alpha\\\\ 843C> {\it t_2} &=& {{1}\over{{\it t_1}}}\\\\ 844C> {\it t_3} &=& \rho_\alpha-\rho_\beta\\\\ 845C> {\it t_4} &=& {{1}\over{{\it t_1}^{{{2}\over{3}}}}}\\\\ 846C> {\it t_5} &=& \left| {\it t_1}\right| \\\\ 847C> {\it t_6} &=& {{1}\over{{\it t_5}^{{{1}\over{3}}}}}\\\\ 848C> {\it t_7} &=& {{0.01443307452126544\,{\it t_6} 849C> +2.843543831490385 \times 10^{-6}\,{\it t_4}+0.002568} 850C> \over{0.01763993811759022\,{\it t_4}\,{\it t_6} 851C> +5.411317332115466\,{\it t_6}+0.1816419932959077\,{\it t_4} 852C> +1.0}}+0.001667\\\\ 853C> {\it t_8} &=& \sigma_{\beta\beta}+2.0\,\sigma_{\alpha\beta} 854C> +\sigma_{\alpha\alpha}\\\\ 855C> {\it t_9} &=& {{1}\over{\rho_s^{{{2}\over{3}}}}}\\\\ 856C> {\it t_{10}} &=& \left| \rho_s\right| \\\\ 857C> {\it t_{11}} &=& {{1}\over{{\it t_{10}}^{{{1}\over{3}}}}}\\\\ 858C> {\it t_{12}} &=& {{0.01443307452126544\,{\it t_{11}} 859C> +2.843543831490385 \times 10^{-6}\,{\it t_9}+0.002568} 860C> \over{0.01763993811759022\,{\it t_9}\,{\it t_{11}} 861C> +5.411317332115466\,{\it t_{11}}+0.1816419932959077\,{\it t_9} 862C> +1.0}}+0.001667\\\\ 863C> f &=& {{0.7937005259840999\,{\it t_2}\,{\it t_7}\,{\it t_6} 864C> \,{\it t_8}}\over{\sqrt{0.3149802624737183\,\left(1.0\,{ 865C> \it t_3}\,{\it t_2}+1.0\right)^{{{5}\over{3}}} 866C> +0.3149802624737183\,\left(1.0-1.0\,{\it t_3}\,{ 867C> \it t_2}\right)^{{{5}\over{3}}}} 868C> \,e^{{{8.129082500000001 \times 10^{-4}\,\sqrt{{ 869C> \it t_8}}}\over{{\it t_7}\,{\it t_5}^{1.166666666666667}}}}}}\\\\ 870C> g &=& 0\\\\ 871C> G &=& {{0.7937005259840999\,{\it t_{12}}\,{\it t_{11}} 872C> \,\sigma_{ss}}\over{\rho_s\,e^{{{8.129082500000001 \times 10^{ 873C> -4}\,\sqrt{\sigma_{ss}}}\over{{\it t_{12}}\,{ 874C> \it t_{10}}^{1.166666666666667}}}}}}\\\\ 875C> \f} 876C> 877C> Code generated with Maxima 5.34.0 [2] 878C> driven by autoxc [3]. 879C> 880C> ### References ### 881C> 882C> [1] JP Perdew, Phys.Rev. B 33, 8822 (1986) , DOI: 883C> <a href="https://doi.org/10.1103/PhysRevB.33.8822 "> 884C> 10.1103/PhysRevB.33.8822 </a> 885C> 886C> [2] Maxima, a computer algebra system, 887C> <a href="http://maxima.sourceforge.net/"> 888C> http://maxima.sourceforge.net/</a> 889C> 890C> [3] autoxc, revision 27097 2015-05-08 891C> 892 subroutine nwxcm_c_perdew86_d3(param,tol_rho,ipol,nq,wght, 893 +rho,rgamma,fnc,Amat,Amat2,Amat3, 894 +Cmat,Cmat2,Cmat3) 895c $Id: $ 896#ifdef NWXC_QUAD_PREC 897 implicit real(kind=selected_real_kind(30))(a-h,o-z),integer(i-n) 898 integer, parameter :: rk=selected_real_kind(30) 899#else 900 implicit real(kind=selected_real_kind(15))(a-h,o-z),integer(i-n) 901 integer, parameter :: rk=selected_real_kind(15) 902#endif 903#include "nwxc_param.fh" 904 double precision param(*) !< [Input] Parameters of functional 905 double precision tol_rho !< [Input] The lower limit on the density 906 integer ipol !< [Input] The number of spin channels 907 integer nq !< [Input] The number of points 908 double precision wght !< [Input] The weight of the functional 909 double precision rho(nq,*) !< [Input] The density 910 double precision rgamma(nq,*) !< [Input] The norm of the density 911 !< gradients 912 double precision fnc(nq) !< [Output] The value of the functional 913c 914c Sampling Matrices for the XC Kernel 915c 916 double precision Amat(nq,*) !< [Output] The derivative wrt rho 917 double precision Cmat(nq,*) !< [Output] The derivative wrt rgamma 918c 919c Sampling Matrices for the XC Kernel 920c 921 double precision Amat2(nq,*) !< [Output] The 2nd derivative wrt rho 922 double precision Cmat2(nq,*) !< [Output] The 2nd derivative wrt rgamma 923 !< and possibly rho 924c 925c Sampling Matrices for the XC Kernel 926c 927 double precision Amat3(nq,*) !< [Output] The 3rd derivative wrt rho 928 double precision Cmat3(nq,*) !< [Output] The 3rd derivative wrt rgamma 929 !< and possibly rho 930 integer iq 931 double precision tmp 932 double precision rhoa,rhob 933 double precision gammaaa,gammaab,gammabb 934 double precision taua,taub 935 double precision nwxcm_heaviside 936 external nwxcm_heaviside 937CDIR$ NOVECTOR 938 do iq = 1, nq 939 if (ipol.eq.1) then 940 rhoa = 0.5d0*rho(iq,R_T) 941 gammaaa = 0.25d0*rgamma(iq,G_TT) 942 if (rhoa.gt.tol_rho) then 943 t1 = 3.968502629920499d-1 944 t2 = 1/rhoa 945 t3 = gammaaa**5.0d-1 946 t4 = 6.299605249474366d-1 947 t5 = 1/rhoa**6.666666666666666d-1 948 t6 = 7.937005259840998d-1 949 t7 = abs(rhoa) 950 t8 = 1/t7**3.333333333333333d-1 951 t9 = 1.443307452126544d-2*t6*t8+2.8435438314903855d-6*t4*t5+ 952 1 2.568d-3 953 t10 = 5.411317332115466d+0*t6*t8+8.819969058795109d-3*t5*t8+ 954 1 1.816419932959077d-1*t4*t5+1.0d+0 955 t11 = 1/t10 956 t12 = t11*t9+1.6670000000000001d-3 957 t13 = 1/t12 958 t14 = 1/t7**1.166666666666667d+0 959 t15 = 1.9842513149602492d-1 960 t16 = 1/rhoa**1.3333333333333333d+0 961 t17 = exp(-7.242189178907065d-4*t13*t14*t3) 962 t18 = 1/t7 963 t19 = rhoa**2 964 t20 = 1/t19 965 t21 = 3.149802624737183d-1 966 t22 = 1/rhoa**1.6666666666666669d+0 967 t23 = 1/rhoa**3.333333333333333d-1 968 t24 = -2.939989686265036d-3*t22*t8-1.8037724440384886d+0*t1* 969 1 t18*t23-1.2109466219727179d-1*t21*t22-1.469994843132518d- 970 2 3*t18*t2 971 t25 = 1/t10**2 972 t26 = -4.8110248404218126d-3*t1*t18*t23-1.895695887660257d-6 973 1 *t21*t22 974 t27 = t11*t26-t24*t25*t9 975 t28 = 1/t7**3.166666666666667d+0 976 t29 = 1/t12**2 977 t30 = 7.242189178907065d-4*t14*t27*t29*t3+4.224610354362455d 978 1 -4*rhoa*t13*t28*t3 979 t31 = 3.1748021039363994d+0*gammaaa*t12*t17*t30*wght+3.17480 980 1 21039363994d+0*gammaaa*t17*t27*wght 981 t32 = 1/t7**1.5000000000000002d+0 982 t33 = 9.921256574801247d-2 983 t34 = 1/rhoa**2.3333333333333334d+0 984 t35 = 4.444444444444445d+0*gammaaa*t12*t17*t18*t33*t34*wght 985 t36 = rhoa**3 986 t37 = 1/t36 987 t38 = t30**2 988 t39 = t24**2 989 t40 = 1/t10**3 990 t41 = 1.5749013123685918d-1 991 t42 = 1/rhoa**2.6666666666666666d+0 992 t43 = 2.4499914052208632d-3*t42*t8+2.01824436995453d-1*t41*t 993 1 42+1.9599931241766905d-3*t18*t20+2.4050299253846513d+0*t1 994 2 5*t16*t18 995 t44 = 3.1594931461004283d-6*t41*t42+6.414699787229083d-3*t15 996 1 *t16*t18 997 t45 = -t25*t43*t9+2*t39*t40*t9+t11*t44-2*t24*t25*t26 998 t46 = 1/t7**5.166666666666667d+0 999 t47 = t27**2 1000 t48 = 1/t12**3 1001 t49 = -1.448437835781413d-3*t14*t3*t47*t48-6.688966394407221 1002 1 d-4*t13*t19*t3*t46+7.242189178907065d-4*t14*t29*t3*t45-8. 1003 2 44922070872491d-4*rhoa*t27*t28*t29*t3+2.1123051771812276d 1004 3 -4*t13*t28*t3 1005 t50 = 6.299605249474365d-1*t2*t6*t8*(3.1748021039363994d+0*g 1006 1 ammaaa*t12*t17*t49*wght+3.1748021039363994d+0*gammaaa*t17 1007 2 *t45*wght+3.1748021039363994d+0*gammaaa*t12*t17*t38*wght+ 1008 3 6.349604207872799d+0*gammaaa*t17*t27*t30*wght) 1009 t51 = -6.349604207872799d+0*gammaaa*t12*t17*t30*wght-6.34960 1010 1 4207872799d+0*gammaaa*t17*t27*wght 1011 t52 = 3.1498026247371824d-1*t20*t51*t6*t8 1012 t53 = -3.1748021039363994d+0*gammaaa*t12*t17*t30*wght-3.1748 1013 1 021039363994d+0*gammaaa*t17*t27*wght 1014 t54 = 4.19973683298291d-1*t15*t16*t18*t53 1015 t55 = -1.0582673679787998d+0*gammaaa*t12*t17*t30*wght-1.0582 1016 1 673679787998d+0*gammaaa*t17*t27*wght 1017 t56 = 1.259921049894873d+0*t15*t16*t18*t55 1018 t57 = 1/t7**2.166666666666667d+0 1019 t58 = 1/t3 1020 t59 = 9.052736473633831d-5*t14*t27*t29*t58+5.280762942953069 1021 1 d-5*rhoa*t13*t28*t58 1022 t60 = 6.299605249474365d-1*t2*t6*t8*(3.1748021039363994d+0*g 1023 1 ammaaa*t12*t17*t59*wght-2.8740646802874464d-4*t14*t17*t3* 1024 2 t30*wght+7.937005259840999d-1*t12*t17*t30*wght-2.87406468 1025 3 02874464d-4*t13*t14*t17*t27*t3*wght+7.937005259840999d-1* 1026 4 t17*t27*wght)-2.5d-1*t12*t17*t20*t6*t8*wght+6.03515764908 1027 5 9219d-5*t1*t16*t17*t3*t57*wght+7.185161700718616d-5*t17*t 1028 6 20*t3*t32*wght-3.333333333333333d-1*t12*t15*t16*t17*t18*w 1029 7 ght 1030 t61 = 1.8105472947267662d-4*t14*t27*t29*t58+1.05615258859061 1031 1 38d-4*rhoa*t13*t28*t58 1032 t62 = 1/t7**2.666666666666667d+0 1033 t63 = 1.300907507941046d-8*t13*t17*t2*t62*wght-5.38887127553 1034 1 8962d-5*t17*t2*t32*t58*wght 1035 t64 = 1/rhoa**3.3333333333333337d+0 1036 t65 = -8.671913154270722d-1*gammaaa*t12*t17*t18*t64*wght 1037 t66 = 1/rhoa**4 1038 t67 = 3.1498026247371824d-1*t20*t6*t8*(-9.524406311809198d+0 1039 1 *gammaaa*t12*t17*t49*wght-9.524406311809198d+0*gammaaa*t1 1040 2 7*t45*wght-9.524406311809198d+0*gammaaa*t12*t17*t38*wght- 1041 3 1.9048812623618397d+1*gammaaa*t17*t27*t30*wght) 1042 t68 = 4.19973683298291d-1*t15*t16*t18*(-6.349604207872799d+0 1043 1 *gammaaa*t12*t17*t49*wght-6.349604207872799d+0*gammaaa*t1 1044 2 7*t45*wght-6.349604207872799d+0*gammaaa*t12*t17*t38*wght- 1045 3 1.26992084157456d+1*gammaaa*t17*t27*t30*wght) 1046 t69 = 1.259921049894873d+0*t15*t16*t18*(-1.0582673679787998d 1047 1 +0*gammaaa*t12*t17*t49*wght-1.0582673679787998d+0*gammaaa 1048 2 *t17*t45*wght-1.0582673679787998d+0*gammaaa*t12*t17*t38*w 1049 3 ght-2.1165347359575995d+0*gammaaa*t17*t27*t30*wght) 1050 t70 = -9.799385943626789d-1*t18*t33*t34*t53 1051 t71 = 1.259921049894873d+0*t18*t33*t34*(3.527557893262666d+0 1052 1 *gammaaa*t12*t17*t30*wght+3.527557893262666d+0*gammaaa*t1 1053 2 7*t27*wght) 1054 t72 = 7.874506561842959d-2 1055 t73 = 1/rhoa**3.6666666666666664d+0 1056 t74 = -t25*(-3.266655206961151d-3*t73*t8-5.381984986545413d- 1057 1 1*t72*t73-3.3483215871351796d-3*t18*t37-5.611736492564187 1058 2 d+0*t18*t33*t34)*t9+6*t24*t40*t43*t9-6*t24**3*t9/t10**4+t 1059 3 11*(-8.425315056267809d-6*t72*t73-1.496763283686786d-2*t1 1060 4 8*t33*t34)-3*t24*t25*t44-3*t25*t26*t43+6*t26*t39*t40 1061 t75 = 6.299605249474365d-1*t2*t6*t8*(3.1748021039363994d+0*g 1062 1 ammaaa*t12*t17*(7.242189178907065d-4*t14*t29*t3*t74+1.727 1063 2 982985221866d-3*t13*t3*t36/t7**7.166666666666667d+0+2.534 1064 3 766212617473d-3*rhoa*t28*t3*t47*t48-4.3453135073442384d-3 1065 4 *t14*t27*t3*t45*t48+2.006689918322166d-3*t19*t27*t29*t3*t 1066 5 46-1.003344959161083d-3*rhoa*t13*t3*t46-1.267383106308736 1067 6 5d-3*rhoa*t28*t29*t3*t45-6.336915531543683d-4*t27*t28*t29 1068 7 *t3+4.3453135073442384d-3*t14*t27**3*t3/t12**4)*wght+3.17 1069 8 48021039363994d+0*gammaaa*t17*t74*wght+9.524406311809198d 1070 9 +0*gammaaa*t12*t17*t30*t49*wght+9.524406311809198d+0*gamm 1071 : aaa*t17*t27*t49*wght+9.524406311809198d+0*gammaaa*t17*t30 1072 ; *t45*wght+9.524406311809198d+0*gammaaa*t17*t27*t38*wght+3 1073 < .1748021039363994d+0*gammaaa*t12*t17*t30**3*wght) 1074 t76 = 6.349604207872799d+0*gammaaa*t17*t27*wght 1075 t77 = 6.349604207872799d+0*gammaaa*t12*t17*t30*wght 1076 t78 = 1.5749013123685912d-1*t37*t6*(t77+t76-2*t51)*t8 1077 t79 = 4.19973683298291d-1*t18*t33*t34*(t77+t76-7*t55) 1078 t80 = -2.011719216363073d-4*t15*t17*t3*t34*t57*wght 1079 t81 = -7.185161700718616d-5*t17*t3*t32*t37*wght 1080 t82 = 1.1111111111111112d+0*t12*t17*t18*t33*t34*wght 1081 t83 = 2.5d-1*t12*t17*t37*t6*t8*wght 1082 t84 = -1.5874010519681997d+0*t17*t27*wght 1083 t85 = 5.748129360574893d-4*t13*t14*t17*t27*t3*wght 1084 t86 = 5.748129360574893d-4*t14*t17*t3*t30*wght 1085 t87 = -1.5874010519681997d+0*t12*t17*t30*wght 1086 t88 = 3.1498026247371824d-1*t20*t6*t8*(-6.349604207872799d+0 1087 1 *gammaaa*t12*t17*t59*wght+t87+t86+t85+t84) 1088 t89 = 4.19973683298291d-1*t15*t16*t18*(-3.1748021039363994d+ 1089 1 0*gammaaa*t12*t17*t59*wght+2.8740646802874464d-4*t14*t17* 1090 2 t3*t30*wght-7.937005259840999d-1*t12*t17*t30*wght+2.87406 1091 3 46802874464d-4*t13*t14*t17*t27*t3*wght-7.937005259840999d 1092 4 -1*t17*t27*wght) 1093 t90 = 1.259921049894873d+0*t15*t16*t18*(-1.0582673679787998d 1094 1 +0*gammaaa*t12*t17*t59*wght+9.580215600958156d-5*t14*t17* 1095 2 t3*t30*wght-2.6456684199469993d-1*t12*t17*t30*wght+9.5802 1096 3 15600958156d-5*t13*t14*t17*t27*t3*wght-2.6456684199469993 1097 4 d-1*t17*t27*wght) 1098 t91 = 6.299605249474365d-1*t2*t6*t8*(6.349604207872799d+0*ga 1099 1 mmaaa*t12*t17*t30*t59*wght+6.349604207872799d+0*gammaaa*t 1100 2 17*t27*t59*wght+3.1748021039363994d+0*gammaaa*t12*t17*(-1 1101 3 .8105472947267662d-4*t14*t47*t48*t58-8.361207993009026d-5 1102 4 *t13*t19*t46*t58+9.052736473633831d-5*t14*t29*t45*t58-1.0 1103 5 561525885906138d-4*rhoa*t27*t28*t29*t58+2.640381471476534 1104 6 5d-5*t13*t28*t58)*wght-2.8740646802874464d-4*t14*t17*t3*t 1105 7 49*wght+7.937005259840999d-1*t12*t17*t49*wght-2.874064680 1106 8 2874464d-4*t13*t14*t17*t3*t45*wght+7.937005259840999d-1*t 1107 9 17*t45*wght-2.8740646802874464d-4*t14*t17*t3*t38*wght+7.9 1108 : 37005259840999d-1*t12*t17*t38*wght-5.748129360574893d-4*t 1109 ; 13*t14*t17*t27*t3*t30*wght+1.5874010519681997d+0*t17*t27* 1110 < t30*wght) 1111 t92 = 9.999999999999998d-1*t2*t6*t8*(2.5146490204538424d-5*t 1112 1 14*t17*t20*t3*wght-6.944444444444446d-2*t12*t17*t20*wght) 1113 2 +t91+t90+t89+t88+t83+t82+t81+t80 1114 t93 = -4.023438432726146d-4*t15*t17*t3*t34*t57*wght 1115 t94 = -1.4370323401437232d-4*t17*t3*t32*t37*wght 1116 t95 = 2.2222222222222224d+0*t12*t17*t18*t33*t34*wght 1117 t96 = 5.0d-1*t12*t17*t37*t6*t8*wght 1118 t97 = 3.1498026247371824d-1*t20*t6*t8*(-6.349604207872799d+0 1119 1 *gammaaa*t12*t17*t61*wght+1.1496258721149785d-3*t14*t17*t 1120 2 3*t30*wght-3.1748021039363994d+0*t12*t17*t30*wght+1.14962 1121 3 58721149785d-3*t13*t14*t17*t27*t3*wght-3.1748021039363994 1122 4 d+0*t17*t27*wght) 1123 t98 = 4.19973683298291d-1*t15*t16*t18*(-3.1748021039363994d+ 1124 1 0*gammaaa*t12*t17*t61*wght+t87+t86+t85+t84) 1125 t99 = 1.259921049894873d+0*t15*t16*t18*(-1.0582673679787998d 1126 1 +0*gammaaa*t12*t17*t61*wght+1.9160431201916311d-4*t14*t17 1127 2 *t3*t30*wght-5.291336839893999d-1*t12*t17*t30*wght+1.9160 1128 3 431201916311d-4*t13*t14*t17*t27*t3*wght-5.291336839893999 1129 4 d-1*t17*t27*wght) 1130 t100 = 6.299605249474365d-1*t2*t6*t8*(6.349604207872799d+0*g 1131 1 ammaaa*t12*t17*t30*t61*wght+6.349604207872799d+0*gammaaa* 1132 2 t17*t27*t61*wght+3.1748021039363994d+0*gammaaa*t12*t17*(- 1133 3 3.6210945894535324d-4*t14*t47*t48*t58-1.6722415986018052d 1134 4 -4*t13*t19*t46*t58+1.8105472947267662d-4*t14*t29*t45*t58- 1135 5 2.1123051771812276d-4*rhoa*t27*t28*t29*t58+5.280762942953 1136 6 069d-5*t13*t28*t58)*wght-5.748129360574893d-4*t14*t17*t3* 1137 7 t49*wght+1.5874010519681997d+0*t12*t17*t49*wght-5.7481293 1138 8 60574893d-4*t13*t14*t17*t3*t45*wght+1.5874010519681997d+0 1139 9 *t17*t45*wght-5.748129360574893d-4*t14*t17*t3*t38*wght+1. 1140 : 5874010519681997d+0*t12*t17*t38*wght-1.1496258721149785d- 1141 ; 3*t13*t14*t17*t27*t3*t30*wght+3.1748021039363994d+0*t17*t 1142 < 27*t30*wght) 1143 t101 = 1/t7**3.3333333333333337d+0 1144 t102 = 1/t7**2.333333333333334d+0 1145 t103 = -5.748129360574893d-4*t14*t17*t3*t59*wght 1146 t104 = 1/t3**3 1147 t105 = 1.5874010519681997d+0*t12*t17*t59*wght 1148 t106 = 6.299605249474365d-1*t2*t6*t8*(-1.0777742551077925d-4 1149 1 *t14*t17*t30*t58*wght-1.0777742551077925d-4*t13*t14*t17*t 1150 2 27*t58*wght+2.601815015882092d-8*t102*t13*t17*t30*wght+3. 1151 3 1748021039363994d+0*gammaaa*t12*t17*(-1.1315920592042289d 1152 4 -5*t104*t14*t27*t29-6.600953678691336d-6*rhoa*t104*t13*t2 1153 5 8)*wght+2.601815015882092d-8*t102*t17*t27*t29*wght+t105+t 1154 6 103)-6.504537539705231d-9*t13*t17*t20*t62*wght+2.26318411 1155 7 84084575d-5*t1*t16*t17*t57*t58*wght+2.694435637769481d-5* 1156 8 t17*t20*t32*t58*wght-5.463469177404018d-9*t1*t101*t13*t16 1157 9 *t17*wght 1158 t107 = 6.299605249474365d-1*t2*t6*t8*(-2.8740646802874464d-4 1159 1 *t14*t17*t3*t61*wght+7.937005259840999d-1*t12*t17*t61*wgh 1160 2 t-2.155548510215585d-4*t14*t17*t30*t58*wght-2.15554851021 1161 3 5585d-4*t13*t14*t17*t27*t58*wght+5.203630031764184d-8*t10 1162 4 2*t13*t17*t30*wght+3.1748021039363994d+0*gammaaa*t12*t17* 1163 5 (-2.2631841184084578d-5*t104*t14*t27*t29-1.32019073573826 1164 6 72d-5*rhoa*t104*t13*t28)*wght+5.203630031764184d-8*t102*t 1165 7 17*t27*t29*wght+t105+t103)-1.300907507941046d-8*t13*t17*t 1166 8 20*t62*wght+4.526368236816915d-5*t1*t16*t17*t57*t58*wght+ 1167 9 5.388871275538962d-5*t17*t20*t32*t58*wght-1.0926938354808 1168 : 036d-8*t1*t101*t13*t16*t17*wght 1169 t108 = 1/t7**3.833333333333334d+0 1170 t109 = 1/gammaaa 1171 t110 = 4.878403154778923d-9*t109*t13*t17*t2*t62*wght-1.17767 1172 1 72845961996d-12*t108*t17*t2*t29*t58*wght+6.73608909442370 1173 2 2d-6*t104*t17*t2*t32*wght 1174 t111 = 9.756806309557846d-9*t109*t13*t17*t2*t62*wght-2.35535 1175 1 4569192399d-12*t108*t17*t2*t29*t58*wght+1.347217818884740 1176 2 6d-5*t104*t17*t2*t32*wght 1177 fnc(iq) = 4.000000000000001d+0*gammaaa*t1*t12*t2*exp(-7.2421 1178 1 89178907064d-4*t13*t14*t3)*t8*wght+fnc(iq) 1179 Amat(iq,D1_RA) = -1.0d+0*gammaaa*t12*t17*t20*t6*t8*wght-1.33 1180 1 33333333333333d+0*gammaaa*t12*t15*t16*t17*t18*wght+6.2996 1181 2 05249474365d-1*t2*t31*t6*t8+Amat(iq,D1_RA) 1182 Cmat(iq,D1_GAA) = 5.0d-1*t12*t17*t2*t6*t8*wght-1.43703234014 1183 1 37232d-4*t17*t2*t3*t32*wght+Cmat(iq,D1_GAA) 1184 Cmat(iq,D1_GAB) = 1.0d+0*t12*t17*t2*t6*t8*wght-2.87406468028 1185 1 74464d-4*t17*t2*t3*t32*wght+Cmat(iq,D1_GAB) 1186 Amat2(iq,D2_RA_RA) = 7.222222222222222d-1*gammaaa*t12*t17*t3 1187 1 7*t6*t8*wght+t56+t54+t52+t50+t35+Amat2(iq,D2_RA_RA) 1188 Amat2(iq,D2_RA_RB) = 1.2777777777777777d+0*gammaaa*t12*t17*t 1189 1 37*t6*t8*wght+t56+t54+t52+t50+t35+Amat2(iq,D2_RA_RB) 1190 Cmat2(iq,D2_RA_GAA) = t60+Cmat2(iq,D2_RA_GAA) 1191 Cmat2(iq,D2_RA_GAB) = 6.299605249474365d-1*t2*t6*t8*(3.17480 1192 1 21039363994d+0*gammaaa*t12*t17*t61*wght-5.748129360574893 1193 2 d-4*t14*t17*t3*t30*wght+1.5874010519681997d+0*t12*t17*t30 1194 3 *wght-5.748129360574893d-4*t13*t14*t17*t27*t3*wght+1.5874 1195 4 010519681997d+0*t17*t27*wght)-5.0d-1*t12*t17*t20*t6*t8*wg 1196 5 ht+1.2070315298178438d-4*t1*t16*t17*t3*t57*wght+1.4370323 1197 6 401437232d-4*t17*t20*t3*t32*wght-6.666666666666666d-1*t12 1198 7 *t15*t16*t17*t18*wght+Cmat2(iq,D2_RA_GAB) 1199 Cmat2(iq,D2_RA_GBB) = t60+Cmat2(iq,D2_RA_GBB) 1200 Cmat2(iq,D2_GAA_GAA) = t63+Cmat2(iq,D2_GAA_GAA) 1201 Cmat2(iq,D2_GAA_GAB) = 2.601815015882092d-8*t13*t17*t2*t62*w 1202 1 ght-1.0777742551077923d-4*t17*t2*t32*t58*wght+Cmat2(iq,D2 1203 2 _GAA_GAB) 1204 Cmat2(iq,D2_GAA_GBB) = t63+Cmat2(iq,D2_GAA_GBB) 1205 Cmat2(iq,D2_GAB_GAB) = 5.203630031764184d-8*t13*t17*t2*t62*w 1206 1 ght-2.1555485102155847d-4*t17*t2*t32*t58*wght+Cmat2(iq,D2 1207 2 _GAB_GAB) 1208 Amat3(iq,D3_RA_RA_RA) = 9.999999999999998d-1*t2*t6*t8*(8.333 1209 1 333333333336d-1*gammaaa*t12*t17*t37*wght-5.55555555555555 1210 2 7d-1*gammaaa*t12*t17*t20*t30*wght-5.555555555555557d-1*ga 1211 3 mmaaa*t17*t20*t27*wght)-1.0833333333333335d+0*gammaaa*t12 1212 4 *t17*t6*t66*t8*wght+5.555555555555556d-1*gammaaa*t12*t15* 1213 5 t17*t18*t64*wght-8.749451735381063d-2*t31*t37*t6*t8+t79+t 1214 6 78+t75+t71+t70+t69+t68+t67+t65+Amat3(iq,D3_RA_RA_RA) 1215 Amat3(iq,D3_RA_RA_RB) = -1.9166666666666665d+0*gammaaa*t12*t 1216 1 17*t6*t66*t8*wght-1.8518518518518515d-1*gammaaa*t12*t15*t 1217 2 17*t18*t64*wght+8.749451735381063d-2*t31*t37*t6*t8+t79+t7 1218 3 8+t75+t71+t70+t69+t68+t67+t65+Amat3(iq,D3_RA_RA_RB) 1219 Cmat3(iq,D3_RA_RA_GAA) = t92+Cmat3(iq,D3_RA_RA_GAA) 1220 Cmat3(iq,D3_RA_RA_GAB) = 9.999999999999998d-1*t2*t6*t8*(5.02 1221 1 9298040907685d-5*t14*t17*t20*t3*wght-1.3888888888888892d- 1222 2 1*t12*t17*t20*wght)+t99+t98+t97+t96+t95+t94+t93+t100+Cmat 1223 3 3(iq,D3_RA_RA_GAB) 1224 Cmat3(iq,D3_RA_RA_GBB) = t92+Cmat3(iq,D3_RA_RA_GBB) 1225 Cmat3(iq,D3_RA_RB_GAA) = 9.999999999999998d-1*t2*t6*t8*(6.94 1226 1 4444444444446d-2*t12*t17*t20*wght-2.5146490204538424d-5*t 1227 2 14*t17*t20*t3*wght)+t91+t90+t89+t88+t83+t82+t81+t80+Cmat3 1228 3 (iq,D3_RA_RB_GAA) 1229 Cmat3(iq,D3_RA_RB_GAB) = 9.999999999999998d-1*t2*t6*t8*(1.38 1230 1 88888888888892d-1*t12*t17*t20*wght-5.029298040907685d-5*t 1231 2 14*t17*t20*t3*wght)+t99+t98+t97+t96+t95+t94+t93+t100+Cmat 1232 3 3(iq,D3_RA_RB_GAB) 1233 Cmat3(iq,D3_RA_GAA_GAA) = t106+Cmat3(iq,D3_RA_GAA_GAA) 1234 Cmat3(iq,D3_RA_GAA_GAB) = t107+Cmat3(iq,D3_RA_GAA_GAB) 1235 Cmat3(iq,D3_RA_GAA_GBB) = t106+Cmat3(iq,D3_RA_GAA_GBB) 1236 Cmat3(iq,D3_RA_GAB_GAB) = 6.299605249474365d-1*t2*t6*t8*(-1. 1237 1 1496258721149785d-3*t14*t17*t3*t61*wght+3.174802103936399 1238 2 4d+0*t12*t17*t61*wght-4.31109702043117d-4*t14*t17*t30*t58 1239 3 *wght-4.31109702043117d-4*t13*t14*t17*t27*t58*wght+1.0407 1240 4 260063528369d-7*t102*t13*t17*t30*wght+3.1748021039363994d 1241 5 +0*gammaaa*t12*t17*(-4.5263682368169156d-5*t104*t14*t27*t 1242 6 29-2.6403814714765345d-5*rhoa*t104*t13*t28)*wght+1.040726 1243 7 0063528369d-7*t102*t17*t27*t29*wght)-2.601815015882092d-8 1244 8 *t13*t17*t20*t62*wght+9.05273647363383d-5*t1*t16*t17*t57* 1245 9 t58*wght+1.0777742551077923d-4*t17*t20*t32*t58*wght-2.185 1246 : 387670961607d-8*t1*t101*t13*t16*t17*wght+Cmat3(iq,D3_RA_G 1247 ; AB_GAB) 1248 Cmat3(iq,D3_RA_GAB_GBB) = t107+Cmat3(iq,D3_RA_GAB_GBB) 1249 Cmat3(iq,D3_RA_GBB_GBB) = t106+Cmat3(iq,D3_RA_GBB_GBB) 1250 Cmat3(iq,D3_GAA_GAA_GAA) = t110+Cmat3(iq,D3_GAA_GAA_GAA) 1251 Cmat3(iq,D3_GAA_GAA_GAB) = t111+Cmat3(iq,D3_GAA_GAA_GAB) 1252 Cmat3(iq,D3_GAA_GAA_GBB) = t110+Cmat3(iq,D3_GAA_GAA_GBB) 1253 Cmat3(iq,D3_GAA_GAB_GAB) = 1.9513612619115692d-8*t109*t13*t1 1254 1 7*t2*t62*wght-4.710709138384798d-12*t108*t17*t2*t29*t58*w 1255 2 ght+2.694435637769481d-5*t104*t17*t2*t32*wght+Cmat3(iq,D3 1256 3 _GAA_GAB_GAB) 1257 Cmat3(iq,D3_GAA_GAB_GBB) = t111+Cmat3(iq,D3_GAA_GAB_GBB) 1258 Cmat3(iq,D3_GAB_GAB_GAB) = 3.9027225238231383d-8*t109*t13*t1 1259 1 7*t2*t62*wght-9.421418276769596d-12*t108*t17*t2*t29*t58*w 1260 2 ght+5.388871275538962d-5*t104*t17*t2*t32*wght+Cmat3(iq,D3 1261 3 _GAB_GAB_GAB) 1262 endif ! rhoa.gt.tol_rho 1263 else ! ipol.eq.1 1264 rhoa = rho(iq,R_A) 1265 rhob = rho(iq,R_B) 1266 gammaaa = rgamma(iq,G_AA) 1267 gammaab = rgamma(iq,G_AB) 1268 gammabb = rgamma(iq,G_BB) 1269 if (rhoa.gt.tol_rho.and.rhob.gt.tol_rho) then 1270 t1 = gammabb+2.0d+0*gammaab+gammaaa 1271 t2 = rhob+rhoa 1272 t3 = 1/t2 1273 t4 = rhoa-rhob 1274 t5 = 1.0d+0-1.0d+0*t3*t4 1275 t6 = 1.0d+0*t3*t4+1.0d+0 1276 t7 = (3.149802624737183d-1*t6**1.6666666666666669d+0+3.14980 1277 1 2624737183d-1*t5**1.6666666666666669d+0)**5.0d-1 1278 t8 = 1/t7 1279 t9 = t1**5.0d-1 1280 t10 = 1/t2**6.666666666666666d-1 1281 t11 = abs(t2) 1282 t12 = 1/t11**3.333333333333333d-1 1283 t13 = 1.443307452126544d-2*t12+2.8435438314903855d-6*t10+2.5 1284 1 68d-3 1285 t14 = 1.7639938117590215d-2*t10*t12+5.411317332115466d+0*t12 1286 1 +1.816419932959077d-1*t10+1.0d+0 1287 t15 = 1/t14 1288 t16 = t13*t15+1.6670000000000001d-3 1289 t17 = 1/t16 1290 t18 = 1/t11**1.166666666666667d+0 1291 t19 = exp(-8.129082500000001d-4*t17*t18*t9) 1292 t20 = 1/t2**1.3333333333333333d+0 1293 t21 = 1/t11 1294 t22 = -2.6456684199469993d-1*t1*t16*t19*t20*t21*t8*wght 1295 t23 = t2**2 1296 t24 = 1/t23 1297 t25 = 1.0d+0*t24*t4 1298 t26 = -1.0d+0*t3 1299 t27 = t26+t25 1300 t28 = t5**6.666666666666666d-1 1301 t29 = -1.0d+0*t24*t4 1302 t30 = 1.0d+0*t3 1303 t31 = t30+t29 1304 t32 = t6**6.666666666666666d-1 1305 t33 = 5.249671041228638d-1*t31*t32+5.249671041228638d-1*t27* 1306 1 t28 1307 t34 = 1/t7**3 1308 t35 = -7.937005259840999d-1*t1*t12*t16*t19*t24*t8*wght 1309 t36 = 1/t2**1.6666666666666669d+0 1310 t37 = 1/t2**3.333333333333333d-1 1311 t38 = -1.8037724440384886d+0*t21*t37-1.1759958745060142d-2*t 1312 1 12*t36-1.2109466219727179d-1*t36-5.879979372530072d-3*t21 1313 2 *t3 1314 t39 = 1/t14**2 1315 t40 = -4.8110248404218126d-3*t21*t37-1.895695887660257d-6*t3 1316 1 6 1317 t41 = t15*t40-t13*t38*t39 1318 t42 = 1/t11**3.166666666666667d+0 1319 t43 = 1/t16**2 1320 t44 = 8.129082500000001d-4*t18*t41*t43*t9+9.483929583333336d 1321 1 -4*t17*t2*t42*t9 1322 t45 = 7.937005259840999d-1*t1*t16*t19*t44*wght+7.93700525984 1323 1 0999d-1*t1*t19*t41*wght 1324 t46 = t12*t3*t45*t8 1325 t47 = t30+t25 1326 t48 = t29+t26 1327 t49 = 5.249671041228638d-1*t32*t48+5.249671041228638d-1*t28* 1328 1 t47 1329 t50 = 1/t11**1.5000000000000002d+0 1330 t51 = 7.937005259840999d-1*t12*t16*t19*t3*t8*wght-3.22602852 1331 1 8009071d-4*t19*t3*t50*t8*t9*wght 1332 t52 = 1/t2**2.3333333333333334d+0 1333 t53 = 8.818894733156665d-1*t1*t16*t19*t21*t52*t8*wght 1334 t54 = t33**2 1335 t55 = 1/t7**5 1336 t56 = t2**3 1337 t57 = 1/t56 1338 t58 = 1.5874010519681997d+0*t1*t12*t16*t19*t57*t8*wght 1339 t59 = t44**2 1340 t60 = t38**2 1341 t61 = 1/t14**3 1342 t62 = 1/t2**2.6666666666666666d+0 1343 t63 = 1.9599931241766905d-2*t12*t62+2.01824436995453d-1*t62+ 1344 1 1.5679944993413522d-2*t21*t24+2.4050299253846513d+0*t20*t 1345 2 21 1346 t64 = 3.1594931461004283d-6*t62+6.414699787229083d-3*t20*t21 1347 t65 = t15*t64-t13*t39*t63+2*t13*t60*t61-2*t38*t39*t40 1348 t66 = 1/t11**5.166666666666667d+0 1349 t67 = t41**2 1350 t68 = 1/t16**3 1351 t69 = -1.6258165000000002d-3*t18*t67*t68*t9-3.00324436805555 1352 1 7d-3*t17*t23*t66*t9+8.129082500000001d-4*t18*t43*t65*t9-1 1353 2 .8967859166666673d-3*t2*t41*t42*t43*t9+9.483929583333336d 1354 3 -4*t17*t42*t9 1355 t70 = 7.937005259840999d-1*t1*t16*t19*t69*wght+7.93700525984 1356 1 0999d-1*t1*t19*t65*wght+7.937005259840999d-1*t1*t16*t19*t 1357 2 59*wght+1.5874010519681997d+0*t1*t19*t41*t44*wght 1358 t71 = t12*t3*t70*t8 1359 t72 = -1.5874010519681997d+0*t1*t16*t19*t44*wght-1.587401051 1360 1 9681997d+0*t1*t19*t41*wght 1361 t73 = t12*t24*t72*t8 1362 t74 = -7.937005259840999d-1*t1*t16*t19*t44*wght-7.9370052598 1363 1 40999d-1*t1*t19*t41*wght 1364 t75 = 3.333333333333333d-1*t20*t21*t74*t8 1365 t76 = -2.6456684199469993d-1*t1*t16*t19*t44*wght-2.645668419 1366 1 9469993d-1*t1*t19*t41*wght 1367 t77 = t20*t21*t76*t8 1368 t78 = -3.9685026299204995d-1*t1*t19*t33*t41*wght 1369 t79 = t27**2 1370 t80 = 1/t5**3.333333333333333d-1 1371 t81 = -2.0d+0*t4*t57 1372 t82 = 2.0d+0*t24 1373 t83 = t82+t81 1374 t84 = t31**2 1375 t85 = 1/t6**3.333333333333333d-1 1376 t86 = 2.0d+0*t4*t57 1377 t87 = -2.0d+0*t24 1378 t88 = t87+t86 1379 t89 = 5.249671041228638d-1*t32*t88+3.4997806941524257d-1*t84 1380 1 *t85+5.249671041228638d-1*t28*t83+3.4997806941524257d-1*t 1381 2 79*t80 1382 t90 = -3.9685026299204995d-1*t1*t16*t19*t33*t44*wght 1383 t91 = -3.9685026299204995d-1*t1*t16*t19*t89*wght+t90+t78 1384 t92 = 1.3228342099734997d-1*t1*t16*t19*t49*wght+1.3228342099 1385 1 734997d-1*t1*t16*t19*t33*wght 1386 t93 = 3.9685026299204995d-1*t1*t16*t19*t49*wght+3.9685026299 1387 1 204995d-1*t1*t16*t19*t33*wght 1388 t94 = -5.0d-1*t12*t3*t34*t45*t49 1389 t95 = 3.4997806941524257d-1*t31*t48*t85+3.4997806941524257d- 1390 1 1*t27*t47*t80+1.0499342082457278d+0*t32*t4*t57-1.04993420 1391 2 82457278d+0*t28*t4*t57 1392 t96 = -3.9685026299204995d-1*t1*t16*t19*t95*wght+t90+t78 1393 t97 = t49**2 1394 t98 = t47**2 1395 t99 = t87+t81 1396 t100 = t48**2 1397 t101 = t86+t82 1398 t102 = 5.249671041228638d-1*t28*t99+3.4997806941524257d-1*t8 1399 1 0*t98+3.4997806941524257d-1*t100*t85+5.249671041228638d-1 1400 2 *t101*t32 1401 t103 = -3.9685026299204995d-1*t1*t16*t19*t44*t49*wght-3.9685 1402 1 026299204995d-1*t1*t19*t41*t49*wght-3.9685026299204995d-1 1403 2 *t1*t102*t16*t19*wght 1404 t104 = 1/t11**2.166666666666667d+0 1405 t105 = 1.0753428426696904d-4*t104*t19*t20*t8*t9*wght 1406 t106 = 3.226028528009071d-4*t19*t24*t50*t8*t9*wght 1407 t107 = -2.6456684199469993d-1*t16*t19*t20*t21*t8*wght 1408 t108 = -7.937005259840999d-1*t12*t16*t19*t24*t8*wght 1409 t109 = 1/t9 1410 t110 = 4.0645412500000005d-4*t109*t18*t41*t43+4.741964791666 1411 1 668d-4*t109*t17*t2*t42 1412 t111 = -3.226028528009071d-4*t18*t19*t44*t9*wght-3.226028528 1413 1 009071d-4*t17*t18*t19*t41*t9*wght+7.937005259840999d-1*t1 1414 2 6*t19*t44*wght+7.937005259840999d-1*t19*t41*wght+7.937005 1415 3 259840999d-1*t1*t110*t16*t19*wght 1416 t112 = t111*t12*t3*t8 1417 t113 = 1.6130142640045356d-4*t19*t3*t33*t34*t50*t9*wght-3.96 1418 1 85026299204995d-1*t12*t16*t19*t3*t33*t34*wght+t112+t108+t 1419 2 107+t106+t105 1420 t114 = 2.1506856853393808d-4*t104*t19*t20*t8*t9*wght 1421 t115 = 6.452057056018142d-4*t19*t24*t50*t8*t9*wght 1422 t116 = -5.291336839893999d-1*t16*t19*t20*t21*t8*wght 1423 t117 = -1.5874010519681997d+0*t12*t16*t19*t24*t8*wght 1424 t118 = 8.129082500000001d-4*t109*t18*t41*t43+9.4839295833333 1425 1 36d-4*t109*t17*t2*t42 1426 t119 = -6.452057056018142d-4*t18*t19*t44*t9*wght-6.452057056 1427 1 018142d-4*t17*t18*t19*t41*t9*wght+1.5874010519681997d+0*t 1428 2 16*t19*t44*wght+1.5874010519681997d+0*t19*t41*wght+7.9370 1429 3 05259840999d-1*t1*t118*t16*t19*wght 1430 t120 = t119*t12*t3*t8 1431 t121 = 1.6130142640045356d-4*t19*t3*t34*t49*t50*t9*wght-3.96 1432 1 85026299204995d-1*t12*t16*t19*t3*t34*t49*wght+t112+t108+t 1433 2 107+t106+t105 1434 t122 = 1/t11**2.666666666666667d+0 1435 t123 = 1.311232602576965d-7*t122*t17*t19*t3*t8*wght-4.839042 1436 1 7920136064d-4*t109*t19*t3*t50*t8*wght 1437 t124 = 2.62246520515393d-7*t122*t17*t19*t3*t8*wght-9.6780855 1438 1 84027213d-4*t109*t19*t3*t50*t8*wght 1439 t125 = -3.468765261708288d+0*t1*t16*t19*t21*t8*wght/t2**3.33 1440 1 33333333333337d+0 1441 t126 = 1/t7**7 1442 t127 = 1/t2**4 1443 t128 = -4.762203155904599d+0*t1*t12*t127*t16*t19*t8*wght 1444 t129 = t12*t24*t8*(-2.3811015779522995d+0*t1*t16*t19*t69*wgh 1445 1 t-2.3811015779522995d+0*t1*t19*t65*wght-2.381101577952299 1446 2 5d+0*t1*t16*t19*t59*wght-4.762203155904599d+0*t1*t19*t41* 1447 3 t44*wght) 1448 t130 = 3.333333333333333d-1*t20*t21*t8*(-1.5874010519681997d 1449 1 +0*t1*t16*t19*t69*wght-1.5874010519681997d+0*t1*t19*t65*w 1450 2 ght-1.5874010519681997d+0*t1*t16*t19*t59*wght-3.174802103 1451 3 9363994d+0*t1*t19*t41*t44*wght) 1452 t131 = t20*t21*t8*(-2.6456684199469993d-1*t1*t16*t19*t69*wgh 1453 1 t-2.6456684199469993d-1*t1*t19*t65*wght-2.645668419946999 1454 2 3d-1*t1*t16*t19*t59*wght-5.291336839893999d-1*t1*t19*t41* 1455 3 t44*wght) 1456 t132 = -7.777777777777778d-1*t21*t52*t74*t8 1457 t133 = t21*t52*t8*(8.818894733156665d-1*t1*t16*t19*t44*wght+ 1458 1 8.818894733156665d-1*t1*t19*t41*wght) 1459 t134 = 2.6456684199469993d-1*t1*t19*t33*t41*wght 1460 t135 = 2.6456684199469993d-1*t1*t16*t19*t33*t44*wght 1461 t136 = 3.9685026299204995d-1*t1*t19*t33*t41*wght 1462 t137 = 3.9685026299204995d-1*t1*t16*t19*t89*wght 1463 t138 = 3.9685026299204995d-1*t1*t16*t19*t33*t44*wght 1464 t139 = 3.333333333333333d-1*(t138+t137+t136)*t20*t21*t34 1465 t140 = 1.1905507889761498d+0*t1*t19*t33*t41*wght 1466 t141 = 1.1905507889761498d+0*t1*t16*t19*t33*t44*wght 1467 t142 = 5.952753944880749d-1*t1*t19*t41*t54*wght 1468 t143 = 5.952753944880749d-1*t1*t16*t19*t44*t54*wght 1469 t144 = 1/t2**3.6666666666666664d+0 1470 t145 = -3*t38*t39*t64+6*t13*t38*t61*t63-3*t39*t40*t63+6*t40* 1471 1 t60*t61-t13*t39*(-5.357314539416287d-2*t21*t57-5.61173649 1472 2 2564187d+0*t21*t52-5.226648331137841d-2*t12*t144-5.381984 1473 3 986545413d-1*t144)+t15*(-1.496763283686786d-2*t21*t52-8.4 1474 4 25315056267809d-6*t144)-6*t13*t38**3/t14**4 1475 t146 = t12*t3*t8*(7.937005259840999d-1*t1*t16*t19*(5.6903577 1476 1 50000001d-3*t2*t42*t67*t68*t9-4.8774495000000007d-3*t18*t 1477 2 41*t65*t68*t9+9.009733104166671d-3*t23*t41*t43*t66*t9-9.0 1478 3 09733104166671d-3*t17*t2*t66*t9-2.8451788750000007d-3*t2* 1479 4 t42*t43*t65*t9+1.5516762568287046d-2*t17*t56*t9/t11**7.16 1480 5 6666666666667d+0-2.8451788750000007d-3*t41*t42*t43*t9+8.1 1481 6 29082500000001d-4*t145*t18*t43*t9+4.8774495000000007d-3*t 1482 7 18*t41**3*t9/t16**4)*wght+2.3811015779522995d+0*t1*t16*t1 1483 8 9*t44*t69*wght+2.3811015779522995d+0*t1*t19*t41*t69*wght+ 1484 9 2.3811015779522995d+0*t1*t19*t44*t65*wght+2.3811015779522 1485 : 995d+0*t1*t19*t41*t59*wght+7.937005259840999d-1*t1*t16*t1 1486 ; 9*t44**3*wght+7.937005259840999d-1*t1*t145*t19*wght) 1487 t147 = -3.9685026299204995d-1*t1*t16*t19*t33*t59*wght 1488 t148 = -3.9685026299204995d-1*t1*t19*t33*t65*wght 1489 t149 = 1/t5**1.3333333333333333d+0 1490 t150 = 6.0d+0*t127*t4 1491 t151 = -6.0d+0*t57 1492 t152 = 1/t6**1.3333333333333333d+0 1493 t153 = -6.0d+0*t127*t4 1494 t154 = 6.0d+0*t57 1495 t155 = -3.9685026299204995d-1*t1*t16*t19*t33*t69*wght 1496 t156 = -7.937005259840999d-1*t1*t19*t33*t41*t44*wght 1497 t157 = -t33*t70 1498 t158 = 1.5874010519681997d+0*t1*t19*t41*wght 1499 t159 = 1.5874010519681997d+0*t1*t16*t19*t44*wght 1500 t160 = t12*t57*(-2*t72+t159+t158)*t8 1501 t161 = 3.333333333333333d-1*t21*t52*(-7*t76+t159+t158)*t8 1502 t162 = t33*t45 1503 t163 = -7.937005259840999d-1*t1*t16*t19*t49*wght 1504 t164 = -5.0d-1*t20*t21*t34*t49*t76 1505 t165 = 7.937005259840999d-1*t1*t16*t19*t95*wght 1506 t166 = -t49*t70 1507 t167 = -t49*t72 1508 t168 = -t49*t74 1509 t169 = 7.5d-1*t12*t3*t45*t55*t97 1510 t170 = -5.0d-1*t102*t12*t3*t34*t45 1511 t171 = 3.9685026299204995d-1*t1*t19*t41*t49*wght 1512 t172 = 3.9685026299204995d-1*t1*t102*t16*t19*wght 1513 t173 = 3.9685026299204995d-1*t1*t16*t19*t44*t49*wght 1514 t174 = t45*t49 1515 t175 = 5.0d-1*t12*(t174+t167)*t24*t34 1516 t176 = 1.6666666666666666d-1*(t174+t168)*t20*t21*t34 1517 t177 = -3.5844761422323007d-4*t104*t19*t52*t8*t9*wght 1518 t178 = -6.452057056018142d-4*t19*t50*t57*t8*t9*wght 1519 t179 = 8.818894733156665d-1*t16*t19*t21*t52*t8*wght 1520 t180 = 1.5874010519681997d+0*t12*t16*t19*t57*t8*wght 1521 t181 = -1.5874010519681997d+0*t19*t41*wght 1522 t182 = 6.452057056018142d-4*t17*t18*t19*t41*t9*wght 1523 t183 = 6.452057056018142d-4*t18*t19*t44*t9*wght 1524 t184 = -1.5874010519681997d+0*t16*t19*t44*wght 1525 t185 = t12*t24*t8*(-1.5874010519681997d+0*t1*t110*t16*t19*wg 1526 1 ht+t184+t183+t182+t181) 1527 t186 = 3.333333333333333d-1*t20*t21*t8*(3.226028528009071d-4 1528 1 *t18*t19*t44*t9*wght+3.226028528009071d-4*t17*t18*t19*t41 1529 2 *t9*wght-7.937005259840999d-1*t16*t19*t44*wght-7.93700525 1530 3 9840999d-1*t19*t41*wght-7.937005259840999d-1*t1*t110*t16* 1531 4 t19*wght) 1532 t187 = t20*t21*t8*(1.0753428426696904d-4*t18*t19*t44*t9*wght 1533 1 +1.0753428426696904d-4*t17*t18*t19*t41*t9*wght-2.64566841 1534 2 99469993d-1*t16*t19*t44*wght-2.6456684199469993d-1*t19*t4 1535 3 1*wght-2.6456684199469993d-1*t1*t110*t16*t19*wght) 1536 t188 = -3.9685026299204995d-1*t19*t33*t41*wght 1537 t189 = 1.6130142640045356d-4*t17*t18*t19*t33*t41*t9*wght 1538 t190 = 1.6130142640045356d-4*t18*t19*t33*t44*t9*wght 1539 t191 = -3.9685026299204995d-1*t1*t110*t16*t19*t33*wght 1540 t192 = -3.9685026299204995d-1*t16*t19*t33*t44*wght 1541 t193 = t12*t3*t8*(-3.226028528009071d-4*t18*t19*t69*t9*wght- 1542 1 3.226028528009071d-4*t17*t18*t19*t65*t9*wght-3.2260285280 1543 2 09071d-4*t18*t19*t59*t9*wght-6.452057056018142d-4*t17*t18 1544 3 *t19*t41*t44*t9*wght+7.937005259840999d-1*t16*t19*t69*wgh 1545 4 t+7.937005259840999d-1*t1*t16*t19*(-8.129082500000001d-4* 1546 5 t109*t18*t67*t68-1.5016221840277785d-3*t109*t17*t23*t66+4 1547 6 .0645412500000005d-4*t109*t18*t43*t65-9.483929583333336d- 1548 7 4*t109*t2*t41*t42*t43+4.741964791666668d-4*t109*t17*t42)* 1549 8 wght+7.937005259840999d-1*t19*t65*wght+7.937005259840999d 1550 9 -1*t16*t19*t59*wght+1.5874010519681997d+0*t19*t41*t44*wgh 1551 : t+1.5874010519681997d+0*t1*t110*t16*t19*t44*wght+1.587401 1552 ; 0519681997d+0*t1*t110*t19*t41*wght) 1553 t194 = t12*t3*t34*(1.6130142640045356d-4*t18*t19*t89*t9*wght 1554 1 -3.9685026299204995d-1*t16*t19*t89*wght+t192+t191+t190+t1 1555 2 89+t188)-2.4195213960068032d-4*t19*t3*t50*t54*t55*t9*wght 1556 3 -3.226028528009071d-4*t19*t24*t33*t34*t50*t9*wght-1.07534 1557 4 28426696904d-4*t104*t19*t20*t33*t34*t9*wght+5.95275394488 1558 5 0749d-1*t12*t16*t19*t3*t54*t55*wght+7.937005259840999d-1* 1559 6 t12*t16*t19*t24*t33*t34*wght+2.6456684199469993d-1*t16*t1 1560 7 9*t20*t21*t33*t34*wght-5.0d-1*t111*t12*t3*t33*t34+t193+t1 1561 8 87+t186+t185+t180+t179+t178+t177 1562 t195 = -7.168952284464601d-4*t104*t19*t52*t8*t9*wght 1563 t196 = -1.2904114112036286d-3*t19*t50*t57*t8*t9*wght 1564 t197 = 1.763778946631333d+0*t16*t19*t21*t52*t8*wght 1565 t198 = 3.1748021039363994d+0*t12*t16*t19*t57*t8*wght 1566 t199 = t12*t24*t8*(1.2904114112036286d-3*t18*t19*t44*t9*wght 1567 1 +1.2904114112036286d-3*t17*t18*t19*t41*t9*wght-3.17480210 1568 2 39363994d+0*t16*t19*t44*wght-3.1748021039363994d+0*t19*t4 1569 3 1*wght-1.5874010519681997d+0*t1*t118*t16*t19*wght) 1570 t200 = 3.333333333333333d-1*t20*t21*t8*(-7.937005259840999d- 1571 1 1*t1*t118*t16*t19*wght+t184+t183+t182+t181) 1572 t201 = t20*t21*t8*(2.1506856853393808d-4*t18*t19*t44*t9*wght 1573 1 +2.1506856853393808d-4*t17*t18*t19*t41*t9*wght-5.29133683 1574 2 9893999d-1*t16*t19*t44*wght-5.291336839893999d-1*t19*t41* 1575 3 wght-2.6456684199469993d-1*t1*t118*t16*t19*wght) 1576 t202 = -7.937005259840999d-1*t19*t33*t41*wght 1577 t203 = 3.226028528009071d-4*t17*t18*t19*t33*t41*t9*wght 1578 t204 = 3.226028528009071d-4*t18*t19*t33*t44*t9*wght 1579 t205 = -3.9685026299204995d-1*t1*t118*t16*t19*t33*wght 1580 t206 = -7.937005259840999d-1*t16*t19*t33*t44*wght 1581 t207 = t12*t3*t8*(-6.452057056018142d-4*t18*t19*t69*t9*wght- 1582 1 6.452057056018142d-4*t17*t18*t19*t65*t9*wght-6.4520570560 1583 2 18142d-4*t18*t19*t59*t9*wght-1.2904114112036286d-3*t17*t1 1584 3 8*t19*t41*t44*t9*wght+1.5874010519681997d+0*t16*t19*t69*w 1585 4 ght+7.937005259840999d-1*t1*t16*t19*(-1.6258165000000002d 1586 5 -3*t109*t18*t67*t68-3.003244368055557d-3*t109*t17*t23*t66 1587 6 +8.129082500000001d-4*t109*t18*t43*t65-1.8967859166666673 1588 7 d-3*t109*t2*t41*t42*t43+9.483929583333336d-4*t109*t17*t42 1589 8 )*wght+1.5874010519681997d+0*t19*t65*wght+1.5874010519681 1590 9 997d+0*t16*t19*t59*wght+3.1748021039363994d+0*t19*t41*t44 1591 : *wght+1.5874010519681997d+0*t1*t118*t16*t19*t44*wght+1.58 1592 ; 74010519681997d+0*t1*t118*t19*t41*wght) 1593 t208 = -5.0d-1*t111*t12*t3*t34*t49 1594 t209 = t12*t3*t34*(1.6130142640045356d-4*t18*t19*t9*t95*wght 1595 1 -3.9685026299204995d-1*t16*t19*t95*wght+t192+t191+t190+t1 1596 2 89+t188)+t20*t21*t34*(-5.376714213348451d-5*t18*t19*t49*t 1597 3 9*wght-5.376714213348451d-5*t18*t19*t33*t9*wght+1.3228342 1598 4 099734997d-1*t16*t19*t49*wght+1.3228342099734997d-1*t16*t 1599 5 19*t33*wght)+t12*t24*t34*(-1.6130142640045356d-4*t18*t19* 1600 6 t49*t9*wght-1.6130142640045356d-4*t18*t19*t33*t9*wght+3.9 1601 7 685026299204995d-1*t16*t19*t49*wght+3.9685026299204995d-1 1602 8 *t16*t19*t33*wght)-2.4195213960068032d-4*t19*t3*t33*t49*t 1603 9 50*t55*t9*wght+5.952753944880749d-1*t12*t16*t19*t3*t33*t4 1604 : 9*t55*wght+t208+t193+t187+t186+t185+t180+t179+t178+t177 1605 t210 = -5.0d-1*t119*t12*t3*t34*t49 1606 t211 = t12*t3*t34*(1.6130142640045356d-4*t18*t19*t44*t49*t9* 1607 1 wght+1.6130142640045356d-4*t17*t18*t19*t41*t49*t9*wght+1. 1608 2 6130142640045356d-4*t102*t18*t19*t9*wght-3.96850262992049 1609 3 95d-1*t16*t19*t44*t49*wght-3.9685026299204995d-1*t19*t41* 1610 4 t49*wght-3.9685026299204995d-1*t1*t110*t16*t19*t49*wght-3 1611 5 .9685026299204995d-1*t102*t16*t19*wght)-2.419521396006803 1612 6 2d-4*t19*t3*t50*t55*t9*t97*wght+5.952753944880749d-1*t12* 1613 7 t16*t19*t3*t55*t97*wght-3.226028528009071d-4*t19*t24*t34* 1614 8 t49*t50*t9*wght-1.0753428426696904d-4*t104*t19*t20*t34*t4 1615 9 9*t9*wght+7.937005259840999d-1*t12*t16*t19*t24*t34*t49*wg 1616 : ht+2.6456684199469993d-1*t16*t19*t20*t21*t34*t49*wght+t20 1617 ; 8+t193+t187+t186+t185+t180+t179+t178+t177 1618 t212 = 1/t11**3.3333333333333337d+0 1619 t213 = -4.370775341923217d-8*t17*t19*t20*t212*t8*wght 1620 t214 = -1.311232602576965d-7*t122*t17*t19*t24*t8*wght 1621 t215 = 1.6130142640045356d-4*t104*t109*t19*t20*t8*wght 1622 t216 = 4.8390427920136064d-4*t109*t19*t24*t50*t8*wght 1623 t217 = 1/t11**2.333333333333334d+0 1624 t218 = -6.452057056018142d-4*t110*t18*t19*t9*wght 1625 t219 = 1/t9**3 1626 t220 = 1.5874010519681997d+0*t110*t16*t19*wght 1627 t221 = t12*t3*t8*(1.311232602576965d-7*t17*t19*t217*t44*wght 1628 1 -4.8390427920136064d-4*t109*t18*t19*t44*wght+7.9370052598 1629 2 40999d-1*t1*t16*t19*(-2.0322706250000003d-4*t18*t219*t41* 1630 3 t43-2.370982395833334d-4*t17*t2*t219*t42)*wght+1.31123260 1631 4 2576965d-7*t19*t217*t41*t43*wght-4.8390427920136064d-4*t1 1632 5 09*t17*t18*t19*t41*wght+t220+t218) 1633 t222 = 2.4195213960068032d-4*t109*t19*t3*t33*t34*t50*wght-6. 1634 1 556163012884825d-8*t122*t17*t19*t3*t33*t34*wght+t221+t216 1635 2 +t215+t214+t213 1636 t223 = -8.741550683846434d-8*t17*t19*t20*t212*t8*wght 1637 t224 = -2.62246520515393d-7*t122*t17*t19*t24*t8*wght 1638 t225 = 3.226028528009071d-4*t104*t109*t19*t20*t8*wght 1639 t226 = 9.678085584027213d-4*t109*t19*t24*t50*t8*wght 1640 t227 = t12*t3*t8*(-3.226028528009071d-4*t118*t18*t19*t9*wght 1641 1 +2.62246520515393d-7*t17*t19*t217*t44*wght-9.678085584027 1642 2 213d-4*t109*t18*t19*t44*wght+7.937005259840999d-1*t1*t16* 1643 3 t19*(-4.0645412500000005d-4*t18*t219*t41*t43-4.7419647916 1644 4 66668d-4*t17*t2*t219*t42)*wght+2.62246520515393d-7*t19*t2 1645 5 17*t41*t43*wght-9.678085584027213d-4*t109*t17*t18*t19*t41 1646 6 *wght+7.937005259840999d-1*t118*t16*t19*wght+t220+t218) 1647 t228 = 4.8390427920136064d-4*t109*t19*t3*t33*t34*t50*wght-1. 1648 1 311232602576965d-7*t122*t17*t19*t3*t33*t34*wght+t227+t226 1649 2 +t225+t224+t223 1650 t229 = -1.748310136769287d-7*t17*t19*t20*t212*t8*wght 1651 t230 = -5.24493041030786d-7*t122*t17*t19*t24*t8*wght 1652 t231 = 6.452057056018142d-4*t104*t109*t19*t20*t8*wght 1653 t232 = 1.9356171168054426d-3*t109*t19*t24*t50*t8*wght 1654 t233 = t12*t3*t8*(-1.2904114112036286d-3*t118*t18*t19*t9*wgh 1655 1 t+5.24493041030786d-7*t17*t19*t217*t44*wght-1.93561711680 1656 2 54426d-3*t109*t18*t19*t44*wght+7.937005259840999d-1*t1*t1 1657 3 6*t19*(-8.129082500000001d-4*t18*t219*t41*t43-9.483929583 1658 4 333336d-4*t17*t2*t219*t42)*wght+5.24493041030786d-7*t19*t 1659 5 217*t41*t43*wght-1.9356171168054426d-3*t109*t17*t18*t19*t 1660 6 41*wght+3.1748021039363994d+0*t118*t16*t19*wght) 1661 t234 = 2.4195213960068032d-4*t109*t19*t3*t34*t49*t50*wght-6. 1662 1 556163012884825d-8*t122*t17*t19*t3*t34*t49*wght+t221+t216 1663 2 +t215+t214+t213 1664 t235 = 4.8390427920136064d-4*t109*t19*t3*t34*t49*t50*wght-1. 1665 1 311232602576965d-7*t122*t17*t19*t3*t34*t49*wght+t227+t226 1666 2 +t225+t224+t223 1667 t236 = 1/t11**3.833333333333334d+0 1668 t237 = 1/t1 1669 t238 = 2.4195213960068032d-4*t19*t219*t3*t50*t8*wght-5.32955 1670 1 9001518931d-11*t109*t19*t236*t3*t43*t8*wght+1.96684890386 1671 2 54477d-7*t122*t17*t19*t237*t3*t8*wght 1672 t239 = 4.8390427920136064d-4*t19*t219*t3*t50*t8*wght-1.06591 1673 1 18003037862d-10*t109*t19*t236*t3*t43*t8*wght+3.9336978077 1674 2 308954d-7*t122*t17*t19*t237*t3*t8*wght 1675 t240 = 9.678085584027213d-4*t19*t219*t3*t50*t8*wght-2.131823 1676 1 6006075725d-10*t109*t19*t236*t3*t43*t8*wght+7.86739561546 1677 2 1791d-7*t122*t17*t19*t237*t3*t8*wght 1678 fnc(iq) = 7.937005259840999d-1*t1*t12*t16*t19*t3*t8*wght+fnc 1679 1 (iq) 1680 Amat(iq,D1_RA) = -3.9685026299204995d-1*t1*t12*t16*t19*t3*t3 1681 1 3*t34*wght+t46+t35+t22+Amat(iq,D1_RA) 1682 Amat(iq,D1_RB) = -3.9685026299204995d-1*t1*t12*t16*t19*t3*t3 1683 1 4*t49*wght+t46+t35+t22+Amat(iq,D1_RB) 1684 Cmat(iq,D1_GAA) = t51+Cmat(iq,D1_GAA) 1685 Cmat(iq,D1_GAB) = -6.452057056018142d-4*t19*t3*t50*t8*t9*wgh 1686 1 t+1.5874010519681997d+0*t12*t16*t19*t3*t8*wght+Cmat(iq,D1 1687 2 _GAB) 1688 Cmat(iq,D1_GBB) = t51+Cmat(iq,D1_GBB) 1689 Amat2(iq,D2_RA_RA) = 5.952753944880749d-1*t1*t12*t16*t19*t3* 1690 1 t54*t55*wght+7.937005259840999d-1*t1*t12*t16*t19*t24*t33* 1691 2 t34*wght+2.6456684199469993d-1*t1*t16*t19*t20*t21*t33*t34 1692 3 *wght+t12*t3*t34*t91+t77+t75+t73+t71+t58+t53-5.0d-1*t12*t 1693 4 3*t33*t34*t45+Amat2(iq,D2_RA_RA) 1694 Amat2(iq,D2_RA_RB) = 5.952753944880749d-1*t1*t12*t16*t19*t3* 1695 1 t33*t49*t55*wght+t12*t3*t34*t96+t94+t12*t24*t34*t93+t20*t 1696 2 21*t34*t92+t77+t75+t73+t71+t58+t53+Amat2(iq,D2_RA_RB) 1697 Amat2(iq,D2_RB_RB) = 5.952753944880749d-1*t1*t12*t16*t19*t3* 1698 1 t55*t97*wght+7.937005259840999d-1*t1*t12*t16*t19*t24*t34* 1699 2 t49*wght+2.6456684199469993d-1*t1*t16*t19*t20*t21*t34*t49 1700 3 *wght+t94+t77+t75+t73+t71+t58+t53+t103*t12*t3*t34+Amat2(i 1701 4 q,D2_RB_RB) 1702 Cmat2(iq,D2_RA_GAA) = t113+Cmat2(iq,D2_RA_GAA) 1703 Cmat2(iq,D2_RA_GAB) = 3.226028528009071d-4*t19*t3*t33*t34*t5 1704 1 0*t9*wght-7.937005259840999d-1*t12*t16*t19*t3*t33*t34*wgh 1705 2 t+t120+t117+t116+t115+t114+Cmat2(iq,D2_RA_GAB) 1706 Cmat2(iq,D2_RA_GBB) = t113+Cmat2(iq,D2_RA_GBB) 1707 Cmat2(iq,D2_RB_GAA) = t121+Cmat2(iq,D2_RB_GAA) 1708 Cmat2(iq,D2_RB_GAB) = 3.226028528009071d-4*t19*t3*t34*t49*t5 1709 1 0*t9*wght-7.937005259840999d-1*t12*t16*t19*t3*t34*t49*wgh 1710 2 t+t120+t117+t116+t115+t114+Cmat2(iq,D2_RB_GAB) 1711 Cmat2(iq,D2_RB_GBB) = t121+Cmat2(iq,D2_RB_GBB) 1712 Cmat2(iq,D2_GAA_GAA) = t123+Cmat2(iq,D2_GAA_GAA) 1713 Cmat2(iq,D2_GAA_GAB) = t124+Cmat2(iq,D2_GAA_GAB) 1714 Cmat2(iq,D2_GAA_GBB) = t123+Cmat2(iq,D2_GAA_GBB) 1715 Cmat2(iq,D2_GAB_GAB) = -1.9356171168054426d-3*t109*t19*t3*t5 1716 1 0*t8*wght+5.24493041030786d-7*t122*t17*t19*t3*t8*wght+Cma 1717 2 t2(iq,D2_GAB_GAB) 1718 Cmat2(iq,D2_GAB_GBB) = t124+Cmat2(iq,D2_GAB_GBB) 1719 Cmat2(iq,D2_GBB_GBB) = t123+Cmat2(iq,D2_GBB_GBB) 1720 Amat3(iq,D3_RA_RA_RA) = t12*t3*t34*(-7.937005259840999d-1*t1 1721 1 *t16*t19*t44*t89*wght-7.937005259840999d-1*t1*t19*t41*t89 1722 2 *wght-3.9685026299204995d-1*t1*t16*t19*(1.049934208245727 1723 3 8d+0*t31*t85*t88+1.0499342082457278d+0*t27*t80*t83+5.2496 1724 4 71041228638d-1*(t154+t153)*t32-1.1665935647174752d-1*t152 1725 5 *t31**3+5.249671041228638d-1*(t151+t150)*t28-1.1665935647 1726 6 174752d-1*t149*t27**3)*wght+t157+t156+t155+t148+t147)+t12 1727 7 *t3*t55*(1.1905507889761498d+0*t1*t16*t19*t33*t89*wght+t1 1728 8 43+t142)+t12*t24*t34*(1.1905507889761498d+0*t1*t16*t19*t8 1729 9 9*wght+t141+t140)+t20*t21*t34*(2.6456684199469993d-1*t1*t 1730 : 16*t19*t89*wght+t135+t134)-2.3811015779522995d+0*t1*t12*t 1731 ; 16*t19*t33*t34*t57*wght-1.7858261834642247d+0*t1*t12*t16* 1732 < t19*t24*t54*t55*wght-5.952753944880749d-1*t1*t16*t19*t20* 1733 = t21*t54*t55*wght-1.3228342099734997d+0*t1*t16*t19*t21*t33 1734 > *t34*t52*wght-1.4881884862201872d+0*t1*t12*t126*t16*t19*t 1735 ? 3*t33**3*wght-1.5d+0*t12*t3*t33*t55*t91-5.0d-1*t12*t3*t34 1736 @ *t45*t89-5.0d-1*t20*t21*t33*t34*t76+1.6666666666666666d-1 1737 1 *t20*t21*t34*(t162-t33*t74)+5.0d-1*t12*t24*t34*(t162-t33* 1738 2 t72)+7.5d-1*t12*t3*t45*t54*t55+t161+t160+t146+t139+t133+t 1739 3 132+t131+t130+t129+t128+t125+Amat3(iq,D3_RA_RA_RA) 1740 Amat3(iq,D3_RA_RA_RB) = t12*t3*t34*(-3.9685026299204995d-1*t 1741 1 1*t16*t19*t44*t95*wght-3.9685026299204995d-1*t1*t19*t41*t 1742 2 95*wght-3.9685026299204995d-1*t1*t16*t19*t44*t89*wght-3.9 1743 3 685026299204995d-1*t1*t19*t41*t89*wght-3.9685026299204995 1744 4 d-1*t1*t16*t19*(3.4997806941524257d-1*t48*t85*t88+1.39991 1745 5 22776609702d+0*t31*t4*t57*t85-1.1665935647174752d-1*t152* 1746 6 t48*t84+3.4997806941524257d-1*t47*t80*t83-1.3999122776609 1747 7 702d+0*t27*t4*t57*t80-1.1665935647174752d-1*t149*t47*t79+ 1748 8 5.249671041228638d-1*t32*(2.0d+0*t57+t153)+5.249671041228 1749 9 638d-1*t28*(t150-2.0d+0*t57))*wght+t156+t155+t148+t147)+t 1750 : 12*t3*t55*(1.1905507889761498d+0*t1*t16*t19*t33*t95*wght+ 1751 ; t143+t142)+t20*t21*t34*(2.6456684199469993d-1*t1*t16*t19* 1752 < t95*wght+t135+t134)+t20*t21*t55*(-1.9842513149602495d-1*t 1753 = 1*t16*t19*t54*wght-3.968502629920499d-1*t1*t16*t19*t33*t4 1754 > 9*wght)+t12*t24*t55*(-5.952753944880749d-1*t1*t16*t19*t54 1755 ? *wght-1.1905507889761498d+0*t1*t16*t19*t33*t49*wght)+t21* 1756 @ t34*t52*(-4.4094473665783324d-1*t1*t16*t19*t49*wght-8.818 1757 1 894733156665d-1*t1*t16*t19*t33*wght)+t12*t34*t57*(t163-1. 1758 2 5874010519681997d+0*t1*t16*t19*t33*wght)-1.48818848622018 1759 3 72d+0*t1*t12*t126*t16*t19*t3*t49*t54*wght+5.0d-1*t12*t3*t 1760 4 34*(-t45*t95+t166+t157)-1.5d+0*t12*t3*t49*t55*t91+7.5d-1* 1761 5 t12*t3*t33*t45*t49*t55+5.0d-1*t12*(t167+t162)*t24*t34+t12 1762 6 *(t165+t141+t140+t137)*t24*t34+1.6666666666666666d-1*(t16 1763 7 8+t162)*t20*t21*t34+t164+t161+t160+t146+t139+t133+t132+t1 1764 8 31+t130+t129+t128+t125+Amat3(iq,D3_RA_RA_RB) 1765 Amat3(iq,D3_RA_RB_RB) = t12*t3*t34*(-3.9685026299204995d-1*t 1766 1 1*t16*t19*(3.4997806941524257d-1*t27*t80*t99-1.1665935647 1767 2 174752d-1*t149*t27*t98+1.3999122776609702d+0*t4*t48*t57*t 1768 3 85+3.4997806941524257d-1*t101*t31*t85-1.3999122776609702d 1769 4 +0*t4*t47*t57*t80-1.0499342082457278d+0*t32*t57+1.0499342 1770 5 082457278d+0*t28*t57-3.149802624737183d+0*t127*t32*t4+3.1 1771 6 49802624737183d+0*t127*t28*t4-1.1665935647174752d-1*t100* 1772 7 t152*t31)*wght-7.937005259840999d-1*t1*t16*t19*t44*t95*wg 1773 8 ht-7.937005259840999d-1*t1*t19*t41*t95*wght+t166+t156+t15 1774 9 5+t148+t147)+t12*t3*t55*(5.952753944880749d-1*t1*t16*t19* 1775 : t49*t95*wght+5.952753944880749d-1*t1*t16*t19*t33*t44*t49* 1776 ; wght+5.952753944880749d-1*t1*t19*t33*t41*t49*wght+5.95275 1777 < 3944880749d-1*t1*t102*t16*t19*t33*wght)+3.333333333333333 1778 = d-1*t20*t21*t34*(3.9685026299204995d-1*t1*t16*t19*t95*wgh 1779 > t+t138+t136)+t20*t21*t34*(1.3228342099734997d-1*t1*t16*t1 1780 ? 9*t95*wght+1.3228342099734997d-1*t1*t16*t19*t44*t49*wght+ 1781 @ 1.3228342099734997d-1*t1*t19*t41*t49*wght+1.3228342099734 1782 1 997d-1*t1*t16*t19*t33*t44*wght+1.3228342099734997d-1*t1*t 1783 2 19*t33*t41*wght+1.3228342099734997d-1*t1*t102*t16*t19*wgh 1784 3 t)+3.333333333333333d-1*t21*t34*t52*(-3.9685026299204995d 1785 4 -1*t1*t16*t19*t49*wght-3.9685026299204995d-1*t1*t16*t19*t 1786 5 33*wght-7*t92)+t12*t24*t34*(7.937005259840999d-1*t1*t16*t 1787 6 19*t33*t44*wght+7.937005259840999d-1*t1*t19*t33*t41*wght+ 1788 7 t173+t172+t171+t165)-1.4881884862201872d+0*t1*t12*t126*t1 1789 8 6*t19*t3*t33*t97*wght-5.952753944880749d-1*t1*t12*t16*t19 1790 9 *t24*t33*t49*t55*wght-1.9842513149602495d-1*t1*t16*t19*t2 1791 : 0*t21*t33*t49*t55*wght-4.4094473665783324d-1*t1*t16*t19*t 1792 ; 21*t34*t49*t52*wght-1.5d+0*t12*t3*t49*t55*t96-1.5d+0*t12* 1793 < t24*t49*t55*t93+t12*t34*t57*(t163-2*t93)-1.5d+0*t20*t21*t 1794 = 49*t55*t92+t176+t175+t170+t169+t164+t161+t160+t146+t133+t 1795 > 132+t131+t130+t129+t128+t125+Amat3(iq,D3_RA_RB_RB) 1796 Amat3(iq,D3_RB_RB_RB) = t12*t3*t34*(-3.9685026299204995d-1*t 1797 1 1*t16*t19*(1.0499342082457278d+0*t47*t80*t99+1.0499342082 1798 2 457278d+0*t101*t48*t85-1.1665935647174752d-1*t152*t48**3- 1799 3 1.1665935647174752d-1*t149*t47**3+5.249671041228638d-1*(t 1800 4 153+t151)*t32+5.249671041228638d-1*(t154+t150)*t28)*wght- 1801 5 3.9685026299204995d-1*t1*t16*t19*t49*t69*wght-3.968502629 1802 6 9204995d-1*t1*t19*t49*t65*wght-3.9685026299204995d-1*t1*t 1803 7 16*t19*t49*t59*wght-7.937005259840999d-1*t1*t19*t41*t44*t 1804 8 49*wght-7.937005259840999d-1*t1*t102*t16*t19*t44*wght-7.9 1805 9 37005259840999d-1*t1*t102*t19*t41*wght+t166)+t12*t3*t55*( 1806 : 5.952753944880749d-1*t1*t16*t19*t44*t97*wght+5.9527539448 1807 ; 80749d-1*t1*t19*t41*t97*wght+1.1905507889761498d+0*t1*t10 1808 < 2*t16*t19*t49*wght)+t12*t24*t34*(1.1905507889761498d+0*t1 1809 = *t16*t19*t44*t49*wght+1.1905507889761498d+0*t1*t19*t41*t4 1810 > 9*wght+1.1905507889761498d+0*t1*t102*t16*t19*wght)+t20*t2 1811 ? 1*t34*(2.6456684199469993d-1*t1*t16*t19*t44*t49*wght+2.64 1812 @ 56684199469993d-1*t1*t19*t41*t49*wght+2.6456684199469993d 1813 1 -1*t1*t102*t16*t19*wght)-1.7858261834642247d+0*t1*t12*t16 1814 2 *t19*t24*t55*t97*wght-5.952753944880749d-1*t1*t16*t19*t20 1815 3 *t21*t55*t97*wght-2.3811015779522995d+0*t1*t12*t16*t19*t3 1816 4 4*t49*t57*wght-1.3228342099734997d+0*t1*t16*t19*t21*t34*t 1817 5 49*t52*wght-1.4881884862201872d+0*t1*t12*t126*t16*t19*t3* 1818 6 t49**3*wght-1.5d+0*t103*t12*t3*t49*t55+3.333333333333333d 1819 7 -1*(t173+t172+t171)*t20*t21*t34+t176+t175+t170+t169+t164+ 1820 8 t161+t160+t146+t133+t132+t131+t130+t129+t128+t125+Amat3(i 1821 9 q,D3_RB_RB_RB) 1822 Cmat3(iq,D3_RA_RA_GAA) = t194+Cmat3(iq,D3_RA_RA_GAA) 1823 Cmat3(iq,D3_RA_RA_GAB) = t12*t3*t34*(3.226028528009071d-4*t1 1824 1 8*t19*t89*t9*wght-7.937005259840999d-1*t16*t19*t89*wght+t 1825 2 206+t205+t204+t203+t202)-4.8390427920136064d-4*t19*t3*t50 1826 3 *t54*t55*t9*wght-6.452057056018142d-4*t19*t24*t33*t34*t50 1827 4 *t9*wght-2.1506856853393808d-4*t104*t19*t20*t33*t34*t9*wg 1828 5 ht+1.1905507889761498d+0*t12*t16*t19*t3*t54*t55*wght+1.58 1829 6 74010519681997d+0*t12*t16*t19*t24*t33*t34*wght+5.29133683 1830 7 9893999d-1*t16*t19*t20*t21*t33*t34*wght-5.0d-1*t119*t12*t 1831 8 3*t33*t34+t207+t201+t200+t199+t198+t197+t196+t195+Cmat3(i 1832 9 q,D3_RA_RA_GAB) 1833 Cmat3(iq,D3_RA_RA_GBB) = t194+Cmat3(iq,D3_RA_RA_GBB) 1834 Cmat3(iq,D3_RA_RB_GAA) = t209+Cmat3(iq,D3_RA_RB_GAA) 1835 Cmat3(iq,D3_RA_RB_GAB) = t12*t3*t34*(3.226028528009071d-4*t1 1836 1 8*t19*t9*t95*wght-7.937005259840999d-1*t16*t19*t95*wght+t 1837 2 206+t205+t204+t203+t202)+t20*t21*t34*(-1.0753428426696904 1838 3 d-4*t18*t19*t49*t9*wght-1.0753428426696904d-4*t18*t19*t33 1839 4 *t9*wght+2.6456684199469993d-1*t16*t19*t49*wght+2.6456684 1840 5 199469993d-1*t16*t19*t33*wght)+t12*t24*t34*(-3.2260285280 1841 6 09071d-4*t18*t19*t49*t9*wght-3.226028528009071d-4*t18*t19 1842 7 *t33*t9*wght+7.937005259840999d-1*t16*t19*t49*wght+7.9370 1843 8 05259840999d-1*t16*t19*t33*wght)-4.8390427920136064d-4*t1 1844 9 9*t3*t33*t49*t50*t55*t9*wght+1.1905507889761498d+0*t12*t1 1845 : 6*t19*t3*t33*t49*t55*wght+t210+t207+t201+t200+t199+t198+t 1846 ; 197+t196+t195+Cmat3(iq,D3_RA_RB_GAB) 1847 Cmat3(iq,D3_RA_RB_GBB) = t209+Cmat3(iq,D3_RA_RB_GBB) 1848 Cmat3(iq,D3_RB_RB_GAA) = t211+Cmat3(iq,D3_RB_RB_GAA) 1849 Cmat3(iq,D3_RB_RB_GAB) = t12*t3*t34*(3.226028528009071d-4*t1 1850 1 8*t19*t44*t49*t9*wght+3.226028528009071d-4*t17*t18*t19*t4 1851 2 1*t49*t9*wght+3.226028528009071d-4*t102*t18*t19*t9*wght-7 1852 3 .937005259840999d-1*t16*t19*t44*t49*wght-7.93700525984099 1853 4 9d-1*t19*t41*t49*wght-3.9685026299204995d-1*t1*t118*t16*t 1854 5 19*t49*wght-7.937005259840999d-1*t102*t16*t19*wght)-4.839 1855 6 0427920136064d-4*t19*t3*t50*t55*t9*t97*wght+1.19055078897 1856 7 61498d+0*t12*t16*t19*t3*t55*t97*wght-6.452057056018142d-4 1857 8 *t19*t24*t34*t49*t50*t9*wght-2.1506856853393808d-4*t104*t 1858 9 19*t20*t34*t49*t9*wght+1.5874010519681997d+0*t12*t16*t19* 1859 : t24*t34*t49*wght+5.291336839893999d-1*t16*t19*t20*t21*t34 1860 ; *t49*wght+t210+t207+t201+t200+t199+t198+t197+t196+t195+Cm 1861 < at3(iq,D3_RB_RB_GAB) 1862 Cmat3(iq,D3_RB_RB_GBB) = t211+Cmat3(iq,D3_RB_RB_GBB) 1863 Cmat3(iq,D3_RA_GAA_GAA) = t222+Cmat3(iq,D3_RA_GAA_GAA) 1864 Cmat3(iq,D3_RA_GAA_GAB) = t228+Cmat3(iq,D3_RA_GAA_GAB) 1865 Cmat3(iq,D3_RA_GAA_GBB) = t222+Cmat3(iq,D3_RA_GAA_GBB) 1866 Cmat3(iq,D3_RA_GAB_GAB) = 9.678085584027213d-4*t109*t19*t3*t 1867 1 33*t34*t50*wght-2.62246520515393d-7*t122*t17*t19*t3*t33*t 1868 2 34*wght+t233+t232+t231+t230+t229+Cmat3(iq,D3_RA_GAB_GAB) 1869 Cmat3(iq,D3_RA_GAB_GBB) = t228+Cmat3(iq,D3_RA_GAB_GBB) 1870 Cmat3(iq,D3_RA_GBB_GBB) = t222+Cmat3(iq,D3_RA_GBB_GBB) 1871 Cmat3(iq,D3_RB_GAA_GAA) = t234+Cmat3(iq,D3_RB_GAA_GAA) 1872 Cmat3(iq,D3_RB_GAA_GAB) = t235+Cmat3(iq,D3_RB_GAA_GAB) 1873 Cmat3(iq,D3_RB_GAA_GBB) = t234+Cmat3(iq,D3_RB_GAA_GBB) 1874 Cmat3(iq,D3_RB_GAB_GAB) = 9.678085584027213d-4*t109*t19*t3*t 1875 1 34*t49*t50*wght-2.62246520515393d-7*t122*t17*t19*t3*t34*t 1876 2 49*wght+t233+t232+t231+t230+t229+Cmat3(iq,D3_RB_GAB_GAB) 1877 Cmat3(iq,D3_RB_GAB_GBB) = t235+Cmat3(iq,D3_RB_GAB_GBB) 1878 Cmat3(iq,D3_RB_GBB_GBB) = t234+Cmat3(iq,D3_RB_GBB_GBB) 1879 Cmat3(iq,D3_GAA_GAA_GAA) = t238+Cmat3(iq,D3_GAA_GAA_GAA) 1880 Cmat3(iq,D3_GAA_GAA_GAB) = t239+Cmat3(iq,D3_GAA_GAA_GAB) 1881 Cmat3(iq,D3_GAA_GAA_GBB) = t238+Cmat3(iq,D3_GAA_GAA_GBB) 1882 Cmat3(iq,D3_GAA_GAB_GAB) = t240+Cmat3(iq,D3_GAA_GAB_GAB) 1883 Cmat3(iq,D3_GAA_GAB_GBB) = t239+Cmat3(iq,D3_GAA_GAB_GBB) 1884 Cmat3(iq,D3_GAA_GBB_GBB) = t238+Cmat3(iq,D3_GAA_GBB_GBB) 1885 Cmat3(iq,D3_GAB_GAB_GAB) = 1.9356171168054426d-3*t19*t219*t3 1886 1 *t50*t8*wght-4.263647201215145d-10*t109*t19*t236*t3*t43*t 1887 2 8*wght+1.573479123092358d-6*t122*t17*t19*t237*t3*t8*wght+ 1888 3 Cmat3(iq,D3_GAB_GAB_GAB) 1889 Cmat3(iq,D3_GAB_GAB_GBB) = t240+Cmat3(iq,D3_GAB_GAB_GBB) 1890 Cmat3(iq,D3_GAB_GBB_GBB) = t239+Cmat3(iq,D3_GAB_GBB_GBB) 1891 Cmat3(iq,D3_GBB_GBB_GBB) = t238+Cmat3(iq,D3_GBB_GBB_GBB) 1892 elseif (rhoa.gt.tol_rho.and.rhob.le.tol_rho) then 1893 t1 = 1/rhoa 1894 t2 = gammaaa**5.0d-1 1895 t3 = 1/rhoa**6.666666666666666d-1 1896 t4 = abs(rhoa) 1897 t5 = 1/t4**3.333333333333333d-1 1898 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 1899 1 -3 1900 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 1901 1 16419932959077d-1*t3+1.0d+0 1902 t8 = 1/t7 1903 t9 = t6*t8+1.6670000000000001d-3 1904 t10 = 1/t9 1905 t11 = 1/t4**1.166666666666667d+0 1906 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 1907 t13 = 1/rhoa**1.3333333333333333d+0 1908 t14 = 1/t4 1909 t15 = rhoa**2 1910 t16 = 1/t15 1911 t17 = 1/rhoa**1.6666666666666669d+0 1912 t18 = 1/rhoa**3.333333333333333d-1 1913 t19 = -1.1759958745060142d-2*t17*t5-1.8037724440384886d+0*t1 1914 1 4*t18-1.2109466219727179d-1*t17-5.879979372530072d-3*t1*t 1915 2 14 1916 t20 = 1/t7**2 1917 t21 = -4.8110248404218126d-3*t14*t18-1.895695887660257d-6*t1 1918 1 7 1919 t22 = t21*t8-t19*t20*t6 1920 t23 = 1/t4**3.166666666666667d+0 1921 t24 = 1/t9**2 1922 t25 = 8.129082500000001d-4*t11*t2*t22*t24+9.483929583333336d 1923 1 -4*rhoa*t10*t2*t23 1924 t26 = 1/t4**1.5000000000000002d+0 1925 t27 = 1/rhoa**2.3333333333333334d+0 1926 t28 = rhoa**3 1927 t29 = 1/t28 1928 t30 = t25**2 1929 t31 = t19**2 1930 t32 = 1/t7**3 1931 t33 = 1/rhoa**2.6666666666666666d+0 1932 t34 = 1.9599931241766905d-2*t33*t5+2.01824436995453d-1*t33+1 1933 1 .5679944993413522d-2*t14*t16+2.4050299253846513d+0*t13*t1 1934 2 4 1935 t35 = 3.1594931461004283d-6*t33+6.414699787229083d-3*t13*t14 1936 t36 = t35*t8-t20*t34*t6+2*t31*t32*t6-2*t19*t20*t21 1937 t37 = 1/t4**5.166666666666667d+0 1938 t38 = t22**2 1939 t39 = 1/t9**3 1940 t40 = -1.6258165000000002d-3*t11*t2*t38*t39-3.00324436805555 1941 1 7d-3*t10*t15*t2*t37+8.129082500000001d-4*t11*t2*t24*t36-1 1942 2 .8967859166666673d-3*rhoa*t2*t22*t23*t24+9.48392958333333 1943 3 6d-4*t10*t2*t23 1944 t41 = -1.5874010519681997d+0*gammaaa*t12*t25*t9*wght-1.58740 1945 1 10519681997d+0*gammaaa*t12*t22*wght 1946 t42 = -7.937005259840999d-1*gammaaa*t12*t25*t9*wght-7.937005 1947 1 259840999d-1*gammaaa*t12*t22*wght 1948 t43 = -2.6456684199469993d-1*gammaaa*t12*t25*t9*wght-2.64566 1949 1 84199469993d-1*gammaaa*t12*t22*wght 1950 t44 = 1/t4**2.166666666666667d+0 1951 t45 = 1/t2 1952 t46 = 4.0645412500000005d-4*t11*t22*t24*t45+4.74196479166666 1953 1 8d-4*rhoa*t10*t23*t45 1954 t47 = 1/t4**2.666666666666667d+0 1955 t48 = 1/rhoa**3.6666666666666664d+0 1956 t49 = (-8.425315056267809d-6*t48-1.496763283686786d-2*t14*t2 1957 1 7)*t8-6*t19**3*t6/t7**4-t20*(-5.226648331137841d-2*t48*t5 1958 2 -5.381984986545413d-1*t48-5.357314539416287d-2*t14*t29-5. 1959 3 611736492564187d+0*t14*t27)*t6+6*t19*t32*t34*t6-3*t19*t20 1960 4 *t35-3*t20*t21*t34+6*t21*t31*t32 1961 t50 = 1.5874010519681997d+0*gammaaa*t12*t22*wght 1962 t51 = 1.5874010519681997d+0*gammaaa*t12*t25*t9*wght 1963 t52 = 1/t4**2.333333333333334d+0 1964 t53 = 1/t2**3 1965 fnc(iq) = 7.937005259840999d-1*gammaaa*t1*t12*t5*t9*wght+fnc 1966 1 (iq) 1967 Amat(iq,D1_RA) = t1*t5*(7.937005259840999d-1*gammaaa*t12*t25 1968 1 *t9*wght+7.937005259840999d-1*gammaaa*t12*t22*wght)-7.937 1969 2 005259840999d-1*gammaaa*t12*t16*t5*t9*wght-2.645668419946 1970 3 9993d-1*gammaaa*t12*t13*t14*t9*wght+Amat(iq,D1_RA) 1971 Cmat(iq,D1_GAA) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 1972 1 26028528009071d-4*t1*t12*t2*t26*wght+Cmat(iq,D1_GAA) 1973 Amat2(iq,D2_RA_RA) = t1*t5*(7.937005259840999d-1*gammaaa*t12 1974 1 *t40*t9*wght+7.937005259840999d-1*gammaaa*t12*t30*t9*wght 1975 2 +7.937005259840999d-1*gammaaa*t12*t36*wght+1.587401051968 1976 3 1997d+0*gammaaa*t12*t22*t25*wght)+1.5874010519681997d+0*g 1977 4 ammaaa*t12*t29*t5*t9*wght+8.818894733156665d-1*gammaaa*t1 1978 5 2*t14*t27*t9*wght+t16*t41*t5+t13*t14*t43+3.33333333333333 1979 6 3d-1*t13*t14*t42+Amat2(iq,D2_RA_RA) 1980 Cmat2(iq,D2_RA_GAA) = t1*t5*(7.937005259840999d-1*gammaaa*t1 1981 1 2*t46*t9*wght+7.937005259840999d-1*t12*t25*t9*wght-3.2260 1982 2 28528009071d-4*t11*t12*t2*t25*wght-3.226028528009071d-4*t 1983 3 10*t11*t12*t2*t22*wght+7.937005259840999d-1*t12*t22*wght) 1984 4 -7.937005259840999d-1*t12*t16*t5*t9*wght-2.64566841994699 1985 5 93d-1*t12*t13*t14*t9*wght+1.0753428426696904d-4*t12*t13*t 1986 6 2*t44*wght+3.226028528009071d-4*t12*t16*t2*t26*wght+Cmat2 1987 7 (iq,D2_RA_GAA) 1988 Cmat2(iq,D2_GAA_GAA) = 1.311232602576965d-7*t1*t10*t12*t47*w 1989 1 ght-4.8390427920136064d-4*t1*t12*t26*t45*wght+Cmat2(iq,D2 1990 2 _GAA_GAA) 1991 Amat3(iq,D3_RA_RA_RA) = t1*t5*(7.937005259840999d-1*gammaaa* 1992 1 t12*(4.8774495000000007d-3*t11*t2*t22**3/t9**4+8.12908250 1993 2 0000001d-4*t11*t2*t24*t49+1.5516762568287046d-2*t10*t2*t2 1994 3 8/t4**7.166666666666667d+0+5.690357750000001d-3*rhoa*t2*t 1995 4 23*t38*t39-4.8774495000000007d-3*t11*t2*t22*t36*t39+9.009 1996 5 733104166671d-3*t15*t2*t22*t24*t37-9.009733104166671d-3*r 1997 6 hoa*t10*t2*t37-2.8451788750000007d-3*rhoa*t2*t23*t24*t36- 1998 7 2.8451788750000007d-3*t2*t22*t23*t24)*t9*wght+2.381101577 1999 8 9522995d+0*gammaaa*t12*t25*t40*t9*wght+7.937005259840999d 2000 9 -1*gammaaa*t12*t25**3*t9*wght+7.937005259840999d-1*gammaa 2001 : a*t12*t49*wght+2.3811015779522995d+0*gammaaa*t12*t22*t40* 2002 ; wght+2.3811015779522995d+0*gammaaa*t12*t25*t36*wght+2.381 2003 < 1015779522995d+0*gammaaa*t12*t22*t30*wght)+t13*t14*(-2.64 2004 = 56684199469993d-1*gammaaa*t12*t40*t9*wght-2.6456684199469 2005 > 993d-1*gammaaa*t12*t30*t9*wght-2.6456684199469993d-1*gamm 2006 ? aaa*t12*t36*wght-5.291336839893999d-1*gammaaa*t12*t22*t25 2007 @ *wght)+3.333333333333333d-1*t13*t14*(-1.5874010519681997d 2008 1 +0*gammaaa*t12*t40*t9*wght-1.5874010519681997d+0*gammaaa* 2009 2 t12*t30*t9*wght-1.5874010519681997d+0*gammaaa*t12*t36*wgh 2010 3 t-3.1748021039363994d+0*gammaaa*t12*t22*t25*wght)+t16*t5* 2011 4 (-2.3811015779522995d+0*gammaaa*t12*t40*t9*wght-2.3811015 2012 5 779522995d+0*gammaaa*t12*t30*t9*wght-2.3811015779522995d+ 2013 6 0*gammaaa*t12*t36*wght-4.762203155904599d+0*gammaaa*t12*t 2014 7 22*t25*wght)+t14*t27*(8.818894733156665d-1*gammaaa*t12*t2 2015 8 5*t9*wght+8.818894733156665d-1*gammaaa*t12*t22*wght)-4.76 2016 9 2203155904599d+0*gammaaa*t12*t5*t9*wght/rhoa**4-3.4687652 2017 : 61708288d+0*gammaaa*t12*t14*t9*wght/rhoa**3.3333333333333 2018 ; 337d+0+3.333333333333333d-1*t14*t27*(t51+t50-7*t43)+t29*t 2019 < 5*(t51+t50-2*t41)-7.777777777777778d-1*t14*t27*t42+Amat3( 2020 = iq,D3_RA_RA_RA) 2021 Cmat3(iq,D3_RA_RA_GAA) = t1*t5*(1.5874010519681997d+0*gammaa 2022 1 a*t12*t25*t46*t9*wght+7.937005259840999d-1*gammaaa*t12*(- 2023 2 8.129082500000001d-4*t11*t38*t39*t45-1.5016221840277785d- 2024 3 3*t10*t15*t37*t45+4.0645412500000005d-4*t11*t24*t36*t45-9 2025 4 .483929583333336d-4*rhoa*t22*t23*t24*t45+4.74196479166666 2026 5 8d-4*t10*t23*t45)*t9*wght+7.937005259840999d-1*t12*t40*t9 2027 6 *wght+7.937005259840999d-1*t12*t30*t9*wght+1.587401051968 2028 7 1997d+0*gammaaa*t12*t22*t46*wght-3.226028528009071d-4*t11 2029 8 *t12*t2*t40*wght-3.226028528009071d-4*t10*t11*t12*t2*t36* 2030 9 wght+7.937005259840999d-1*t12*t36*wght-3.226028528009071d 2031 : -4*t11*t12*t2*t30*wght-6.452057056018142d-4*t10*t11*t12*t 2032 ; 2*t22*t25*wght+1.5874010519681997d+0*t12*t22*t25*wght)+t1 2033 < 3*t14*(-2.6456684199469993d-1*gammaaa*t12*t46*t9*wght-2.6 2034 = 456684199469993d-1*t12*t25*t9*wght+1.0753428426696904d-4* 2035 > t11*t12*t2*t25*wght+1.0753428426696904d-4*t10*t11*t12*t2* 2036 ? t22*wght-2.6456684199469993d-1*t12*t22*wght)+3.3333333333 2037 @ 33333d-1*t13*t14*(-7.937005259840999d-1*gammaaa*t12*t46*t 2038 1 9*wght-7.937005259840999d-1*t12*t25*t9*wght+3.22602852800 2039 2 9071d-4*t11*t12*t2*t25*wght+3.226028528009071d-4*t10*t11* 2040 3 t12*t2*t22*wght-7.937005259840999d-1*t12*t22*wght)+t16*t5 2041 4 *(-1.5874010519681997d+0*gammaaa*t12*t46*t9*wght-1.587401 2042 5 0519681997d+0*t12*t25*t9*wght+6.452057056018142d-4*t11*t1 2043 6 2*t2*t25*wght+6.452057056018142d-4*t10*t11*t12*t2*t22*wgh 2044 7 t-1.5874010519681997d+0*t12*t22*wght)+1.5874010519681997d 2045 8 +0*t12*t29*t5*t9*wght+8.818894733156665d-1*t12*t14*t27*t9 2046 9 *wght-3.5844761422323007d-4*t12*t2*t27*t44*wght-6.4520570 2047 : 56018142d-4*t12*t2*t26*t29*wght+Cmat3(iq,D3_RA_RA_GAA) 2048 Cmat3(iq,D3_RA_GAA_GAA) = t1*t5*(7.937005259840999d-1*gammaa 2049 1 a*t12*(-2.0322706250000003d-4*t11*t22*t24*t53-2.370982395 2050 2 833334d-4*rhoa*t10*t23*t53)*t9*wght+1.5874010519681997d+0 2051 3 *t12*t46*t9*wght+1.311232602576965d-7*t10*t12*t25*t52*wgh 2052 4 t+1.311232602576965d-7*t12*t22*t24*t52*wght-6.45205705601 2053 5 8142d-4*t11*t12*t2*t46*wght-4.8390427920136064d-4*t11*t12 2054 6 *t25*t45*wght-4.8390427920136064d-4*t10*t11*t12*t22*t45*w 2055 7 ght)-1.311232602576965d-7*t10*t12*t16*t47*wght+1.61301426 2056 8 40045356d-4*t12*t13*t44*t45*wght+4.8390427920136064d-4*t1 2057 9 2*t16*t26*t45*wght-4.370775341923217d-8*t10*t12*t13*wght/ 2058 : t4**3.3333333333333337d+0+Cmat3(iq,D3_RA_GAA_GAA) 2059 Cmat3(iq,D3_GAA_GAA_GAA) = 2.4195213960068032d-4*t1*t12*t26* 2060 1 t53*wght+1.9668489038654477d-7*t1*t10*t12*t47*wght/gammaa 2061 2 a-5.329559001518931d-11*t1*t12*t24*t45*wght/t4**3.8333333 2062 3 33333334d+0+Cmat3(iq,D3_GAA_GAA_GAA) 2063 elseif (rhoa.le.tol_rho.and.rhob.gt.tol_rho) then 2064 t1 = 1/rhob 2065 t2 = gammabb**5.0d-1 2066 t3 = 1/rhob**6.666666666666666d-1 2067 t4 = abs(rhob) 2068 t5 = 1/t4**3.333333333333333d-1 2069 t6 = 1.443307452126544d-2*t5+2.8435438314903855d-6*t3+2.568d 2070 1 -3 2071 t7 = 1.7639938117590215d-2*t3*t5+5.411317332115466d+0*t5+1.8 2072 1 16419932959077d-1*t3+1.0d+0 2073 t8 = 1/t7 2074 t9 = t6*t8+1.6670000000000001d-3 2075 t10 = 1/t9 2076 t11 = 1/t4**1.166666666666667d+0 2077 t12 = exp(-8.129082500000001d-4*t10*t11*t2) 2078 t13 = 1/rhob**1.3333333333333333d+0 2079 t14 = 1/t4 2080 t15 = rhob**2 2081 t16 = 1/t15 2082 t17 = 1/rhob**1.6666666666666669d+0 2083 t18 = 1/rhob**3.333333333333333d-1 2084 t19 = -1.1759958745060142d-2*t17*t5-1.8037724440384886d+0*t1 2085 1 4*t18-1.2109466219727179d-1*t17-5.879979372530072d-3*t1*t 2086 2 14 2087 t20 = 1/t7**2 2088 t21 = -4.8110248404218126d-3*t14*t18-1.895695887660257d-6*t1 2089 1 7 2090 t22 = t21*t8-t19*t20*t6 2091 t23 = 1/t4**3.166666666666667d+0 2092 t24 = 1/t9**2 2093 t25 = 8.129082500000001d-4*t11*t2*t22*t24+9.483929583333336d 2094 1 -4*rhob*t10*t2*t23 2095 t26 = 1/t4**1.5000000000000002d+0 2096 t27 = 1/rhob**2.3333333333333334d+0 2097 t28 = rhob**3 2098 t29 = 1/t28 2099 t30 = t25**2 2100 t31 = t19**2 2101 t32 = 1/t7**3 2102 t33 = 1/rhob**2.6666666666666666d+0 2103 t34 = 1.9599931241766905d-2*t33*t5+2.01824436995453d-1*t33+1 2104 1 .5679944993413522d-2*t14*t16+2.4050299253846513d+0*t13*t1 2105 2 4 2106 t35 = 3.1594931461004283d-6*t33+6.414699787229083d-3*t13*t14 2107 t36 = t35*t8-t20*t34*t6+2*t31*t32*t6-2*t19*t20*t21 2108 t37 = 1/t4**5.166666666666667d+0 2109 t38 = t22**2 2110 t39 = 1/t9**3 2111 t40 = -1.6258165000000002d-3*t11*t2*t38*t39-3.00324436805555 2112 1 7d-3*t10*t15*t2*t37+8.129082500000001d-4*t11*t2*t24*t36-1 2113 2 .8967859166666673d-3*rhob*t2*t22*t23*t24+9.48392958333333 2114 3 6d-4*t10*t2*t23 2115 t41 = -1.5874010519681997d+0*gammabb*t12*t25*t9*wght-1.58740 2116 1 10519681997d+0*gammabb*t12*t22*wght 2117 t42 = -7.937005259840999d-1*gammabb*t12*t25*t9*wght-7.937005 2118 1 259840999d-1*gammabb*t12*t22*wght 2119 t43 = -2.6456684199469993d-1*gammabb*t12*t25*t9*wght-2.64566 2120 1 84199469993d-1*gammabb*t12*t22*wght 2121 t44 = 1/t4**2.166666666666667d+0 2122 t45 = 1/t2 2123 t46 = 4.0645412500000005d-4*t11*t22*t24*t45+4.74196479166666 2124 1 8d-4*rhob*t10*t23*t45 2125 t47 = 1/t4**2.666666666666667d+0 2126 t48 = 1/rhob**3.6666666666666664d+0 2127 t49 = (-8.425315056267809d-6*t48-1.496763283686786d-2*t14*t2 2128 1 7)*t8-6*t19**3*t6/t7**4-t20*(-5.226648331137841d-2*t48*t5 2129 2 -5.381984986545413d-1*t48-5.357314539416287d-2*t14*t29-5. 2130 3 611736492564187d+0*t14*t27)*t6+6*t19*t32*t34*t6-3*t19*t20 2131 4 *t35-3*t20*t21*t34+6*t21*t31*t32 2132 t50 = 1.5874010519681997d+0*gammabb*t12*t22*wght 2133 t51 = 1.5874010519681997d+0*gammabb*t12*t25*t9*wght 2134 t52 = 1/t4**2.333333333333334d+0 2135 t53 = 1/t2**3 2136 fnc(iq) = 7.937005259840999d-1*gammabb*t1*t12*t5*t9*wght+fnc 2137 1 (iq) 2138 Amat(iq,D1_RB) = t1*t5*(7.937005259840999d-1*gammabb*t12*t25 2139 1 *t9*wght+7.937005259840999d-1*gammabb*t12*t22*wght)-7.937 2140 2 005259840999d-1*gammabb*t12*t16*t5*t9*wght-2.645668419946 2141 3 9993d-1*gammabb*t12*t13*t14*t9*wght+Amat(iq,D1_RB) 2142 Cmat(iq,D1_GBB) = 7.937005259840999d-1*t1*t12*t5*t9*wght-3.2 2143 1 26028528009071d-4*t1*t12*t2*t26*wght+Cmat(iq,D1_GBB) 2144 Amat2(iq,D2_RB_RB) = t1*t5*(7.937005259840999d-1*gammabb*t12 2145 1 *t40*t9*wght+7.937005259840999d-1*gammabb*t12*t30*t9*wght 2146 2 +7.937005259840999d-1*gammabb*t12*t36*wght+1.587401051968 2147 3 1997d+0*gammabb*t12*t22*t25*wght)+1.5874010519681997d+0*g 2148 4 ammabb*t12*t29*t5*t9*wght+8.818894733156665d-1*gammabb*t1 2149 5 2*t14*t27*t9*wght+t16*t41*t5+t13*t14*t43+3.33333333333333 2150 6 3d-1*t13*t14*t42+Amat2(iq,D2_RB_RB) 2151 Cmat2(iq,D2_RB_GBB) = t1*t5*(7.937005259840999d-1*gammabb*t1 2152 1 2*t46*t9*wght+7.937005259840999d-1*t12*t25*t9*wght-3.2260 2153 2 28528009071d-4*t11*t12*t2*t25*wght-3.226028528009071d-4*t 2154 3 10*t11*t12*t2*t22*wght+7.937005259840999d-1*t12*t22*wght) 2155 4 -7.937005259840999d-1*t12*t16*t5*t9*wght-2.64566841994699 2156 5 93d-1*t12*t13*t14*t9*wght+1.0753428426696904d-4*t12*t13*t 2157 6 2*t44*wght+3.226028528009071d-4*t12*t16*t2*t26*wght+Cmat2 2158 7 (iq,D2_RB_GBB) 2159 Cmat2(iq,D2_GBB_GBB) = 1.311232602576965d-7*t1*t10*t12*t47*w 2160 1 ght-4.8390427920136064d-4*t1*t12*t26*t45*wght+Cmat2(iq,D2 2161 2 _GBB_GBB) 2162 Amat3(iq,D3_RB_RB_RB) = t1*t5*(7.937005259840999d-1*gammabb* 2163 1 t12*(4.8774495000000007d-3*t11*t2*t22**3/t9**4+8.12908250 2164 2 0000001d-4*t11*t2*t24*t49+1.5516762568287046d-2*t10*t2*t2 2165 3 8/t4**7.166666666666667d+0+5.690357750000001d-3*rhob*t2*t 2166 4 23*t38*t39-4.8774495000000007d-3*t11*t2*t22*t36*t39+9.009 2167 5 733104166671d-3*t15*t2*t22*t24*t37-9.009733104166671d-3*r 2168 6 hob*t10*t2*t37-2.8451788750000007d-3*rhob*t2*t23*t24*t36- 2169 7 2.8451788750000007d-3*t2*t22*t23*t24)*t9*wght+2.381101577 2170 8 9522995d+0*gammabb*t12*t25*t40*t9*wght+7.937005259840999d 2171 9 -1*gammabb*t12*t25**3*t9*wght+7.937005259840999d-1*gammab 2172 : b*t12*t49*wght+2.3811015779522995d+0*gammabb*t12*t22*t40* 2173 ; wght+2.3811015779522995d+0*gammabb*t12*t25*t36*wght+2.381 2174 < 1015779522995d+0*gammabb*t12*t22*t30*wght)+t13*t14*(-2.64 2175 = 56684199469993d-1*gammabb*t12*t40*t9*wght-2.6456684199469 2176 > 993d-1*gammabb*t12*t30*t9*wght-2.6456684199469993d-1*gamm 2177 ? abb*t12*t36*wght-5.291336839893999d-1*gammabb*t12*t22*t25 2178 @ *wght)+3.333333333333333d-1*t13*t14*(-1.5874010519681997d 2179 1 +0*gammabb*t12*t40*t9*wght-1.5874010519681997d+0*gammabb* 2180 2 t12*t30*t9*wght-1.5874010519681997d+0*gammabb*t12*t36*wgh 2181 3 t-3.1748021039363994d+0*gammabb*t12*t22*t25*wght)+t16*t5* 2182 4 (-2.3811015779522995d+0*gammabb*t12*t40*t9*wght-2.3811015 2183 5 779522995d+0*gammabb*t12*t30*t9*wght-2.3811015779522995d+ 2184 6 0*gammabb*t12*t36*wght-4.762203155904599d+0*gammabb*t12*t 2185 7 22*t25*wght)+t14*t27*(8.818894733156665d-1*gammabb*t12*t2 2186 8 5*t9*wght+8.818894733156665d-1*gammabb*t12*t22*wght)-4.76 2187 9 2203155904599d+0*gammabb*t12*t5*t9*wght/rhob**4-3.4687652 2188 : 61708288d+0*gammabb*t12*t14*t9*wght/rhob**3.3333333333333 2189 ; 337d+0+3.333333333333333d-1*t14*t27*(t51+t50-7*t43)+t29*t 2190 < 5*(t51+t50-2*t41)-7.777777777777778d-1*t14*t27*t42+Amat3( 2191 = iq,D3_RB_RB_RB) 2192 Cmat3(iq,D3_RB_RB_GBB) = t1*t5*(1.5874010519681997d+0*gammab 2193 1 b*t12*t25*t46*t9*wght+7.937005259840999d-1*gammabb*t12*(- 2194 2 8.129082500000001d-4*t11*t38*t39*t45-1.5016221840277785d- 2195 3 3*t10*t15*t37*t45+4.0645412500000005d-4*t11*t24*t36*t45-9 2196 4 .483929583333336d-4*rhob*t22*t23*t24*t45+4.74196479166666 2197 5 8d-4*t10*t23*t45)*t9*wght+7.937005259840999d-1*t12*t40*t9 2198 6 *wght+7.937005259840999d-1*t12*t30*t9*wght+1.587401051968 2199 7 1997d+0*gammabb*t12*t22*t46*wght-3.226028528009071d-4*t11 2200 8 *t12*t2*t40*wght-3.226028528009071d-4*t10*t11*t12*t2*t36* 2201 9 wght+7.937005259840999d-1*t12*t36*wght-3.226028528009071d 2202 : -4*t11*t12*t2*t30*wght-6.452057056018142d-4*t10*t11*t12*t 2203 ; 2*t22*t25*wght+1.5874010519681997d+0*t12*t22*t25*wght)+t1 2204 < 3*t14*(-2.6456684199469993d-1*gammabb*t12*t46*t9*wght-2.6 2205 = 456684199469993d-1*t12*t25*t9*wght+1.0753428426696904d-4* 2206 > t11*t12*t2*t25*wght+1.0753428426696904d-4*t10*t11*t12*t2* 2207 ? t22*wght-2.6456684199469993d-1*t12*t22*wght)+3.3333333333 2208 @ 33333d-1*t13*t14*(-7.937005259840999d-1*gammabb*t12*t46*t 2209 1 9*wght-7.937005259840999d-1*t12*t25*t9*wght+3.22602852800 2210 2 9071d-4*t11*t12*t2*t25*wght+3.226028528009071d-4*t10*t11* 2211 3 t12*t2*t22*wght-7.937005259840999d-1*t12*t22*wght)+t16*t5 2212 4 *(-1.5874010519681997d+0*gammabb*t12*t46*t9*wght-1.587401 2213 5 0519681997d+0*t12*t25*t9*wght+6.452057056018142d-4*t11*t1 2214 6 2*t2*t25*wght+6.452057056018142d-4*t10*t11*t12*t2*t22*wgh 2215 7 t-1.5874010519681997d+0*t12*t22*wght)+1.5874010519681997d 2216 8 +0*t12*t29*t5*t9*wght+8.818894733156665d-1*t12*t14*t27*t9 2217 9 *wght-3.5844761422323007d-4*t12*t2*t27*t44*wght-6.4520570 2218 : 56018142d-4*t12*t2*t26*t29*wght+Cmat3(iq,D3_RB_RB_GBB) 2219 Cmat3(iq,D3_RB_GBB_GBB) = t1*t5*(7.937005259840999d-1*gammab 2220 1 b*t12*(-2.0322706250000003d-4*t11*t22*t24*t53-2.370982395 2221 2 833334d-4*rhob*t10*t23*t53)*t9*wght+1.5874010519681997d+0 2222 3 *t12*t46*t9*wght+1.311232602576965d-7*t10*t12*t25*t52*wgh 2223 4 t+1.311232602576965d-7*t12*t22*t24*t52*wght-6.45205705601 2224 5 8142d-4*t11*t12*t2*t46*wght-4.8390427920136064d-4*t11*t12 2225 6 *t25*t45*wght-4.8390427920136064d-4*t10*t11*t12*t22*t45*w 2226 7 ght)-1.311232602576965d-7*t10*t12*t16*t47*wght+1.61301426 2227 8 40045356d-4*t12*t13*t44*t45*wght+4.8390427920136064d-4*t1 2228 9 2*t16*t26*t45*wght-4.370775341923217d-8*t10*t12*t13*wght/ 2229 : t4**3.3333333333333337d+0+Cmat3(iq,D3_RB_GBB_GBB) 2230 Cmat3(iq,D3_GBB_GBB_GBB) = 2.4195213960068032d-4*t1*t12*t26* 2231 1 t53*wght+1.9668489038654477d-7*t1*t10*t12*t47*wght/gammab 2232 2 b-5.329559001518931d-11*t1*t12*t24*t45*wght/t4**3.8333333 2233 3 33333334d+0+Cmat3(iq,D3_GBB_GBB_GBB) 2234 endif ! rhoa.gt.tol_rho.and.rhob.gt.tol_rho 2235 endif ! ipol.eq.1 2236 enddo ! iq 2237 end 2238C> @} 2239