cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 1 | /* |
Cristy | 7ce65e7 | 2015-12-12 18:03:16 -0500 | [diff] [blame] | 2 | Copyright 1999-2016 ImageMagick Studio LLC, a non-profit organization |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 3 | dedicated to making software imaging solutions freely available. |
| 4 | |
| 5 | You may not use this file except in compliance with the License. |
| 6 | obtain a copy of the License at |
| 7 | |
| 8 | http://www.imagemagick.org/script/license.php |
| 9 | |
| 10 | Unless required by applicable law or agreed to in writing, software |
| 11 | distributed under the License is distributed on an "AS IS" BASIS, |
| 12 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. |
| 13 | See the License for the specific language governing permissions and |
| 14 | limitations under the License. |
| 15 | |
| 16 | MagickCore quantum inline methods. |
| 17 | */ |
| 18 | #ifndef _MAGICKCORE_QUANTUM_PRIVATE_H |
| 19 | #define _MAGICKCORE_QUANTUM_PRIVATE_H |
| 20 | |
cristy | 6398ec7 | 2013-11-28 02:00:27 +0000 | [diff] [blame] | 21 | #include "MagickCore/cache.h" |
| 22 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 23 | #if defined(__cplusplus) || defined(c_plusplus) |
| 24 | extern "C" { |
| 25 | #endif |
| 26 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 27 | typedef struct _QuantumState |
| 28 | { |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 29 | double |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 30 | inverse_scale; |
| 31 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 32 | unsigned int |
cristy | 6b825f9 | 2010-06-04 02:01:10 +0000 | [diff] [blame] | 33 | pixel; |
| 34 | |
| 35 | size_t |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 36 | bits; |
| 37 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 38 | const unsigned int |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 39 | *mask; |
| 40 | } QuantumState; |
| 41 | |
| 42 | struct _QuantumInfo |
| 43 | { |
cristy | bb50337 | 2010-05-27 20:51:26 +0000 | [diff] [blame] | 44 | size_t |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 45 | depth, |
| 46 | quantum; |
| 47 | |
| 48 | QuantumFormatType |
| 49 | format; |
| 50 | |
| 51 | double |
| 52 | minimum, |
| 53 | maximum, |
| 54 | scale; |
| 55 | |
| 56 | size_t |
| 57 | pad; |
| 58 | |
| 59 | MagickBooleanType |
| 60 | min_is_white, |
| 61 | pack; |
| 62 | |
| 63 | QuantumAlphaType |
| 64 | alpha_type; |
| 65 | |
cristy | bb50337 | 2010-05-27 20:51:26 +0000 | [diff] [blame] | 66 | size_t |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 67 | number_threads; |
| 68 | |
| 69 | unsigned char |
| 70 | **pixels; |
| 71 | |
| 72 | size_t |
| 73 | extent; |
| 74 | |
cristy | 32c6843 | 2012-01-12 15:06:28 +0000 | [diff] [blame] | 75 | EndianType |
| 76 | endian; |
| 77 | |
| 78 | QuantumState |
| 79 | state; |
| 80 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 81 | SemaphoreInfo |
| 82 | *semaphore; |
| 83 | |
cristy | bb50337 | 2010-05-27 20:51:26 +0000 | [diff] [blame] | 84 | size_t |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 85 | signature; |
| 86 | }; |
| 87 | |
cristy | 7963022 | 2012-01-14 01:38:37 +0000 | [diff] [blame] | 88 | extern MagickPrivate void |
| 89 | ResetQuantumState(QuantumInfo *); |
| 90 | |
cristy | bb50337 | 2010-05-27 20:51:26 +0000 | [diff] [blame] | 91 | static inline MagickSizeType GetQuantumRange(const size_t depth) |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 92 | { |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 93 | MagickSizeType |
| 94 | one; |
| 95 | |
| 96 | one=1; |
| 97 | return((MagickSizeType) ((one << (depth-1))+((one << (depth-1))-1))); |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 98 | } |
| 99 | |
| 100 | static inline float HalfToSinglePrecision(const unsigned short half) |
cristy | 5ab5fc8 | 2010-01-05 20:52:48 +0000 | [diff] [blame] | 101 | { |
cristy | 8a6a980 | 2010-04-14 13:22:18 +0000 | [diff] [blame] | 102 | #define ExponentBias (127-15) |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 103 | #define ExponentMask 0x7c00 |
| 104 | #define ExponentShift 23 |
cristy | bccb64f | 2010-01-10 21:16:36 +0000 | [diff] [blame] | 105 | #define SignBitShift 31 |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 106 | #define SignificandShift 13 |
| 107 | #define SignificandMask 0x00000400 |
cristy | 5ab5fc8 | 2010-01-05 20:52:48 +0000 | [diff] [blame] | 108 | |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 109 | typedef union _SinglePrecision |
| 110 | { |
| 111 | unsigned int |
| 112 | fixed_point; |
| 113 | |
| 114 | float |
| 115 | single_precision; |
| 116 | } SinglePrecision; |
| 117 | |
| 118 | register unsigned int |
| 119 | exponent, |
| 120 | significand, |
| 121 | sign_bit; |
| 122 | |
| 123 | SinglePrecision |
| 124 | map; |
| 125 | |
| 126 | unsigned int |
| 127 | value; |
| 128 | |
| 129 | /* |
cristy | bccb64f | 2010-01-10 21:16:36 +0000 | [diff] [blame] | 130 | The IEEE 754 standard specifies half precision as having: |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 131 | |
| 132 | Sign bit: 1 bit |
| 133 | Exponent width: 5 bits |
| 134 | Significand precision: 11 (10 explicitly stored) |
| 135 | */ |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 136 | sign_bit=(unsigned int) ((half >> 15) & 0x00000001); |
| 137 | exponent=(unsigned int) ((half >> 10) & 0x0000001f); |
| 138 | significand=(unsigned int) (half & 0x000003ff); |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 139 | if (exponent == 0) |
| 140 | { |
cristy | 92caafa | 2010-05-10 12:56:29 +0000 | [diff] [blame] | 141 | if (significand == 0) |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 142 | value=sign_bit << SignBitShift; |
cristy | 92caafa | 2010-05-10 12:56:29 +0000 | [diff] [blame] | 143 | else |
| 144 | { |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 145 | while ((significand & SignificandMask) == 0) |
| 146 | { |
| 147 | significand<<=1; |
| 148 | exponent--; |
| 149 | } |
cristy | 92caafa | 2010-05-10 12:56:29 +0000 | [diff] [blame] | 150 | exponent++; |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 151 | significand&=(~SignificandMask); |
| 152 | exponent+=ExponentBias; |
| 153 | value=(sign_bit << SignBitShift) | (exponent << ExponentShift) | |
| 154 | (significand << SignificandShift); |
cristy | 92caafa | 2010-05-10 12:56:29 +0000 | [diff] [blame] | 155 | } |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 156 | } |
| 157 | else |
| 158 | if (exponent == SignBitShift) |
| 159 | { |
| 160 | value=(sign_bit << SignBitShift) | 0x7f800000; |
| 161 | if (significand != 0) |
| 162 | value|=(significand << SignificandShift); |
| 163 | } |
| 164 | else |
| 165 | { |
| 166 | exponent+=ExponentBias; |
| 167 | significand<<=SignificandShift; |
| 168 | value=(sign_bit << SignBitShift) | (exponent << ExponentShift) | |
| 169 | significand; |
| 170 | } |
| 171 | map.fixed_point=value; |
| 172 | return(map.single_precision); |
| 173 | } |
cristy | 5ab5fc8 | 2010-01-05 20:52:48 +0000 | [diff] [blame] | 174 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 175 | static inline unsigned char *PopCharPixel(const unsigned char pixel, |
| 176 | unsigned char *pixels) |
| 177 | { |
| 178 | *pixels++=pixel; |
| 179 | return(pixels); |
| 180 | } |
| 181 | |
| 182 | static inline unsigned char *PopLongPixel(const EndianType endian, |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 183 | const unsigned int pixel,unsigned char *pixels) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 184 | { |
| 185 | register unsigned int |
| 186 | quantum; |
| 187 | |
| 188 | quantum=(unsigned int) pixel; |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 189 | if (endian == LSBEndian) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 190 | { |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 191 | *pixels++=(unsigned char) (quantum); |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 192 | *pixels++=(unsigned char) (quantum >> 8); |
| 193 | *pixels++=(unsigned char) (quantum >> 16); |
| 194 | *pixels++=(unsigned char) (quantum >> 24); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 195 | return(pixels); |
| 196 | } |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 197 | *pixels++=(unsigned char) (quantum >> 24); |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 198 | *pixels++=(unsigned char) (quantum >> 16); |
| 199 | *pixels++=(unsigned char) (quantum >> 8); |
| 200 | *pixels++=(unsigned char) (quantum); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 201 | return(pixels); |
| 202 | } |
| 203 | |
| 204 | static inline unsigned char *PopShortPixel(const EndianType endian, |
| 205 | const unsigned short pixel,unsigned char *pixels) |
| 206 | { |
| 207 | register unsigned int |
| 208 | quantum; |
| 209 | |
| 210 | quantum=pixel; |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 211 | if (endian == LSBEndian) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 212 | { |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 213 | *pixels++=(unsigned char) (quantum); |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 214 | *pixels++=(unsigned char) (quantum >> 8); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 215 | return(pixels); |
| 216 | } |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 217 | *pixels++=(unsigned char) (quantum >> 8); |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 218 | *pixels++=(unsigned char) (quantum); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 219 | return(pixels); |
| 220 | } |
| 221 | |
| 222 | static inline const unsigned char *PushCharPixel(const unsigned char *pixels, |
| 223 | unsigned char *pixel) |
| 224 | { |
| 225 | *pixel=(*pixels++); |
| 226 | return(pixels); |
| 227 | } |
| 228 | |
| 229 | static inline const unsigned char *PushLongPixel(const EndianType endian, |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 230 | const unsigned char *pixels,unsigned int *pixel) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 231 | { |
Cristy | a1a36e7 | 2015-09-27 13:19:14 -0400 | [diff] [blame] | 232 | register unsigned int |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 233 | quantum; |
| 234 | |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 235 | if (endian == LSBEndian) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 236 | { |
Cristy | a1a36e7 | 2015-09-27 13:19:14 -0400 | [diff] [blame] | 237 | quantum=((unsigned int) *pixels++); |
| 238 | quantum|=((unsigned int) *pixels++ << 8); |
| 239 | quantum|=((unsigned int) *pixels++ << 16); |
| 240 | quantum|=((unsigned int) *pixels++ << 24); |
| 241 | *pixel=quantum; |
cristy | 70d7eae | 2012-11-17 19:40:21 +0000 | [diff] [blame] | 242 | return(pixels); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 243 | } |
Cristy | a1a36e7 | 2015-09-27 13:19:14 -0400 | [diff] [blame] | 244 | quantum=((unsigned int) *pixels++ << 24); |
| 245 | quantum|=((unsigned int) *pixels++ << 16); |
| 246 | quantum|=((unsigned int) *pixels++ << 8); |
| 247 | quantum|=((unsigned int) *pixels++); |
| 248 | *pixel=quantum; |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 249 | return(pixels); |
| 250 | } |
| 251 | |
| 252 | static inline const unsigned char *PushShortPixel(const EndianType endian, |
| 253 | const unsigned char *pixels,unsigned short *pixel) |
| 254 | { |
| 255 | register unsigned int |
| 256 | quantum; |
| 257 | |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 258 | if (endian == LSBEndian) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 259 | { |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 260 | quantum=(unsigned int) *pixels++; |
| 261 | quantum|=(unsigned int) (*pixels++ << 8); |
cristy | 70d7eae | 2012-11-17 19:40:21 +0000 | [diff] [blame] | 262 | *pixel=(unsigned short) (quantum & 0xffff); |
| 263 | return(pixels); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 264 | } |
cristy | 1715f79 | 2012-11-22 16:56:12 +0000 | [diff] [blame] | 265 | quantum=(unsigned int) (*pixels++ << 8); |
| 266 | quantum|=(unsigned int) *pixels++; |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 267 | *pixel=(unsigned short) (quantum & 0xffff); |
| 268 | return(pixels); |
| 269 | } |
| 270 | |
| 271 | static inline Quantum ScaleAnyToQuantum(const QuantumAny quantum, |
| 272 | const QuantumAny range) |
| 273 | { |
cristy | da950bb | 2015-07-18 23:17:28 +0000 | [diff] [blame] | 274 | if (quantum > range) |
| 275 | return(QuantumRange); |
cristy | b3963e6 | 2010-05-09 21:36:40 +0000 | [diff] [blame] | 276 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | a19f1d7 | 2012-08-07 18:24:38 +0000 | [diff] [blame] | 277 | return((Quantum) (((double) QuantumRange*quantum)/range+0.5)); |
cristy | b3963e6 | 2010-05-09 21:36:40 +0000 | [diff] [blame] | 278 | #else |
cristy | a19f1d7 | 2012-08-07 18:24:38 +0000 | [diff] [blame] | 279 | return((Quantum) (((double) QuantumRange*quantum)/range)); |
cristy | b3963e6 | 2010-05-09 21:36:40 +0000 | [diff] [blame] | 280 | #endif |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 281 | } |
| 282 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 283 | static inline QuantumAny ScaleQuantumToAny(const Quantum quantum, |
| 284 | const QuantumAny range) |
| 285 | { |
cristy | a19f1d7 | 2012-08-07 18:24:38 +0000 | [diff] [blame] | 286 | return((QuantumAny) (((double) range*quantum)/QuantumRange+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 287 | } |
| 288 | |
| 289 | #if (MAGICKCORE_QUANTUM_DEPTH == 8) |
| 290 | static inline Quantum ScaleCharToQuantum(const unsigned char value) |
| 291 | { |
| 292 | return((Quantum) value); |
| 293 | } |
| 294 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 295 | static inline Quantum ScaleLongToQuantum(const unsigned int value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 296 | { |
| 297 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 298 | return((Quantum) ((value+8421504UL)/16843009UL)); |
| 299 | #else |
| 300 | return((Quantum) (value/16843009.0)); |
| 301 | #endif |
| 302 | } |
| 303 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 304 | static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value) |
| 305 | { |
| 306 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 307 | return((Quantum) ((value+MagickULLConstant(551911719039))/ |
| 308 | MagickULLConstant(1103823438079))); |
| 309 | #else |
| 310 | return((Quantum) (value/1103823438079.0)); |
| 311 | #endif |
| 312 | } |
| 313 | |
| 314 | static inline Quantum ScaleMapToQuantum(const MagickRealType value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 315 | { |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 316 | if (value <= 0.0) |
| 317 | return((Quantum) 0); |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 318 | if (value >= MaxMap) |
cristy | 6e963d8 | 2012-06-19 15:23:24 +0000 | [diff] [blame] | 319 | return(QuantumRange); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 320 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 321 | return((Quantum) (value+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 322 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 323 | return((Quantum) value); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 324 | #endif |
| 325 | } |
| 326 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 327 | static inline unsigned int ScaleQuantumToLong(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 328 | { |
| 329 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 330 | return((unsigned int) (16843009UL*quantum)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 331 | #else |
| 332 | if (quantum <= 0.0) |
| 333 | return(0UL); |
| 334 | if ((16843009.0*quantum) >= 4294967295.0) |
| 335 | return(4294967295UL); |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 336 | return((unsigned int) (16843009.0*quantum+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 337 | #endif |
| 338 | } |
| 339 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 340 | static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum) |
| 341 | { |
| 342 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 343 | return((MagickSizeType) (MagickULLConstant(551911719039)*quantum)); |
| 344 | #else |
| 345 | if (quantum <= 0.0) |
| 346 | return(0UL); |
| 347 | if ((551911719039.0*quantum) >= 18446744073709551615.0) |
| 348 | return(MagickULLConstant(18446744073709551615)); |
| 349 | return((MagickSizeType) (1103823438079.0*quantum+0.5)); |
| 350 | #endif |
| 351 | } |
| 352 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 353 | static inline unsigned int ScaleQuantumToMap(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 354 | { |
| 355 | if (quantum >= (Quantum) MaxMap) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 356 | return((unsigned int) MaxMap); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 357 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 358 | return((unsigned int) quantum); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 359 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 360 | if (quantum < 0.0) |
| 361 | return(0UL); |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 362 | return((unsigned int) (quantum+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 363 | #endif |
| 364 | } |
| 365 | |
| 366 | static inline unsigned short ScaleQuantumToShort(const Quantum quantum) |
| 367 | { |
| 368 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 369 | return((unsigned short) (257UL*quantum)); |
| 370 | #else |
| 371 | if (quantum <= 0.0) |
| 372 | return(0); |
cristy | cd817db | 2010-05-09 03:05:41 +0000 | [diff] [blame] | 373 | if ((257.0*quantum) >= 65535.0) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 374 | return(65535); |
cristy | cd817db | 2010-05-09 03:05:41 +0000 | [diff] [blame] | 375 | return((unsigned short) (257.0*quantum+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 376 | #endif |
| 377 | } |
| 378 | |
| 379 | static inline Quantum ScaleShortToQuantum(const unsigned short value) |
| 380 | { |
| 381 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | e90d740 | 2010-03-14 18:21:29 +0000 | [diff] [blame] | 382 | return((Quantum) ((value+128U)/257U)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 383 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 384 | return((Quantum) (value/257.0)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 385 | #endif |
| 386 | } |
| 387 | #elif (MAGICKCORE_QUANTUM_DEPTH == 16) |
| 388 | static inline Quantum ScaleCharToQuantum(const unsigned char value) |
| 389 | { |
| 390 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | e90d740 | 2010-03-14 18:21:29 +0000 | [diff] [blame] | 391 | return((Quantum) (257U*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 392 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 393 | return((Quantum) (257.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 394 | #endif |
| 395 | } |
| 396 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 397 | static inline Quantum ScaleLongToQuantum(const unsigned int value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 398 | { |
| 399 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 400 | return((Quantum) ((value+MagickULLConstant(32768))/ |
| 401 | MagickULLConstant(65537))); |
| 402 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 403 | return((Quantum) (value/65537.0)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 404 | #endif |
| 405 | } |
| 406 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 407 | static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value) |
| 408 | { |
| 409 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 410 | return((Quantum) ((value+MagickULLConstant(8421376))/ |
| 411 | MagickULLConstant(16842752))); |
| 412 | #else |
| 413 | return((Quantum) (value/16842752.0)); |
| 414 | #endif |
| 415 | } |
| 416 | |
| 417 | static inline Quantum ScaleMapToQuantum(const MagickRealType value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 418 | { |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 419 | if (value <= 0.0) |
| 420 | return((Quantum) 0); |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 421 | if (value >= MaxMap) |
cristy | 6e963d8 | 2012-06-19 15:23:24 +0000 | [diff] [blame] | 422 | return(QuantumRange); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 423 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 424 | return((Quantum) (value+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 425 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 426 | return((Quantum) value); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 427 | #endif |
| 428 | } |
| 429 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 430 | static inline unsigned int ScaleQuantumToLong(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 431 | { |
| 432 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 433 | return((unsigned int) (65537UL*quantum)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 434 | #else |
| 435 | if (quantum <= 0.0) |
| 436 | return(0UL); |
cristy | cd817db | 2010-05-09 03:05:41 +0000 | [diff] [blame] | 437 | if ((65537.0*quantum) >= 4294967295.0) |
cristy | 9f48ca6 | 2010-11-25 03:06:31 +0000 | [diff] [blame] | 438 | return(4294967295U); |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 439 | return((unsigned int) (65537.0*quantum+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 440 | #endif |
| 441 | } |
| 442 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 443 | static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum) |
| 444 | { |
| 445 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 446 | return((MagickSizeType) (MagickULLConstant(16842752)*quantum)); |
| 447 | #else |
| 448 | if (quantum <= 0.0) |
| 449 | return(0UL); |
| 450 | if ((65537.0*quantum) >= 18446744073709551615.0) |
| 451 | return(MagickULLConstant(18446744073709551615)); |
| 452 | return((MagickSizeType) (16842752.0*quantum+0.5)); |
| 453 | #endif |
| 454 | } |
| 455 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 456 | static inline unsigned int ScaleQuantumToMap(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 457 | { |
| 458 | if (quantum >= (Quantum) MaxMap) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 459 | return((unsigned int) MaxMap); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 460 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 461 | return((unsigned int) quantum); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 462 | #else |
cristy | a830757 | 2010-03-31 14:24:19 +0000 | [diff] [blame] | 463 | if (quantum < 0.0) |
| 464 | return(0UL); |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 465 | return((unsigned int) (quantum+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 466 | #endif |
| 467 | } |
| 468 | |
| 469 | static inline unsigned short ScaleQuantumToShort(const Quantum quantum) |
| 470 | { |
| 471 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 472 | return((unsigned short) quantum); |
| 473 | #else |
| 474 | if (quantum <= 0.0) |
| 475 | return(0); |
| 476 | if (quantum >= 65535.0) |
| 477 | return(65535); |
| 478 | return((unsigned short) (quantum+0.5)); |
| 479 | #endif |
| 480 | } |
| 481 | |
| 482 | static inline Quantum ScaleShortToQuantum(const unsigned short value) |
| 483 | { |
| 484 | return((Quantum) value); |
| 485 | } |
| 486 | #elif (MAGICKCORE_QUANTUM_DEPTH == 32) |
| 487 | static inline Quantum ScaleCharToQuantum(const unsigned char value) |
| 488 | { |
| 489 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 490 | return((Quantum) (16843009UL*value)); |
| 491 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 492 | return((Quantum) (16843009.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 493 | #endif |
| 494 | } |
| 495 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 496 | static inline Quantum ScaleLongToQuantum(const unsigned int value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 497 | { |
| 498 | return((Quantum) value); |
| 499 | } |
| 500 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 501 | static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value) |
| 502 | { |
| 503 | return((Quantum) value); |
| 504 | } |
| 505 | |
| 506 | static inline Quantum ScaleMapToQuantum(const MagickRealType value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 507 | { |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 508 | if (value <= 0.0) |
| 509 | return((Quantum) 0); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 510 | if (value >= (Quantum) MaxMap) |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 511 | return(QuantumRange); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 512 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 513 | return((Quantum) (65537.0*value+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 514 | #else |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 515 | return((Quantum) (65537.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 516 | #endif |
| 517 | } |
| 518 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 519 | static inline unsigned int ScaleQuantumToLong(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 520 | { |
cristy | 0402ad3 | 2010-05-11 13:57:55 +0000 | [diff] [blame] | 521 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 522 | return((unsigned int) quantum); |
cristy | 0402ad3 | 2010-05-11 13:57:55 +0000 | [diff] [blame] | 523 | #else |
cristy | da97dba | 2014-05-24 22:59:47 +0000 | [diff] [blame] | 524 | if (quantum <= 0.0) |
| 525 | return(0); |
| 526 | if ((quantum) >= 4294967295.0) |
| 527 | return(4294967295); |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 528 | return((unsigned int) (quantum+0.5)); |
cristy | 0402ad3 | 2010-05-11 13:57:55 +0000 | [diff] [blame] | 529 | #endif |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 530 | } |
| 531 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 532 | static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum) |
| 533 | { |
| 534 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 535 | return((MagickSizeType) quantum); |
| 536 | #else |
| 537 | return((MagickSizeType) (quantum+0.5)); |
| 538 | #endif |
| 539 | } |
| 540 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 541 | static inline unsigned int ScaleQuantumToMap(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 542 | { |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 543 | if (quantum < 0.0) |
| 544 | return(0UL); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 545 | if ((quantum/65537) >= (Quantum) MaxMap) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 546 | return((unsigned int) MaxMap); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 547 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 548 | return((unsigned int) ((quantum+MagickULLConstant(32768))/ |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 549 | MagickULLConstant(65537))); |
| 550 | #else |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 551 | return((unsigned int) (quantum/65537.0+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 552 | #endif |
| 553 | } |
| 554 | |
| 555 | static inline unsigned short ScaleQuantumToShort(const Quantum quantum) |
| 556 | { |
| 557 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 558 | return((unsigned short) ((quantum+MagickULLConstant(32768))/ |
| 559 | MagickULLConstant(65537))); |
| 560 | #else |
| 561 | if (quantum <= 0.0) |
| 562 | return(0); |
cristy | cd817db | 2010-05-09 03:05:41 +0000 | [diff] [blame] | 563 | if ((quantum/65537.0) >= 65535.0) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 564 | return(65535); |
cristy | cd817db | 2010-05-09 03:05:41 +0000 | [diff] [blame] | 565 | return((unsigned short) (quantum/65537.0+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 566 | #endif |
| 567 | } |
| 568 | |
| 569 | static inline Quantum ScaleShortToQuantum(const unsigned short value) |
| 570 | { |
| 571 | #if !defined(MAGICKCORE_HDRI_SUPPORT) |
| 572 | return((Quantum) (65537UL*value)); |
| 573 | #else |
cristy | 41f82a6 | 2010-05-10 13:05:18 +0000 | [diff] [blame] | 574 | return((Quantum) (65537.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 575 | #endif |
| 576 | } |
| 577 | #elif (MAGICKCORE_QUANTUM_DEPTH == 64) |
| 578 | static inline Quantum ScaleCharToQuantum(const unsigned char value) |
| 579 | { |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 580 | return((Quantum) (72340172838076673.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 581 | } |
| 582 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 583 | static inline Quantum ScaleLongToQuantum(const unsigned int value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 584 | { |
cristy | ff2c025 | 2010-05-11 14:00:48 +0000 | [diff] [blame] | 585 | return((Quantum) (4294967297.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 586 | } |
| 587 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 588 | static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value) |
| 589 | { |
| 590 | return((Quantum) (18446744073709551615.0*value)); |
| 591 | } |
| 592 | |
| 593 | static inline Quantum ScaleMapToQuantum(const MagickRealType value) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 594 | { |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 595 | if (value <= 0.0) |
| 596 | return((Quantum) 0); |
cristy | a45da3b | 2010-03-31 02:06:45 +0000 | [diff] [blame] | 597 | if (value >= MaxMap) |
| 598 | return(QuantumRange); |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 599 | return((Quantum) (281479271743489.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 600 | } |
| 601 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 602 | static inline unsigned int ScaleQuantumToLong(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 603 | { |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 604 | return((unsigned int) (quantum/4294967297.0+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 605 | } |
| 606 | |
cristy | 8b8290b | 2014-05-26 12:55:32 +0000 | [diff] [blame] | 607 | static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum) |
| 608 | { |
| 609 | return((MagickSizeType) (quantum/18446744073709551615.0+0.5)); |
| 610 | } |
| 611 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 612 | static inline unsigned int ScaleQuantumToMap(const Quantum quantum) |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 613 | { |
cristy | 41f82a6 | 2010-05-10 13:05:18 +0000 | [diff] [blame] | 614 | if (quantum <= 0.0) |
cristy | dd2fd9f | 2010-05-10 01:43:37 +0000 | [diff] [blame] | 615 | return(0UL); |
cristy | 4351352 | 2010-05-11 13:46:36 +0000 | [diff] [blame] | 616 | if ((quantum/281479271743489.0) >= MaxMap) |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 617 | return((unsigned int) MaxMap); |
| 618 | return((unsigned int) (quantum/281479271743489.0+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 619 | } |
| 620 | |
| 621 | static inline unsigned short ScaleQuantumToShort(const Quantum quantum) |
| 622 | { |
cristy | 5f54f2b | 2010-05-11 14:29:17 +0000 | [diff] [blame] | 623 | if (quantum <= 0.0) |
| 624 | return(0); |
cristy | f4ad8d9 | 2010-05-12 17:58:41 +0000 | [diff] [blame] | 625 | if ((quantum/281479271743489.0) >= 65535.0) |
cristy | 5f54f2b | 2010-05-11 14:29:17 +0000 | [diff] [blame] | 626 | return(65535); |
cristy | 4351352 | 2010-05-11 13:46:36 +0000 | [diff] [blame] | 627 | return((unsigned short) (quantum/281479271743489.0+0.5)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 628 | } |
| 629 | |
| 630 | static inline Quantum ScaleShortToQuantum(const unsigned short value) |
| 631 | { |
cristy | b3963e6 | 2010-05-09 21:36:40 +0000 | [diff] [blame] | 632 | return((Quantum) (281479271743489.0*value)); |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 633 | } |
| 634 | #endif |
| 635 | |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 636 | static inline unsigned short SinglePrecisionToHalf(const float value) |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 637 | { |
| 638 | typedef union _SinglePrecision |
| 639 | { |
| 640 | unsigned int |
| 641 | fixed_point; |
| 642 | |
| 643 | float |
| 644 | single_precision; |
| 645 | } SinglePrecision; |
| 646 | |
| 647 | register int |
| 648 | exponent; |
| 649 | |
| 650 | register unsigned int |
| 651 | significand, |
| 652 | sign_bit; |
| 653 | |
| 654 | SinglePrecision |
| 655 | map; |
| 656 | |
| 657 | unsigned short |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 658 | half; |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 659 | |
| 660 | /* |
cristy | bccb64f | 2010-01-10 21:16:36 +0000 | [diff] [blame] | 661 | The IEEE 754 standard specifies half precision as having: |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 662 | |
| 663 | Sign bit: 1 bit |
| 664 | Exponent width: 5 bits |
| 665 | Significand precision: 11 (10 explicitly stored) |
| 666 | */ |
| 667 | map.single_precision=value; |
| 668 | sign_bit=(map.fixed_point >> 16) & 0x00008000; |
| 669 | exponent=(int) ((map.fixed_point >> ExponentShift) & 0x000000ff)-ExponentBias; |
| 670 | significand=map.fixed_point & 0x007fffff; |
| 671 | if (exponent <= 0) |
| 672 | { |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 673 | int |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 674 | shift; |
| 675 | |
| 676 | if (exponent < -10) |
| 677 | return((unsigned short) sign_bit); |
| 678 | significand=significand | 0x00800000; |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 679 | shift=(int) (14-exponent); |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 680 | significand=(unsigned int) ((significand+((1 << (shift-1))-1)+ |
| 681 | ((significand >> shift) & 0x01)) >> shift); |
| 682 | return((unsigned short) (sign_bit | significand)); |
| 683 | } |
| 684 | else |
| 685 | if (exponent == (0xff-ExponentBias)) |
| 686 | { |
| 687 | if (significand == 0) |
| 688 | return((unsigned short) (sign_bit | ExponentMask)); |
| 689 | else |
| 690 | { |
| 691 | significand>>=SignificandShift; |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 692 | half=(unsigned short) (sign_bit | significand | |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 693 | (significand == 0) | ExponentMask); |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 694 | return(half); |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 695 | } |
| 696 | } |
| 697 | significand=significand+((significand >> SignificandShift) & 0x01)+0x00000fff; |
| 698 | if ((significand & 0x00800000) != 0) |
| 699 | { |
| 700 | significand=0; |
| 701 | exponent++; |
| 702 | } |
| 703 | if (exponent > 30) |
| 704 | { |
| 705 | float |
| 706 | alpha; |
| 707 | |
cristy | 4cb162a | 2010-05-30 03:04:47 +0000 | [diff] [blame] | 708 | register int |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 709 | i; |
| 710 | |
| 711 | /* |
| 712 | Float overflow. |
| 713 | */ |
| 714 | alpha=1.0e10; |
| 715 | for (i=0; i < 10; i++) |
| 716 | alpha*=alpha; |
| 717 | return((unsigned short) (sign_bit | ExponentMask)); |
| 718 | } |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 719 | half=(unsigned short) (sign_bit | (exponent << 10) | |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 720 | (significand >> SignificandShift)); |
cristy | 2a4d01c | 2010-01-10 21:14:51 +0000 | [diff] [blame] | 721 | return(half); |
cristy | 24fb7dc | 2010-01-10 20:15:51 +0000 | [diff] [blame] | 722 | } |
| 723 | |
cristy | 3ed852e | 2009-09-05 21:47:34 +0000 | [diff] [blame] | 724 | #if defined(__cplusplus) || defined(c_plusplus) |
| 725 | } |
| 726 | #endif |
| 727 | |
| 728 | #endif |