numbertheory.h 2.4 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576
  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_NUMBERTHEORY_H_
  28. #define LIBMPDEC_NUMBERTHEORY_H_
  29. #include "mpdecimal.h"
  30. #include "constants.h"
  31. /* Internal header file: all symbols have local scope in the DSO */
  32. MPD_PRAGMA(MPD_HIDE_SYMBOLS_START)
  33. /* transform parameters */
  34. struct fnt_params {
  35. int modnum;
  36. mpd_uint_t modulus;
  37. mpd_uint_t kernel;
  38. mpd_uint_t wtable[];
  39. };
  40. mpd_uint_t _mpd_getkernel(mpd_uint_t n, int sign, int modnum);
  41. struct fnt_params *_mpd_init_fnt_params(mpd_size_t n, int sign, int modnum);
  42. void _mpd_init_w3table(mpd_uint_t w3table[3], int sign, int modnum);
  43. #ifdef PPRO
  44. static inline void
  45. ppro_setmodulus(int modnum, mpd_uint_t *umod, double *dmod, uint32_t dinvmod[3])
  46. {
  47. *dmod = *umod = mpd_moduli[modnum];
  48. dinvmod[0] = mpd_invmoduli[modnum][0];
  49. dinvmod[1] = mpd_invmoduli[modnum][1];
  50. dinvmod[2] = mpd_invmoduli[modnum][2];
  51. }
  52. #else
  53. static inline void
  54. std_setmodulus(int modnum, mpd_uint_t *umod)
  55. {
  56. *umod = mpd_moduli[modnum];
  57. }
  58. #endif
  59. MPD_PRAGMA(MPD_HIDE_SYMBOLS_END) /* restore previous scope rules */
  60. #endif /* LIBMPDEC_NUMBERTHEORY_H_ */