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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
|
#ifndef MINC_BASIC_HEADER_FILE
#define MINC_BASIC_HEADER_FILE
/* ----------------------------- MNI Header -----------------------------------
@NAME : minc_basic.h
@INPUT :
@OUTPUT :
@RETURNS :
@DESCRIPTION: Constants and macros for private use by MINC routines.
@METHOD :
@GLOBALS :
@CALLS :
@CREATED : August 28, 1992 (Peter Neelin)
@MODIFIED :
* $Log: minc_basic.h,v $
* Revision 6.8 2010-03-27 15:09:28 rotor
* * back to 1000000
*
* Revision 6.7 2010-03-02 12:23:14 rotor
* * ported HDF calls to 1.8.x
* * Makefile.am: updated for minccmp
*
* Revision 6.6 2008/04/11 05:15:00 rotor
* * rewrote error code (Claude) to remove global defs that were
* causing build problems with DYLIB on OSX
*
* Revision 6.5 2008/01/12 01:05:37 rotor
* * initial commits from Steve to remove warnings
*
* Revision 6.4 2007/08/09 17:05:25 rotor
* * added some fixes of Claudes for chunking and internal compression
*
* Revision 6.3 2004/04/27 15:48:15 bert
* Minor changes
*
* Revision 6.2 2001/04/17 18:40:13 neelin
* Modifications to work with NetCDF 3.x
* In particular, changed NC_LONG to NC_INT (and corresponding longs to ints).
* Changed NC_UNSPECIFIED to NC_NAT.
* A few fixes to the configure script.
*
* Revision 6.1 1999/10/19 14:45:08 neelin
* Fixed Log subsitutions for CVS
*
* Revision 6.0 1997/09/12 13:24:54 neelin
* Release of minc version 0.6
*
* Revision 5.0 1997/08/21 13:25:53 neelin
* Release of minc version 0.5
*
* Revision 4.0 1997/05/07 20:07:52 neelin
* Release of minc version 0.4
*
* Revision 3.0 1995/05/15 19:33:12 neelin
* Release of minc version 0.3
*
* Revision 2.0 1994/09/28 10:38:01 neelin
* Release of minc version 0.2
*
* Revision 1.8 94/09/28 10:37:26 neelin
* Pre-release
*
* Revision 1.7 93/10/28 10:18:23 neelin
* Added FILLVALUE_EPSILON for doing fillvalue checking in icv's.
*
* Revision 1.6 93/08/11 12:06:37 neelin
* Added RCS logging in source.
*
@COPYRIGHT :
Copyright 1993 Peter Neelin, McConnell Brain Imaging Centre,
Montreal Neurological Institute, McGill University.
Permission to use, copy, modify, and distribute this
software and its documentation for any purpose and without
fee is hereby granted, provided that the above copyright
notice appear in all copies. The author and McGill University
make no representations about the suitability of this
software for any purpose. It is provided "as is" without
express or implied warranty.
@RCSID : $Header: /private-cvsroot/minc/libsrc/minc_basic.h,v 6.8 2010-03-27 15:09:28 rotor Exp $ MINC (MNI)
---------------------------------------------------------------------------- */
#include <minc.h>
/* --------- MINC specific constants ------------- */
/* Maximum buffer size for conversions. Should not be a power of 2 - this
can cause poor performance on some systems (e.g. SGI) due to caching-
related inefficiencies */
#define MI_MAX_VAR_BUFFER_SIZE 1000000
/* Possible values for sign of a value */
#define MI_PRIV_DEFSIGN 0
#define MI_PRIV_SIGNED 1
#define MI_PRIV_UNSIGNED 2
/* Operations for MI_varaccess */
#define MI_PRIV_GET 10
#define MI_PRIV_PUT 11
/* Suffix for dimension width variable names */
#define MI_WIDTH_SUFFIX "-width"
/* Epsilon for detecting fillvalues */
#define FILLVALUE_EPSILON (10.0 * FLT_EPSILON)
/* NetCDF routine name variable (for error logging) */
extern char *cdf_routine_name ; /* defined in globdef.c */
#define MI_NC_ROUTINE_VAR cdf_routine_name
/* Macros for logging errors. All routines should start with MI_SAVE_ROUTINE
and exit with MI_RETURN (which includes MI_RETURN_ERROR and
MI_CHK_ERROR). All the macros except MI_CHK_ERROR are single line
commands. MI_CHK_ERROR is in a block and so should not be followed by
a ';' */
#define MI_SAVE_ROUTINE_NAME(name) MI_save_routine_name(name)
#define MI_RETURN(value) \
return( MI_return() ? (value) : (value) )
#define MI_RETURN_ERROR(error) \
return( MI_return_error() ? (error) : (error) )
#define MI_LOG_PKG_ERROR2(p1,p2) MI_log_pkg_error2(p1, p2)
#define MI_LOG_PKG_ERROR3(p1,p2,p3) MI_log_pkg_error3(p1, p2, p3)
#define MI_LOG_SYS_ERROR1(p1) MI_log_sys_error1(p1)
#define MI_CHK_ERR(expr) {if ((expr)<0) MI_RETURN_ERROR(MI_ERROR);}
/* Macros for converting data types. These macros are compound statements,
so don't put a semi-colon after them. dvalue should be a double, type
is an int NetCDF type, sign is one of MI_PRIV_UNSIGNED and
MI_PRIV_SIGNED and ptr is a void pointer */
#define MI_TO_DOUBLE(dvalue, type, sign, ptr) \
switch (type) { \
case NC_BYTE : \
case NC_CHAR: \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = (double) *((unsigned char *) ptr); break; \
case MI_PRIV_SIGNED : \
dvalue = (double) *((signed char *) ptr); break; \
} \
break; \
case NC_SHORT : \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = (double) *((unsigned short *) ptr); break; \
case MI_PRIV_SIGNED : \
dvalue = (double) *((signed short *) ptr); break; \
} \
break; \
case NC_INT : \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = (double) *((unsigned int *) ptr); break; \
case MI_PRIV_SIGNED : \
dvalue = (double) *((signed int *) ptr); break; \
} \
break; \
case NC_FLOAT : \
dvalue = (double) *((float *) ptr); \
break; \
case NC_DOUBLE : \
dvalue = (double) *((double *) ptr); \
break; \
case NC_NAT : \
MI_LOG_PKG_ERROR2(MI_ERR_NONNUMERIC, \
"Attempt to convert NC_NAT value to double"); \
dvalue = 0; \
break; \
}
#define MI_FROM_DOUBLE(dvalue, type, sign, ptr) \
switch (type) { \
case NC_BYTE : \
case NC_CHAR : \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = MAX(0, dvalue); \
dvalue = MIN(UCHAR_MAX, dvalue); \
*((unsigned char *) ptr) = ROUND(dvalue); \
break; \
case MI_PRIV_SIGNED : \
dvalue = MAX(SCHAR_MIN, dvalue); \
dvalue = MIN(SCHAR_MAX, dvalue); \
*((signed char *) ptr) = ROUND(dvalue); \
break; \
} \
break; \
case NC_SHORT : \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = MAX(0, dvalue); \
dvalue = MIN(USHRT_MAX, dvalue); \
*((unsigned short *) ptr) = ROUND(dvalue); \
break; \
case MI_PRIV_SIGNED : \
dvalue = MAX(SHRT_MIN, dvalue); \
dvalue = MIN(SHRT_MAX, dvalue); \
*((signed short *) ptr) = ROUND(dvalue); \
break; \
} \
break; \
case NC_INT : \
switch (sign) { \
case MI_PRIV_UNSIGNED : \
dvalue = MAX(0, dvalue); \
dvalue = MIN(UINT_MAX, dvalue); \
*((unsigned int *) ptr) = ROUND(dvalue); \
break; \
case MI_PRIV_SIGNED : \
dvalue = MAX(INT_MIN, dvalue); \
dvalue = MIN(INT_MAX, dvalue); \
*((signed int *) ptr) = ROUND(dvalue); \
break; \
} \
break; \
case NC_FLOAT : \
dvalue = MAX(-FLT_MAX,dvalue); \
*((float *) ptr) = MIN(FLT_MAX,dvalue); \
break; \
case NC_DOUBLE : \
*((double *) ptr) = dvalue; \
break; \
case NC_NAT : \
MI_LOG_PKG_ERROR2(MI_ERR_NONNUMERIC, \
"Attempt to convert to NC_NAT from double"); \
dvalue = 0; \
break; \
}
/**/
#define _(x) x /* For future gettext */
#endif
|
