sandy
/
PL286
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
1 Commit
1 Branch
0 Tags
18 GiB
Branch: master
Branches Tags
${ item.name }
Create tag ${ searchTerm }
Create branch ${ searchTerm }
from 'master'
${ noResults }
PL286/external/ImageMagick/MagickCore/quantum-private.h

728 lines
17 KiB
Raw Permalink Blame History

/*
Copyright 1999-2016 ImageMagick Studio LLC, a non-profit organization
dedicated to making software imaging solutions freely available.
You may not use this file except in compliance with the License.
obtain a copy of the License at
http://www.imagemagick.org/script/license.php
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
MagickCore quantum inline methods.
*/
#ifndef MAGICKCORE_QUANTUM_PRIVATE_H
#define MAGICKCORE_QUANTUM_PRIVATE_H
#include "MagickCore/cache.h"
#if defined(__cplusplus) || defined(c_plusplus)
extern "C" {
#endif
typedef struct _QuantumState
{
double
inverse_scale;
unsigned int
pixel;
size_t
bits;
const unsigned int
*mask;
} QuantumState;
struct _QuantumInfo
{
size_t
depth,
quantum;
QuantumFormatType
format;
double
minimum,
maximum,
scale;
size_t
pad;
MagickBooleanType
min_is_white,
pack;
QuantumAlphaType
alpha_type;
size_t
number_threads;
unsigned char
**pixels;
size_t
extent;
EndianType
endian;
QuantumState
state;
SemaphoreInfo
*semaphore;
size_t
signature;
};
extern MagickPrivate void
ResetQuantumState(QuantumInfo *);
static inline MagickSizeType GetQuantumRange(const size_t depth)
{
MagickSizeType
one;
one=1;
return((MagickSizeType) ((one << (depth-1))+((one << (depth-1))-1)));
}
static inline float HalfToSinglePrecision(const unsigned short half)
{
#define ExponentBias (127-15)
#define ExponentMask 0x7c00
#define ExponentShift 23
#define SignBitShift 31
#define SignificandShift 13
#define SignificandMask 0x00000400
typedef union _SinglePrecision
{
unsigned int
fixed_point;
float
single_precision;
} SinglePrecision;
register unsigned int
exponent,
significand,
sign_bit;
SinglePrecision
map;
unsigned int
value;
/*
The IEEE 754 standard specifies half precision as having:
Sign bit: 1 bit
Exponent width: 5 bits
Significand precision: 11 (10 explicitly stored)
*/
sign_bit=(unsigned int) ((half >> 15) & 0x00000001);
exponent=(unsigned int) ((half >> 10) & 0x0000001f);
significand=(unsigned int) (half & 0x000003ff);
if (exponent == 0)
{
if (significand == 0)
value=sign_bit << SignBitShift;
else
{
while ((significand & SignificandMask) == 0)
{
significand<<=1;
exponent--;
}
exponent++;
significand&=(~SignificandMask);
exponent+=ExponentBias;
value=(sign_bit << SignBitShift) | (exponent << ExponentShift) |
(significand << SignificandShift);
}
}
else
if (exponent == SignBitShift)
{
value=(sign_bit << SignBitShift) | 0x7f800000;
if (significand != 0)
value|=(significand << SignificandShift);
}
else
{
exponent+=ExponentBias;
significand<<=SignificandShift;
value=(sign_bit << SignBitShift) | (exponent << ExponentShift) |
significand;
}
map.fixed_point=value;
return(map.single_precision);
}
static inline unsigned char *PopCharPixel(const unsigned char pixel,
unsigned char *pixels)
{
*pixels++=pixel;
return(pixels);
}
static inline unsigned char *PopLongPixel(const EndianType endian,
const unsigned int pixel,unsigned char *pixels)
{
register unsigned int
quantum;
quantum=(unsigned int) pixel;
if (endian == LSBEndian)
{
*pixels++=(unsigned char) (quantum);
*pixels++=(unsigned char) (quantum >> 8);
*pixels++=(unsigned char) (quantum >> 16);
*pixels++=(unsigned char) (quantum >> 24);
return(pixels);
}
*pixels++=(unsigned char) (quantum >> 24);
*pixels++=(unsigned char) (quantum >> 16);
*pixels++=(unsigned char) (quantum >> 8);
*pixels++=(unsigned char) (quantum);
return(pixels);
}
static inline unsigned char *PopShortPixel(const EndianType endian,
const unsigned short pixel,unsigned char *pixels)
{
register unsigned int
quantum;
quantum=pixel;
if (endian == LSBEndian)
{
*pixels++=(unsigned char) (quantum);
*pixels++=(unsigned char) (quantum >> 8);
return(pixels);
}
*pixels++=(unsigned char) (quantum >> 8);
*pixels++=(unsigned char) (quantum);
return(pixels);
}
static inline const unsigned char *PushCharPixel(const unsigned char *pixels,
unsigned char *pixel)
{
*pixel=(*pixels++);
return(pixels);
}
static inline const unsigned char *PushLongPixel(const EndianType endian,
const unsigned char *pixels,unsigned int *pixel)
{
register unsigned int
quantum;
if (endian == LSBEndian)
{
quantum=((unsigned int) *pixels++);
quantum|=((unsigned int) *pixels++ << 8);
quantum|=((unsigned int) *pixels++ << 16);
quantum|=((unsigned int) *pixels++ << 24);
*pixel=quantum;
return(pixels);
}
quantum=((unsigned int) *pixels++ << 24);
quantum|=((unsigned int) *pixels++ << 16);
quantum|=((unsigned int) *pixels++ << 8);
quantum|=((unsigned int) *pixels++);
*pixel=quantum;
return(pixels);
}
static inline const unsigned char *PushShortPixel(const EndianType endian,
const unsigned char *pixels,unsigned short *pixel)
{
register unsigned int
quantum;
if (endian == LSBEndian)
{
quantum=(unsigned int) *pixels++;
quantum|=(unsigned int) (*pixels++ << 8);
*pixel=(unsigned short) (quantum & 0xffff);
return(pixels);
}
quantum=(unsigned int) (*pixels++ << 8);
quantum|=(unsigned int) *pixels++;
*pixel=(unsigned short) (quantum & 0xffff);
return(pixels);
}
static inline Quantum ScaleAnyToQuantum(const QuantumAny quantum,
const QuantumAny range)
{
if (quantum > range)
return(QuantumRange);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (((double) QuantumRange*quantum)/range+0.5));
#else
return((Quantum) (((double) QuantumRange*quantum)/range));
#endif
}
static inline QuantumAny ScaleQuantumToAny(const Quantum quantum,
const QuantumAny range)
{
return((QuantumAny) (((double) range*quantum)/QuantumRange+0.5));
}
#if (MAGICKCORE_QUANTUM_DEPTH == 8)
static inline Quantum ScaleCharToQuantum(const unsigned char value)
{
return((Quantum) value);
}
static inline Quantum ScaleLongToQuantum(const unsigned int value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) ((value+8421504UL)/16843009UL));
#else
return((Quantum) (value/16843009.0));
#endif
}
static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) ((value+MagickULLConstant(551911719039))/
MagickULLConstant(1103823438079)));
#else
return((Quantum) (value/1103823438079.0));
#endif
}
static inline Quantum ScaleMapToQuantum(const MagickRealType value)
{
if (value <= 0.0)
return((Quantum) 0);
if (value >= MaxMap)
return(QuantumRange);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (value+0.5));
#else
return((Quantum) value);
#endif
}
static inline unsigned int ScaleQuantumToLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) (16843009UL*quantum));
#else
if (quantum <= 0.0)
return(0UL);
if ((16843009.0*quantum) >= 4294967295.0)
return(4294967295UL);
return((unsigned int) (16843009.0*quantum+0.5));
#endif
}
static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((MagickSizeType) (MagickULLConstant(551911719039)*quantum));
#else
if (quantum <= 0.0)
return(0UL);
if ((551911719039.0*quantum) >= 18446744073709551615.0)
return(MagickULLConstant(18446744073709551615));
return((MagickSizeType) (1103823438079.0*quantum+0.5));
#endif
}
static inline unsigned int ScaleQuantumToMap(const Quantum quantum)
{
if (quantum >= (Quantum) MaxMap)
return((unsigned int) MaxMap);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) quantum);
#else
if (quantum < 0.0)
return(0UL);
return((unsigned int) (quantum+0.5));
#endif
}
static inline unsigned short ScaleQuantumToShort(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned short) (257UL*quantum));
#else
if (quantum <= 0.0)
return(0);
if ((257.0*quantum) >= 65535.0)
return(65535);
return((unsigned short) (257.0*quantum+0.5));
#endif
}
static inline Quantum ScaleShortToQuantum(const unsigned short value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) ((value+128U)/257U));
#else
return((Quantum) (value/257.0));
#endif
}
#elif (MAGICKCORE_QUANTUM_DEPTH == 16)
static inline Quantum ScaleCharToQuantum(const unsigned char value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (257U*value));
#else
return((Quantum) (257.0*value));
#endif
}
static inline Quantum ScaleLongToQuantum(const unsigned int value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) ((value+MagickULLConstant(32768))/
MagickULLConstant(65537)));
#else
return((Quantum) (value/65537.0));
#endif
}
static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) ((value+MagickULLConstant(8421376))/
MagickULLConstant(16842752)));
#else
return((Quantum) (value/16842752.0));
#endif
}
static inline Quantum ScaleMapToQuantum(const MagickRealType value)
{
if (value <= 0.0)
return((Quantum) 0);
if (value >= MaxMap)
return(QuantumRange);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (value+0.5));
#else
return((Quantum) value);
#endif
}
static inline unsigned int ScaleQuantumToLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) (65537UL*quantum));
#else
if (quantum <= 0.0)
return(0UL);
if ((65537.0*quantum) >= 4294967295.0)
return(4294967295U);
return((unsigned int) (65537.0*quantum+0.5));
#endif
}
static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((MagickSizeType) (MagickULLConstant(16842752)*quantum));
#else
if (quantum <= 0.0)
return(0UL);
if ((65537.0*quantum) >= 18446744073709551615.0)
return(MagickULLConstant(18446744073709551615));
return((MagickSizeType) (16842752.0*quantum+0.5));
#endif
}
static inline unsigned int ScaleQuantumToMap(const Quantum quantum)
{
if (quantum >= (Quantum) MaxMap)
return((unsigned int) MaxMap);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) quantum);
#else
if (quantum < 0.0)
return(0UL);
return((unsigned int) (quantum+0.5));
#endif
}
static inline unsigned short ScaleQuantumToShort(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned short) quantum);
#else
if (quantum <= 0.0)
return(0);
if (quantum >= 65535.0)
return(65535);
return((unsigned short) (quantum+0.5));
#endif
}
static inline Quantum ScaleShortToQuantum(const unsigned short value)
{
return((Quantum) value);
}
#elif (MAGICKCORE_QUANTUM_DEPTH == 32)
static inline Quantum ScaleCharToQuantum(const unsigned char value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (16843009UL*value));
#else
return((Quantum) (16843009.0*value));
#endif
}
static inline Quantum ScaleLongToQuantum(const unsigned int value)
{
return((Quantum) value);
}
static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value)
{
return((Quantum) value);
}
static inline Quantum ScaleMapToQuantum(const MagickRealType value)
{
if (value <= 0.0)
return((Quantum) 0);
if (value >= (Quantum) MaxMap)
return(QuantumRange);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (65537.0*value+0.5));
#else
return((Quantum) (65537.0*value));
#endif
}
static inline unsigned int ScaleQuantumToLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) quantum);
#else
if (quantum <= 0.0)
return(0);
if ((quantum) >= 4294967295.0)
return(4294967295);
return((unsigned int) (quantum+0.5));
#endif
}
static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((MagickSizeType) quantum);
#else
return((MagickSizeType) (quantum+0.5));
#endif
}
static inline unsigned int ScaleQuantumToMap(const Quantum quantum)
{
if (quantum < 0.0)
return(0UL);
if ((quantum/65537) >= (Quantum) MaxMap)
return((unsigned int) MaxMap);
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned int) ((quantum+MagickULLConstant(32768))/
MagickULLConstant(65537)));
#else
return((unsigned int) (quantum/65537.0+0.5));
#endif
}
static inline unsigned short ScaleQuantumToShort(const Quantum quantum)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((unsigned short) ((quantum+MagickULLConstant(32768))/
MagickULLConstant(65537)));
#else
if (quantum <= 0.0)
return(0);
if ((quantum/65537.0) >= 65535.0)
return(65535);
return((unsigned short) (quantum/65537.0+0.5));
#endif
}
static inline Quantum ScaleShortToQuantum(const unsigned short value)
{
#if !defined(MAGICKCORE_HDRI_SUPPORT)
return((Quantum) (65537UL*value));
#else
return((Quantum) (65537.0*value));
#endif
}
#elif (MAGICKCORE_QUANTUM_DEPTH == 64)
static inline Quantum ScaleCharToQuantum(const unsigned char value)
{
return((Quantum) (72340172838076673.0*value));
}
static inline Quantum ScaleLongToQuantum(const unsigned int value)
{
return((Quantum) (4294967297.0*value));
}
static inline Quantum ScaleLongLongToQuantum(const MagickSizeType value)
{
return((Quantum) (18446744073709551615.0*value));
}
static inline Quantum ScaleMapToQuantum(const MagickRealType value)
{
if (value <= 0.0)
return((Quantum) 0);
if (value >= MaxMap)
return(QuantumRange);
return((Quantum) (281479271743489.0*value));
}
static inline unsigned int ScaleQuantumToLong(const Quantum quantum)
{
return((unsigned int) (quantum/4294967297.0+0.5));
}
static inline MagickSizeType ScaleQuantumToLongLong(const Quantum quantum)
{
return((MagickSizeType) (quantum/18446744073709551615.0+0.5));
}
static inline unsigned int ScaleQuantumToMap(const Quantum quantum)
{
if (quantum <= 0.0)
return(0UL);
if ((quantum/281479271743489.0) >= MaxMap)
return((unsigned int) MaxMap);
return((unsigned int) (quantum/281479271743489.0+0.5));
}
static inline unsigned short ScaleQuantumToShort(const Quantum quantum)
{
if (quantum <= 0.0)
return(0);
if ((quantum/281479271743489.0) >= 65535.0)
return(65535);
return((unsigned short) (quantum/281479271743489.0+0.5));
}
static inline Quantum ScaleShortToQuantum(const unsigned short value)
{
return((Quantum) (281479271743489.0*value));
}
#endif
static inline unsigned short SinglePrecisionToHalf(const float value)
{
typedef union _SinglePrecision
{
unsigned int
fixed_point;
float
single_precision;
} SinglePrecision;
register int
exponent;
register unsigned int
significand,
sign_bit;
SinglePrecision
map;
unsigned short
half;
/*
The IEEE 754 standard specifies half precision as having:
Sign bit: 1 bit
Exponent width: 5 bits
Significand precision: 11 (10 explicitly stored)
*/
map.single_precision=value;
sign_bit=(map.fixed_point >> 16) & 0x00008000;
exponent=(int) ((map.fixed_point >> ExponentShift) & 0x000000ff)-ExponentBias;
significand=map.fixed_point & 0x007fffff;
if (exponent <= 0)
{
int
shift;
if (exponent < -10)
return((unsigned short) sign_bit);
significand=significand | 0x00800000;
shift=(int) (14-exponent);
significand=(unsigned int) ((significand+((1 << (shift-1))-1)+
((significand >> shift) & 0x01)) >> shift);
return((unsigned short) (sign_bit | significand));
}
else
if (exponent == (0xff-ExponentBias))
{
if (significand == 0)
return((unsigned short) (sign_bit | ExponentMask));
else
{
significand>>=SignificandShift;
half=(unsigned short) (sign_bit | significand |
(significand == 0) | ExponentMask);
return(half);
}
}
significand=significand+((significand >> SignificandShift) & 0x01)+0x00000fff;
if ((significand & 0x00800000) != 0)
{
significand=0;
exponent++;
}
if (exponent > 30)
{
float
alpha;
register int
i;
/*
Float overflow.
*/
alpha=1.0e10;
for (i=0; i < 10; i++)
alpha*=alpha;
return((unsigned short) (sign_bit | ExponentMask));
}
half=(unsigned short) (sign_bit | (exponent << 10) |
(significand >> SignificandShift));
return(half);
}
#if defined(__cplusplus) || defined(c_plusplus)
}
#endif
#endif