Appendix A: Summary of the implementation method

In the following appendices a set of algorithms is presented as they were used in varies applications of the methods described in the present paper. The algorithms are built using a vectorial storage of symmetric and upper or lower triangular arrays, implicitly using the knowledge of symmetry or zero elements. This reduces storage needs and in some cases facilitates implementation. All algorithms are given in Fortran 77.

As a first step, the abscissae for all stars to which a combined and/or constrained solution is to be applied are extracted from the intermediate astrometric data file, where they are sorted on star number and for each star on orbit number. The abscissa file thus created is then sorted on orbit number and within each orbit on star number.