($Id: cc523e423547ec1e335b86a2ead55e4b2e6b3040 $) using OpenMP (TM) compiled with gcc 12.2.0 < Millepede II-P starting ... Tue Nov 22 11:07:54 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 mp2param.txt ! start values (from truth) 7 mp2tst.bin ! binary data file 8 Cfiles ! following bin files are Cfiles 9 10 * Selection of variable global parameters 11 *entries 10 ! lower limit on number of entries/parameter 12 *entries 25 5 ! default limit on number of entries/parameter 13 entries 50 ! higher limit on number of entries/parameter 14 15 * Initial parameter values, presigmas 16 *Parameter ! define parameter attributes (start of list) 17 *205 -0.01 0. ! start value 18 *206 0.01 0. ! start value 19 *215 0.01 -1. ! fix parameter at value 20 *216 0. -1. ! fix parameter at value 21 22 * Handling of outliers, tails etc 23 *hugecut 50.0 !cut factor in iteration 0 24 chisqcut 30.0 6.0 ! cut factor in iterations 1 and 2 25 *outlierdownweighting 2 ! number of internal iterations (> 1) 26 *dwfractioncut 0.2 ! 0 < value < 0.5 27 *presigma 0.01 ! default value for presigma 28 *regularisation 1.0 ! regularisation factor 29 *regularisation 1.0 0.01 ! regularisation factor, pre-sigma 30 31 * Solution methods 32 method fullMINRES 3 0.01 ! (approximate) MINRES, full storage 33 method sparseMINRES 3 0.01 ! (approximate) MINRES, sparse storage 34 *mrestol 1.0D-8 ! epsilon for MINRES convergence 35 *bandwidth 0 ! width of MINRES precond. band matrix 36 method diagonalization 3 0.001 ! diagonalization (-> millepede.eve) 37 method decomposition 3 0.001 ! Cholesky decomposition 38 method inversion 3 0.001 ! Gauss matrix inversion 39 * last method is applied 40 41 * Additional output. monitoring 42 printcounts ! print number of entries 43 *monitorresiduals ! poor man's DMR (-> millepede.mon) 44 *printrecord 1 2 ! debug printout for records 45 *printrecord -1 -1 ! debug printout for bad data records 46 47 end ! optional for end-of-data end-statement after 47 text lines ------------------------- lun 11 0 1 binary files opened Processing text files ... File mp2str.txt end-statement after 47 text lines File mp2con.txt end-of-file after 24 text lines File mp2param.txt end-of-file after 701 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 = 14 /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% === 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.99726E+05 30.0 0 00000 FMS 1 2 0.99112E+05 0.61E+03 0 0 6.0 0 00000 F 00000 S 2 3 0.99051E+05 0.86E-15 0 0 1.0 24 00000 F Data rejected in last loop: 0 (rank deficit/NaN) 0 (Ndf=0) 0 (huge) 24 (large) Sum(Chi^2)/Sum(Ndf) = 99050.817948768963 / ( 99992 - 645 ) = 0.99701871167492695 _____________________________Iteration-end________________________________ 0 h 0 min 0.7 sec total elapsed 0 h 0 min 0.5 sec Result of fit for global parameters =================================== label parameter presigma differ error ----------- ---------------------------------------------------- 50 0.13336E-01 0.0000 -0.53782E-04 0.33879E-03 51 -0.42513E-03 0.0000 -0.55126E-04 0.29960E-03 52 0.40542E-02 0.0000 0.24205E-04 0.28379E-03 53 0.10450E-01 0.0000 -0.15978E-03 0.28508E-03 54 0.97931E-04 0.0000 0.97931E-04 0.23374E-03 55 -0.74712E-02 0.0000 -0.35117E-03 0.26657E-03 56 0.11643E-02 0.0000 -0.25747E-04 0.27588E-03 57 0.14816E-01 0.0000 -0.32412E-03 0.27975E-03 58 0.10151E-01 0.0000 -0.31890E-03 0.29362E-03 59 0.10525E-01 0.0000 0.41482E-03 0.35736E-03 ... (further printout suppressed, but see log file) < Millepede II-P ending ... Tue Nov 22 11:07:55 2022 Peak dynamic memory allocation: 0.102013 GB