constants.h 3.4 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889
  1. /*
  2. * Copyright (c) 2008-2020 Stefan Krah. All rights reserved.
  3. *
  4. * Redistribution and use in source and binary forms, with or without
  5. * modification, are permitted provided that the following conditions
  6. * are met:
  7. *
  8. * 1. Redistributions of source code must retain the above copyright
  9. * notice, this list of conditions and the following disclaimer.
  10. *
  11. * 2. Redistributions in binary form must reproduce the above copyright
  12. * notice, this list of conditions and the following disclaimer in the
  13. * documentation and/or other materials provided with the distribution.
  14. *
  15. * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
  16. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  17. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18. * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  19. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  20. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  21. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  22. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  23. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  24. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  25. * SUCH DAMAGE.
  26. */
  27. #ifndef LIBMPDEC_CONSTANTS_H_
  28. #define LIBMPDEC_CONSTANTS_H_
  29. #include "mpdecimal.h"
  30. #include <stdint.h>
  31. /* Internal header file: all symbols have local scope in the DSO */
  32. MPD_PRAGMA(MPD_HIDE_SYMBOLS_START)
  33. /* choice of optimized functions */
  34. #if defined(CONFIG_64)
  35. /* x64 */
  36. #define MULMOD(a, b) x64_mulmod(a, b, umod)
  37. #define MULMOD2C(a0, a1, w) x64_mulmod2c(a0, a1, w, umod)
  38. #define MULMOD2(a0, b0, a1, b1) x64_mulmod2(a0, b0, a1, b1, umod)
  39. #define POWMOD(base, exp) x64_powmod(base, exp, umod)
  40. #define SETMODULUS(modnum) std_setmodulus(modnum, &umod)
  41. #define SIZE3_NTT(x0, x1, x2, w3table) std_size3_ntt(x0, x1, x2, w3table, umod)
  42. #elif defined(PPRO)
  43. /* PentiumPro (or later) gcc inline asm */
  44. #define MULMOD(a, b) ppro_mulmod(a, b, &dmod, dinvmod)
  45. #define MULMOD2C(a0, a1, w) ppro_mulmod2c(a0, a1, w, &dmod, dinvmod)
  46. #define MULMOD2(a0, b0, a1, b1) ppro_mulmod2(a0, b0, a1, b1, &dmod, dinvmod)
  47. #define POWMOD(base, exp) ppro_powmod(base, exp, &dmod, dinvmod)
  48. #define SETMODULUS(modnum) ppro_setmodulus(modnum, &umod, &dmod, dinvmod)
  49. #define SIZE3_NTT(x0, x1, x2, w3table) ppro_size3_ntt(x0, x1, x2, w3table, umod, &dmod, dinvmod)
  50. #else
  51. /* ANSI C99 */
  52. #define MULMOD(a, b) std_mulmod(a, b, umod)
  53. #define MULMOD2C(a0, a1, w) std_mulmod2c(a0, a1, w, umod)
  54. #define MULMOD2(a0, b0, a1, b1) std_mulmod2(a0, b0, a1, b1, umod)
  55. #define POWMOD(base, exp) std_powmod(base, exp, umod)
  56. #define SETMODULUS(modnum) std_setmodulus(modnum, &umod)
  57. #define SIZE3_NTT(x0, x1, x2, w3table) std_size3_ntt(x0, x1, x2, w3table, umod)
  58. #endif
  59. /* PentiumPro (or later) gcc inline asm */
  60. extern const float MPD_TWO63;
  61. extern const uint32_t mpd_invmoduli[3][3];
  62. enum {P1, P2, P3};
  63. extern const mpd_uint_t mpd_moduli[];
  64. extern const mpd_uint_t mpd_roots[];
  65. extern const mpd_size_t mpd_bits[];
  66. extern const mpd_uint_t mpd_pow10[];
  67. extern const mpd_uint_t INV_P1_MOD_P2;
  68. extern const mpd_uint_t INV_P1P2_MOD_P3;
  69. extern const mpd_uint_t LH_P1P2;
  70. extern const mpd_uint_t UH_P1P2;
  71. MPD_PRAGMA(MPD_HIDE_SYMBOLS_END) /* restore previous scope rules */
  72. #endif /* LIBMPDEC_CONSTANTS_H_ */