This routine replaces the supplied Hermitian image with one in which the real terms have been inverted (i.e. replaced by their reciprocal) and the imaginary terms have been set to zero.

To form 1/H where H is a general complex Hermitian image, the numerator and denominator are both multiplied by the complex conjugate of H (" H$\ast$" ). The denominator then becomes purely real and equal to the modulus squared of H. The reciprocal of the denominator is then taken using this routine and the resulting array is multiplied by H$\ast$. Thus the steps are:

1) Call KPG1_HCONx to form H$\ast$, the complex conjugate of H. 2) Call KPG1_HMLTx to multiply H by H$\ast$. This gives a purely real image holding the squared modulus of H. 3) Call KPG1_HRCPx to take the reciprocal of H.H$\ast$. 4) Call KPG1_HMLTx to multiply the reciprocal of H.H$\ast$ by H$\ast$.

CALL KPG1_HRCPx( M, N, BADVAL, D, STATUS )

M = INTEGER (Given)

N = INTEGER (Given)

BADVAL = ? (Given)

D( M, N ) = ? (Given and Returned)

STATUS = INTEGER (Given and Returned)

• There is a routine for the double precision and real data types: replace " x" in the routine names by D or R as appropriate. The BADVAL and D arguments must have the data type specified.