OLD | NEW |
(Empty) | |
| 1 /* |
| 2 * Copyright (c) 2008-2010 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 |
| 28 |
| 29 #ifndef CONSTANTS_H |
| 30 #define CONSTANTS_H |
| 31 |
| 32 |
| 33 #include "mpdecimal.h" |
| 34 |
| 35 |
| 36 /* choice of optimized functions */ |
| 37 #if defined(CONFIG_64) |
| 38 /* x64 */ |
| 39 #define MULMOD(a, b) x64_mulmod(a, b, umod) |
| 40 #define MULMOD2C(a0, a1, w) x64_mulmod2c(a0, a1, w, umod) |
| 41 #define MULMOD2(a0, b0, a1, b1) x64_mulmod2(a0, b0, a1, b1, umod) |
| 42 #define POWMOD(base, exp) x64_powmod(base, exp, umod) |
| 43 #define SETMODULUS(modnum) std_setmodulus(modnum, &umod) |
| 44 #define SIZE3_NTT(x0, x1, x2, w3table) std_size3_ntt(x0, x1, x2, w3table, umod
) |
| 45 #elif defined(PPRO) |
| 46 /* PentiumPro (or later) gcc inline asm */ |
| 47 #define MULMOD(a, b) ppro_mulmod(a, b, &dmod, dinvmod) |
| 48 #define MULMOD2C(a0, a1, w) ppro_mulmod2c(a0, a1, w, &dmod, dinvmod) |
| 49 #define MULMOD2(a0, b0, a1, b1) ppro_mulmod2(a0, b0, a1, b1, &dmod, dinvmod) |
| 50 #define POWMOD(base, exp) ppro_powmod(base, exp, &dmod, dinvmod) |
| 51 #define SETMODULUS(modnum) ppro_setmodulus(modnum, &umod, &dmod, dinvmod) |
| 52 #define SIZE3_NTT(x0, x1, x2, w3table) ppro_size3_ntt(x0, x1, x2, w3table, umo
d, &dmod, dinvmod) |
| 53 #else |
| 54 /* ANSI C99 */ |
| 55 #define MULMOD(a, b) std_mulmod(a, b, umod) |
| 56 #define MULMOD2C(a0, a1, w) std_mulmod2c(a0, a1, w, umod) |
| 57 #define MULMOD2(a0, b0, a1, b1) std_mulmod2(a0, b0, a1, b1, umod) |
| 58 #define POWMOD(base, exp) std_powmod(base, exp, umod) |
| 59 #define SETMODULUS(modnum) std_setmodulus(modnum, &umod) |
| 60 #define SIZE3_NTT(x0, x1, x2, w3table) std_size3_ntt(x0, x1, x2, w3table, umod
) |
| 61 #endif |
| 62 |
| 63 /* PentiumPro (or later) gcc inline asm */ |
| 64 extern const float MPD_TWO63; |
| 65 extern const uint32_t mpd_invmoduli[3][3]; |
| 66 |
| 67 enum {P1, P2, P3}; |
| 68 |
| 69 extern const mpd_uint_t mpd_moduli[]; |
| 70 extern const mpd_uint_t mpd_roots[]; |
| 71 extern const mpd_size_t mpd_bits[]; |
| 72 extern const mpd_uint_t mpd_pow10[]; |
| 73 |
| 74 extern const mpd_uint_t INV_P1_MOD_P2; |
| 75 extern const mpd_uint_t INV_P1P2_MOD_P3; |
| 76 extern const mpd_uint_t LH_P1P2; |
| 77 extern const mpd_uint_t UH_P1P2; |
| 78 |
| 79 |
| 80 #endif /* CONSTANTS_H */ |
| 81 |
| 82 |
| 83 |
OLD | NEW |