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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include "mwcs.h"
# MW_SLTERMD -- Set the Lterm, double precision version. Since all floating
# data is stored as double internally, we merely copy the data in.
procedure mw_sltermd (mw, ltm, ltv, ndim)
pointer mw #I pointer to MWCS descriptor
double ltm[ndim,ndim] #I linear transformation matrix
double ltv[ndim] #I translation vector
int ndim #I dimensionality of system
pointer mw_allocd()
errchk syserrs, mw_allocd
begin
# The dimensionality of the data must match that of the current Lterm.
if (ndim != MI_NDIM(mw))
call syserrs (SYS_MWNDIM, "mw_sltermd")
# Copy in the data. Cobber the old data if the Lterm has been set,
# otherwise allocate space in the global data area.
if (MI_LTM(mw) == NULL)
MI_LTM(mw) = mw_allocd (mw, ndim*ndim)
call amovd (ltm, D(mw,MI_LTM(mw)), ndim*ndim)
if (MI_LTV(mw) == NULL)
MI_LTV(mw) = mw_allocd (mw, ndim)
call amovd (ltv, D(mw,MI_LTV(mw)), ndim)
end
|
