SLA_INVF

Invert Linear Model

ACTION:

Invert a linear model of the type produced by the sla_FITXY routine.

CALL:

CALL sla_INVF (FWDS,BKWDS,J)

GIVEN:

RETURNED:

NOTES:

(1)

The models relate two sets of $\left[x,y\right]$ coordinates as follows. Naming the six elements of FWDS $a,b,c,d,e$ & $f$, where two sets of coordinates $\left[x_1,y_1\right]$ and $\left[x_2,y_2\right]$ are related thus:

$x_2=a+b{x}_1+c{y}_1$
$y_2=d+e{x}_1+f{y}_1$

The present routine generates a new set of coefficients $p,q,r,s,t$ & $u$ (the array BKWDS) such that:

$x_1=p+q{x}_2+r{y}_2$
$y_1=s+t{x}_2+u{y}_2$

(2)

Two successive calls to this routine will deliver a set of coefficients equal to the starting values.

(3)

To comply with the ANSI Fortran 77 standard, FWDS and BKWDS must not be the same array. The routine is, in fact, coded so as to work properly with many Fortran compilers even if this rule is violated, something that is not, however, recommended.

(4)

See also sla_FITXY, sla_PXY, sla_XY2XY, sla_DCMPF.