blob: 09b894e3c813c441850d0f1a02063814a59ba11c (
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
51
52
53
54
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <math/gsurfit.h>
$if (datatype == r)
include "gsurfitdef.h"
$else
include "dgsurfitdef.h"
$endif
# GSSCOEFF -- Procedure to set a particular coefficient.
# If the requested coefficient is undefined then the coefficient is not set.
$if (datatype == r)
procedure gsscoeff (sf, xorder, yorder, coeff)
$else
procedure dgsscoeff (sf, xorder, yorder, coeff)
$endif
pointer sf # pointer to the surface fitting descriptor
int xorder # X order of desired coefficent
int yorder # Y order of desired coefficent
PIXEL coeff # Coefficient value
int i, n, maxorder, xincr
begin
if ((xorder > GS_XORDER(sf)) || (yorder > GS_YORDER(sf)))
return
switch (GS_XTERMS(sf)) {
case GS_XNONE:
if (yorder == 1)
n = xorder
else if (xorder == 1)
n = GS_NXCOEFF(sf) + yorder - 1
else
return
case GS_XHALF:
maxorder = max (GS_XORDER(sf) + 1, GS_YORDER(sf) + 1)
if ((xorder + yorder) > maxorder)
return
n = xorder
xincr = GS_XORDER(sf)
do i = 2, yorder {
n = n + xincr
if ((i + GS_XORDER(sf) + 1) > maxorder)
xincr = xincr - 1
}
case GS_XFULL:
n = xorder + (yorder - 1) * GS_NXCOEFF(sf)
}
COEFF(GS_COEFF(sf) + n - 1) = coeff
end
|
