APFEL error at the end of the fit
Created by: wilsonmr
Hello,
Having an issue running fits on the cluster at Edinburgh, the fits appear to be finishing however they are not outputting the LHAPDF grids, so we don't get the results.
@nhartland suggests it is a problem concerning APFEL.
I have pasted below the final bit of output from the fit, I have the full outputs available for a few different configurations of fits:
**** Producing T0 Predictions with Set 170206-003
- Final Positivity Test
- Passed all points for POSF2U
- Passed all points for POSF2DW
- Passed all points for POSF2S
- Passed all points for POSFLL
- Passed all points for POSDYU
- Passed all points for POSDYD
- Passed all points for POSDYS
- Writing fitinfo file...
- Computing arclengths...
- Writing sumrules file...
- Writing preproc file...
- Writing params file...
- Writing out LHAPDF grid: 180321-edi-002
- Solving DGLAP for LHAPDF grid...
Thanks for using LHAPDF 6.2.1. Please make sure to cite the paper:
Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
Checking APFEL v3.0.2 ...
WARNING: FONLL-C is a VFN scheme
... setting VFNS PDF evolution
WARNING: FONLL-C is a NNLO scheme
... setting NNLO perturbative order
Welcome to
_/_/_/ _/_/_/_/ _/_/_/_/ _/_/_/_/ _/
_/ _/ _/ _/ _/ _/ _/
_/_/_/_/ _/_/_/_/ _/_/_/ _/_/_/ _/
_/ _/ _/ _/ _/ _/
_/ _/ _/ _/ _/_/_/_/ _/_/_/_/
_____v3.0.2 A PDF Evolution Library, arXiv:1310.1394
Authors: V. Bertone, S. Carrazza, J. Rojo
Report of the evolution parameters:
QCD evolution
Space-like evolution (PDFs)
Unpolarized evolution
Evolution scheme: VFNS at N2LO
Solution of the DGLAP equation: 'exactalpha' with maximum 6 active flavours
Solution of the coupling equations: 'exact' with maximum 6 active flavours
Coupling reference value:
- AlphaQCD( 1.4142 GeV) = 0.350000
Pole heavy quark masses:
- Mc = 1.4142 GeV
- Mb = 4.5000 GeV
- Mt = 175.0000 GeV
The matching thresholds coincide with the physical masses
muR / muF = 1.0000
Allowed evolution range [ 1.0000 : 10000.0000 ] GeV
The internal subgrids will be locked
Fast evolution enabled
Initialization of the evolution completed in 4.256 s
Report of the electroweak parameters:
Mass of the Z = 91.188 GeV
Mass of the W = 80.385 GeV
Mass of the proton = 0.9383 GeV
sin^2(thetaW) = 0.2313
GFermi = 1.16638E-05
| 0.9743 0.2254 0.0036 |
CKM = | 0.2252 0.9734 0.0414 |
| 0.0089 0.0405 0.9991 |
Z propagator correction = 0.00000
Report of the DIS parameters:
Computation in the FONLL-C mass scheme
Electromagnetic (EM) process
Scattering electron - proton
muR / Q = 1.0000
muF / Q = 1.0000
Target Mass corrections disabled
FONLL damping factor for charm enabled with suppression power = 2
FONLL damping factor for bottom enabled with suppression power = 2
FONLL damping factor for top enabled with suppression power = 2
Intrinsic charm disabled
Initialization of the DIS module completed in 45.390 s
Check ... succeded
WARNING: if there are external grids they cannot be locked
... unlocking subgrids
Welcome to
_/_/_/ _/_/_/_/ _/_/_/_/ _/_/_/_/ _/
_/ _/ _/ _/ _/ _/ _/
_/_/_/_/ _/_/_/_/ _/_/_/ _/_/_/ _/
_/ _/ _/ _/ _/ _/
_/ _/ _/ _/ _/_/_/_/ _/_/_/_/
_____v3.0.2 A PDF Evolution Library, arXiv:1310.1394
Authors: V. Bertone, S. Carrazza, J. Rojo
Report of the evolution parameters:
QCD evolution
Space-like evolution (PDFs)
Unpolarized evolution
Evolution scheme: VFNS at N2LO
Solution of the DGLAP equation: 'truncated' with maximum 5 active flavours
- value of the truncation parameter epsilon = 1.000E-02
Solution of the coupling equations: 'expanded' with maximum 5 active flavours
Coupling reference value:
- AlphaQCD( 91.2000 GeV) = 0.118000
Pole heavy quark masses:
- Mc = 1.5100 GeV
- Mb = 4.9200 GeV
- Mt = 172.5000 GeV
The matching thresholds coincide with the physical masses
muR / muF = 1.0000
Allowed evolution range [ 1.6500 : 100000.0000 ] GeV
Initialization of the evolution completed in 917.283 s
In odeintns.f:
stepsize underflow in rkqsns