LALSimulation  5.4.0.1-fe68b98
SEOBNRv2_opt_3derivs.h
Go to the documentation of this file.
1 REAL8 tmp12=sigmaKerr->data[2]*sigmaKerr->data[2];
3 REAL8 tmp21=coeffs->KK*eta;
5 REAL8 tmp14=1/x->data[0];
6 REAL8 tmp18=x->data[0]*x->data[0];
7 REAL8 tmp5=pow(x->data[0],-5.);
8 REAL8 tmp25=(1.0/(x->data[0]*x->data[0]));
9 REAL8 tmp30=pow(x->data[0],-6.);
10 REAL8 tmp33=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
11 REAL8 tmp35=1./(x->data[0]*x->data[0]*x->data[0]);
14 REAL8 tmp40=coeffs->k5l*tmp39;
15 REAL8 tmp41=coeffs->k5+tmp40;
22 REAL8 tmp44=1.*coeffs->k4*tmp33;
23 REAL8 tmp45=1.*coeffs->k3*tmp35;
24 REAL8 tmp46=1.*coeffs->k2*tmp25;
25 REAL8 tmp47=1.*coeffs->k1*tmp14;
28 REAL8 tmp52=coeffs->k0*eta;
33 REAL8 tmp7=1.*e3z;
35 REAL8 tmp9=0.+p->data[1];
42 REAL8 tmp16=2.*tmp13*x->data[0];
47 REAL8 tmp79=-3.*eta;
48 REAL8 tmp31=-(coeffs->k5l*tmp30);
49 REAL8 tmp32=-4.*coeffs->k4*tmp5;
50 REAL8 tmp34=-3.*coeffs->k3*tmp33;
51 REAL8 tmp36=-2.*coeffs->k2*tmp35;
52 REAL8 tmp37=-(coeffs->k1*tmp25);
56 REAL8 tmp85=0.+p->data[0];
62 REAL8 tmp91=6.*eta*tmp25;
69 REAL8 tmp20=4.*tmp19*x->data[0];
71 REAL8 tmp56=-2.*tmp12*tmp29*tmp55*x->data[0];
76 REAL8 tmp101=(x->data[0]*x->data[0]*x->data[0]);
78 REAL8 tmp81=1.*p->data[0];
90 REAL8 tmp97=-12.*eta*tmp35;
93 REAL8 tmp133=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
106 REAL8 tmp159=pow(x->data[0],6.);
107 REAL8 tmp162=eta*eta;
110 REAL8 tmp171=103.*eta;
113 REAL8 tmp176=3.*eta;
115 REAL8 tmp181=-16.*eta;
118 REAL8 tmp188=-47.*eta;
123 REAL8 tmp156=(eta*eta*eta);
124 REAL8 tmp165=pow(x->data[0],7.);
126 REAL8 tmp168=pow(x->data[0],5.);
127 REAL8 tmp213=-109.*eta;
130 REAL8 tmp220=-6.*eta;
133 REAL8 tmp227=16.*eta;
138 REAL8 tmp128=-3.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp33;
139 REAL8 tmp129=-36.*eta*sigmaKerr->data[2]*tmp101*tmp29*tmp43*tmp50*tmp75*tmp88*tmp94;
140 REAL8 tmp130=-30.*eta*sigmaStar->data[2]*tmp101*tmp29*tmp43*tmp50*tmp75*tmp88*tmp94;
141 REAL8 tmp131=-36.*sigmaKerr->data[2]*tmp101*tmp29*tmp55*tmp75*tmp88*tmp98*tmp99;
142 REAL8 tmp132=-30.*sigmaStar->data[2]*tmp101*tmp29*tmp55*tmp75*tmp88*tmp98*tmp99;
143 REAL8 tmp134=72.*sigmaKerr->data[2]*tmp102*tmp133*tmp29*tmp55*tmp88*tmp94;
144 REAL8 tmp135=60.*sigmaStar->data[2]*tmp102*tmp133*tmp29*tmp55*tmp88*tmp94;
145 REAL8 tmp136=-108.*sigmaKerr->data[2]*tmp18*tmp29*tmp55*tmp75*tmp88*tmp94;
146 REAL8 tmp137=-90.*sigmaStar->data[2]*tmp18*tmp29*tmp55*tmp75*tmp88*tmp94;
147 REAL8 tmp138=-36.*sigmaKerr->data[2]*tmp101*tmp55*tmp59*tmp75*tmp88*tmp94;
148 REAL8 tmp139=-30.*sigmaStar->data[2]*tmp101*tmp55*tmp59*tmp75*tmp88*tmp94;
149 REAL8 tmp141=3.*sigmaKerr->data[2]*tmp140*x->data[0];
150 REAL8 tmp142=4.*sigmaStar->data[2]*tmp140*x->data[0];
151 REAL8 tmp144=3.*sigmaKerr->data[2]*tmp143;
152 REAL8 tmp145=4.*sigmaStar->data[2]*tmp143;
154 REAL8 tmp147=0.08333333333333333*eta*tmp14*tmp146;
155 REAL8 tmp148=-8.*sigmaKerr->data[2];
156 REAL8 tmp149=14.*sigmaStar->data[2];
157 REAL8 tmp150=-36.*sigmaKerr->data[2]*tmp101*tmp29*tmp55*tmp75*tmp88*tmp94;
158 REAL8 tmp151=-30.*sigmaStar->data[2]*tmp101*tmp29*tmp55*tmp75*tmp88*tmp94;
159 REAL8 tmp152=3.*sigmaKerr->data[2]*tmp143*x->data[0];
160 REAL8 tmp153=4.*sigmaStar->data[2]*tmp143*x->data[0];
162 REAL8 tmp155=-0.08333333333333333*eta*tmp154*tmp25;
168 REAL8 tmp174=-2.*tmp140*tmp173*x->data[0];
171 REAL8 tmp180=-2.*eta*tmp177*tmp179*x->data[0];
182 REAL8 tmp198=0.013888888888888888*sigmaStar->data[2]*tmp197*tmp25;
183 REAL8 tmp199=27.*eta;
187 REAL8 tmp203=-2.*tmp143*tmp173*x->data[0];
191 REAL8 tmp207=-0.027777777777777776*sigmaStar->data[2]*tmp206*tmp35;
197 REAL8 tmp216=4.*tmp140*tmp215*x->data[0];
200 REAL8 tmp219=-90.*eta*tmp179*x->data[0];
211 REAL8 tmp237=0.006944444444444444*sigmaKerr->data[2]*tmp236*tmp25;
213 REAL8 tmp239=-112.*eta*tmp238;
215 REAL8 tmp241=4.*tmp143*tmp215*x->data[0];
219 REAL8 tmp245=-0.013888888888888888*sigmaKerr->data[2]*tmp244*tmp35;
221 REAL8 tmp259=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp35;
222 REAL8 tmp260=0.08333333333333333*eta*tmp14*tmp154;
223 REAL8 tmp261=0.013888888888888888*sigmaStar->data[2]*tmp206*tmp25;
224 REAL8 tmp262=0.006944444444444444*sigmaKerr->data[2]*tmp244*tmp25;
225 REAL8 tmp263=0.+sigmaStar->data[2]+tmp259+tmp260+tmp261+tmp262;
239 REAL8 tmp257=1./sqrt(tmp256);
240 REAL8 tmp272=1./sqrt(tmp74);
241 REAL8 tmp275=1./sqrt(tmp249);
242 REAL8 tmp282=1./sqrt(tmp252);
243 REAL8 tmp277=2.*tmp29*tmp55*x->data[0];
259 REAL8 tmp300=-4.*tmp19*x->data[0];
260 REAL8 tmp301=2.*coeffs->k2;
261 REAL8 tmp302=3.*coeffs->k3;
262 REAL8 tmp303=4.*coeffs->k4;
270 REAL8 tmp311=coeffs->k1+tmp310;
300 REAL8 tmp413=1./sqrt(tmp412);
305 REAL8 tmp418=1.*tmp75*x->data[0];
321 REAL8 tmp402=1./sqrt(tmp296);
328 REAL8 tmp368=2.*eta*tmp29*tmp43*tmp50*x->data[0];
329 REAL8 tmp369=-5.*coeffs->k5l*tmp25;
347 REAL8 tmp387=2.*tmp55*tmp59*x->data[0];
349 REAL8 tmp1=s1Vec->data[2]*s1Vec->data[2];
350 REAL8 tmp3=s2Vec->data[2]*s2Vec->data[2];
352 REAL8 tmp78=1./sqrt(tmp77);
356 REAL8 tmp115=1./sqrt(tmp114);
358 REAL8 tmp506=-66.*sigmaKerr->data[2]*tmp101*tmp29*tmp55*tmp75*tmp87*tmp94;
359 REAL8 tmp507=-52.*sigmaStar->data[2]*tmp101*tmp29*tmp55*tmp75*tmp87*tmp94;
361 REAL8 tmp509=0.08333333333333333*eta*tmp14*tmp508;
368 REAL8 tmp517=0.013888888888888888*sigmaStar->data[2]*tmp25*tmp516;
375 REAL8 tmp524=0.006944444444444444*sigmaKerr->data[2]*tmp25*tmp523;
378 REAL8 tmp398=pow(tmp74,-2.5);
384 REAL8 tmp488=1.*coeffs->dheffSSv2*eta*tmp33*tmp4;
395 REAL8 tmp499=1./sqrt(tmp498);
396 REAL8 tmp556=6.*sigmaKerr->data[2]*tmp101*tmp11*tmp70*tmp74*tmp8;
397 REAL8 tmp557=8.*sigmaStar->data[2]*tmp101*tmp11*tmp70*tmp74*tmp8;
399 REAL8 tmp559=0.08333333333333333*eta*tmp14*tmp558;
404 REAL8 tmp564=0.006944444444444444*sigmaKerr->data[2]*tmp25*tmp563;
409 REAL8 tmp569=0.013888888888888888*sigmaStar->data[2]*tmp25*tmp568;
416 REAL8 g0=1.*eta*tmp499*(1.5*tmp293*tmp33-tmp246*tmp263*tmp35-4.*coeffs->dheffSSv2*eta*tmp4*tmp5-tmp257*tmp275*tmp282*tmp295*tmp297*tmp321*tmp324*tmp361*tmp396*tmp61-0.25*tmp140*tmp269*tmp275*tmp282*tmp295*tmp297*tmp321*tmp324*tmp361*tmp65-0.25*tmp140*tmp275*tmp282*tmp295*tmp297*tmp321*tmp357*tmp358*tmp361*tmp65+0.25*tmp257*tmp275*tmp282*tmp295*tmp321*tmp324*tmp361*tmp402*tmp407*tmp65-tmp17*tmp266*tmp61*tmp65+0.5*tmp257*tmp275*tmp282*tmp295*tmp297*tmp324*tmp361*tmp65*(0.+tmp17*(-8.*tmp18-4.*tmp19+1.*tmp12*tmp388)+tmp298*tmp61+tmp319*tmp69)+1.*tmp11*tmp17*tmp70+1.*e3z*tmp17*tmp246*tmp70+1.*tmp11*tmp248*tmp255*tmp257*tmp266*tmp70+1.*tmp11*tmp18*tmp255*tmp257*tmp266*tmp272*tmp70+1.*tmp266*tmp69*tmp70-eta*tmp120*tmp25*tmp253*tmp43*tmp431*tmp433*tmp50*tmp73*tmp75-2.*tmp102*tmp14*tmp253*tmp431*tmp433*tmp73*tmp76-tmp124*tmp25*tmp253*tmp431*tmp433*tmp59*tmp75*tmp76-tmp25*tmp253*tmp344*tmp411*tmp431*tmp73*tmp75*tmp76+0.5*tmp25*tmp282*tmp288*tmp431*tmp433*tmp73*tmp75*tmp76-2.*tmp253*tmp35*tmp431*tmp433*tmp73*tmp75*tmp76+0.5*tmp115*tmp78*(tmp100+tmp103+tmp104+tmp105+tmp107+tmp108+tmp109-4.*eta*tmp35*tmp80*tmp83+tmp95)-1.5*tmp257*tmp275*tmp282*tmp297*tmp321*tmp324*tmp361*tmp398*tmp65*x->data[0]-tmp11*tmp17*tmp61*tmp65*x->data[0]-tmp11*tmp248*tmp255*tmp257*tmp266*tmp61*tmp65*x->data[0]+1.*e3z*tmp11*tmp246*tmp248*tmp255*tmp257*tmp70*x->data[0]-0.5*tmp11*tmp140*tmp248*tmp255*tmp266*tmp269*tmp70*x->data[0]+1.*tmp11*tmp69*tmp70*x->data[0]+1.*tmp11*tmp248*tmp257*tmp266*tmp70*x->data[0]*(-0.5*tmp275*tmp279+0.5*tmp248*tmp282*tmp288+1.*tmp253*tmp272*x->data[0])+0.5*tmp257*tmp275*tmp282*tmp295*tmp297*tmp321*tmp324*tmp65*(1.*eta*tmp18*tmp29*tmp351*tmp43*tmp50+1.*tmp18*tmp351*tmp55*tmp59-tmp106*tmp133*tmp29*tmp327*tmp330*tmp55*tmp61*tmp65+1.*eta*tmp106*tmp133*tmp29*tmp327*tmp330*tmp43*tmp50*tmp70+1.*tmp106*tmp133*tmp246*tmp29*tmp326*tmp327*tmp55*tmp70+4.*tmp101*tmp106*tmp29*tmp327*tmp330*tmp55*tmp70+1.*tmp106*tmp133*tmp327*tmp330*tmp55*tmp59*tmp70+4.*tmp106*tmp168*tmp29*tmp330*tmp55*tmp70*tmp74+2.*tmp29*tmp351*tmp55*x->data[0]+1.*tmp18*tmp29*tmp55*(tmp246*tmp326*tmp337*tmp74+tmp330*tmp344*tmp74-eta*tmp18*tmp29*tmp330*tmp43*tmp50*tmp88*tmp94-tmp18*tmp246*tmp29*tmp326*tmp55*tmp88*tmp94-tmp18*tmp330*tmp55*tmp59*tmp88*tmp94-tmp18*tmp29*tmp330*tmp55*tmp88*tmp98*tmp99+2.*tmp330*tmp337*x->data[0]-2.*tmp29*tmp330*tmp55*tmp88*tmp94*x->data[0]))-0.25*tmp257*tmp279*tmp282*tmp295*tmp297*tmp321*tmp324*tmp361*tmp65*(1.0/sqrt(tmp249*tmp249*tmp249))-0.25*tmp257*tmp275*tmp288*tmp295*tmp297*tmp321*tmp324*tmp361*tmp65*(1.0/sqrt(tmp252*tmp252*tmp252))+1.*tmp25*tmp253*tmp433*tmp73*tmp75*tmp76*(0.5*tmp402*tmp407*tmp429+tmp297*(-0.5*tmp11*tmp248*tmp253*tmp323*tmp330*tmp413*tmp416-0.5*tmp11*tmp18*tmp253*tmp272*tmp323*tmp330*tmp413*tmp416+0.5*tmp275*tmp279*tmp427-0.5*tmp11*tmp246*tmp248*tmp253*tmp323*tmp326*tmp413*tmp416*x->data[0]-0.25*tmp11*tmp140*tmp248*tmp253*tmp257*tmp330*tmp413*tmp416*x->data[0]-0.25*tmp11*tmp248*tmp282*tmp288*tmp323*tmp330*tmp413*tmp416*x->data[0]-0.5*tmp11*tmp248*tmp253*tmp323*tmp330*tmp413*x->data[0]*(tmp297*tmp388+0.5*tmp317*tmp402*tmp407-2.*eta*tmp18*tmp29*tmp43*tmp50-2.*tmp18*tmp55*tmp59-4.*tmp29*tmp55*x->data[0])+tmp250*(1.*tmp11*tmp248*tmp253*tmp330*tmp423*tmp425+1.*tmp11*tmp18*tmp253*tmp272*tmp330*tmp423*tmp425+1.*tmp11*tmp140*tmp248*tmp253*tmp257*tmp330*tmp423*x->data[0]+1.*tmp11*tmp246*tmp248*tmp253*tmp326*tmp423*tmp425*x->data[0]+0.5*tmp11*tmp248*tmp282*tmp288*tmp330*tmp423*tmp425*x->data[0]+1.*tmp11*tmp248*tmp253*tmp330*tmp425*x->data[0]*(-2.*tmp102*tmp18-0.5*eta*tmp120*tmp19*tmp25*tmp421*tmp43*tmp50*tmp70*tmp73+1.*tmp75-0.5*tmp124*tmp19*tmp25*tmp421*tmp59*tmp70*tmp76-0.5*tmp19*tmp25*tmp421*tmp61*tmp65*tmp73*tmp76+1.*tmp14*tmp421*tmp70*tmp73*tmp76-tmp19*tmp35*tmp421*tmp70*tmp73*tmp76+0.5*tmp19*tmp25*tmp70*tmp73*tmp76*(tmp19*tmp388-4.*eta*tmp101*tmp29*tmp43*tmp50-12.*tmp18*tmp29*tmp55-4.*tmp101*tmp55*tmp59+2.*tmp317*x->data[0])))+0.25*tmp11*tmp248*tmp253*tmp323*tmp330*tmp416*(4.*tmp101*tmp158*tmp163*tmp94+2.*eta*tmp133*tmp158*tmp43*tmp50*tmp55*tmp94+2.*tmp133*tmp163*tmp29*tmp59*tmp94+tmp133*tmp158*tmp163*tmp98*tmp99)*x->data[0]*(1.0/sqrt(tmp412*tmp412*tmp412))))-0.5*tmp118*(-(eta*tmp120*tmp25*tmp43*tmp50*tmp64*tmp73*tmp75)-2.*tmp102*tmp14*tmp64*tmp73*tmp76-tmp124*tmp25*tmp59*tmp64*tmp75*tmp76+tmp25*tmp61*tmp73*tmp75*tmp76-2.*tmp35*tmp64*tmp73*tmp75*tmp76)*(1.0/sqrt(tmp77*tmp77*tmp77)));
REAL8 tmp265
REAL8 tmp236
REAL8 tmp312
REAL8 tmp292
REAL8 g4
REAL8 tmp149
REAL8 tmp494
REAL8 tmp241
REAL8 tmp167
REAL8 tmp279
REAL8 tmp424
REAL8 tmp488
REAL8 tmp384
REAL8 tmp421
REAL8 tmp557
REAL8 tmp140
REAL8 tmp542
REAL8 tmp422
REAL8 tmp79
REAL8 tmp105
REAL8 tmp131
REAL8 tmp560
REAL8 tmp210
REAL8 tmp382
REAL8 tmp63
REAL8 tmp269
REAL8 tmp97
REAL8 tmp224
REAL8 tmp415
REAL8 tmp73
REAL8 tmp253
REAL8 tmp33
REAL8 tmp173
REAL8 tmp296
REAL8 tmp226
REAL8 tmp387
REAL8 g3
REAL8 tmp569
REAL8 tmp564
REAL8 tmp205
REAL8 tmp47
REAL8 tmp359
REAL8 tmp511
REAL8 tmp176
REAL8 tmp90
REAL8 tmp263
REAL8 tmp566
REAL8 tmp48
REAL8 tmp244
REAL8 tmp109
REAL8 tmp319
REAL8 tmp518
REAL8 tmp229
REAL8 tmp155
REAL8 tmp203
REAL8 tmp107
REAL8 tmp78
REAL8 tmp190
REAL8 tmp403
REAL8 tmp322
REAL8 tmp493
REAL8 tmp170
REAL8 tmp169
REAL8 tmp128
REAL8 tmp172
REAL8 tmp58
REAL8 tmp297
REAL8 tmp3
REAL8 tmp368
REAL8 tmp318
REAL8 tmp287
REAL8 tmp41
REAL8 tmp248
REAL8 tmp301
REAL8 tmp381
REAL8 tmp106
REAL8 tmp199
REAL8 tmp521
REAL8 tmp50
REAL8 tmp114
REAL8 tmp558
REAL8 tmp144
REAL8 tmp396
REAL8 tmp44
REAL8 tmp249
REAL8 tmp277
REAL8 tmp118
REAL8 tmp74
REAL8 tmp181
REAL8 tmp92
REAL8 tmp412
REAL8 tmp522
REAL8 tmp299
REAL8 tmp220
REAL8 tmp180
REAL8 tmp223
REAL8 tmp89
REAL8 tmp404
REAL8 tmp17
REAL8 tmp543
REAL8 tmp261
REAL8 tmp18
REAL8 tmp208
REAL8 tmp320
REAL8 tmp166
REAL8 tmp559
REAL8 tmp562
REAL8 tmp233
REAL8 tmp187
REAL8 tmp213
REAL8 tmp29
REAL8 tmp69
REAL8 tmp283
REAL8 tmp519
REAL8 tmp82
REAL8 tmp19
REAL8 tmp197
REAL8 tmp161
REAL8 tmp11
REAL8 tmp94
REAL8 tmp42
REAL8 tmp27
REAL8 tmp57
REAL8 tmp14
REAL8 tmp36
REAL8 tmp130
REAL8 tmp512
REAL8 tmp43
REAL8 tmp300
REAL8 tmp83
REAL8 tmp76
REAL8 tmp491
REAL8 tmp370
REAL8 g1
REAL8 tmp55
REAL8 tmp192
REAL8 tmp225
REAL8 tmp70
REAL8 tmp251
REAL8 tmp257
REAL8 tmp24
REAL8 tmp369
REAL8 tmp234
REAL8 tmp417
REAL8 tmp427
REAL8 tmp111
REAL8 tmp108
REAL8 tmp5
REAL8 tmp147
REAL8 tmp416
REAL8 tmp21
REAL8 tmp413
REAL8 tmp88
REAL8 tmp13
REAL8 tmp146
REAL8 tmp316
REAL8 tmp503
REAL8 tmp32
REAL8 tmp174
REAL8 tmp4
REAL8 tmp40
REAL8 tmp282
REAL8 tmp314
REAL8 tmp113
REAL8 tmp37
REAL8 tmp216
REAL8 tmp499
REAL8 tmp581
REAL8 tmp56
REAL8 tmp212
REAL8 tmp60
REAL8 tmp103
REAL8 tmp81
REAL8 tmp10
REAL8 tmp506
REAL8 tmp191
REAL8 tmp95
REAL8 tmp425
REAL8 tmp145
REAL8 tmp77
REAL8 tmp372
REAL8 tmp313
REAL8 tmp183
REAL8 tmp329
REAL8 tmp152
REAL8 tmp304
REAL8 tmp28
REAL8 tmp104
REAL8 tmp65
REAL8 tmp515
REAL8 tmp376
REAL8 tmp93
REAL8 tmp134
REAL8 tmp101
REAL8 tmp120
REAL8 tmp317
REAL8 tmp204
REAL8 tmp182
REAL8 tmp496
REAL8 tmp583
REAL8 tmp8
REAL8 tmp517
REAL8 tmp570
REAL8 tmp508
REAL8 tmp507
REAL8 tmp374
REAL8 tmp75
REAL8 tmp242
REAL8 tmp177
REAL8 tmp431
REAL8 tmp163
REAL8 tmp67
REAL8 tmp165
REAL8 tmp214
REAL8 tmp151
REAL8 tmp185
REAL8 tmp385
REAL8 tmp164
REAL8 tmp286
REAL8 tmp360
REAL8 tmp510
REAL8 tmp207
REAL8 tmp206
REAL8 tmp315
REAL8 tmp39
REAL8 tmp262
REAL8 tmp309
REAL8 tmp153
REAL8 tmp137
REAL8 tmp358
REAL8 tmp198
REAL8 tmp26
REAL8 tmp520
REAL8 tmp217
REAL8 tmp337
REAL8 tmp293
REAL8 tmp260
REAL8 tmp7
REAL8 tmp87
REAL8 tmp509
REAL8 tmp135
REAL8 tmp115
REAL8 tmp194
REAL8 tmp256
REAL8 tmp189
REAL8 tmp228
REAL8 tmp85
REAL8 g0
REAL8 tmp278
REAL8 tmp310
REAL8 tmp568
REAL8 tmp62
REAL8 tmp154
REAL8 tmp284
REAL8 tmp302
REAL8 tmp371
REAL8 tmp160
REAL8 tmp49
REAL8 tmp235
REAL8 tmp428
REAL8 tmp490
REAL8 tmp159
REAL8 tmp186
REAL8 tmp157
REAL8 tmp426
REAL8 tmp99
REAL8 tmp361
REAL8 tmp68
REAL8 tmp16
REAL8 tmp343
REAL8 tmp124
REAL8 tmp305
REAL8 tmp196
REAL8 tmp373
REAL8 tmp1
REAL8 tmp398
REAL8 tmp245
REAL8 tmp222
REAL8 tmp380
REAL8 tmp326
REAL8 tmp383
REAL8 tmp96
REAL8 tmp523
REAL8 tmp288
REAL8 tmp266
REAL8 tmp252
REAL8 tmp139
REAL8 tmp61
REAL8 tmp327
REAL8 tmp52
REAL8 tmp516
REAL8 tmp22
REAL8 tmp23
REAL8 tmp406
REAL8 tmp232
REAL8 tmp295
REAL8 tmp168
REAL8 tmp188
REAL8 tmp259
REAL8 tmp498
REAL8 tmp561
REAL8 tmp357
REAL8 tmp64
REAL8 tmp31
REAL8 tmp30
REAL8 tmp433
REAL8 tmp238
REAL8 tmp215
REAL8 tmp377
REAL8 tmp556
REAL8 tmp240
REAL8 tmp378
REAL8 tmp20
REAL8 tmp243
REAL8 tmp141
REAL8 tmp429
REAL8 tmp272
REAL8 tmp414
REAL8 tmp321
REAL8 tmp230
REAL8 tmp102
REAL8 tmp514
REAL8 tmp45
REAL8 tmp201
REAL8 tmp9
REAL8 tmp80
REAL8 tmp350
REAL8 tmp582
REAL8 tmp351
REAL8 tmp524
REAL8 tmp405
REAL8 tmp237
REAL8 tmp218
REAL8 tmp375
REAL8 tmp419
REAL8 tmp525
REAL8 tmp175
REAL8 tmp513
REAL8 tmp239
REAL8 tmp324
REAL8 tmp303
REAL8 tmp86
REAL8 tmp112
REAL8 tmp308
REAL8 tmp158
REAL8 tmp574
REAL8 tmp418
REAL8 tmp489
REAL8 tmp179
REAL8 g2
REAL8 tmp171
REAL8 tmp195
REAL8 tmp162
REAL8 tmp565
REAL8 tmp12
REAL8 tmp178
REAL8 tmp307
REAL8 tmp46
REAL8 tmp563
REAL8 tmp148
REAL8 tmp492
REAL8 tmp34
REAL8 tmp330
REAL8 tmp402
REAL8 tmp54
REAL8 tmp221
REAL8 tmp388
REAL8 tmp98
REAL8 tmp143
REAL8 tmp15
REAL8 tmp142
REAL8 tmp202
REAL8 tmp231
REAL8 tmp255
REAL8 tmp91
REAL8 tmp136
REAL8 tmp349
REAL8 tmp344
REAL8 tmp285
REAL8 tmp193
REAL8 tmp150
REAL8 tmp100
REAL8 tmp298
REAL8 tmp567
REAL8 tmp51
REAL8 tmp367
REAL8 tmp411
REAL8 tmp325
REAL8 tmp306
REAL8 tmp254
REAL8 tmp200
REAL8 tmp227
REAL8 tmp495
REAL8 tmp430
REAL8 tmp184
REAL8 tmp59
REAL8 tmp209
REAL8 tmp246
REAL8 tmp129
REAL8 tmp25
REAL8 tmp138
REAL8 tmp219
REAL8 tmp156
REAL8 tmp276
REAL8 tmp38
REAL8 tmp132
REAL8 g5
REAL8 tmp275
REAL8 tmp211
REAL8 tmp250
REAL8 tmp379
REAL8 tmp311
REAL8 tmp323
REAL8 tmp497
REAL8 tmp133
REAL8 tmp35
REAL8 tmp386
REAL8 tmp420
REAL8 tmp423
REAL8 tmp53
REAL8 tmp407
double REAL8
list x
list p