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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <mach.h>
# IMPUTD -- Put an image header parameter of type double.
procedure imputd (im, key, dval)
pointer im # image descriptor
char key[ARB] # parameter to be set
double dval # double precision value
int junk, i
pointer sp, sval
int dtoc(), strlen()
begin
call smark (sp)
call salloc (sval, SZ_FNAME, TY_CHAR)
# Reduce the precision of the encoded value if necessary to fit in
# the FITS value field. Start with NDIGITS_DP-1 as the precision
# estimate NDIGITS_DP is only approximate, and if we make up half a
# digit of precision the result can be 1.00000000000000001 instead
# of 1.0.
for (i=NDIGITS_DP-1; i >= NDIGITS_RP; i=i-1) {
junk = dtoc (dval, Memc[sval], SZ_FNAME, i, 'g', SZ_FNAME)
if (strlen (Memc[sval]) < 20)
break
}
# Write the new value to the header.
call impstr (im, key, Memc[sval])
call sfree (sp)
end
|
