blob: 1a6b4c1c15e2157ead4e329c97a839e311bafd99 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
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
|
