include/tinyhtm/common.h

Go to the documentation of this file.

#ifndef HTM_COMMON_H
#define HTM_COMMON_H

#include <stddef.h>
#include <stdint.h>
#include <math.h>

#include "config.h"


#if HAVE_ATTRIBUTE_UNUSED
#   define HTM_UNUSED __attribute__ ((unused))
#else
#   define HTM_UNUSED
#endif

#if __STDC_VERSION__ >= 199901L || __GNUC__
/* "inline" is a keyword */
#   define HTM_INLINE static inline
#else
#   define HTM_INLINE static HTM_UNUSED
#endif

#if HAVE_ATTRIBUTE_ALIGNED
#   define HTM_ALIGNED(x) __attribute__ ((aligned(x)))
#else
#   define HTM_ALIGNED(x)
#endif

#if __STDC_VERSION__ >= 199901L
#   define HTM_ISNAN(x) isnan(x)
#   define HTM_ISSPECIAL(x) (!isfinite(x))
#else
#   define HTM_ISNAN(x) ((x) != (x))
#   define HTM_ISSPECIAL(x) ((x) != (x) || ((x) != 0.0 && (x) == 2*(x)))
#endif

/* M_PI may not be available in <math.h> */
#ifndef M_PI
#   define M_PI 3.141592653589793238462643
#endif


#ifdef __cplusplus
extern "C" {
#endif

struct htm_range {
    int64_t min;  
    int64_t max;  
};

HTM_INLINE int htm_popcount(uint64_t x)
{
    const uint64_t c1 = UINT64_C(0x5555555555555555); /* -1/3 */
    const uint64_t c2 = UINT64_C(0x3333333333333333); /* -1/5 */
    const uint64_t c4 = UINT64_C(0x0f0f0f0f0f0f0f0f); /* -1/17 */
    const uint64_t f  = UINT64_C(0x0101010101010101); /* -1/255 */
    x = x - ((x >> 1)  & c1);
    x = (x & c2) + ((x >> 2)  & c2);
    x = (x + (x >> 4)) & c4;
    return (int) ((x * f) >> 56);
}

enum htm_errcode {
    HTM_OK = 0,     
    HTM_ENOMEM,     
    HTM_ENULLPTR,   
    HTM_ENANINF,    
    HTM_EZERONORM,  
    HTM_ELAT,       
    HTM_EANG,       
    HTM_EHEMIS,     
    HTM_ELEN,       
    HTM_EDEGEN,     
    HTM_EID,        
    HTM_ELEVEL,     
    HTM_EIO,        
    HTM_EMMAN,      
    HTM_EINV,       
    HTM_ETREE,      
    HTM_NUM_CODES
};


const char * htm_errmsg(enum htm_errcode err);


#define HTM_DEG_PER_RAD    57.2957795130823208767981548141

#define HTM_RAD_PER_DEG    0.0174532925199432957692369076849

#define HTM_ARCSEC_PER_DEG 3600.0


HTM_INLINE double htm_angred(double angle_deg)
{
    double angle = fmod(angle_deg, 360.0);
    if (angle < 0.0) {
        angle += 360.0;
        if (angle == 360.0) {
            angle = 0.0;
        }
    }
    return angle;
}

HTM_INLINE double htm_clamp(double x, double min, double max)
{
    return (x < min) ? min : ((x > max) ? max : x);
}

#ifdef __cplusplus
}
#endif

#endif /* HTM_COMMON_H */