mpmod.f90 File Reference

Data structures. More...

Go to the source code of this file.

Modules
module	mpmod
	Parameters, variables, dynamic arrays.

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

Detailed Description

Data structures.

Author

Claus Kleinwort, DESY, 2012 (Claus.nosp@m..Kle.nosp@m.inwor.nosp@m.t@de.nosp@m.sy.de)

Copyright

Copyright (c) 2012 - 2018 Deutsches Elektronen-Synchroton, Member of the Helmholtz Association, (DESY), HAMBURG, GERMANY

This library is free software; you can redistribute it and/or modify it under the terms of the GNU Library General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public License for more details.

You should have received a copy of the GNU Library General Public License along with this program (see the file COPYING.LIB for more details); if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

Definition in file mpmod.f90.