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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
/*============================================================================
WCSLIB 7.12 - an implementation of the FITS WCS standard.
Copyright (C) 1995-2022, Mark Calabretta
This file is part of WCSLIB.
WCSLIB is free software: you can redistribute it and/or modify it under the
terms of the GNU Lesser General Public License as published by the Free
Software Foundation, either version 3 of the License, or (at your option)
any later version.
WCSLIB is distributed in the hope that it will be useful, but WITHOUT ANY
WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
more details.
You should have received a copy of the GNU Lesser General Public License
along with WCSLIB. If not, see http://www.gnu.org/licenses.
Author: Mark Calabretta, Australia Telescope National Facility, CSIRO.
http://www.atnf.csiro.au/people/Mark.Calabretta
$Id: log.c,v 7.12 2022/09/09 04:57:58 mcalabre Exp $
*===========================================================================*/
#include <math.h>
#include "log.h"
// Map status return value to message.
const char *log_errmsg[] = {
"Success",
"",
"Invalid log-coordinate reference value",
"One or more of the x coordinates were invalid",
"One or more of the world coordinates were invalid"};
//----------------------------------------------------------------------------
int logx2s(
double crval,
int nx,
int sx,
int slogc,
const double x[],
double logc[],
int stat[])
{
register int ix;
register int *statp;
register const double *xp;
register double *logcp;
if (crval <= 0.0) {
return LOGERR_BAD_LOG_REF_VAL;
}
xp = x;
logcp = logc;
statp = stat;
for (ix = 0; ix < nx; ix++, xp += sx, logcp += slogc, statp++) {
*logcp = crval * exp((*xp) / crval);
*statp = 0;
}
return 0;
}
//----------------------------------------------------------------------------
int logs2x(
double crval,
int nlogc,
int slogc,
int sx,
const double logc[],
double x[],
int stat[])
{
int status;
register int ilogc;
register int *statp;
register const double *logcp;
register double *xp;
if (crval <= 0.0) {
return LOGERR_BAD_LOG_REF_VAL;
}
xp = x;
logcp = logc;
statp = stat;
status = 0;
for (ilogc = 0; ilogc < nlogc; ilogc++, logcp += slogc, xp += sx, statp++) {
if (*logcp > 0.0) {
*xp = crval * log(*logcp / crval);
*statp = 0;
} else {
status = LOGERR_BAD_WORLD;
*statp = 1;
}
}
return status;
}
|
