NormalsEliminate

This program eliminates parameters from a system of inputfileNormalEquations. To just remove (cutting out) parameters use NormalsReorder.

The remainingParameters allows the selection of parameters that will remain, all others will be eliminated. The order of remaining parameters can be modified via the parameter selection. Block size of the output normal matrix can be adjusted with outBlockSize. If it is set to zero, the outputfileNormalEquation is written to a single block file.

For example the normal equations are divided into two groups of parameters $\hat{\M x}_1$ and $\hat{\M x}_2$ according to \[ \begin{pmatrix} \M N_{11} & \M N_{12} \\ \M N_{21} & \M N_{22} \end{pmatrix} \begin{pmatrix} \hat{\M x}_1 \\ \hat{\M x}_2 \end{pmatrix} = \begin{pmatrix} \M n_1 \\ \M n_2 \end{pmatrix}. \]and $\hat{\M x}_2$ shall be eliminated, the reduced system of normal equations is given by \[ \bar{\M N}\hat{\M x} = \bar{\M n} \qquad\text{with}\qquad \bar{\M N}=\M N_{11}-\M N_{12}\M N_{22}^{-1}\M N_{12}^T \qquad\text{and}\qquad\bar{\M n} = \M n_1 - \M N_{12}\M N_{22}^{-1}\M n_2. \] See also NormalsReorder.

NameTypeAnnotation
outputfileNormalEquation
filename
inputfileNormalEquation
filename
remainingParameters
parameterSelectorparameter order/selection of output normal equations
outBlockSize
uintblock size for distributing the normal equations, 0: one block
This program is parallelized.