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