|
| integer(mpi) | ictest =0 |
| | test mode '-t'
|
| |
| integer(mpi) | metsol =0 |
| | solution method (1: inversion, 2: diagonalization, 3: MINRES-QLP)
|
| |
| integer(mpi) | matsto =2 |
| | (global) matrix storage mode (1: full, 2: sparse)
|
| |
| integer(mpi) | mprint =1 |
| | print flag (0: minimal, 1: normal, >1: more)
|
| |
| integer(mpi) | mdebug =0 |
| | debug flag (number of records to print)
|
| |
| integer(mpi) | mdebg2 =10 |
| | number of measurements for record debug printout
|
| |
| integer(mpi) | mreqenf =25 |
| | required number of entries (for variable global parameter from binary Files)
|
| |
| integer(mpi) | mreqena =10 |
| | required number of entries (for variable global parameter from Accepted local fits)
|
| |
| integer(mpi) | mitera =1 |
| | number of iterations
|
| |
| integer(mpi) | nloopn =0 |
| | number of data reading, fitting loops
|
| |
| integer(mpi) | mbandw =0 |
| | band width of preconditioner matrix
|
| |
| integer(mpi) | lprecm =0 |
| | additional flag for preconditioner (band) matrix (>0: preserve rank by skyline matrix)
|
| |
| integer(mpi) | lunkno =0 |
| | flag for unkown keywords
|
| |
| integer(mpi) | lhuber =0 |
| | Huber down-weighting flag.
|
| |
| real(mps) | chicut =0.0 |
| | cut in terms of 3-sigma cut, first iteration
|
| |
| real(mps) | chirem =0.0 |
| | cut in terms of 3-sigma cut, other iterations, approaching 1.
|
| |
| real(mps) | chhuge =50.0 |
| | cut in terms of 3-sigma for unreasonable data, all iterations
|
| |
| integer(mpi) | nrecpr =0 |
| | record number with printout
|
| |
| integer(mpi) | nrecp2 =0 |
| | record number with printout
|
| |
| integer(mpi) | nrec1 =0 |
| | record number with largest residual
|
| |
| integer(mpi) | nrec2 =0 |
| | record number with largest chi^2/Ndf
|
| |
| real(mps) | value1 =0.0 |
| | largest residual
|
| |
| real(mps) | value2 =0.0 |
| | largest chi^2/Ndf
|
| |
| real(mps) | dwcut =0.0 |
| | down-weight fraction cut
|
| |
| integer(mpi) | isubit =0 |
| | subito flag '-s'
|
| |
| real(mps) | wolfc1 =0.0 |
| | C_1 of strong Wolfe condition.
|
| |
| real(mps) | wolfc2 =0.0 |
| | C_2 of strong Wolfe condition.
|
| |
| real(mpd) | mrestl =1.0E-06 |
| | tolerance criterion for MINRES-QLP
|
| |
| real(mpd) | mrtcnd =1.0E+07 |
| | transition (QR -> QLP) (matrix) condition for MINRES-QLP
|
| |
| integer(mpi) | mrmode =0 |
| | MINRES-QLP mode (0: QR+QLP, 1: only QR, 2: only QLP factorization)
|
| |
| integer(mpi) | nofeas =0 |
| | flag for skipping making parameters feasible
|
| |
| integer(mpi) | nhistp =0 |
| | flag for histogram printout
|
| |
| real(mps) | delfun =0.0 |
| | expected function change
|
| |
| real(mps) | actfun =0.0 |
| | actual function change
|
| |
| real(mps) | angras =0.0 |
| | angle between gradient and search direction
|
| |
| integer(mpi) | iterat =0 |
| | iterations in solution
|
| |
| integer(mpi) | nregul =0 |
| | regularization flag
|
| |
| real(mps) | regula =1.0 |
| | regularization parameter, add regula * norm(global par.) to objective function
|
| |
| real(mps) | regpre =0.0 |
| | default presigma
|
| |
| integer(mpi) | matrit =0 |
| | matrix calculation up to iteration MATRIT
|
| |
| integer(mpi) | icalcm =0 |
| | calculation mode (for XLOOPN) , >0: calculate matrix
|
| |
| integer(mpi), dimension(2) | nbndr =0 |
| | number of records with bordered band matrix for local fit (upper/left, lower/right)
|
| |
| integer(mpi) | nbdrx =0 |
| | max border size for local fit
|
| |
| integer(mpi) | nbndx =0 |
| | max band width for local fit
|
| |
| integer(mpi) | nrecer =0 |
| | record with error (rank deficit or Not-a-Number) for printout
|
| |
| integer(mpi) | nrec3 = huge(nrec3) |
| | (1.) record number with error
|
| |
| integer(mpi) | mreqpe =1 |
| | min number of pair entries
|
| |
| integer(mpi) | mhispe =0 |
| | upper bound for pair entry histogrammimg
|
| |
| integer(mpi) | msngpe =-1 |
| | upper bound for pair entry single precision storage
|
| |
| integer(mpi) | mcmprs =0 |
| | compression flag for sparsity (column indices)
|
| |
| integer(mpi) | mextnd =0 |
| | flag for extended storage (both 'halves' of sym. mat. for improved access patterns)
|
| |
| integer(mpi) | mthrd =1 |
| | number of (OpenMP) threads
|
| |
| integer(mpi) | mxrec =0 |
| | max number of records
|
| |
| integer(mpi) | matmon =0 |
| | record interval for monitoring of (sparse) matrix construction
|
| |
| integer(mpi) | lfitnp =huge(lfitnp) |
| | local fit: number of iteration to calculate pulls
|
| |
| integer(mpi) | lfitbb =1 |
| | local fit: check for bordered band matrix (if >0)
|
| |
| integer(mpi) | mnrsel =0 |
| | number of MINRES error labels in LBMNRS (calc err, corr with SOLGLO)
|
| |
| integer(mpi) | ncache =-1 |
| | buffer size for caching (default 100MB per thread)
|
| |
| real(mps), dimension(3) | fcache = (/ 0.8, 0., 0. /) |
| | read cache, average fill level; write cache; dynamic size
|
| |
| integer(mpi) | mthrdr =1 |
| | number of threads for reading binary files
|
| |
| integer(mpi) | mnrsit =0 |
| | total number of MINRES internal iterations
|
| |
| integer(mpi) | iforce =0 |
| | switch to SUBITO for (global) rank defects if zero
|
| |
| integer(mpi) | igcorr =0 |
| | flag for output of global correlations for inversion, =0: none
|
| |
| integer(mpi) | memdbg =0 |
| | debug flag for memory management
|
| |
| real(mps) | prange =0.0 |
| | range (-PRANGE..PRANGE) for histograms of pulls, norm. residuals
|
| |
| integer(mpi) | lsearch =2 |
| | iterations (solutions) with line search: >2: all, =2: all with (next) Chi2 cut scaling factor =1., =1: last, <1: none
|
| |
| integer(mpi) | ipcntr =0 |
| | flag for output of global parameter counts (entries), =0: none, =1: local fits, >1: binary files
|
| |
| integer(mpi) | iwcons =0 |
| | flag for weighting of constraints (>0: weighting with globalParCounts, else: none)
|
| |
| integer(mpi) | icelim =1 |
| | flag for using elimination (instead of multipliers) for constraints
|
| |
| integer(mpi) | icheck =0 |
| | flag for checking input only (no solution determined)
|
| |
| integer(mpi) | iteren =0 |
| | entries cut is iterated for parameters with less entries (if > mreqenf)
|
| |
| integer(mpi) | iskpec =0 |
| | flag for skipping empty constraints (no variable parameters)
|
| |
| integer(mpi) | imonit =0 |
| | flag for monitoring residuals per local fit cycle (=0: none, <0: all, bit 0: first, bit 1: last)
|
| |
| integer(mpi) | measbins =100 |
| | number of bins per measurement for monitoring
|
| |
| integer(mpi) | imonmd =0 |
| | monitoring mode: 0:residuals (normalized to average error), 1:pulls
|
| |
| integer(mpi) | iscerr =0 |
| | flag for scaling of errors
|
| |
| real(mpd), dimension(2) | dscerr = (/ 1.0, 1.0 /) |
| | scaling factors for errors of 'global' and 'local' measurement
|
| |
| integer(mpi) | keepopen =1 |
| | flag for keeping binary files open
|
| |
| integer(mpi) | lunmon |
| | unit for monitoring output file
|
| |
| integer(mpi) | lunlog |
| | unit for logfile
|
| |
| integer(mpi) | lvllog |
| | log level
|
| |
| integer(mpi) | ntgb |
| | total number of global parameters
|
| |
| integer(mpi) | nvgb |
| | number of variable global parameters
|
| |
| integer(mpi) | nagb |
| | number of all parameters (global par. + Lagrange mult.)
|
| |
| integer(mpi) | nfgb |
| | number of fit parameters
|
| |
| integer(mpi) | ncgb |
| | number of constraints
|
| |
| integer(mpi) | ncgbe |
| | number of empty constraints (no variable parameters)
|
| |
| integer(mpi) | ncblck |
| | number of (disjoint) constraint blocks
|
| |
| integer(mpi) | mszcon |
| | (integrated block) matrix size for constraint matrix
|
| |
| integer(mpi) | mszprd |
| | (integrated block) matrix size for (constraint) product matrix
|
| |
| integer(mpi), dimension(2) | nprecond |
| | number of constraints, matrix size for preconditioner
|
| |
| integer(mpi) | nagbn |
| | max number of global paramters per record
|
| |
| integer(mpi) | nalcn |
| | max number of local paramters per record
|
| |
| integer(mpi) | naeqn |
| | max number of equations (measurements) per record
|
| |
| integer(mpi) | nrec |
| | number of records read
|
| |
| integer(mpi) | nrecd |
| | number of records read containing doubles
|
| |
| real(mps) | dflim |
| | convergence limit
|
| |
| integer(mpi), dimension(0:3) | nrejec |
| | rejected events
|
| |
| real(mps), dimension(0:8) | times |
| | cpu time counters
|
| |
| real(mps) | stepl |
| | step length (line search)
|
| |
| character(len=74) | textl |
| | name of current MP 'module' (step)
|
| |
| logical | newite |
| | flag for new iteration
|
| |
| integer(mpi) | ndfsum |
| | sum(ndf)
|
| |
| integer(mpi) | iitera |
| | MINRES iterations.
|
| |
| integer(mpi) | istopa |
| | MINRES istop (convergence)
|
| |
| integer(mpi) | lsinfo |
| | line search: returned information
|
| |
| real | rstart |
| | cpu start time for solution iterations
|
| |
| real(mps) | deltim |
| | cpu time difference
|
| |
| integer(mpi) | npresg |
| | number of pre-sigmas
|
| |
| integer(mpi) | nrecal |
| | number of records
|
| |
| integer(mpi) | ndefec =0 |
| | rank deficit for global matrix (from inversion)
|
| |
| integer(mpi) | nmiss1 =0 |
| | rank deficit for constraints
|
| |
| integer(mpi) | nalow =0 |
| | (sum of) global parameters with too few accepted entries
|
| |
| integer(mpi) | lcalcm |
| | last calclation mode
|
| |
| integer(mpi) | nspc =1 |
| | number of precision for sparse global matrix (1=D, 2=D+F)
|
| |
| integer(mpi) | nencdb |
| | encoding info (number bits for column counter)
|
| |
| integer(mpi) | nummeas |
| | number of measurement groups for monitoring
|
| |
| real(mpd), parameter | measbinsize =0.1 |
| | bins size for monitoring
|
| |
| integer(mpi), dimension(100) | lbmnrs |
| | MINRES error labels.
|
| |
| real(mpd) | fvalue |
| | function value (chi2 sum) solution
|
| |
| real(mpd) | flines |
| | function value line search
|
| |
| real(mpd) | sumndf |
| | weighted sum(ndf)
|
| |
| integer(mpi) | numreadbuffer |
| | number of buffers (records) in (read) block
|
| |
| integer(mpi) | numblocks |
| | number of (read) blocks
|
| |
| integer(mpi) | sumrecords |
| | sum of records
|
| |
| integer(mpi) | skippedrecords |
| | number of skipped records (buffer too small)
|
| |
| integer(mpi) | minrecordsinblock |
| | min. records in block
|
| |
| integer(mpi) | maxrecordsinblock |
| | max. records in block
|
| |
| integer(mpi), parameter | nexp20 =1048576 |
| |
| real(mpd) | accuratedsum =0.0_mpd |
| | fractional part of sum
|
| |
| integer(mpi) | accuratensum =0 |
| | sum mod 2**20
|
| |
| integer(mpi) | accuratenexp =0 |
| | sum / 2**20
|
| |
| integer(mpi) | lenglobalvec |
| | length of global vector 'b' (A*x=b)
|
| |
| real(mpd), dimension(:), allocatable | globalparameter |
| | global parameters (start values + sum(x_i))
|
| |
| real(mpd), dimension(:), allocatable | globalparcopy |
| | copy of global parameters
|
| |
| real(mpd), dimension(:), allocatable | globalcorrections |
| | correction x_i (from A*x_i=b_i in iteration i)
|
| |
| real(mpd), dimension(:), allocatable | globalparstart |
| | start value for global parameters
|
| |
| real(mpd), dimension(:), allocatable | globalparpresigma |
| | pre-sigma for global parameters
|
| |
| real(mpd), dimension(:), allocatable | globalparpreweight |
| | weight from pre-sigma
|
| |
| integer(mpi), dimension(:), allocatable | globalparcounts |
| | global parameters counts (from binary files)
|
| |
| integer(mpi), dimension(:), allocatable | globalparcons |
| | global parameters (number of) constraints
|
| |
| real(mpd), dimension(:), allocatable | globalmatd |
| | global matrix 'A' (double, full or sparse)
|
| |
| real(mps), dimension(:), allocatable | globalmatf |
| | global matrix 'A' (float part for compressed sparse)
|
| |
| real(mpd), dimension(:), allocatable | globalvector |
| | global vector 'x' (in A*x=b)
|
| |
| integer(mpi), dimension(:), allocatable | globalcounter |
| | global counter (entries in 'x')
|
| |
| real(mpd), dimension(:), allocatable | vecxav |
| | vector x for AVPROD (A*x=b)
|
| |
| real(mpd), dimension(:), allocatable | vecbav |
| | vector b for AVPROD (A*x=b)
|
| |
| real(mpd), dimension(:), allocatable | matprecond |
| | preconditioner (band) matrix
|
| |
| integer(mpi), dimension(:), allocatable | indprecond |
| | preconditioner pointer array
|
| |
| real(mpd), dimension(:), allocatable | workspaced |
| | (general) workspace (D)
|
| |
| real(mpd), dimension(:), allocatable | workspacediag |
| | diagonal of global matrix (for global corr.)
|
| |
| real(mpd), dimension(:), allocatable | workspacelinesearch |
| | workspace line search
|
| |
| real(mpd), dimension(:), allocatable | workspacediagonalization |
| | workspace diag.
|
| |
| real(mpd), dimension(:), allocatable | workspaceeigenvalues |
| | workspace eigen values
|
| |
| real(mpd), dimension(:), allocatable | workspaceeigenvectors |
| | workspace eigen vectors
|
| |
| integer(mpi), dimension(:), allocatable | workspacei |
| | (general) workspace (I)
|
| |
| real(mpd), dimension(:), allocatable | matconsproduct |
| | product matrix of constraints
|
| |
| real(mpd), dimension(:), allocatable | vecconsresiduals |
| | residuals of constraints
|
| |
| real(mpd), dimension(:), allocatable | vecconssolution |
| | solution for constraint elimination
|
| |
| integer(mpi), dimension(:), allocatable | vecconsstart |
| | start of constraint in listConstraints (unsorted input)
|
| |
| integer(mpi), dimension(:,:), allocatable | matconssort |
| | keys and index for sorting
|
| |
| integer(mpi), dimension(:,:), allocatable | matconsblocks |
| | start of constraint blocks, parameter range
|
| |
| integer(mpi), dimension(:), allocatable | measindex |
| | mapping of 1. global label to measurement index
|
| |
| integer(mpi), dimension(:), allocatable | meashists |
| | measurement histograms (100 bins per thread)
|
| |
| real(mpd), dimension(:), allocatable | measres |
| | average measurement error
|
| |
| integer(mpi), dimension(:,:), allocatable | globalparlabelindex |
| | global parameters label, total -> var. index
|
| |
| integer(mpi), dimension(:), allocatable | globalparhashtable |
| | global parameters hash table
|
| |
| integer(mpi), dimension(:), allocatable | globalparvartototal |
| | global parameters variable -> total index
|
| |
| integer(mpi), dimension(-7:0) | globalparheader = 0 |
| | global parameters (mapping) header
|
| |
| integer(mpi), dimension(:), allocatable | sparsematrixcompression |
| | compression info (per row)
|
| |
| integer(mpi), dimension(:), allocatable | sparsematrixcolumns |
| | (compressed) list of columns for sparse matrix
|
| |
| integer(mpl), dimension(:,:), allocatable | sparsematrixoffsets |
| | row offsets for column list, sparse matrix elements
|
| |
| integer(mpi), dimension(:,:), allocatable | readbufferinfo |
| | buffer management (per thread)
|
| |
| integer(mpi), dimension(:), allocatable | readbufferpointer |
| | pointer to used buffers
|
| |
| integer(mpi), dimension(:), allocatable | readbufferdatai |
| | integer data
|
| |
| real(mpr4), dimension(:), allocatable | readbufferdataf |
| | float data
|
| |
| real(mpr8), dimension(:), allocatable | readbufferdatad |
| | double data
|
| |
| integer(mpi), dimension(:), allocatable | globalindexusage |
| | indices of global par in record
|
| |
| integer(mpi), dimension(:), allocatable | backindexusage |
| | list of global par in record
|
| |
| integer(mpi), dimension(:), allocatable | appearancecounter |
| | appearance statistics for global par (first/last file,record)
|
| |
| integer(mpi), dimension(:), allocatable | paircounter |
| | number of paired parameters (in equations)
|
| |
| real(mpd), dimension(:), allocatable | blvec |
| | local fit vector 'b' (in A*x=b), replaced by 'x'
|
| |
| real(mpd), dimension(:), allocatable | clmat |
| | local fit matrix 'A' (in A*x=b)
|
| |
| integer(mpi), dimension(:), allocatable | ibandh |
| | local fit 'band width histogram' (band size autodetection)
|
| |
| real(mpd), dimension(:), allocatable | vbnd |
| | local fit band part of 'A'
|
| |
| real(mpd), dimension(:), allocatable | vbdr |
| | local fit border part of 'A'
|
| |
| real(mpd), dimension(:), allocatable | aux |
| | local fit 'solutions for border rows'
|
| |
| real(mpd), dimension(:), allocatable | vbk |
| | local fit 'matrix for border solution'
|
| |
| real(mpd), dimension(:), allocatable | vzru |
| | local fit 'border solution'
|
| |
| real(mpd), dimension(:), allocatable | scdiag |
| | local fit workspace (D)
|
| |
| integer(mpi), dimension(:), allocatable | scflag |
| | local fit workspace (I)
|
| |
| real(mpd), dimension(:), allocatable | localcorrections |
| | local fit corrections (to residuals)
|
| |
| real(mpd), dimension(:), allocatable | localglobalmatrix |
| | matrix correlating local and global par
|
| |
| integer(mpi), dimension(:,:), allocatable | writebufferinfo |
| | write buffer management (per thread)
|
| |
| real(mps), dimension(:,:), allocatable | writebufferdata |
| | write buffer data (largest residual, Chi2/ndf, per thread)
|
| |
| integer(mpi), dimension(:), allocatable | writebufferindices |
| | write buffer for indices
|
| |
| real(mpd), dimension(:), allocatable | writebufferupdates |
| | write buffer for update matrices
|
| |
| integer(mpi), dimension(-6:6) | writebufferheader = 0 |
| | write buffer header (-6..-1: updates, 1..6: indices)
|
| |
| integer(mpi) | lenparameters =0 |
| | list items from steering file
|
| |
| type(listitem), dimension(:), allocatable | listparameters |
| | list of parameters from steering file
|
| |
| integer(mpi) | lenpresigmas =0 |
| | length of list of pre-sigmas from steering file
|
| |
| type(listitem), dimension(:), allocatable | listpresigmas |
| | list of pre-sgmas from steering file
|
| |
| integer(mpi) | lenconstraints =0 |
| | length of list of constraints from steering file
|
| |
| type(listitem), dimension(:), allocatable | listconstraints |
| | list of constraints from steering file
|
| |
| integer(mpi) | nummeasurements =0 |
| | number of (external) measurements from steering file
|
| |
| integer(mpi) | lenmeasurements =0 |
| | length of list of (external) measurements from steering file
|
| |
| type(listitem), dimension(:), allocatable | listmeasurements |
| | list of (external) measurements from steering file
|
| |
| integer(mpi), dimension(:), allocatable | mfd |
| | file mode: cbinary =1, text =2, fbinary=3
|
| |
| integer(mpi), dimension(:), allocatable | lfd |
| | length of file name
|
| |
| integer(mpi), dimension(:), allocatable | nfd |
| | index (line) in (steering) file
|
| |
| integer(mpi), dimension(:,:), allocatable | kfd |
| | (1,.)= number of records in file, (2,..)= file order
|
| |
| integer(mpi), dimension(:), allocatable | ifd |
| | file: integrated record numbers (=offset)
|
| |
| integer(mpi), dimension(:), allocatable | jfd |
| | file: number of accepted records
|
| |
| integer(mpi), dimension(:), allocatable | dfd |
| | file: ndf sum
|
| |
| integer(mpi), dimension(:), allocatable | xfd |
| | file: max. record size
|
| |
| real(mps), dimension(:), allocatable | cfd |
| | file: chi2 sum
|
| |
| real(mps), dimension(:), allocatable | ofd |
| | file: option
|
| |
| real(mps), dimension(:), allocatable | wfd |
| | binary file: weight
|
| |
| integer(mpi), dimension(:,:), allocatable | sfd |
| | offset (1,..), length (2,..) of binary file name in tfd
|
| |
| integer(mpi), dimension(:), allocatable | yfd |
| | binary file: modification date
|
| |
| character(len=1024) | filnam |
| | name of steering file
|
| |
| integer(mpi) | nfnam |
| | length of sterring file name
|
| |
| character, dimension(:), allocatable | tfd |
| | file names (concatenation)
|
| |
| integer(mpi) | ifile |
| | current file (index)
|
| |
| integer(mpi) | nfiles |
| | number of files
|
| |
| integer(mpi) | nfilb |
| | number of binary files
|
| |
| integer(mpi) | nfilf |
| | number of Fortran binary files
|
| |
| integer(mpi) | nfilc |
| | number of C binary files
|
| |
| integer(mpi) | nfilw |
| | number of weighted binary files
|
| |
| integer(mpi) | ndimbuf =10000 |
| | default read buffer size (I/F words, half record length)
|
| |
1.9.8