1 files changed, 50 insertions, 0 deletions
diff --git a/math/iminterp/msirestore.x b/math/iminterp/msirestore.x
new file mode 100644
index 00000000..1a6b4c1c
--- /dev/null
+++ b/math/iminterp/msirestore.x
@@ -0,0 +1,50 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include "im2interpdef.h"
+include	<math/iminterp.h>
+
+# MSIRESTORE -- Procedure to restore the interpolant stored by MSISAVE
+# for use by MSIEVAL, MSIVECTOR, MSIDER and MSIGRL.
+
+procedure msirestore (msi, interpolant)
+
+pointer	msi			# interpolant descriptor
+real	interpolant[ARB]	# array containing the interpolant
+
+int	interp_type, npix
+
+begin
+	interp_type = nint (MSI_SAVETYPE(interpolant))
+	if (interp_type < 1 || interp_type > II_NTYPES)
+	    call error (0, "MSIRESTORE: Unknown interpolant type.")
+
+	# allocate the interpolant descriptor structure and restore
+	# interpolant parameters
+	call malloc (msi, LEN_MSISTRUCT, TY_STRUCT)
+	MSI_TYPE(msi) = interp_type
+	MSI_NSINC(msi) = nint (MSI_SAVENSINC(interpolant))
+	MSI_NXINCR(msi) = nint (MSI_SAVENXINCR(interpolant))
+	MSI_NYINCR(msi) = nint (MSI_SAVENYINCR(interpolant))
+	MSI_XSHIFT(msi) = MSI_SAVEXSHIFT(interpolant)
+	MSI_YSHIFT(msi) = MSI_SAVEYSHIFT(interpolant)
+	MSI_XPIXFRAC(msi) = MSI_SAVEXPIXFRAC(interpolant)
+	MSI_YPIXFRAC(msi) = MSI_SAVEYPIXFRAC(interpolant)
+	MSI_NXCOEFF(msi) = nint (MSI_SAVENXCOEFF(interpolant))
+	MSI_NYCOEFF(msi) = nint (MSI_SAVENYCOEFF(interpolant))
+	MSI_FSTPNT(msi) = nint (MSI_SAVEFSTPNT(interpolant))
+	MSI_BADVAL(msi) = MSI_SAVEBADVAL(interpolant)
+
+	# allocate space for and restore coefficients
+	call malloc (MSI_COEFF(msi), MSI_NXCOEFF(msi) * MSI_NYCOEFF(msi),
+	    TY_REAL)
+	call amovr (interpolant[1+MSI_SAVECOEFF], COEFF(MSI_COEFF(msi)),
+	    MSI_NXCOEFF(msi) * MSI_NYCOEFF(msi))
+
+	# allocate space for and restore the look-up table
+	if (MSI_NXINCR(msi) > 0 && MSI_NYINCR(msi) > 0) {
+	    npix = (2.0 * MSI_NSINC(msi) + 1) ** 2 * MSI_NXINCR(msi) *
+		MSI_NYINCR(msi)
+	    call amovr (interpolant[1+MSI_SAVECOEFF+MSI_NXCOEFF(msi) *
+	        MSI_NYCOEFF(msi)], LTABLE(MSI_LTABLE(msi)), npix)
+	}
+end