($Id: cc523e423547ec1e335b86a2ead55e4b2e6b3040 $) using OpenMP (TM) compiled with gcc 12.2.0 < Millepede II-P starting ... Tue Nov 22 11:00:40 2022 ________________________________FILETC/X__________________________________ Command line options: --------------------- 1 mp2str.txt --------------------- Listing of steering file: mp2str.txt ------------------------- 1 * *** Default test steering file *** 2 3 * Additinal files 4 fortranfiles ! following bin files are fortran 5 mp2con.txt ! constraints text file 6 mp2tst.bin ! binary data file 7 Cfiles ! following bin files are Cfiles 8 9 * Selection of variable global parameters 10 *entries 10 ! lower limit on number of entries/parameter 11 *entries 25 5 ! default limit on number of entries/parameter 12 entries 50 ! higher limit on number of entries/parameter 13 14 * Initial parameter values, presigmas 15 *Parameter ! define parameter attributes (start of list) 16 *205 -0.01 0. ! start value 17 *206 0.01 0. ! start value 18 *215 0.01 -1. ! fix parameter at value 19 *216 0. -1. ! fix parameter at value 20 21 * Handling of outliers, tails etc 22 *hugecut 50.0 !cut factor in iteration 0 23 chisqcut 30.0 6.0 ! cut factor in iterations 1 and 2 24 *outlierdownweighting 2 ! number of internal iterations (> 1) 25 *dwfractioncut 0.2 ! 0 < value < 0.5 26 *presigma 0.01 ! default value for presigma 27 *regularisation 1.0 ! regularisation factor 28 *regularisation 1.0 0.01 ! regularisation factor, pre-sigma 29 30 * Solution methods 31 method fullMINRES 3 0.01 ! (approximate) MINRES, full storage 32 method sparseMINRES 3 0.01 ! (approximate) MINRES, sparse storage 33 *mrestol 1.0D-8 ! epsilon for MINRES convergence 34 *bandwidth 0 ! width of MINRES precond. band matrix 35 method diagonalization 3 0.001 ! diagonalization (-> millepede.eve) 36 method decomposition 3 0.001 ! Cholesky decomposition 37 method inversion 3 0.001 ! Gauss matrix inversion 38 * last method is applied 39 40 * Additional output. monitoring 41 printcounts ! print number of entries 42 *monitorresiduals ! poor man's DMR (-> millepede.mon) 43 *printrecord 1 2 ! debug printout for records 44 *printrecord -1 -1 ! debug printout for bad data records 45 46 end ! optional for end-of-data end-statement after 46 text lines ------------------------- lun 11 0 1 binary files opened Processing text files ... File mp2str.txt end-statement after 46 text lines File mp2con.txt end-of-file after 24 text lines ... end of text file processing. Solution method and matrix-storage mode: METSOL = 1: matrix inversion with 3 iterations MATSTO = 1: full symmetric matrix, (n*n+n)/2 elements Chi square cut equiv 3 st.dev applied ... in first iteration with factor 30.0000000 in second iteration with factor 6.00000000 (reduced by sqrt in next iterations) Iterations (solutions) with line search: All with Chi square cut scaling factor <= 1. ______________________________FILETC/X-end________________________________ Number of processors available: 8 Maximum number of OpenMP threads: 1 Number of threads for processing: 1 Number of threads for reading: 1 _________________________________LOOP1____________________________________ Read all binary data files: PEREAD: file 1 read the first time, found 10000 records Record 1 Record 10 Record 100 Record 1000 Record 10000 PEREAD: file 1 with max record size 370 Read cache usage (#blocks, #records, min,max records/block 1 10000 10000 10000 SORT22 (quicksort): maxlevel used/available = 12 /64 LOOP1: 0 is number of pre-sigmas Warning: no pre-sigmas defined Default pre-sigma = 0.00E+00 (if no individual pre-sigma defined) Pre-sigma factor is 1.00000000 No regularization will be done NTGB = 700 = total number of parameters NVGB = 649 = number of variable parameters NREC = 10000 = number of records MREQENF = 50 = required number of entries (eqns in binary files) MREQENA = 10 = required number of entries (from accepted fits) NTGB = 700 = total number of parameters NVGB = 649 = number of variable parameters _______________________________LOOP1-end__________________________________ 0 h 0 min 0.1 sec total elapsed 0 h 0 min 0.1 sec _________________________________LOOP2____________________________________ PRPCON: 4 constraints accepted SORT2I (quicksort): maxlevel used/available = 2 /64 SORT2I (quicksort): maxlevel used/available = 2 /64 PRPCON: constraints split into 4 (disjoint) groups, groups combined into 4 (non overlapping) blocks max group size (cons., par.) 1 41 max block size (cons., par.) 1 41 Record 1 Record 10 Record 100 Record 1000 Record 10000 Read cache usage (#blocks, #records, min,max records/block 1 10000 10000 10000 NTGB = 700 = total number of parameters (all parameters, appearing in binary files) NVGB = 649 = number of variable parameters (appearing in fit matrix/vectors) NAGB = 649 = number of all parameters (including Lagrange multiplier or reduced) NTPGRP = 700 = total number of parameter groups NVPGRP = 649 = number of variable parameter groups NFGB = 645 = number of fit parameters NOFF = 210276 = max number of off-diagonal elements NCGB = 4 = number of constraints NAGBN = 18 = max number of global parameters in an event NALCN = 28 = max number of local parameters in an event NAEQN = 38 = max number of equations in an event NCACHE = 25000000 = number of words for caching cache splitting 57.1 % 3.0 % 39.9 % Solution method and matrix-storage mode: METSOL = 1: matrix inversion with 3 iterations MATSTO = 1: full symmetric matrix, (n*n+n)/2 elements Convergence assumed, if expected dF < 0.1000E-02 Constraints handled by elimination Rank of product matrix of constraints is 4 for 4 constraint equations QL decomposition of constraints matrix largest |eigenvalue| of L: -2.2360679774997898 smallest |eigenvalue| of L: -2.2360679774997898 Size of global matrix: 1 MB _______________________________LOOP2-end__________________________________ 0 h 0 min 0.2 sec total elapsed 0 h 0 min 0.1 sec _______________________________Iteration__________________________________ Solution algorithm: ======================================== solution method: matrix inversion convergence limit at Delta F= 0.10E-02 maximum number of iterations= 3 using pre-sigmas: no regularization: no Chi square cut equiv 3 st.dev applied ... in first iteration with factor 0.30E+02 ... in second iteration with factor 0.60E+01 (reduced by sqrt in next iterations) Checking feasibility of parameters: parameters are feasible (i.e. satisfy constraints) Reading files and accumulating vectors/matrices ... Read cache usage (#blocks, #records, min,max records/block 1 10000 10000 10000 Write cache usage (#flush,#overrun,,peak(levels)) 1 0, 22.3% 21.0% 22.3% 21.0% Data rejected in initial loop: 0 (rank deficit/NaN) 0 (Ndf=0) 0 (huge) 1 (large) === local fits have bordered band matrix structure === NBNDR = 10000 = number of records (upper/left border) NBDRX = 0 = max border size NBNDX = 4 = max band width it fc fcn_value dfcn_exp slpr costh iit st ls step cutf rejects hhmmss FMS -- -- ----------- -------- ---- ----- --- -- -- ----- ---- ------- ------ --- 0 1 0.19854E+07 30.0 1 00000 FMS 1 2 0.11643E+06 0.19E+07 0 0 6.0 14 00000 F 00000 S 2 3 0.11188E+06 0.13E+02 0 0 2.4 117 00000 F 00000 S 3 4 0.10807E+06 0.16E+03 0 0 1.6 182 00000 F 00000 S 4 5 0.10450E+06 0.17E+03 0.428 0.141 0 0 1 1.000 1.0 276 00000 F Data rejected in last loop: 0 (rank deficit/NaN) 0 (Ndf=0) 0 (huge) 276 (large) Sum(Chi^2)/Sum(Ndf) = 104495.47920560479 / ( 99992 - 645 ) = 1.0518231975359578 WarningWarningWarningWarningWarningWarningWarningWarningWar arningWarningWarningWarningWarningWarningWarningWarningWarn rningWarningWarningWarningWarningWarningWarningWarningWarni ningWarningWarningWarningWarningWarningWarningWarningWarnin ingWarningWarningWarningWarningWarningWarningWarningWarning ngWarningWarningWarningWarningWarningWarningWarningWarningW gWarningWarningWarningWarningWarningWarningWarningWarningWa Fraction of rejects = 2.76 % (should be far below 1 %) => please provide correct mille data WarningWarningWarningWarningWarningWarningWarningWarningWar arningWarningWarningWarningWarningWarningWarningWarningWarn rningWarningWarningWarningWarningWarningWarningWarningWarni ningWarningWarningWarningWarningWarningWarningWarningWarnin ingWarningWarningWarningWarningWarningWarningWarningWarning ngWarningWarningWarningWarningWarningWarningWarningWarningW gWarningWarningWarningWarningWarningWarningWarningWarningWa _____________________________Iteration-end________________________________ 0 h 0 min 0.9 sec total elapsed 0 h 0 min 0.8 sec Result of fit for global parameters =================================== label parameter presigma differ error ----------- ---------------------------------------------------- 50 0.15275E-01 0.0000 0.15275E-01 0.33879E-03 51 -0.43139E-03 0.0000 -0.43139E-03 0.29960E-03 52 0.39162E-02 0.0000 0.39162E-02 0.28379E-03 53 0.10414E-01 0.0000 0.10414E-01 0.28508E-03 54 0.19977E-03 0.0000 0.19977E-03 0.23375E-03 55 -0.71871E-02 0.0000 -0.71871E-02 0.26658E-03 56 0.15601E-02 0.0000 0.15601E-02 0.27588E-03 57 0.15350E-01 0.0000 0.15350E-01 0.27976E-03 58 0.10681E-01 0.0000 0.10681E-01 0.29363E-03 59 0.10742E-01 0.0000 0.10742E-01 0.35736E-03 ... (further printout suppressed, but see log file) < Millepede II-P ending ... Tue Nov 22 11:00:41 2022 Peak dynamic memory allocation: 0.102007 GB