mille.f90 File Reference

Write Millepede-II F-binary record. More...

Go to the source code of this file.

Functions/Subroutines
subroutine	mille (nlc, derlc, ngl, dergl, label, rmeas, sigma)
	Add data block to record. Called from user code.

---

Detailed Description

Write Millepede-II F-binary record.

Definition in file mille.f90.

---

Function/Subroutine Documentation

subroutine mille	(	integer, intent(in)	nlc,
		real, dimension(nlc), intent(in)	derlc,
		integer, intent(in)	ngl,
		real, dimension(ngl), intent(in)	dergl,
		integer, dimension(ngl), intent(in)	label,
		real, intent(in)	rmeas,
		real, intent(in)	sigma
	)

Add data block to record. Called from user code.

    CALL MILLE(...)       ! measured value, derivatives (one set)
    CALL ENDLE            ! complete, write record (many sets)
(or CALL KILLE            ! stop record)

The data transmitted by MILLE calls are collected in two arrays, a real array and an integer array, of same length. The collected data are written at the ENDLE call. The content of the arrays:

    real array              integer array
1   0.0                     error count (this record)
2   RMEAS, measured value   0                            JA
3   local derivative        index of local derivative
4   local derivative        index of local derivative
5    ...
6   SIGMA, error (>0)       0                            JB
    global derivative       label of global derivative
    global derivative       label of global derivative   IST
    RMEAS, measured value   0
    local derivative        index of local derivative
    local derivative        index of local derivative
    ...
    SIGMA, error            0
    global derivative       label of global derivative
    global derivative       label of global derivative
    ...
NR  global derivative       label of global derivative

The 0's in the integer array allow to recognize the start of a new set, the measured value and the error. The local and the global derivatives are inbetween, with a positive value in the integer array, the index of the local derivative or the label of the global derivative.

If more than one output unit is needed: duplicate this subroutine change the entry names to e.g. AMILLE, AENDLE, AKILLE and change the value of LUN and evtl. the dimension parameter in the parameter statements.

Parameters:

[in]	NLC	number of local derivatives
[in]	DERLC	local derivatives
[in]	NGL	number of global derivatives
[in]	DERGL	global derivatives
[in]	LABEL	labels for global derivatives
[in]	RMEAS	measurement
[in]	SIGMA	error of measurement

Definition at line 65 of file mille.f90.

Referenced by mptest(), and mptst2().