Changeset 7cfef0d for libcfa/src/math.hfa
- Timestamp:
- Aug 24, 2020, 8:34:24 PM (4 years ago)
- Branches:
- ADT, arm-eh, ast-experimental, enum, forall-pointer-decay, jacob/cs343-translation, master, new-ast-unique-expr, pthread-emulation, qualifiedEnum
- Children:
- 0254454
- Parents:
- 97cba9f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
libcfa/src/math.hfa
r97cba9f r7cfef0d 10 10 // Created On : Mon Apr 18 23:37:04 2016 11 11 // Last Modified By : Peter A. Buhr 12 // Last Modified On : Tue Feb 4 10:27:11202013 // Update Count : 1 1712 // Last Modified On : Mon Aug 24 08:56:20 2020 13 // Update Count : 126 14 14 // 15 15 … … 19 19 #include <complex.h> 20 20 21 //--------------------------------------- 22 23 #include "common.hfa" 24 21 25 //---------------------- General ---------------------- 22 26 23 static inline float ?%?( float x, float y ) { return fmodf( x, y ); } 24 static inline float fmod( float x, float y ) { return fmodf( x, y ); } 25 static inline double ?%?( double x, double y ) { return fmod( x, y ); } 26 // extern "C" { double fmod( double, double ); } 27 static inline long double ?%?( long double x, long double y ) { return fmodl( x, y ); } 28 static inline long double fmod( long double x, long double y ) { return fmodl( x, y ); } 29 30 static inline float remainder( float x, float y ) { return remainderf( x, y ); } 31 // extern "C" { double remainder( double, double ); } 32 static inline long double remainder( long double x, long double y ) { return remainderl( x, y ); } 33 34 static inline float remquo( float x, float y, int * quo ) { return remquof( x, y, quo ); } 35 // extern "C" { double remquo( double x, double y, int * quo ); } 36 static inline long double remquo( long double x, long double y, int * quo ) { return remquol( x, y, quo ); } 37 static inline [ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; } 38 static inline [ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; } 39 static inline [ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; } 40 41 static inline [ float, float ] div( float x, float y ) { y = modff( x / y, &x ); return [ x, y ]; } 42 static inline [ double, double ] div( double x, double y ) { y = modf( x / y, &x ); return [ x, y ]; } 43 static inline [ long double, long double ] div( long double x, long double y ) { y = modfl( x / y, &x ); return [ x, y ]; } 44 45 static inline float fma( float x, float y, float z ) { return fmaf( x, y, z ); } 46 // extern "C" { double fma( double, double, double ); } 47 static inline long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); } 48 49 static inline float fdim( float x, float y ) { return fdimf( x, y ); } 50 // extern "C" { double fdim( double, double ); } 51 static inline long double fdim( long double x, long double y ) { return fdiml( x, y ); } 52 53 static inline float nan( const char tag[] ) { return nanf( tag ); } 54 // extern "C" { double nan( const char [] ); } 55 static inline long double nan( const char tag[] ) { return nanl( tag ); } 27 static inline { 28 float ?%?( float x, float y ) { return fmodf( x, y ); } 29 float fmod( float x, float y ) { return fmodf( x, y ); } 30 double ?%?( double x, double y ) { return fmod( x, y ); } 31 // extern "C" { double fmod( double, double ); } 32 long double ?%?( long double x, long double y ) { return fmodl( x, y ); } 33 long double fmod( long double x, long double y ) { return fmodl( x, y ); } 34 35 float remainder( float x, float y ) { return remainderf( x, y ); } 36 // extern "C" { double remainder( double, double ); } 37 long double remainder( long double x, long double y ) { return remainderl( x, y ); } 38 39 float remquo( float x, float y, int * quo ) { return remquof( x, y, quo ); } 40 // extern "C" { double remquo( double x, double y, int * quo ); } 41 long double remquo( long double x, long double y, int * quo ) { return remquol( x, y, quo ); } 42 [ int, float ] remquo( float x, float y ) { int quo; x = remquof( x, y, &quo ); return [ quo, x ]; } 43 [ int, double ] remquo( double x, double y ) { int quo; x = remquo( x, y, &quo ); return [ quo, x ]; } 44 [ int, long double ] remquo( long double x, long double y ) { int quo; x = remquol( x, y, &quo ); return [ quo, x ]; } 45 46 [ float, float ] div( float x, float y ) { y = modff( x / y, &x ); return [ x, y ]; } 47 [ double, double ] div( double x, double y ) { y = modf( x / y, &x ); return [ x, y ]; } 48 [ long double, long double ] div( long double x, long double y ) { y = modfl( x / y, &x ); return [ x, y ]; } 49 50 float fma( float x, float y, float z ) { return fmaf( x, y, z ); } 51 // extern "C" { double fma( double, double, double ); } 52 long double fma( long double x, long double y, long double z ) { return fmal( x, y, z ); } 53 54 float fdim( float x, float y ) { return fdimf( x, y ); } 55 // extern "C" { double fdim( double, double ); } 56 long double fdim( long double x, long double y ) { return fdiml( x, y ); } 57 58 float nan( const char tag[] ) { return nanf( tag ); } 59 // extern "C" { double nan( const char [] ); } 60 long double nan( const char tag[] ) { return nanl( tag ); } 61 } // distribution 56 62 57 63 //---------------------- Exponential ---------------------- 58 64 59 static inline float exp( float x ) { return expf( x ); } 60 // extern "C" { double exp( double ); } 61 static inline long double exp( long double x ) { return expl( x ); } 62 static inline float _Complex exp( float _Complex x ) { return cexpf( x ); } 63 static inline double _Complex exp( double _Complex x ) { return cexp( x ); } 64 static inline long double _Complex exp( long double _Complex x ) { return cexpl( x ); } 65 66 static inline float exp2( float x ) { return exp2f( x ); } 67 // extern "C" { double exp2( double ); } 68 static inline long double exp2( long double x ) { return exp2l( x ); } 69 //static inline float _Complex exp2( float _Complex x ) { return cexp2f( x ); } 70 //static inline double _Complex exp2( double _Complex x ) { return cexp2( x ); } 71 //static inline long double _Complex exp2( long double _Complex x ) { return cexp2l( x ); } 72 73 static inline float expm1( float x ) { return expm1f( x ); } 74 // extern "C" { double expm1( double ); } 75 static inline long double expm1( long double x ) { return expm1l( x ); } 76 77 static inline float pow( float x, float y ) { return powf( x, y ); } 78 // extern "C" { double pow( double, double ); } 79 static inline long double pow( long double x, long double y ) { return powl( x, y ); } 80 static inline float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); } 81 static inline double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); } 82 static inline long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); } 65 static inline { 66 float exp( float x ) { return expf( x ); } 67 // extern "C" { double exp( double ); } 68 long double exp( long double x ) { return expl( x ); } 69 float _Complex exp( float _Complex x ) { return cexpf( x ); } 70 double _Complex exp( double _Complex x ) { return cexp( x ); } 71 long double _Complex exp( long double _Complex x ) { return cexpl( x ); } 72 73 float exp2( float x ) { return exp2f( x ); } 74 // extern "C" { double exp2( double ); } 75 long double exp2( long double x ) { return exp2l( x ); } 76 //float _Complex exp2( float _Complex x ) { return cexp2f( x ); } 77 //double _Complex exp2( double _Complex x ) { return cexp2( x ); } 78 //long double _Complex exp2( long double _Complex x ) { return cexp2l( x ); } 79 80 float expm1( float x ) { return expm1f( x ); } 81 // extern "C" { double expm1( double ); } 82 long double expm1( long double x ) { return expm1l( x ); } 83 84 float pow( float x, float y ) { return powf( x, y ); } 85 // extern "C" { double pow( double, double ); } 86 long double pow( long double x, long double y ) { return powl( x, y ); } 87 float _Complex pow( float _Complex x, float _Complex y ) { return cpowf( x, y ); } 88 double _Complex pow( double _Complex x, double _Complex y ) { return cpow( x, y ); } 89 long double _Complex pow( long double _Complex x, long double _Complex y ) { return cpowl( x, y ); } 90 } // distribution 83 91 84 92 //---------------------- Logarithm ---------------------- 85 93 86 static inline float log( float x ) { return logf( x ); } 87 // extern "C" { double log( double ); } 88 static inline long double log( long double x ) { return logl( x ); } 89 static inline float _Complex log( float _Complex x ) { return clogf( x ); } 90 static inline double _Complex log( double _Complex x ) { return clog( x ); } 91 static inline long double _Complex log( long double _Complex x ) { return clogl( x ); } 92 93 static inline float log2( float x ) { return log2f( x ); } 94 // extern "C" { double log2( double ); } 95 static inline long double log2( long double x ) { return log2l( x ); } 96 // static inline float _Complex log2( float _Complex x ) { return clog2f( x ); } 97 // static inline double _Complex log2( double _Complex x ) { return clog2( x ); } 98 // static inline long double _Complex log2( long double _Complex x ) { return clog2l( x ); } 99 100 static inline float log10( float x ) { return log10f( x ); } 101 // extern "C" { double log10( double ); } 102 static inline long double log10( long double x ) { return log10l( x ); } 103 // static inline float _Complex log10( float _Complex x ) { return clog10f( x ); } 104 // static inline double _Complex log10( double _Complex x ) { return clog10( x ); } 105 // static inline long double _Complex log10( long double _Complex x ) { return clog10l( x ); } 106 107 static inline float log1p( float x ) { return log1pf( x ); } 108 // extern "C" { double log1p( double ); } 109 static inline long double log1p( long double x ) { return log1pl( x ); } 110 111 static inline int ilogb( float x ) { return ilogbf( x ); } 112 // extern "C" { int ilogb( double ); } 113 static inline int ilogb( long double x ) { return ilogbl( x ); } 114 115 static inline float logb( float x ) { return logbf( x ); } 116 // extern "C" { double logb( double ); } 117 static inline long double logb( long double x ) { return logbl( x ); } 118 119 static inline float sqrt( float x ) { return sqrtf( x ); } 120 // extern "C" { double sqrt( double ); } 121 static inline long double sqrt( long double x ) { return sqrtl( x ); } 122 static inline float _Complex sqrt( float _Complex x ) { return csqrtf( x ); } 123 static inline double _Complex sqrt( double _Complex x ) { return csqrt( x ); } 124 static inline long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); } 125 126 static inline float cbrt( float x ) { return cbrtf( x ); } 127 // extern "C" { double cbrt( double ); } 128 static inline long double cbrt( long double x ) { return cbrtl( x ); } 129 130 static inline float hypot( float x, float y ) { return hypotf( x, y ); } 131 // extern "C" { double hypot( double, double ); } 132 static inline long double hypot( long double x, long double y ) { return hypotl( x, y ); } 94 static inline { 95 float log( float x ) { return logf( x ); } 96 // extern "C" { double log( double ); } 97 long double log( long double x ) { return logl( x ); } 98 float _Complex log( float _Complex x ) { return clogf( x ); } 99 double _Complex log( double _Complex x ) { return clog( x ); } 100 long double _Complex log( long double _Complex x ) { return clogl( x ); } 101 102 float log2( float x ) { return log2f( x ); } 103 // extern "C" { double log2( double ); } 104 long double log2( long double x ) { return log2l( x ); } 105 // float _Complex log2( float _Complex x ) { return clog2f( x ); } 106 // double _Complex log2( double _Complex x ) { return clog2( x ); } 107 // long double _Complex log2( long double _Complex x ) { return clog2l( x ); } 108 109 float log10( float x ) { return log10f( x ); } 110 // extern "C" { double log10( double ); } 111 long double log10( long double x ) { return log10l( x ); } 112 // float _Complex log10( float _Complex x ) { return clog10f( x ); } 113 // double _Complex log10( double _Complex x ) { return clog10( x ); } 114 // long double _Complex log10( long double _Complex x ) { return clog10l( x ); } 115 116 float log1p( float x ) { return log1pf( x ); } 117 // extern "C" { double log1p( double ); } 118 long double log1p( long double x ) { return log1pl( x ); } 119 120 int ilogb( float x ) { return ilogbf( x ); } 121 // extern "C" { int ilogb( double ); } 122 int ilogb( long double x ) { return ilogbl( x ); } 123 124 float logb( float x ) { return logbf( x ); } 125 // extern "C" { double logb( double ); } 126 long double logb( long double x ) { return logbl( x ); } 127 128 float sqrt( float x ) { return sqrtf( x ); } 129 // extern "C" { double sqrt( double ); } 130 long double sqrt( long double x ) { return sqrtl( x ); } 131 float _Complex sqrt( float _Complex x ) { return csqrtf( x ); } 132 double _Complex sqrt( double _Complex x ) { return csqrt( x ); } 133 long double _Complex sqrt( long double _Complex x ) { return csqrtl( x ); } 134 135 float cbrt( float x ) { return cbrtf( x ); } 136 // extern "C" { double cbrt( double ); } 137 long double cbrt( long double x ) { return cbrtl( x ); } 138 139 float hypot( float x, float y ) { return hypotf( x, y ); } 140 // extern "C" { double hypot( double, double ); } 141 long double hypot( long double x, long double y ) { return hypotl( x, y ); } 142 } // distribution 133 143 134 144 //---------------------- Trigonometric ---------------------- 135 145 136 static inline float sin( float x ) { return sinf( x ); } 137 // extern "C" { double sin( double ); } 138 static inline long double sin( long double x ) { return sinl( x ); } 139 static inline float _Complex sin( float _Complex x ) { return csinf( x ); } 140 static inline double _Complex sin( double _Complex x ) { return csin( x ); } 141 static inline long double _Complex sin( long double _Complex x ) { return csinl( x ); } 142 143 static inline float cos( float x ) { return cosf( x ); } 144 // extern "C" { double cos( double ); } 145 static inline long double cos( long double x ) { return cosl( x ); } 146 static inline float _Complex cos( float _Complex x ) { return ccosf( x ); } 147 static inline double _Complex cos( double _Complex x ) { return ccos( x ); } 148 static inline long double _Complex cos( long double _Complex x ) { return ccosl( x ); } 149 150 static inline float tan( float x ) { return tanf( x ); } 151 // extern "C" { double tan( double ); } 152 static inline long double tan( long double x ) { return tanl( x ); } 153 static inline float _Complex tan( float _Complex x ) { return ctanf( x ); } 154 static inline double _Complex tan( double _Complex x ) { return ctan( x ); } 155 static inline long double _Complex tan( long double _Complex x ) { return ctanl( x ); } 156 157 static inline float asin( float x ) { return asinf( x ); } 158 // extern "C" { double asin( double ); } 159 static inline long double asin( long double x ) { return asinl( x ); } 160 static inline float _Complex asin( float _Complex x ) { return casinf( x ); } 161 static inline double _Complex asin( double _Complex x ) { return casin( x ); } 162 static inline long double _Complex asin( long double _Complex x ) { return casinl( x ); } 163 164 static inline float acos( float x ) { return acosf( x ); } 165 // extern "C" { double acos( double ); } 166 static inline long double acos( long double x ) { return acosl( x ); } 167 static inline float _Complex acos( float _Complex x ) { return cacosf( x ); } 168 static inline double _Complex acos( double _Complex x ) { return cacos( x ); } 169 static inline long double _Complex acos( long double _Complex x ) { return cacosl( x ); } 170 171 static inline float atan( float x ) { return atanf( x ); } 172 // extern "C" { double atan( double ); } 173 static inline long double atan( long double x ) { return atanl( x ); } 174 static inline float _Complex atan( float _Complex x ) { return catanf( x ); } 175 static inline double _Complex atan( double _Complex x ) { return catan( x ); } 176 static inline long double _Complex atan( long double _Complex x ) { return catanl( x ); } 177 178 static inline float atan2( float x, float y ) { return atan2f( x, y ); } 179 // extern "C" { double atan2( double, double ); } 180 static inline long double atan2( long double x, long double y ) { return atan2l( x, y ); } 181 182 // alternative name for atan2 183 static inline float atan( float x, float y ) { return atan2f( x, y ); } 184 static inline double atan( double x, double y ) { return atan2( x, y ); } 185 static inline long double atan( long double x, long double y ) { return atan2l( x, y ); } 146 static inline { 147 float sin( float x ) { return sinf( x ); } 148 // extern "C" { double sin( double ); } 149 long double sin( long double x ) { return sinl( x ); } 150 float _Complex sin( float _Complex x ) { return csinf( x ); } 151 double _Complex sin( double _Complex x ) { return csin( x ); } 152 long double _Complex sin( long double _Complex x ) { return csinl( x ); } 153 154 float cos( float x ) { return cosf( x ); } 155 // extern "C" { double cos( double ); } 156 long double cos( long double x ) { return cosl( x ); } 157 float _Complex cos( float _Complex x ) { return ccosf( x ); } 158 double _Complex cos( double _Complex x ) { return ccos( x ); } 159 long double _Complex cos( long double _Complex x ) { return ccosl( x ); } 160 161 float tan( float x ) { return tanf( x ); } 162 // extern "C" { double tan( double ); } 163 long double tan( long double x ) { return tanl( x ); } 164 float _Complex tan( float _Complex x ) { return ctanf( x ); } 165 double _Complex tan( double _Complex x ) { return ctan( x ); } 166 long double _Complex tan( long double _Complex x ) { return ctanl( x ); } 167 168 float asin( float x ) { return asinf( x ); } 169 // extern "C" { double asin( double ); } 170 long double asin( long double x ) { return asinl( x ); } 171 float _Complex asin( float _Complex x ) { return casinf( x ); } 172 double _Complex asin( double _Complex x ) { return casin( x ); } 173 long double _Complex asin( long double _Complex x ) { return casinl( x ); } 174 175 float acos( float x ) { return acosf( x ); } 176 // extern "C" { double acos( double ); } 177 long double acos( long double x ) { return acosl( x ); } 178 float _Complex acos( float _Complex x ) { return cacosf( x ); } 179 double _Complex acos( double _Complex x ) { return cacos( x ); } 180 long double _Complex acos( long double _Complex x ) { return cacosl( x ); } 181 182 float atan( float x ) { return atanf( x ); } 183 // extern "C" { double atan( double ); } 184 long double atan( long double x ) { return atanl( x ); } 185 float _Complex atan( float _Complex x ) { return catanf( x ); } 186 double _Complex atan( double _Complex x ) { return catan( x ); } 187 long double _Complex atan( long double _Complex x ) { return catanl( x ); } 188 189 float atan2( float x, float y ) { return atan2f( x, y ); } 190 // extern "C" { double atan2( double, double ); } 191 long double atan2( long double x, long double y ) { return atan2l( x, y ); } 192 193 // alternative name for atan2 194 float atan( float x, float y ) { return atan2f( x, y ); } 195 double atan( double x, double y ) { return atan2( x, y ); } 196 long double atan( long double x, long double y ) { return atan2l( x, y ); } 197 } // distribution 186 198 187 199 //---------------------- Hyperbolic ---------------------- 188 200 189 static inline float sinh( float x ) { return sinhf( x ); } 190 // extern "C" { double sinh( double ); } 191 static inline long double sinh( long double x ) { return sinhl( x ); } 192 static inline float _Complex sinh( float _Complex x ) { return csinhf( x ); } 193 static inline double _Complex sinh( double _Complex x ) { return csinh( x ); } 194 static inline long double _Complex sinh( long double _Complex x ) { return csinhl( x ); } 195 196 static inline float cosh( float x ) { return coshf( x ); } 197 // extern "C" { double cosh( double ); } 198 static inline long double cosh( long double x ) { return coshl( x ); } 199 static inline float _Complex cosh( float _Complex x ) { return ccoshf( x ); } 200 static inline double _Complex cosh( double _Complex x ) { return ccosh( x ); } 201 static inline long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); } 202 203 static inline float tanh( float x ) { return tanhf( x ); } 204 // extern "C" { double tanh( double ); } 205 static inline long double tanh( long double x ) { return tanhl( x ); } 206 static inline float _Complex tanh( float _Complex x ) { return ctanhf( x ); } 207 static inline double _Complex tanh( double _Complex x ) { return ctanh( x ); } 208 static inline long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); } 209 210 static inline float asinh( float x ) { return asinhf( x ); } 211 // extern "C" { double asinh( double ); } 212 static inline long double asinh( long double x ) { return asinhl( x ); } 213 static inline float _Complex asinh( float _Complex x ) { return casinhf( x ); } 214 static inline double _Complex asinh( double _Complex x ) { return casinh( x ); } 215 static inline long double _Complex asinh( long double _Complex x ) { return casinhl( x ); } 216 217 static inline float acosh( float x ) { return acoshf( x ); } 218 // extern "C" { double acosh( double ); } 219 static inline long double acosh( long double x ) { return acoshl( x ); } 220 static inline float _Complex acosh( float _Complex x ) { return cacoshf( x ); } 221 static inline double _Complex acosh( double _Complex x ) { return cacosh( x ); } 222 static inline long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); } 223 224 static inline float atanh( float x ) { return atanhf( x ); } 225 // extern "C" { double atanh( double ); } 226 static inline long double atanh( long double x ) { return atanhl( x ); } 227 static inline float _Complex atanh( float _Complex x ) { return catanhf( x ); } 228 static inline double _Complex atanh( double _Complex x ) { return catanh( x ); } 229 static inline long double _Complex atanh( long double _Complex x ) { return catanhl( x ); } 201 static inline { 202 float sinh( float x ) { return sinhf( x ); } 203 // extern "C" { double sinh( double ); } 204 long double sinh( long double x ) { return sinhl( x ); } 205 float _Complex sinh( float _Complex x ) { return csinhf( x ); } 206 double _Complex sinh( double _Complex x ) { return csinh( x ); } 207 long double _Complex sinh( long double _Complex x ) { return csinhl( x ); } 208 209 float cosh( float x ) { return coshf( x ); } 210 // extern "C" { double cosh( double ); } 211 long double cosh( long double x ) { return coshl( x ); } 212 float _Complex cosh( float _Complex x ) { return ccoshf( x ); } 213 double _Complex cosh( double _Complex x ) { return ccosh( x ); } 214 long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); } 215 216 float tanh( float x ) { return tanhf( x ); } 217 // extern "C" { double tanh( double ); } 218 long double tanh( long double x ) { return tanhl( x ); } 219 float _Complex tanh( float _Complex x ) { return ctanhf( x ); } 220 double _Complex tanh( double _Complex x ) { return ctanh( x ); } 221 long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); } 222 223 float asinh( float x ) { return asinhf( x ); } 224 // extern "C" { double asinh( double ); } 225 long double asinh( long double x ) { return asinhl( x ); } 226 float _Complex asinh( float _Complex x ) { return casinhf( x ); } 227 double _Complex asinh( double _Complex x ) { return casinh( x ); } 228 long double _Complex asinh( long double _Complex x ) { return casinhl( x ); } 229 230 float acosh( float x ) { return acoshf( x ); } 231 // extern "C" { double acosh( double ); } 232 long double acosh( long double x ) { return acoshl( x ); } 233 float _Complex acosh( float _Complex x ) { return cacoshf( x ); } 234 double _Complex acosh( double _Complex x ) { return cacosh( x ); } 235 long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); } 236 237 float atanh( float x ) { return atanhf( x ); } 238 // extern "C" { double atanh( double ); } 239 long double atanh( long double x ) { return atanhl( x ); } 240 float _Complex atanh( float _Complex x ) { return catanhf( x ); } 241 double _Complex atanh( double _Complex x ) { return catanh( x ); } 242 long double _Complex atanh( long double _Complex x ) { return catanhl( x ); } 243 } // distribution 230 244 231 245 //---------------------- Error / Gamma ---------------------- 232 246 233 static inline float erf( float x ) { return erff( x ); } 234 // extern "C" { double erf( double ); } 235 static inline long double erf( long double x ) { return erfl( x ); } 236 // float _Complex erf( float _Complex ); 237 // double _Complex erf( double _Complex ); 238 // long double _Complex erf( long double _Complex ); 239 240 static inline float erfc( float x ) { return erfcf( x ); } 241 // extern "C" { double erfc( double ); } 242 static inline long double erfc( long double x ) { return erfcl( x ); } 243 // float _Complex erfc( float _Complex ); 244 // double _Complex erfc( double _Complex ); 245 // long double _Complex erfc( long double _Complex ); 246 247 static inline float lgamma( float x ) { return lgammaf( x ); } 248 // extern "C" { double lgamma( double ); } 249 static inline long double lgamma( long double x ) { return lgammal( x ); } 250 static inline float lgamma( float x, int * sign ) { return lgammaf_r( x, sign ); } 251 static inline double lgamma( double x, int * sign ) { return lgamma_r( x, sign ); } 252 static inline long double lgamma( long double x, int * sign ) { return lgammal_r( x, sign ); } 253 254 static inline float tgamma( float x ) { return tgammaf( x ); } 255 // extern "C" { double tgamma( double ); } 256 static inline long double tgamma( long double x ) { return tgammal( x ); } 247 static inline { 248 float erf( float x ) { return erff( x ); } 249 // extern "C" { double erf( double ); } 250 long double erf( long double x ) { return erfl( x ); } 251 // float _Complex erf( float _Complex ); 252 // double _Complex erf( double _Complex ); 253 // long double _Complex erf( long double _Complex ); 254 255 float erfc( float x ) { return erfcf( x ); } 256 // extern "C" { double erfc( double ); } 257 long double erfc( long double x ) { return erfcl( x ); } 258 // float _Complex erfc( float _Complex ); 259 // double _Complex erfc( double _Complex ); 260 // long double _Complex erfc( long double _Complex ); 261 262 float lgamma( float x ) { return lgammaf( x ); } 263 // extern "C" { double lgamma( double ); } 264 long double lgamma( long double x ) { return lgammal( x ); } 265 float lgamma( float x, int * sign ) { return lgammaf_r( x, sign ); } 266 double lgamma( double x, int * sign ) { return lgamma_r( x, sign ); } 267 long double lgamma( long double x, int * sign ) { return lgammal_r( x, sign ); } 268 269 float tgamma( float x ) { return tgammaf( x ); } 270 // extern "C" { double tgamma( double ); } 271 long double tgamma( long double x ) { return tgammal( x ); } 272 } // distribution 257 273 258 274 //---------------------- Nearest Integer ---------------------- 259 275 260 static inline float floor( float x ) { return floorf( x ); } 261 // extern "C" { double floor( double ); } 262 static inline long double floor( long double x ) { return floorl( x ); } 263 264 static inline float ceil( float x ) { return ceilf( x ); } 265 // extern "C" { double ceil( double ); } 266 static inline long double ceil( long double x ) { return ceill( x ); } 267 268 static inline float trunc( float x ) { return truncf( x ); } 269 // extern "C" { double trunc( double ); } 270 static inline long double trunc( long double x ) { return truncl( x ); } 271 272 static inline float rint( float x ) { return rintf( x ); } 273 // extern "C" { double rint( double x ); } 274 static inline long double rint( long double x ) { return rintl( x ); } 275 static inline long int rint( float x ) { return lrintf( x ); } 276 static inline long int rint( double x ) { return lrint( x ); } 277 static inline long int rint( long double x ) { return lrintl( x ); } 278 static inline long long int rint( float x ) { return llrintf( x ); } 279 static inline long long int rint( double x ) { return llrint( x ); } 280 static inline long long int rint( long double x ) { return llrintl( x ); } 281 282 static inline long int lrint( float x ) { return lrintf( x ); } 283 // extern "C" { long int lrint( double ); } 284 static inline long int lrint( long double x ) { return lrintl( x ); } 285 static inline long long int llrint( float x ) { return llrintf( x ); } 286 // extern "C" { long long int llrint( double ); } 287 static inline long long int llrint( long double x ) { return llrintl( x ); } 288 289 static inline float nearbyint( float x ) { return nearbyintf( x ); } 290 // extern "C" { double nearbyint( double ); } 291 static inline long double nearbyint( long double x ) { return nearbyintl( x ); } 292 293 static inline float round( float x ) { return roundf( x ); } 294 // extern "C" { double round( double x ); } 295 static inline long double round( long double x ) { return roundl( x ); } 296 static inline long int round( float x ) { return lroundf( x ); } 297 static inline long int round( double x ) { return lround( x ); } 298 static inline long int round( long double x ) { return lroundl( x ); } 299 static inline long long int round( float x ) { return llroundf( x ); } 300 static inline long long int round( double x ) { return llround( x ); } 301 static inline long long int round( long double x ) { return llroundl( x ); } 302 303 static inline long int lround( float x ) { return lroundf( x ); } 304 // extern "C" { long int lround( double ); } 305 static inline long int lround( long double x ) { return lroundl( x ); } 306 static inline long long int llround( float x ) { return llroundf( x ); } 307 // extern "C" { long long int llround( double ); } 308 static inline long long int llround( long double x ) { return llroundl( x ); } 276 static inline { 277 signed char floor( signed char n, signed char align ) { return n / align * align; } 278 unsigned char floor( unsigned char n, unsigned char align ) { return n / align * align; } 279 short int floor( short int n, short int align ) { return n / align * align; } 280 unsigned short int floor( unsigned short int n, unsigned short int align ) { return n / align * align; } 281 int floor( int n, int align ) { return n / align * align; } 282 unsigned int floor( unsigned int n, unsigned int align ) { return n / align * align; } 283 long int floor( long int n, long int align ) { return n / align * align; } 284 unsigned long int floor( unsigned long int n, unsigned long int align ) { return n / align * align; } 285 long long int floor( long long int n, long long int align ) { return n / align * align; } 286 unsigned long long int floor( unsigned long long int n, unsigned long long int align ) { return n / align * align; } 287 288 // forall( otype T | { T ?/?( T, T ); T ?*?( T, T ); } ) 289 // T floor( T n, T align ) { return n / align * align; } 290 291 signed char ceiling_div( signed char n, char align ) { return (n + (align - 1)) / align; } 292 unsigned char ceiling_div( unsigned char n, unsigned char align ) { return (n + (align - 1)) / align; } 293 short int ceiling_div( short int n, short int align ) { return (n + (align - 1)) / align; } 294 unsigned short int ceiling_div( unsigned short int n, unsigned short int align ) { return (n + (align - 1)) / align; } 295 int ceiling_div( int n, int align ) { return (n + (align - 1)) / align; } 296 unsigned int ceiling_div( unsigned int n, unsigned int align ) { return (n + (align - 1)) / align; } 297 long int ceiling_div( long int n, long int align ) { return (n + (align - 1)) / align; } 298 unsigned long int ceiling_div( unsigned long int n, unsigned long int align ) { return (n + (align - 1)) / align; } 299 long long int ceiling_div( long long int n, long long int align ) { return (n + (align - 1)) / align; } 300 unsigned long long int ceiling_div( unsigned long long int n, unsigned long long int align ) { return (n + (align - 1)) / align; } 301 302 // forall( otype T | { T ?+?( T, T ); T ?-?( T, T ); T ?%?( T, T ); } ) 303 // T ceiling_div( T n, T align ) { verify( is_pow2( align ) );return (n + (align - 1)) / align; } 304 305 // gcc notices the div/mod pair and saves both so only one div. 306 signed char ceiling( signed char n, signed char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 307 unsigned char ceiling( unsigned char n, unsigned char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 308 short int ceiling( short int n, short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 309 unsigned short int ceiling( unsigned short int n, unsigned short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 310 int ceiling( int n, int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 311 unsigned int ceiling( unsigned int n, unsigned int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 312 long int ceiling( long int n, long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 313 unsigned long int ceiling( unsigned long int n, unsigned long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0) , align); } 314 long long int ceiling( long long int n, long long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 315 unsigned long long int ceiling( unsigned long long int n, unsigned long long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 316 317 // forall( otype T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T ); T ?/?( T, T ); } ) 318 // T ceiling( T n, T align ) { return return floor( n + (n % align != 0 ? align - 1 : 0), align ); *} 319 320 float floor( float x ) { return floorf( x ); } 321 // extern "C" { double floor( double ); } 322 long double floor( long double x ) { return floorl( x ); } 323 324 float ceil( float x ) { return ceilf( x ); } 325 // extern "C" { double ceil( double ); } 326 long double ceil( long double x ) { return ceill( x ); } 327 328 float trunc( float x ) { return truncf( x ); } 329 // extern "C" { double trunc( double ); } 330 long double trunc( long double x ) { return truncl( x ); } 331 332 float rint( float x ) { return rintf( x ); } 333 // extern "C" { double rint( double x ); } 334 long double rint( long double x ) { return rintl( x ); } 335 long int rint( float x ) { return lrintf( x ); } 336 long int rint( double x ) { return lrint( x ); } 337 long int rint( long double x ) { return lrintl( x ); } 338 long long int rint( float x ) { return llrintf( x ); } 339 long long int rint( double x ) { return llrint( x ); } 340 long long int rint( long double x ) { return llrintl( x ); } 341 342 long int lrint( float x ) { return lrintf( x ); } 343 // extern "C" { long int lrint( double ); } 344 long int lrint( long double x ) { return lrintl( x ); } 345 long long int llrint( float x ) { return llrintf( x ); } 346 // extern "C" { long long int llrint( double ); } 347 long long int llrint( long double x ) { return llrintl( x ); } 348 349 float nearbyint( float x ) { return nearbyintf( x ); } 350 // extern "C" { double nearbyint( double ); } 351 long double nearbyint( long double x ) { return nearbyintl( x ); } 352 353 float round( float x ) { return roundf( x ); } 354 // extern "C" { double round( double x ); } 355 long double round( long double x ) { return roundl( x ); } 356 long int round( float x ) { return lroundf( x ); } 357 long int round( double x ) { return lround( x ); } 358 long int round( long double x ) { return lroundl( x ); } 359 long long int round( float x ) { return llroundf( x ); } 360 long long int round( double x ) { return llround( x ); } 361 long long int round( long double x ) { return llroundl( x ); } 362 363 long int lround( float x ) { return lroundf( x ); } 364 // extern "C" { long int lround( double ); } 365 long int lround( long double x ) { return lroundl( x ); } 366 long long int llround( float x ) { return llroundf( x ); } 367 // extern "C" { long long int llround( double ); } 368 long long int llround( long double x ) { return llroundl( x ); } 369 } // distribution 309 370 310 371 //---------------------- Manipulation ---------------------- 311 372 312 static inline float copysign( float x, float y ) { return copysignf( x, y ); } 313 // extern "C" { double copysign( double, double ); } 314 static inline long double copysign( long double x, long double y ) { return copysignl( x, y ); } 315 316 static inline float frexp( float x, int * ip ) { return frexpf( x, ip ); } 317 // extern "C" { double frexp( double, int * ); } 318 static inline long double frexp( long double x, int * ip ) { return frexpl( x, ip ); } 319 320 static inline float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); } 321 // extern "C" { double ldexp( double, int ); } 322 static inline long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); } 323 324 static inline [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; } 325 static inline float modf( float x, float * i ) { return modff( x, i ); } 326 static inline [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; } 327 // extern "C" { double modf( double, double * ); } 328 static inline [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; } 329 static inline long double modf( long double x, long double * i ) { return modfl( x, i ); } 330 331 static inline float nextafter( float x, float y ) { return nextafterf( x, y ); } 332 // extern "C" { double nextafter( double, double ); } 333 static inline long double nextafter( long double x, long double y ) { return nextafterl( x, y ); } 334 335 static inline float nexttoward( float x, long double y ) { return nexttowardf( x, y ); } 336 // extern "C" { double nexttoward( double, long double ); } 337 static inline long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); } 338 339 static inline float scalbn( float x, int exp ) { return scalbnf( x, exp ); } 340 // extern "C" { double scalbn( double, int ); } 341 static inline long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); } 342 static inline float scalbn( float x, long int exp ) { return scalblnf( x, exp ); } 343 static inline double scalbn( double x, long int exp ) { return scalbln( x, exp ); } 344 static inline long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); } 345 346 static inline float scalbln( float x, long int exp ) { return scalblnf( x, exp ); } 347 // extern "C" { double scalbln( double, long int ); } 348 static inline long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); } 373 static inline { 374 float copysign( float x, float y ) { return copysignf( x, y ); } 375 // extern "C" { double copysign( double, double ); } 376 long double copysign( long double x, long double y ) { return copysignl( x, y ); } 377 378 float frexp( float x, int * ip ) { return frexpf( x, ip ); } 379 // extern "C" { double frexp( double, int * ); } 380 long double frexp( long double x, int * ip ) { return frexpl( x, ip ); } 381 382 float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); } 383 // extern "C" { double ldexp( double, int ); } 384 long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); } 385 386 [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; } 387 float modf( float x, float * i ) { return modff( x, i ); } 388 [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; } 389 // extern "C" { double modf( double, double * ); } 390 [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; } 391 long double modf( long double x, long double * i ) { return modfl( x, i ); } 392 393 float nextafter( float x, float y ) { return nextafterf( x, y ); } 394 // extern "C" { double nextafter( double, double ); } 395 long double nextafter( long double x, long double y ) { return nextafterl( x, y ); } 396 397 float nexttoward( float x, long double y ) { return nexttowardf( x, y ); } 398 // extern "C" { double nexttoward( double, long double ); } 399 long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); } 400 401 float scalbn( float x, int exp ) { return scalbnf( x, exp ); } 402 // extern "C" { double scalbn( double, int ); } 403 long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); } 404 float scalbn( float x, long int exp ) { return scalblnf( x, exp ); } 405 double scalbn( double x, long int exp ) { return scalbln( x, exp ); } 406 long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); } 407 408 float scalbln( float x, long int exp ) { return scalblnf( x, exp ); } 409 // extern "C" { double scalbln( double, long int ); } 410 long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); } 411 } // distribution 349 412 350 413 //--------------------------------------- 351 414 352 #include "common.hfa" 353 354 //--------------------------------------- 355 356 forall( otype T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T );T ?*?( T, T ); } ) 357 T lerp( T x, T y, T a ) { return x * ((T){1} - a) + y * a; } 358 359 forall( otype T | { void ?{}( T &, zero_t ); void ?{}( T &, one_t ); int ?<?( T, T ); } ) 360 T step( T edge, T x ) { return x < edge ? (T){0} : (T){1}; } 361 362 forall( otype T | { void ?{}( T &, int ); T clamp( T, T, T ); T ?-?( T, T ); T ?*?( T, T ); T ?/?( T, T ); } ) 363 T smoothstep( T edge0, T edge1, T x ) { T t = clamp( (x - edge0) / (edge1 - edge0), (T){0}, (T){1} ); return t * t * ((T){3} - (T){2} * t); } 415 static inline { 416 forall( otype T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T );T ?*?( T, T ); } ) 417 T lerp( T x, T y, T a ) { return x * ((T){1} - a) + y * a; } 418 419 forall( otype T | { void ?{}( T &, zero_t ); void ?{}( T &, one_t ); int ?<?( T, T ); } ) 420 T step( T edge, T x ) { return x < edge ? (T){0} : (T){1}; } 421 422 forall( otype T | { void ?{}( T &, int ); T clamp( T, T, T ); T ?-?( T, T ); T ?*?( T, T ); T ?/?( T, T ); } ) 423 T smoothstep( T edge0, T edge1, T x ) { T t = clamp( (x - edge0) / (edge1 - edge0), (T){0}, (T){1} ); return t * t * ((T){3} - (T){2} * t); } 424 } // distribution 364 425 365 426 // Local Variables: //
Note: See TracChangeset
for help on using the changeset viewer.