Changeset 95bda0a for libcfa/src/math.hfa
- Timestamp:
- Oct 8, 2022, 4:29:32 PM (19 months ago)
- Branches:
- ADT, ast-experimental, master
- Children:
- 301071a, affb51b
- Parents:
- b2ddaf3
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
libcfa/src/math.hfa
rb2ddaf3 r95bda0a 10 10 // Created On : Mon Apr 18 23:37:04 2016 11 11 // Last Modified By : Peter A. Buhr 12 // Last Modified On : Thu Apr 15 11:47:56 202113 // Update Count : 13 212 // Last Modified On : Sat Oct 8 08:40:42 2022 13 // Update Count : 136 14 14 // 15 15 … … 26 26 //---------------------- General ---------------------- 27 27 28 static inline {28 static inline __attribute__((always_inline)) { 29 29 float ?%?( float x, float y ) { return fmodf( x, y ); } 30 30 float fmod( float x, float y ) { return fmodf( x, y ); } … … 64 64 //---------------------- Exponential ---------------------- 65 65 66 static inline {66 static inline __attribute__((always_inline)) { 67 67 float exp( float x ) { return expf( x ); } 68 68 // extern "C" { double exp( double ); } … … 93 93 //---------------------- Logarithm ---------------------- 94 94 95 static inline {95 static inline __attribute__((always_inline)) { 96 96 float log( float x ) { return logf( x ); } 97 97 // extern "C" { double log( double ); } … … 148 148 } // distribution 149 149 150 //---------------------- Trigonometric ---------------------- 151 152 static inline { 153 float sin( float x ) { return sinf( x ); } 154 // extern "C" { double sin( double ); } 155 long double sin( long double x ) { return sinl( x ); } 156 float _Complex sin( float _Complex x ) { return csinf( x ); } 157 double _Complex sin( double _Complex x ) { return csin( x ); } 158 long double _Complex sin( long double _Complex x ) { return csinl( x ); } 159 160 float cos( float x ) { return cosf( x ); } 161 // extern "C" { double cos( double ); } 162 long double cos( long double x ) { return cosl( x ); } 163 float _Complex cos( float _Complex x ) { return ccosf( x ); } 164 double _Complex cos( double _Complex x ) { return ccos( x ); } 165 long double _Complex cos( long double _Complex x ) { return ccosl( x ); } 166 167 float tan( float x ) { return tanf( x ); } 168 // extern "C" { double tan( double ); } 169 long double tan( long double x ) { return tanl( x ); } 170 float _Complex tan( float _Complex x ) { return ctanf( x ); } 171 double _Complex tan( double _Complex x ) { return ctan( x ); } 172 long double _Complex tan( long double _Complex x ) { return ctanl( x ); } 173 174 float asin( float x ) { return asinf( x ); } 175 // extern "C" { double asin( double ); } 176 long double asin( long double x ) { return asinl( x ); } 177 float _Complex asin( float _Complex x ) { return casinf( x ); } 178 double _Complex asin( double _Complex x ) { return casin( x ); } 179 long double _Complex asin( long double _Complex x ) { return casinl( x ); } 180 181 float acos( float x ) { return acosf( x ); } 182 // extern "C" { double acos( double ); } 183 long double acos( long double x ) { return acosl( x ); } 184 float _Complex acos( float _Complex x ) { return cacosf( x ); } 185 double _Complex acos( double _Complex x ) { return cacos( x ); } 186 long double _Complex acos( long double _Complex x ) { return cacosl( x ); } 187 188 float atan( float x ) { return atanf( x ); } 189 // extern "C" { double atan( double ); } 190 long double atan( long double x ) { return atanl( x ); } 191 float _Complex atan( float _Complex x ) { return catanf( x ); } 192 double _Complex atan( double _Complex x ) { return catan( x ); } 193 long double _Complex atan( long double _Complex x ) { return catanl( x ); } 194 195 float atan2( float x, float y ) { return atan2f( x, y ); } 196 // extern "C" { double atan2( double, double ); } 197 long double atan2( long double x, long double y ) { return atan2l( x, y ); } 198 199 // alternative name for atan2 200 float atan( float x, float y ) { return atan2f( x, y ); } 201 double atan( double x, double y ) { return atan2( x, y ); } 202 long double atan( long double x, long double y ) { return atan2l( x, y ); } 203 } // distribution 204 205 //---------------------- Hyperbolic ---------------------- 206 207 static inline { 208 float sinh( float x ) { return sinhf( x ); } 209 // extern "C" { double sinh( double ); } 210 long double sinh( long double x ) { return sinhl( x ); } 211 float _Complex sinh( float _Complex x ) { return csinhf( x ); } 212 double _Complex sinh( double _Complex x ) { return csinh( x ); } 213 long double _Complex sinh( long double _Complex x ) { return csinhl( x ); } 214 215 float cosh( float x ) { return coshf( x ); } 216 // extern "C" { double cosh( double ); } 217 long double cosh( long double x ) { return coshl( x ); } 218 float _Complex cosh( float _Complex x ) { return ccoshf( x ); } 219 double _Complex cosh( double _Complex x ) { return ccosh( x ); } 220 long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); } 221 222 float tanh( float x ) { return tanhf( x ); } 223 // extern "C" { double tanh( double ); } 224 long double tanh( long double x ) { return tanhl( x ); } 225 float _Complex tanh( float _Complex x ) { return ctanhf( x ); } 226 double _Complex tanh( double _Complex x ) { return ctanh( x ); } 227 long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); } 228 229 float asinh( float x ) { return asinhf( x ); } 230 // extern "C" { double asinh( double ); } 231 long double asinh( long double x ) { return asinhl( x ); } 232 float _Complex asinh( float _Complex x ) { return casinhf( x ); } 233 double _Complex asinh( double _Complex x ) { return casinh( x ); } 234 long double _Complex asinh( long double _Complex x ) { return casinhl( x ); } 235 236 float acosh( float x ) { return acoshf( x ); } 237 // extern "C" { double acosh( double ); } 238 long double acosh( long double x ) { return acoshl( x ); } 239 float _Complex acosh( float _Complex x ) { return cacoshf( x ); } 240 double _Complex acosh( double _Complex x ) { return cacosh( x ); } 241 long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); } 242 243 float atanh( float x ) { return atanhf( x ); } 244 // extern "C" { double atanh( double ); } 245 long double atanh( long double x ) { return atanhl( x ); } 246 float _Complex atanh( float _Complex x ) { return catanhf( x ); } 247 double _Complex atanh( double _Complex x ) { return catanh( x ); } 248 long double _Complex atanh( long double _Complex x ) { return catanhl( x ); } 249 } // distribution 250 251 //---------------------- Error / Gamma ---------------------- 252 253 static inline { 254 float erf( float x ) { return erff( x ); } 255 // extern "C" { double erf( double ); } 256 long double erf( long double x ) { return erfl( x ); } 257 // float _Complex erf( float _Complex ); 258 // double _Complex erf( double _Complex ); 259 // long double _Complex erf( long double _Complex ); 260 261 float erfc( float x ) { return erfcf( x ); } 262 // extern "C" { double erfc( double ); } 263 long double erfc( long double x ) { return erfcl( x ); } 264 // float _Complex erfc( float _Complex ); 265 // double _Complex erfc( double _Complex ); 266 // long double _Complex erfc( long double _Complex ); 267 268 float lgamma( float x ) { return lgammaf( x ); } 269 // extern "C" { double lgamma( double ); } 270 long double lgamma( long double x ) { return lgammal( x ); } 271 float lgamma( float x, int * sign ) { return lgammaf_r( x, sign ); } 272 double lgamma( double x, int * sign ) { return lgamma_r( x, sign ); } 273 long double lgamma( long double x, int * sign ) { return lgammal_r( x, sign ); } 274 275 float tgamma( float x ) { return tgammaf( x ); } 276 // extern "C" { double tgamma( double ); } 277 long double tgamma( long double x ) { return tgammal( x ); } 278 } // distribution 279 280 //---------------------- Nearest Integer ---------------------- 281 282 static inline { 283 signed char floor( signed char n, signed char align ) { return n / align * align; } 284 unsigned char floor( unsigned char n, unsigned char align ) { return n / align * align; } 285 short int floor( short int n, short int align ) { return n / align * align; } 286 unsigned short int floor( unsigned short int n, unsigned short int align ) { return n / align * align; } 287 int floor( int n, int align ) { return n / align * align; } 288 unsigned int floor( unsigned int n, unsigned int align ) { return n / align * align; } 289 long int floor( long int n, long int align ) { return n / align * align; } 290 unsigned long int floor( unsigned long int n, unsigned long int align ) { return n / align * align; } 291 long long int floor( long long int n, long long int align ) { return n / align * align; } 292 unsigned long long int floor( unsigned long long int n, unsigned long long int align ) { return n / align * align; } 293 294 // forall( T | { T ?/?( T, T ); T ?*?( T, T ); } ) 295 // T floor( T n, T align ) { return n / align * align; } 296 297 signed char ceiling_div( signed char n, char align ) { return (n + (align - 1)) / align; } 298 unsigned char ceiling_div( unsigned char n, unsigned char align ) { return (n + (align - 1)) / align; } 299 short int ceiling_div( short int n, short int align ) { return (n + (align - 1)) / align; } 300 unsigned short int ceiling_div( unsigned short int n, unsigned short int align ) { return (n + (align - 1)) / align; } 301 int ceiling_div( int n, int align ) { return (n + (align - 1)) / align; } 302 unsigned int ceiling_div( unsigned int n, unsigned int align ) { return (n + (align - 1)) / align; } 303 long int ceiling_div( long int n, long int align ) { return (n + (align - 1)) / align; } 304 unsigned long int ceiling_div( unsigned long int n, unsigned long int align ) { return (n + (align - 1)) / align; } 305 long long int ceiling_div( long long int n, long long int align ) { return (n + (align - 1)) / align; } 306 unsigned long long int ceiling_div( unsigned long long int n, unsigned long long int align ) { return (n + (align - 1)) / align; } 307 308 // forall( T | { T ?+?( T, T ); T ?-?( T, T ); T ?%?( T, T ); } ) 309 // T ceiling_div( T n, T align ) { verify( is_pow2( align ) );return (n + (align - 1)) / align; } 310 311 // gcc notices the div/mod pair and saves both so only one div. 312 signed char ceiling( signed char n, signed char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 313 unsigned char ceiling( unsigned char n, unsigned char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 314 short int ceiling( short int n, short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 315 unsigned short int ceiling( unsigned short int n, unsigned short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 316 int ceiling( int n, int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 317 unsigned int ceiling( unsigned int n, unsigned int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 318 long int ceiling( long int n, long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 319 unsigned long int ceiling( unsigned long int n, unsigned long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0) , align); } 320 long long int ceiling( long long int n, long long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 321 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 ); } 322 323 // forall( T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T ); T ?/?( T, T ); } ) 324 // T ceiling( T n, T align ) { return return floor( n + (n % align != 0 ? align - 1 : 0), align ); *} 325 326 float floor( float x ) { return floorf( x ); } 327 // extern "C" { double floor( double ); } 328 long double floor( long double x ) { return floorl( x ); } 329 330 float ceil( float x ) { return ceilf( x ); } 331 // extern "C" { double ceil( double ); } 332 long double ceil( long double x ) { return ceill( x ); } 333 334 float trunc( float x ) { return truncf( x ); } 335 // extern "C" { double trunc( double ); } 336 long double trunc( long double x ) { return truncl( x ); } 337 338 float rint( float x ) { return rintf( x ); } 339 // extern "C" { double rint( double x ); } 340 long double rint( long double x ) { return rintl( x ); } 341 long int rint( float x ) { return lrintf( x ); } 342 long int rint( double x ) { return lrint( x ); } 343 long int rint( long double x ) { return lrintl( x ); } 344 long long int rint( float x ) { return llrintf( x ); } 345 long long int rint( double x ) { return llrint( x ); } 346 long long int rint( long double x ) { return llrintl( x ); } 347 348 long int lrint( float x ) { return lrintf( x ); } 349 // extern "C" { long int lrint( double ); } 350 long int lrint( long double x ) { return lrintl( x ); } 351 long long int llrint( float x ) { return llrintf( x ); } 352 // extern "C" { long long int llrint( double ); } 353 long long int llrint( long double x ) { return llrintl( x ); } 354 355 float nearbyint( float x ) { return nearbyintf( x ); } 356 // extern "C" { double nearbyint( double ); } 357 long double nearbyint( long double x ) { return nearbyintl( x ); } 358 359 float round( float x ) { return roundf( x ); } 360 // extern "C" { double round( double x ); } 361 long double round( long double x ) { return roundl( x ); } 362 long int round( float x ) { return lroundf( x ); } 363 long int round( double x ) { return lround( x ); } 364 long int round( long double x ) { return lroundl( x ); } 365 long long int round( float x ) { return llroundf( x ); } 366 long long int round( double x ) { return llround( x ); } 367 long long int round( long double x ) { return llroundl( x ); } 368 369 long int lround( float x ) { return lroundf( x ); } 370 // extern "C" { long int lround( double ); } 371 long int lround( long double x ) { return lroundl( x ); } 372 long long int llround( float x ) { return llroundf( x ); } 373 // extern "C" { long long int llround( double ); } 374 long long int llround( long double x ) { return llroundl( x ); } 375 } // distribution 376 377 //---------------------- Manipulation ---------------------- 378 379 static inline { 380 float copysign( float x, float y ) { return copysignf( x, y ); } 381 // extern "C" { double copysign( double, double ); } 382 long double copysign( long double x, long double y ) { return copysignl( x, y ); } 383 384 float frexp( float x, int * ip ) { return frexpf( x, ip ); } 385 // extern "C" { double frexp( double, int * ); } 386 long double frexp( long double x, int * ip ) { return frexpl( x, ip ); } 387 388 float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); } 389 // extern "C" { double ldexp( double, int ); } 390 long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); } 391 392 [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; } 393 float modf( float x, float * i ) { return modff( x, i ); } 394 [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; } 395 // extern "C" { double modf( double, double * ); } 396 [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; } 397 long double modf( long double x, long double * i ) { return modfl( x, i ); } 398 399 float nextafter( float x, float y ) { return nextafterf( x, y ); } 400 // extern "C" { double nextafter( double, double ); } 401 long double nextafter( long double x, long double y ) { return nextafterl( x, y ); } 402 403 float nexttoward( float x, long double y ) { return nexttowardf( x, y ); } 404 // extern "C" { double nexttoward( double, long double ); } 405 long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); } 406 407 float scalbn( float x, int exp ) { return scalbnf( x, exp ); } 408 // extern "C" { double scalbn( double, int ); } 409 long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); } 410 float scalbn( float x, long int exp ) { return scalblnf( x, exp ); } 411 double scalbn( double x, long int exp ) { return scalbln( x, exp ); } 412 long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); } 413 414 float scalbln( float x, long int exp ) { return scalblnf( x, exp ); } 415 // extern "C" { double scalbln( double, long int ); } 416 long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); } 417 } // distribution 418 419 //--------------------------------------- 420 421 static inline { 422 forall( T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T );T ?*?( T, T ); } ) 423 T lerp( T x, T y, T a ) { return x * ((T){1} - a) + y * a; } 424 425 forall( T | { void ?{}( T &, zero_t ); void ?{}( T &, one_t ); int ?<?( T, T ); } ) 426 T step( T edge, T x ) { return x < edge ? (T){0} : (T){1}; } 427 428 forall( T | { void ?{}( T &, int ); T clamp( T, T, T ); T ?-?( T, T ); T ?*?( T, T ); T ?/?( T, T ); } ) 429 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); } 430 } // distribution 431 432 static inline unsigned long long log2_u32_32(unsigned long long val) { 150 static inline unsigned long long log2_u32_32( unsigned long long val ) { 433 151 enum { 434 152 TABLE_BITS = 6, … … 520 238 } 521 239 240 //---------------------- Trigonometric ---------------------- 241 242 static inline __attribute__((always_inline)) { 243 float sin( float x ) { return sinf( x ); } 244 // extern "C" { double sin( double ); } 245 long double sin( long double x ) { return sinl( x ); } 246 float _Complex sin( float _Complex x ) { return csinf( x ); } 247 double _Complex sin( double _Complex x ) { return csin( x ); } 248 long double _Complex sin( long double _Complex x ) { return csinl( x ); } 249 250 float cos( float x ) { return cosf( x ); } 251 // extern "C" { double cos( double ); } 252 long double cos( long double x ) { return cosl( x ); } 253 float _Complex cos( float _Complex x ) { return ccosf( x ); } 254 double _Complex cos( double _Complex x ) { return ccos( x ); } 255 long double _Complex cos( long double _Complex x ) { return ccosl( x ); } 256 257 float tan( float x ) { return tanf( x ); } 258 // extern "C" { double tan( double ); } 259 long double tan( long double x ) { return tanl( x ); } 260 float _Complex tan( float _Complex x ) { return ctanf( x ); } 261 double _Complex tan( double _Complex x ) { return ctan( x ); } 262 long double _Complex tan( long double _Complex x ) { return ctanl( x ); } 263 264 float asin( float x ) { return asinf( x ); } 265 // extern "C" { double asin( double ); } 266 long double asin( long double x ) { return asinl( x ); } 267 float _Complex asin( float _Complex x ) { return casinf( x ); } 268 double _Complex asin( double _Complex x ) { return casin( x ); } 269 long double _Complex asin( long double _Complex x ) { return casinl( x ); } 270 271 float acos( float x ) { return acosf( x ); } 272 // extern "C" { double acos( double ); } 273 long double acos( long double x ) { return acosl( x ); } 274 float _Complex acos( float _Complex x ) { return cacosf( x ); } 275 double _Complex acos( double _Complex x ) { return cacos( x ); } 276 long double _Complex acos( long double _Complex x ) { return cacosl( x ); } 277 278 float atan( float x ) { return atanf( x ); } 279 // extern "C" { double atan( double ); } 280 long double atan( long double x ) { return atanl( x ); } 281 float _Complex atan( float _Complex x ) { return catanf( x ); } 282 double _Complex atan( double _Complex x ) { return catan( x ); } 283 long double _Complex atan( long double _Complex x ) { return catanl( x ); } 284 285 float atan2( float x, float y ) { return atan2f( x, y ); } 286 // extern "C" { double atan2( double, double ); } 287 long double atan2( long double x, long double y ) { return atan2l( x, y ); } 288 289 // alternative name for atan2 290 float atan( float x, float y ) { return atan2f( x, y ); } 291 double atan( double x, double y ) { return atan2( x, y ); } 292 long double atan( long double x, long double y ) { return atan2l( x, y ); } 293 } // distribution 294 295 //---------------------- Hyperbolic ---------------------- 296 297 static inline __attribute__((always_inline)) { 298 float sinh( float x ) { return sinhf( x ); } 299 // extern "C" { double sinh( double ); } 300 long double sinh( long double x ) { return sinhl( x ); } 301 float _Complex sinh( float _Complex x ) { return csinhf( x ); } 302 double _Complex sinh( double _Complex x ) { return csinh( x ); } 303 long double _Complex sinh( long double _Complex x ) { return csinhl( x ); } 304 305 float cosh( float x ) { return coshf( x ); } 306 // extern "C" { double cosh( double ); } 307 long double cosh( long double x ) { return coshl( x ); } 308 float _Complex cosh( float _Complex x ) { return ccoshf( x ); } 309 double _Complex cosh( double _Complex x ) { return ccosh( x ); } 310 long double _Complex cosh( long double _Complex x ) { return ccoshl( x ); } 311 312 float tanh( float x ) { return tanhf( x ); } 313 // extern "C" { double tanh( double ); } 314 long double tanh( long double x ) { return tanhl( x ); } 315 float _Complex tanh( float _Complex x ) { return ctanhf( x ); } 316 double _Complex tanh( double _Complex x ) { return ctanh( x ); } 317 long double _Complex tanh( long double _Complex x ) { return ctanhl( x ); } 318 319 float asinh( float x ) { return asinhf( x ); } 320 // extern "C" { double asinh( double ); } 321 long double asinh( long double x ) { return asinhl( x ); } 322 float _Complex asinh( float _Complex x ) { return casinhf( x ); } 323 double _Complex asinh( double _Complex x ) { return casinh( x ); } 324 long double _Complex asinh( long double _Complex x ) { return casinhl( x ); } 325 326 float acosh( float x ) { return acoshf( x ); } 327 // extern "C" { double acosh( double ); } 328 long double acosh( long double x ) { return acoshl( x ); } 329 float _Complex acosh( float _Complex x ) { return cacoshf( x ); } 330 double _Complex acosh( double _Complex x ) { return cacosh( x ); } 331 long double _Complex acosh( long double _Complex x ) { return cacoshl( x ); } 332 333 float atanh( float x ) { return atanhf( x ); } 334 // extern "C" { double atanh( double ); } 335 long double atanh( long double x ) { return atanhl( x ); } 336 float _Complex atanh( float _Complex x ) { return catanhf( x ); } 337 double _Complex atanh( double _Complex x ) { return catanh( x ); } 338 long double _Complex atanh( long double _Complex x ) { return catanhl( x ); } 339 } // distribution 340 341 //---------------------- Error / Gamma ---------------------- 342 343 static inline __attribute__((always_inline)) { 344 float erf( float x ) { return erff( x ); } 345 // extern "C" { double erf( double ); } 346 long double erf( long double x ) { return erfl( x ); } 347 // float _Complex erf( float _Complex ); 348 // double _Complex erf( double _Complex ); 349 // long double _Complex erf( long double _Complex ); 350 351 float erfc( float x ) { return erfcf( x ); } 352 // extern "C" { double erfc( double ); } 353 long double erfc( long double x ) { return erfcl( x ); } 354 // float _Complex erfc( float _Complex ); 355 // double _Complex erfc( double _Complex ); 356 // long double _Complex erfc( long double _Complex ); 357 358 float lgamma( float x ) { return lgammaf( x ); } 359 // extern "C" { double lgamma( double ); } 360 long double lgamma( long double x ) { return lgammal( x ); } 361 float lgamma( float x, int * sign ) { return lgammaf_r( x, sign ); } 362 double lgamma( double x, int * sign ) { return lgamma_r( x, sign ); } 363 long double lgamma( long double x, int * sign ) { return lgammal_r( x, sign ); } 364 365 float tgamma( float x ) { return tgammaf( x ); } 366 // extern "C" { double tgamma( double ); } 367 long double tgamma( long double x ) { return tgammal( x ); } 368 } // distribution 369 370 //---------------------- Nearest Integer ---------------------- 371 372 inline __attribute__((always_inline)) static { 373 signed char floor( signed char n, signed char align ) { return n / align * align; } 374 unsigned char floor( unsigned char n, unsigned char align ) { return n / align * align; } 375 short int floor( short int n, short int align ) { return n / align * align; } 376 unsigned short int floor( unsigned short int n, unsigned short int align ) { return n / align * align; } 377 int floor( int n, int align ) { return n / align * align; } 378 unsigned int floor( unsigned int n, unsigned int align ) { return n / align * align; } 379 long int floor( long int n, long int align ) { return n / align * align; } 380 unsigned long int floor( unsigned long int n, unsigned long int align ) { return n / align * align; } 381 long long int floor( long long int n, long long int align ) { return n / align * align; } 382 unsigned long long int floor( unsigned long long int n, unsigned long long int align ) { return n / align * align; } 383 384 // forall( T | { T ?/?( T, T ); T ?*?( T, T ); } ) 385 // T floor( T n, T align ) { return n / align * align; } 386 387 signed char ceiling_div( signed char n, char align ) { return (n + (align - 1)) / align; } 388 unsigned char ceiling_div( unsigned char n, unsigned char align ) { return (n + (align - 1)) / align; } 389 short int ceiling_div( short int n, short int align ) { return (n + (align - 1)) / align; } 390 unsigned short int ceiling_div( unsigned short int n, unsigned short int align ) { return (n + (align - 1)) / align; } 391 int ceiling_div( int n, int align ) { return (n + (align - 1)) / align; } 392 unsigned int ceiling_div( unsigned int n, unsigned int align ) { return (n + (align - 1)) / align; } 393 long int ceiling_div( long int n, long int align ) { return (n + (align - 1)) / align; } 394 unsigned long int ceiling_div( unsigned long int n, unsigned long int align ) { return (n + (align - 1)) / align; } 395 long long int ceiling_div( long long int n, long long int align ) { return (n + (align - 1)) / align; } 396 unsigned long long int ceiling_div( unsigned long long int n, unsigned long long int align ) { return (n + (align - 1)) / align; } 397 398 // forall( T | { T ?+?( T, T ); T ?-?( T, T ); T ?%?( T, T ); } ) 399 // T ceiling_div( T n, T align ) { verify( is_pow2( align ) );return (n + (align - 1)) / align; } 400 401 // gcc notices the div/mod pair and saves both so only one div. 402 signed char ceiling( signed char n, signed char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 403 unsigned char ceiling( unsigned char n, unsigned char align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 404 short int ceiling( short int n, short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 405 unsigned short int ceiling( unsigned short int n, unsigned short int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 406 int ceiling( int n, int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 407 unsigned int ceiling( unsigned int n, unsigned int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 408 long int ceiling( long int n, long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 409 unsigned long int ceiling( unsigned long int n, unsigned long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0) , align); } 410 long long int ceiling( long long int n, long long int align ) { return floor( n + (n % align != 0 ? align - 1 : 0), align ); } 411 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 ); } 412 413 // forall( T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T ); T ?/?( T, T ); } ) 414 // T ceiling( T n, T align ) { return return floor( n + (n % align != 0 ? align - 1 : 0), align ); *} 415 416 float floor( float x ) { return floorf( x ); } 417 // extern "C" { double floor( double ); } 418 long double floor( long double x ) { return floorl( x ); } 419 420 float ceil( float x ) { return ceilf( x ); } 421 // extern "C" { double ceil( double ); } 422 long double ceil( long double x ) { return ceill( x ); } 423 424 float trunc( float x ) { return truncf( x ); } 425 // extern "C" { double trunc( double ); } 426 long double trunc( long double x ) { return truncl( x ); } 427 428 float rint( float x ) { return rintf( x ); } 429 // extern "C" { double rint( double x ); } 430 long double rint( long double x ) { return rintl( x ); } 431 long int rint( float x ) { return lrintf( x ); } 432 long int rint( double x ) { return lrint( x ); } 433 long int rint( long double x ) { return lrintl( x ); } 434 long long int rint( float x ) { return llrintf( x ); } 435 long long int rint( double x ) { return llrint( x ); } 436 long long int rint( long double x ) { return llrintl( x ); } 437 438 long int lrint( float x ) { return lrintf( x ); } 439 // extern "C" { long int lrint( double ); } 440 long int lrint( long double x ) { return lrintl( x ); } 441 long long int llrint( float x ) { return llrintf( x ); } 442 // extern "C" { long long int llrint( double ); } 443 long long int llrint( long double x ) { return llrintl( x ); } 444 445 float nearbyint( float x ) { return nearbyintf( x ); } 446 // extern "C" { double nearbyint( double ); } 447 long double nearbyint( long double x ) { return nearbyintl( x ); } 448 449 float round( float x ) { return roundf( x ); } 450 // extern "C" { double round( double x ); } 451 long double round( long double x ) { return roundl( x ); } 452 long int round( float x ) { return lroundf( x ); } 453 long int round( double x ) { return lround( x ); } 454 long int round( long double x ) { return lroundl( x ); } 455 long long int round( float x ) { return llroundf( x ); } 456 long long int round( double x ) { return llround( x ); } 457 long long int round( long double x ) { return llroundl( x ); } 458 459 long int lround( float x ) { return lroundf( x ); } 460 // extern "C" { long int lround( double ); } 461 long int lround( long double x ) { return lroundl( x ); } 462 long long int llround( float x ) { return llroundf( x ); } 463 // extern "C" { long long int llround( double ); } 464 long long int llround( long double x ) { return llroundl( x ); } 465 } // distribution 466 467 //---------------------- Manipulation ---------------------- 468 469 static inline __attribute__((always_inline)) { 470 float copysign( float x, float y ) { return copysignf( x, y ); } 471 // extern "C" { double copysign( double, double ); } 472 long double copysign( long double x, long double y ) { return copysignl( x, y ); } 473 474 float frexp( float x, int * ip ) { return frexpf( x, ip ); } 475 // extern "C" { double frexp( double, int * ); } 476 long double frexp( long double x, int * ip ) { return frexpl( x, ip ); } 477 478 float ldexp( float x, int exp2 ) { return ldexpf( x, exp2 ); } 479 // extern "C" { double ldexp( double, int ); } 480 long double ldexp( long double x, int exp2 ) { return ldexpl( x, exp2 ); } 481 482 [ float, float ] modf( float x ) { float i; x = modff( x, &i ); return [ i, x ]; } 483 float modf( float x, float * i ) { return modff( x, i ); } 484 [ double, double ] modf( double x ) { double i; x = modf( x, &i ); return [ i, x ]; } 485 // extern "C" { double modf( double, double * ); } 486 [ long double, long double ] modf( long double x ) { long double i; x = modfl( x, &i ); return [ i, x ]; } 487 long double modf( long double x, long double * i ) { return modfl( x, i ); } 488 489 float nextafter( float x, float y ) { return nextafterf( x, y ); } 490 // extern "C" { double nextafter( double, double ); } 491 long double nextafter( long double x, long double y ) { return nextafterl( x, y ); } 492 493 float nexttoward( float x, long double y ) { return nexttowardf( x, y ); } 494 // extern "C" { double nexttoward( double, long double ); } 495 long double nexttoward( long double x, long double y ) { return nexttowardl( x, y ); } 496 497 float scalbn( float x, int exp ) { return scalbnf( x, exp ); } 498 // extern "C" { double scalbn( double, int ); } 499 long double scalbn( long double x, int exp ) { return scalbnl( x, exp ); } 500 float scalbn( float x, long int exp ) { return scalblnf( x, exp ); } 501 double scalbn( double x, long int exp ) { return scalbln( x, exp ); } 502 long double scalbn( long double x, long int exp ) { return scalblnl( x, exp ); } 503 504 float scalbln( float x, long int exp ) { return scalblnf( x, exp ); } 505 // extern "C" { double scalbln( double, long int ); } 506 long double scalbln( long double x, long int exp ) { return scalblnl( x, exp ); } 507 } // distribution 508 509 //--------------------------------------- 510 511 static inline __attribute__((always_inline)) { 512 forall( T | { void ?{}( T &, one_t ); T ?+?( T, T ); T ?-?( T, T );T ?*?( T, T ); } ) 513 T lerp( T x, T y, T a ) { return x * ((T){1} - a) + y * a; } 514 515 forall( T | { void ?{}( T &, zero_t ); void ?{}( T &, one_t ); int ?<?( T, T ); } ) 516 T step( T edge, T x ) { return x < edge ? (T){0} : (T){1}; } 517 518 forall( T | { void ?{}( T &, int ); T clamp( T, T, T ); T ?-?( T, T ); T ?*?( T, T ); T ?/?( T, T ); } ) 519 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); } 520 } // distribution 521 522 522 // Local Variables: // 523 523 // mode: c //
Note: See TracChangeset
for help on using the changeset viewer.