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:

FWDS model coefficients

RETURNED:

BKWDS inverse model

J I status: 0 = OK, -1 = no inverse

NOTES:

1.

The models relate two sets of $[\,x,y\,]$ coordinates as follows. Naming the six elements of FWDS a,b,c,d,e & f, where two sets of coordinates [x₁,y₁] and $[x_{2},y_{2}\,]$ are related thus:

x₂ = a + bx₁ + cy₁
y₂ = d + ex₁ + fy₁

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

x₁ = p + qx₂ + ry₂
y₁ = s + tx₂ + uy₂

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 on the VAX and many other systems even if this rule is violated, something that is not, however, recommended.

4.

See also sla_FITXY, sla_PXY, sla_XY2XY, sla_DCMPF.