Lines Matching refs:t
90 real t[7],zk; in pz81a_energy() local
94 t[1] = rhob+rhoa; in pz81a_energy()
95 t[2] = 1/POW(t[1],0.33333333333333); in pz81a_energy()
96 t[3] = 1/POW(t[1],0.16666666666667); in pz81a_energy()
97 t[4] = 1/(0.78762331789974*b1u*t[3]+0.6203504908994*b2u*t[2]+1.0); in pz81a_energy()
98 t[5] = rhoa-1.0*rhob; in pz81a_energy()
99 t[6] = 1/t[1]; in pz81a_energy()
100 …t[1]*((POW(t[5]*t[6]+1.0,1.333333333333333)+POW(1.0-1.0*t[5]*t[6],1.333333333333333)-2.0)*(gp/(0.7… in pz81a_energy()
107 real t[28]; in pz81a_first() local
112 t[1] = rhob+rhoa; in pz81a_first()
113 t[2] = 1/POW(t[1],0.33333333333333); in pz81a_first()
114 t[3] = 1/POW(t[1],0.16666666666667); in pz81a_first()
115 t[4] = 0.78762331789974*b1u*t[3]+0.6203504908994*b2u*t[2]+1.0; in pz81a_first()
116 t[5] = 1/t[4]; in pz81a_first()
117 t[6] = gu*t[5]; in pz81a_first()
118 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81a_first()
119 t[8] = rhoa-1.0*rhob; in pz81a_first()
120 t[9] = 1/t[1]; in pz81a_first()
121 t[10] = 1.0-1.0*t[8]*t[9]; in pz81a_first()
122 t[11] = t[8]*t[9]+1.0; in pz81a_first()
123 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81a_first()
124 t[13] = 0.78762331789974*b1p*t[3]+0.6203504908994*b2p*t[2]+1.0; in pz81a_first()
125 t[14] = gp/t[13]-1.0*t[5]*gu; in pz81a_first()
126 t[15] = t[7]*t[12]*t[14]; in pz81a_first()
127 t[16] = 1/POW(t[1],1.333333333333333); in pz81a_first()
128 t[17] = 1/POW(t[1],1.166666666666667); in pz81a_first()
129 t[18] = -0.13127055298329*b1u*t[17]-0.20678349696647*b2u*t[16]; in pz81a_first()
130 t[19] = 1/POW(t[4],2.0); in pz81a_first()
131 t[20] = -1.0*t[18]*t[19]*gu; in pz81a_first()
132 …t[21] = t[12]*t[7]*(gu*t[18]*t[19]-1.0*(-0.13127055298329*b1p*t[17]-0.20678349696647*b2p*t[16])*gp… in pz81a_first()
133 t[22] = 1/POW(t[1],2.0); in pz81a_first()
134 t[23] = t[8]*t[22]; in pz81a_first()
135 t[24] = -1.0*t[9]; in pz81a_first()
136 t[25] = POW(t[10],0.33333333333333); in pz81a_first()
137 t[26] = -1.0*t[22]*t[8]; in pz81a_first()
138 t[27] = POW(t[11],0.33333333333333); in pz81a_first()
139 …dfdra = t[1]*(t[14]*t[7]*(1.333333333333333*t[27]*(t[9]+t[26])+1.333333333333333*(t[24]+t[23])*t[2… in pz81a_first()
140 …dfdrb = t[1]*(t[14]*t[7]*(1.333333333333333*t[25]*(t[9]+t[23])+1.333333333333333*(t[24]+t[26])*t[2… in pz81a_first()
148 real t[56]; in pz81a_second() local
154 t[1] = rhob+rhoa; in pz81a_second()
155 t[2] = 1/POW(t[1],0.33333333333333); in pz81a_second()
156 t[3] = 1/POW(t[1],0.16666666666667); in pz81a_second()
157 t[4] = 0.78762331789974*b1u*t[3]+0.6203504908994*b2u*t[2]+1.0; in pz81a_second()
158 t[5] = 1/t[4]; in pz81a_second()
159 t[6] = gu*t[5]; in pz81a_second()
160 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81a_second()
161 t[8] = rhoa-1.0*rhob; in pz81a_second()
162 t[9] = 1/t[1]; in pz81a_second()
163 t[10] = 1.0-1.0*t[8]*t[9]; in pz81a_second()
164 t[11] = t[8]*t[9]+1.0; in pz81a_second()
165 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81a_second()
166 t[13] = 0.78762331789974*b1p*t[3]+0.6203504908994*b2p*t[2]+1.0; in pz81a_second()
167 t[14] = gp/t[13]-1.0*t[5]*gu; in pz81a_second()
168 t[15] = t[7]*t[12]*t[14]; in pz81a_second()
169 t[16] = 1/POW(t[1],1.333333333333333); in pz81a_second()
170 t[17] = 1/POW(t[1],1.166666666666667); in pz81a_second()
171 t[18] = -0.13127055298329*b1u*t[17]-0.20678349696647*b2u*t[16]; in pz81a_second()
172 t[19] = 1/POW(t[4],2.0); in pz81a_second()
173 t[20] = -1.0*t[18]*t[19]*gu; in pz81a_second()
174 t[21] = -0.13127055298329*b1p*t[17]-0.20678349696647*b2p*t[16]; in pz81a_second()
175 t[22] = 1/POW(t[13],2.0); in pz81a_second()
176 t[23] = gu*t[18]*t[19]-1.0*t[21]*t[22]*gp; in pz81a_second()
177 t[24] = t[7]*t[12]*t[23]; in pz81a_second()
178 t[25] = 1/POW(t[1],2.0); in pz81a_second()
179 t[26] = t[8]*t[25]; in pz81a_second()
180 t[27] = -1.0*t[9]; in pz81a_second()
181 t[28] = t[27]+t[26]; in pz81a_second()
182 t[29] = POW(t[10],0.33333333333333); in pz81a_second()
183 t[30] = -1.0*t[25]*t[8]; in pz81a_second()
184 t[31] = t[9]+t[30]; in pz81a_second()
185 t[32] = POW(t[11],0.33333333333333); in pz81a_second()
186 t[33] = 1.333333333333333*t[31]*t[32]+1.333333333333333*t[28]*t[29]; in pz81a_second()
187 t[34] = t[7]*t[33]*t[14]; in pz81a_second()
188 t[35] = t[9]+t[26]; in pz81a_second()
189 t[36] = t[27]+t[30]; in pz81a_second()
190 t[37] = 1.333333333333333*t[32]*t[36]+1.333333333333333*t[29]*t[35]; in pz81a_second()
191 t[38] = t[7]*t[37]*t[14]; in pz81a_second()
192 t[39] = -2.0*t[18]*t[19]*gu; in pz81a_second()
193 t[40] = 2.0*t[12]*t[23]*t[7]; in pz81a_second()
194 t[41] = POW(t[18],2.0); in pz81a_second()
195 t[42] = 1/POW(t[4],3.0); in pz81a_second()
196 t[43] = 2.0*t[41]*t[42]*gu; in pz81a_second()
197 t[44] = 1/POW(t[1],2.333333333333334); in pz81a_second()
198 t[45] = 1/POW(t[1],2.166666666666667); in pz81a_second()
199 t[46] = 0.15314897848051*b1u*t[45]+0.27571132928862*b2u*t[44]; in pz81a_second()
200 t[47] = -1.0*t[19]*t[46]*gu; in pz81a_second()
201 …t[48] = t[12]*t[7]*(-2.0*t[41]*t[42]*gu-1.0*t[22]*(0.15314897848051*b1p*t[45]+0.27571132928862*b2p… in pz81a_second()
202 t[49] = 1/POW(t[10],0.66666666666667); in pz81a_second()
203 t[50] = 1/POW(t[1],3.0); in pz81a_second()
204 t[51] = -2.0*t[50]*t[8]; in pz81a_second()
205 t[52] = 2.0*t[25]; in pz81a_second()
206 t[53] = 1/POW(t[11],0.66666666666667); in pz81a_second()
207 t[54] = 2.0*t[50]*t[8]; in pz81a_second()
208 t[55] = -2.0*t[25]; in pz81a_second()
209 dfdra = t[1]*(t[34]+t[24]+t[20])+t[15]+t[6]; in pz81a_second()
210 dfdrb = t[1]*(t[38]+t[24]+t[20])+t[15]+t[6]; in pz81a_second()
211 …t[1]*(t[14]*(1.333333333333333*t[32]*(t[55]+t[54])+0.44444444444444*POW(t[31],2.0)*t[53]+1.3333333… in pz81a_second()
212 …t[1]*(t[14]*t[7]*(2.666666666666667*t[32]*t[50]*t[8]-2.666666666666667*t[29]*t[50]*t[8]+0.44444444… in pz81a_second()
213 …t[1]*(t[14]*(1.333333333333333*t[29]*(t[55]+t[51])+0.44444444444444*POW(t[36],2.0)*t[53]+1.3333333… in pz81a_second()
224 real t[95]; in pz81a_third() local
231 t[1] = rhob+rhoa; in pz81a_third()
232 t[2] = 1/POW(t[1],0.33333333333333); in pz81a_third()
233 t[3] = 1/POW(t[1],0.16666666666667); in pz81a_third()
234 t[4] = 0.78762331789974*b1u*t[3]+0.6203504908994*b2u*t[2]+1.0; in pz81a_third()
235 t[5] = 1/t[4]; in pz81a_third()
236 t[6] = gu*t[5]; in pz81a_third()
237 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81a_third()
238 t[8] = rhoa-1.0*rhob; in pz81a_third()
239 t[9] = 1/t[1]; in pz81a_third()
240 t[10] = 1.0-1.0*t[8]*t[9]; in pz81a_third()
241 t[11] = t[8]*t[9]+1.0; in pz81a_third()
242 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81a_third()
243 t[13] = 0.78762331789974*b1p*t[3]+0.6203504908994*b2p*t[2]+1.0; in pz81a_third()
244 t[14] = gp/t[13]-1.0*t[5]*gu; in pz81a_third()
245 t[15] = t[7]*t[12]*t[14]; in pz81a_third()
246 t[16] = 1/POW(t[1],1.333333333333333); in pz81a_third()
247 t[17] = 1/POW(t[1],1.166666666666667); in pz81a_third()
248 t[18] = -0.13127055298329*b1u*t[17]-0.20678349696647*b2u*t[16]; in pz81a_third()
249 t[19] = 1/POW(t[4],2.0); in pz81a_third()
250 t[20] = -1.0*t[18]*t[19]*gu; in pz81a_third()
251 t[21] = -0.13127055298329*b1p*t[17]-0.20678349696647*b2p*t[16]; in pz81a_third()
252 t[22] = 1/POW(t[13],2.0); in pz81a_third()
253 t[23] = gu*t[18]*t[19]-1.0*t[21]*t[22]*gp; in pz81a_third()
254 t[24] = t[7]*t[12]*t[23]; in pz81a_third()
255 t[25] = 1/POW(t[1],2.0); in pz81a_third()
256 t[26] = t[8]*t[25]; in pz81a_third()
257 t[27] = -1.0*t[9]; in pz81a_third()
258 t[28] = t[27]+t[26]; in pz81a_third()
259 t[29] = POW(t[10],0.33333333333333); in pz81a_third()
260 t[30] = -1.0*t[25]*t[8]; in pz81a_third()
261 t[31] = t[9]+t[30]; in pz81a_third()
262 t[32] = POW(t[11],0.33333333333333); in pz81a_third()
263 t[33] = 1.333333333333333*t[31]*t[32]+1.333333333333333*t[28]*t[29]; in pz81a_third()
264 t[34] = t[7]*t[33]*t[14]; in pz81a_third()
265 t[35] = t[9]+t[26]; in pz81a_third()
266 t[36] = t[27]+t[30]; in pz81a_third()
267 t[37] = 1.333333333333333*t[32]*t[36]+1.333333333333333*t[29]*t[35]; in pz81a_third()
268 t[38] = t[7]*t[37]*t[14]; in pz81a_third()
269 t[39] = -2.0*t[18]*t[19]*gu; in pz81a_third()
270 t[40] = 2.0*t[12]*t[23]*t[7]; in pz81a_third()
271 t[41] = POW(t[18],2.0); in pz81a_third()
272 t[42] = 1/POW(t[4],3.0); in pz81a_third()
273 t[43] = 2.0*t[41]*t[42]*gu; in pz81a_third()
274 t[44] = 1/POW(t[1],2.333333333333334); in pz81a_third()
275 t[45] = 1/POW(t[1],2.166666666666667); in pz81a_third()
276 t[46] = 0.15314897848051*b1u*t[45]+0.27571132928862*b2u*t[44]; in pz81a_third()
277 t[47] = -1.0*t[19]*t[46]*gu; in pz81a_third()
278 t[48] = 1/POW(t[13],3.0); in pz81a_third()
279 t[49] = 0.15314897848051*b1p*t[45]+0.27571132928862*b2p*t[44]; in pz81a_third()
280 t[50] = -2.0*t[41]*t[42]*gu-1.0*t[22]*t[49]*gp+2.0*POW(t[21],2.0)*t[48]*gp+gu*t[46]*t[19]; in pz81a_third()
281 t[51] = t[7]*t[12]*t[50]; in pz81a_third()
282 t[52] = 2.0*t[23]*t[33]*t[7]; in pz81a_third()
283 t[53] = POW(t[28],2.0); in pz81a_third()
284 t[54] = 1/POW(t[10],0.66666666666667); in pz81a_third()
285 t[55] = 1/POW(t[1],3.0); in pz81a_third()
286 t[56] = -2.0*t[55]*t[8]; in pz81a_third()
287 t[57] = 2.0*t[25]; in pz81a_third()
288 t[58] = t[57]+t[56]; in pz81a_third()
289 t[59] = POW(t[31],2.0); in pz81a_third()
290 t[60] = 1/POW(t[11],0.66666666666667); in pz81a_third()
291 t[61] = 2.0*t[55]*t[8]; in pz81a_third()
292 t[62] = -2.0*t[25]; in pz81a_third()
293 t[63] = t[62]+t[61]; in pz81a_third()
294 …t[64] = 1.333333333333333*t[32]*t[63]+0.44444444444444*t[59]*t[60]+1.333333333333333*t[29]*t[58]+0… in pz81a_third()
295 t[65] = t[7]*t[64]*t[14]; in pz81a_third()
296 …t[66] = 2.666666666666667*t[32]*t[55]*t[8]-2.666666666666667*t[29]*t[55]*t[8]+0.44444444444444*t[3… in pz81a_third()
297 t[67] = 2.0*t[23]*t[37]*t[7]; in pz81a_third()
298 t[68] = POW(t[35],2.0); in pz81a_third()
299 t[69] = t[62]+t[56]; in pz81a_third()
300 t[70] = POW(t[36],2.0); in pz81a_third()
301 t[71] = t[57]+t[61]; in pz81a_third()
302 …t[72] = 1.333333333333333*t[32]*t[71]+0.44444444444444*t[60]*t[70]+1.333333333333333*t[29]*t[69]+0… in pz81a_third()
303 t[73] = t[7]*t[72]*t[14]; in pz81a_third()
304 t[74] = 6.0*t[41]*t[42]*gu; in pz81a_third()
305 t[75] = -3.0*t[19]*t[46]*gu; in pz81a_third()
306 t[76] = 3.0*t[12]*t[50]*t[7]; in pz81a_third()
307 t[77] = 2.0*t[14]*t[66]*t[7]; in pz81a_third()
308 t[78] = POW(t[18],3.0); in pz81a_third()
309 t[79] = 1/POW(t[4],4.0); in pz81a_third()
310 t[80] = -6.0*t[78]*t[79]*gu; in pz81a_third()
311 t[81] = 6.0*t[18]*t[42]*t[46]*gu; in pz81a_third()
312 t[82] = 1/POW(t[1],3.333333333333334); in pz81a_third()
313 t[83] = 1/POW(t[1],3.166666666666667); in pz81a_third()
314 t[84] = -0.33182278670776*b1u*t[83]-0.64332643500679*b2u*t[82]; in pz81a_third()
315 t[85] = -1.0*t[19]*t[84]*gu; in pz81a_third()
316 …t[86] = t[12]*t[7]*(6.0*t[78]*t[79]*gu-6.0*t[18]*t[42]*t[46]*gu-1.0*t[22]*(-0.33182278670776*b1p*t… in pz81a_third()
317 t[87] = 2.0*t[23]*t[66]*t[7]; in pz81a_third()
318 t[88] = 1/POW(t[10],1.666666666666667); in pz81a_third()
319 t[89] = 1/POW(t[1],4.0); in pz81a_third()
320 t[90] = 6.0*t[8]*t[89]; in pz81a_third()
321 t[91] = 1/POW(t[11],1.666666666666667); in pz81a_third()
322 t[92] = -6.0*t[8]*t[89]; in pz81a_third()
323 t[93] = -6.0*t[55]; in pz81a_third()
324 t[94] = 6.0*t[55]; in pz81a_third()
325 dfdra = t[1]*(t[34]+t[24]+t[20])+t[15]+t[6]; in pz81a_third()
326 dfdrb = t[1]*(t[38]+t[24]+t[20])+t[15]+t[6]; in pz81a_third()
327 d2fdrara = 2.0*t[14]*t[33]*t[7]+t[1]*(t[65]+t[52]+t[51]+t[47]+t[43])+t[40]+t[39]; in pz81a_third()
328 …d2fdrarb = t[1]*(t[7]*t[66]*t[14]+t[7]*t[33]*t[23]+t[7]*t[37]*t[23]+t[51]+t[47]+t[43])+t[34]+t[38]… in pz81a_third()
329 d2fdrbrb = t[1]*(t[73]+t[67]+t[51]+t[47]+t[43])+2.0*t[14]*t[37]*t[7]+t[40]+t[39]; in pz81a_third()
330 …t[1]*(t[14]*t[7]*(1.333333333333333*t[32]*(t[92]+2.0*t[55])-0.2962962962963*t[36]*t[59]*t[91]+1.33… in pz81a_third()
331 …t[1]*(t[14]*t[7]*(-0.2962962962963*t[31]*t[70]*t[91]-8.0*t[32]*t[8]*t[89]+8.0*t[29]*t[8]*t[89]-0.2… in pz81a_third()
332 …t[1]*(t[14]*t[7]*(1.333333333333333*t[32]*(t[94]+t[92])+1.333333333333333*t[29]*(t[93]+t[90])-0.29… in pz81a_third()
333 …t[1]*(t[14]*t[7]*(1.333333333333333*t[29]*(t[94]+t[90])+1.333333333333333*t[32]*(t[93]+t[92])-0.29… in pz81a_third()
354 real t[6],zk; in pz81b_energy() local
358 t[1] = rhob+rhoa; in pz81b_energy()
359 t[2] = 1/POW(t[1],0.33333333333333); in pz81b_energy()
360 t[3] = LOG(0.6203504908994*t[2]); in pz81b_energy()
361 t[4] = rhoa-1.0*rhob; in pz81b_energy()
362 t[5] = 1/t[1]; in pz81b_energy()
363 …t[1]*(Bu+(POW(t[4]*t[5]+1.0,1.333333333333333)+POW(1.0-1.0*t[4]*t[5],1.333333333333333)-2.0)*(-1.0… in pz81b_energy()
370 real t[27]; in pz81b_first() local
375 t[1] = rhob+rhoa; in pz81b_first()
376 t[2] = 1/POW(t[1],0.33333333333333); in pz81b_first()
377 t[3] = 0.6203504908994*Du*t[2]; in pz81b_first()
378 t[4] = LOG(0.6203504908994*t[2]); in pz81b_first()
379 t[5] = Au*t[4]; in pz81b_first()
380 t[6] = 0.6203504908994*Cu*t[2]*t[4]; in pz81b_first()
381 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81b_first()
382 t[8] = rhoa-1.0*rhob; in pz81b_first()
383 t[9] = 1/t[1]; in pz81b_first()
384 t[10] = 1.0-1.0*t[8]*t[9]; in pz81b_first()
385 t[11] = t[8]*t[9]+1.0; in pz81b_first()
386 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81b_first()
387 …t[13] = -1.0*Bu+Bp-1.0*t[4]*Au+Ap*t[4]-0.6203504908994*Cu*t[2]*t[4]+0.6203504908994*Cp*t[2]*t[4]-0… in pz81b_first()
388 t[14] = t[7]*t[12]*t[13]; in pz81b_first()
389 t[15] = 1/POW(t[1],1.333333333333333); in pz81b_first()
390 t[16] = -0.20678349696647*Cu*t[15]; in pz81b_first()
391 t[17] = -0.20678349696647*Du*t[15]; in pz81b_first()
392 t[18] = -0.33333333333333*Au*t[9]; in pz81b_first()
393 t[19] = -0.20678349696647*Cu*t[15]*t[4]; in pz81b_first()
394 …t[20] = t[7]*t[12]*(0.20678349696647*Cu*t[15]*t[4]-0.20678349696647*Cp*t[15]*t[4]+0.33333333333333… in pz81b_first()
395 t[21] = 1/POW(t[1],2.0); in pz81b_first()
396 t[22] = t[8]*t[21]; in pz81b_first()
397 t[23] = -1.0*t[9]; in pz81b_first()
398 t[24] = POW(t[10],0.33333333333333); in pz81b_first()
399 t[25] = -1.0*t[21]*t[8]; in pz81b_first()
400 t[26] = POW(t[11],0.33333333333333); in pz81b_first()
401 … = Bu+t[1]*(t[13]*t[7]*(1.333333333333333*t[26]*(t[9]+t[25])+1.333333333333333*(t[23]+t[22])*t[24]… in pz81b_first()
402 … = Bu+t[1]*(t[13]*t[7]*(1.333333333333333*t[24]*(t[9]+t[22])+1.333333333333333*(t[23]+t[25])*t[26]… in pz81b_first()
410 real t[54]; in pz81b_second() local
416 t[1] = rhob+rhoa; in pz81b_second()
417 t[2] = 1/POW(t[1],0.33333333333333); in pz81b_second()
418 t[3] = 0.6203504908994*Du*t[2]; in pz81b_second()
419 t[4] = LOG(0.6203504908994*t[2]); in pz81b_second()
420 t[5] = Au*t[4]; in pz81b_second()
421 t[6] = 0.6203504908994*Cu*t[2]*t[4]; in pz81b_second()
422 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81b_second()
423 t[8] = rhoa-1.0*rhob; in pz81b_second()
424 t[9] = 1/t[1]; in pz81b_second()
425 t[10] = 1.0-1.0*t[8]*t[9]; in pz81b_second()
426 t[11] = t[8]*t[9]+1.0; in pz81b_second()
427 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81b_second()
428 …t[13] = -1.0*Bu+Bp-1.0*t[4]*Au+Ap*t[4]-0.6203504908994*Cu*t[2]*t[4]+0.6203504908994*Cp*t[2]*t[4]-0… in pz81b_second()
429 t[14] = t[7]*t[12]*t[13]; in pz81b_second()
430 t[15] = 1/POW(t[1],1.333333333333333); in pz81b_second()
431 t[16] = -0.20678349696647*Cu*t[15]; in pz81b_second()
432 t[17] = -0.20678349696647*Du*t[15]; in pz81b_second()
433 t[18] = -0.33333333333333*Au*t[9]; in pz81b_second()
434 t[19] = -0.20678349696647*Cu*t[15]*t[4]; in pz81b_second()
435 …t[20] = 0.20678349696647*Cu*t[15]*t[4]-0.20678349696647*Cp*t[15]*t[4]+0.33333333333333*Au*t[9]-0.3… in pz81b_second()
436 t[21] = t[7]*t[12]*t[20]; in pz81b_second()
437 t[22] = 1/POW(t[1],2.0); in pz81b_second()
438 t[23] = t[8]*t[22]; in pz81b_second()
439 t[24] = -1.0*t[9]; in pz81b_second()
440 t[25] = t[24]+t[23]; in pz81b_second()
441 t[26] = POW(t[10],0.33333333333333); in pz81b_second()
442 t[27] = -1.0*t[22]*t[8]; in pz81b_second()
443 t[28] = t[9]+t[27]; in pz81b_second()
444 t[29] = POW(t[11],0.33333333333333); in pz81b_second()
445 t[30] = 1.333333333333333*t[28]*t[29]+1.333333333333333*t[25]*t[26]; in pz81b_second()
446 t[31] = t[7]*t[30]*t[13]; in pz81b_second()
447 t[32] = t[9]+t[23]; in pz81b_second()
448 t[33] = t[24]+t[27]; in pz81b_second()
449 t[34] = 1.333333333333333*t[29]*t[33]+1.333333333333333*t[26]*t[32]; in pz81b_second()
450 t[35] = t[7]*t[34]*t[13]; in pz81b_second()
451 t[36] = -0.41356699393293*Cu*t[15]; in pz81b_second()
452 t[37] = -0.41356699393293*Du*t[15]; in pz81b_second()
453 t[38] = -0.66666666666667*Au*t[9]; in pz81b_second()
454 t[39] = -0.41356699393293*Cu*t[15]*t[4]; in pz81b_second()
455 t[40] = 2.0*t[12]*t[20]*t[7]; in pz81b_second()
456 t[41] = 1/POW(t[1],2.333333333333334); in pz81b_second()
457 t[42] = 0.34463916161078*Cu*t[41]; in pz81b_second()
458 t[43] = 0.27571132928862*Du*t[41]; in pz81b_second()
459 t[44] = 0.33333333333333*Au*t[22]; in pz81b_second()
460 t[45] = 0.27571132928862*Cu*t[41]*t[4]; in pz81b_second()
461 …t[46] = t[7]*t[12]*(-0.27571132928862*Cu*t[41]*t[4]+0.27571132928862*Cp*t[41]*t[4]-0.3333333333333… in pz81b_second()
462 t[47] = 1/POW(t[10],0.66666666666667); in pz81b_second()
463 t[48] = 1/POW(t[1],3.0); in pz81b_second()
464 t[49] = -2.0*t[48]*t[8]; in pz81b_second()
465 t[50] = 2.0*t[22]; in pz81b_second()
466 t[51] = 1/POW(t[11],0.66666666666667); in pz81b_second()
467 t[52] = 2.0*t[48]*t[8]; in pz81b_second()
468 t[53] = -2.0*t[22]; in pz81b_second()
469 dfdra = t[1]*(t[31]+t[21]+t[19]+t[18]+t[17]+t[16])+t[14]+t[6]+t[5]+t[3]+Bu; in pz81b_second()
470 dfdrb = t[1]*(t[35]+t[21]+t[19]+t[18]+t[17]+t[16])+t[14]+t[6]+t[5]+t[3]+Bu; in pz81b_second()
471 …t[1]*(t[13]*(1.333333333333333*t[29]*(t[53]+t[52])+0.44444444444444*POW(t[28],2.0)*t[51]+1.3333333… in pz81b_second()
472 …t[1]*(t[13]*t[7]*(2.666666666666667*t[29]*t[48]*t[8]-2.666666666666667*t[26]*t[48]*t[8]+0.44444444… in pz81b_second()
473 …t[1]*(t[13]*(1.333333333333333*t[26]*(t[53]+t[49])+0.44444444444444*POW(t[33],2.0)*t[51]+1.3333333… in pz81b_second()
484 real t[90]; in pz81b_third() local
491 t[1] = rhob+rhoa; in pz81b_third()
492 t[2] = 1/POW(t[1],0.33333333333333); in pz81b_third()
493 t[3] = 0.6203504908994*Du*t[2]; in pz81b_third()
494 t[4] = LOG(0.6203504908994*t[2]); in pz81b_third()
495 t[5] = Au*t[4]; in pz81b_third()
496 t[6] = 0.6203504908994*Cu*t[2]*t[4]; in pz81b_third()
497 t[7] = 1/(2.0*POW(2.0,0.33333333333333)-2.0); in pz81b_third()
498 t[8] = rhoa-1.0*rhob; in pz81b_third()
499 t[9] = 1/t[1]; in pz81b_third()
500 t[10] = 1.0-1.0*t[8]*t[9]; in pz81b_third()
501 t[11] = t[8]*t[9]+1.0; in pz81b_third()
502 t[12] = POW(t[11],1.333333333333333)+POW(t[10],1.333333333333333)-2.0; in pz81b_third()
503 …t[13] = -1.0*Bu+Bp-1.0*t[4]*Au+Ap*t[4]-0.6203504908994*Cu*t[2]*t[4]+0.6203504908994*Cp*t[2]*t[4]-0… in pz81b_third()
504 t[14] = t[7]*t[12]*t[13]; in pz81b_third()
505 t[15] = 1/POW(t[1],1.333333333333333); in pz81b_third()
506 t[16] = -0.20678349696647*Cu*t[15]; in pz81b_third()
507 t[17] = -0.20678349696647*Du*t[15]; in pz81b_third()
508 t[18] = -0.33333333333333*Au*t[9]; in pz81b_third()
509 t[19] = -0.20678349696647*Cu*t[15]*t[4]; in pz81b_third()
510 …t[20] = 0.20678349696647*Cu*t[15]*t[4]-0.20678349696647*Cp*t[15]*t[4]+0.33333333333333*Au*t[9]-0.3… in pz81b_third()
511 t[21] = t[7]*t[12]*t[20]; in pz81b_third()
512 t[22] = 1/POW(t[1],2.0); in pz81b_third()
513 t[23] = t[8]*t[22]; in pz81b_third()
514 t[24] = -1.0*t[9]; in pz81b_third()
515 t[25] = t[24]+t[23]; in pz81b_third()
516 t[26] = POW(t[10],0.33333333333333); in pz81b_third()
517 t[27] = -1.0*t[22]*t[8]; in pz81b_third()
518 t[28] = t[9]+t[27]; in pz81b_third()
519 t[29] = POW(t[11],0.33333333333333); in pz81b_third()
520 t[30] = 1.333333333333333*t[28]*t[29]+1.333333333333333*t[25]*t[26]; in pz81b_third()
521 t[31] = t[7]*t[30]*t[13]; in pz81b_third()
522 t[32] = t[9]+t[23]; in pz81b_third()
523 t[33] = t[24]+t[27]; in pz81b_third()
524 t[34] = 1.333333333333333*t[29]*t[33]+1.333333333333333*t[26]*t[32]; in pz81b_third()
525 t[35] = t[7]*t[34]*t[13]; in pz81b_third()
526 t[36] = -0.41356699393293*Cu*t[15]; in pz81b_third()
527 t[37] = -0.41356699393293*Du*t[15]; in pz81b_third()
528 t[38] = -0.66666666666667*Au*t[9]; in pz81b_third()
529 t[39] = -0.41356699393293*Cu*t[15]*t[4]; in pz81b_third()
530 t[40] = 2.0*t[12]*t[20]*t[7]; in pz81b_third()
531 t[41] = 1/POW(t[1],2.333333333333334); in pz81b_third()
532 t[42] = 0.34463916161078*Cu*t[41]; in pz81b_third()
533 t[43] = 0.27571132928862*Du*t[41]; in pz81b_third()
534 t[44] = 0.33333333333333*Au*t[22]; in pz81b_third()
535 t[45] = 0.27571132928862*Cu*t[41]*t[4]; in pz81b_third()
536 …t[46] = -0.27571132928862*Cu*t[41]*t[4]+0.27571132928862*Cp*t[41]*t[4]-0.33333333333333*Au*t[22]+0… in pz81b_third()
537 t[47] = t[7]*t[12]*t[46]; in pz81b_third()
538 t[48] = 2.0*t[20]*t[30]*t[7]; in pz81b_third()
539 t[49] = POW(t[25],2.0); in pz81b_third()
540 t[50] = 1/POW(t[10],0.66666666666667); in pz81b_third()
541 t[51] = 1/POW(t[1],3.0); in pz81b_third()
542 t[52] = -2.0*t[51]*t[8]; in pz81b_third()
543 t[53] = 2.0*t[22]; in pz81b_third()
544 t[54] = t[53]+t[52]; in pz81b_third()
545 t[55] = POW(t[28],2.0); in pz81b_third()
546 t[56] = 1/POW(t[11],0.66666666666667); in pz81b_third()
547 t[57] = 2.0*t[51]*t[8]; in pz81b_third()
548 t[58] = -2.0*t[22]; in pz81b_third()
549 t[59] = t[58]+t[57]; in pz81b_third()
550 …t[60] = 1.333333333333333*t[29]*t[59]+0.44444444444444*t[55]*t[56]+1.333333333333333*t[26]*t[54]+0… in pz81b_third()
551 t[61] = t[7]*t[60]*t[13]; in pz81b_third()
552 …t[62] = 2.666666666666667*t[29]*t[51]*t[8]-2.666666666666667*t[26]*t[51]*t[8]+0.44444444444444*t[2… in pz81b_third()
553 t[63] = 2.0*t[20]*t[34]*t[7]; in pz81b_third()
554 t[64] = POW(t[32],2.0); in pz81b_third()
555 t[65] = t[58]+t[52]; in pz81b_third()
556 t[66] = POW(t[33],2.0); in pz81b_third()
557 t[67] = t[53]+t[57]; in pz81b_third()
558 …t[68] = 1.333333333333333*t[29]*t[67]+0.44444444444444*t[56]*t[66]+1.333333333333333*t[26]*t[65]+0… in pz81b_third()
559 t[69] = t[7]*t[68]*t[13]; in pz81b_third()
560 t[70] = 1.033917484832333*Cu*t[41]; in pz81b_third()
561 t[71] = 0.82713398786587*Du*t[41]; in pz81b_third()
562 t[72] = Au*t[22]; in pz81b_third()
563 t[73] = 0.82713398786587*Cu*t[41]*t[4]; in pz81b_third()
564 t[74] = 3.0*t[12]*t[46]*t[7]; in pz81b_third()
565 t[75] = 2.0*t[13]*t[62]*t[7]; in pz81b_third()
566 t[76] = 1/POW(t[1],3.333333333333334); in pz81b_third()
567 t[77] = -0.89606182018802*Cu*t[76]; in pz81b_third()
568 t[78] = -0.64332643500679*Du*t[76]; in pz81b_third()
569 t[79] = -0.66666666666667*Au*t[51]; in pz81b_third()
570 t[80] = -0.64332643500679*Cu*t[76]*t[4]; in pz81b_third()
571 …t[81] = t[7]*t[12]*(0.64332643500679*Cu*t[76]*t[4]-0.64332643500679*Cp*t[76]*t[4]+0.66666666666667… in pz81b_third()
572 t[82] = 2.0*t[20]*t[62]*t[7]; in pz81b_third()
573 t[83] = 1/POW(t[10],1.666666666666667); in pz81b_third()
574 t[84] = 1/POW(t[1],4.0); in pz81b_third()
575 t[85] = 6.0*t[8]*t[84]; in pz81b_third()
576 t[86] = 1/POW(t[11],1.666666666666667); in pz81b_third()
577 t[87] = -6.0*t[8]*t[84]; in pz81b_third()
578 t[88] = -6.0*t[51]; in pz81b_third()
579 t[89] = 6.0*t[51]; in pz81b_third()
580 dfdra = t[1]*(t[31]+t[21]+t[19]+t[18]+t[17]+t[16])+t[14]+t[6]+t[5]+t[3]+Bu; in pz81b_third()
581 dfdrb = t[1]*(t[35]+t[21]+t[19]+t[18]+t[17]+t[16])+t[14]+t[6]+t[5]+t[3]+Bu; in pz81b_third()
582 …d2fdrara = 2.0*t[13]*t[30]*t[7]+t[1]*(t[61]+t[48]+t[47]+t[45]+t[44]+t[43]+t[42])+t[40]+t[39]+t[38]… in pz81b_third()
583 …d2fdrarb = t[1]*(t[7]*t[62]*t[13]+t[7]*t[30]*t[20]+t[7]*t[34]*t[20]+t[47]+t[45]+t[44]+t[43]+t[42])… in pz81b_third()
584 …d2fdrbrb = 2.0*t[13]*t[34]*t[7]+t[1]*(t[69]+t[63]+t[47]+t[45]+t[44]+t[43]+t[42])+t[40]+t[39]+t[38]… in pz81b_third()
585 …t[1]*(t[13]*t[7]*(1.333333333333333*t[29]*(t[87]+2.0*t[51])-0.2962962962963*t[33]*t[55]*t[86]+1.33… in pz81b_third()
586 …t[1]*(t[13]*t[7]*(-0.2962962962963*t[28]*t[66]*t[86]-8.0*t[29]*t[8]*t[84]+8.0*t[26]*t[8]*t[84]-0.2… in pz81b_third()
587 …t[1]*(t[13]*t[7]*(1.333333333333333*t[29]*(t[89]+t[87])+1.333333333333333*t[26]*(t[88]+t[85])-0.29… in pz81b_third()
588 …t[1]*(t[13]*t[7]*(1.333333333333333*t[26]*(t[89]+t[85])+1.333333333333333*t[29]*(t[88]+t[87])-0.29… in pz81b_third()