Hex Artifact Content
Not logged in

Artifact 0d8de19c7ee3756e580521fadf3aee0e0c941e4e:


0000: 2f 2a 20 67 6d 70 73 74 61 74 2e 68 20 2a 2f 0a  /* gmpstat.h */.
0010: 0a 2f 2a 0a 43 6f 70 79 72 69 67 68 74 20 31 39  ./*.Copyright 19
0020: 39 39 20 46 72 65 65 20 53 6f 66 74 77 61 72 65  99 Free Software
0030: 20 46 6f 75 6e 64 61 74 69 6f 6e 2c 20 49 6e 63   Foundation, Inc
0040: 2e 0a 0a 54 68 69 73 20 66 69 6c 65 20 69 73 20  ...This file is 
0050: 70 61 72 74 20 6f 66 20 74 68 65 20 47 4e 55 20  part of the GNU 
0060: 4d 50 20 4c 69 62 72 61 72 79 2e 0a 0a 54 68 65  MP Library...The
0070: 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79 20   GNU MP Library 
0080: 69 73 20 66 72 65 65 20 73 6f 66 74 77 61 72 65  is free software
0090: 3b 20 79 6f 75 20 63 61 6e 20 72 65 64 69 73 74  ; you can redist
00a0: 72 69 62 75 74 65 20 69 74 20 61 6e 64 2f 6f 72  ribute it and/or
00b0: 20 6d 6f 64 69 66 79 0a 69 74 20 75 6e 64 65 72   modify.it under
00c0: 20 74 68 65 20 74 65 72 6d 73 20 6f 66 20 74 68   the terms of th
00d0: 65 20 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e  e GNU Lesser Gen
00e0: 65 72 61 6c 20 50 75 62 6c 69 63 20 4c 69 63 65  eral Public Lice
00f0: 6e 73 65 20 61 73 20 70 75 62 6c 69 73 68 65 64  nse as published
0100: 20 62 79 0a 74 68 65 20 46 72 65 65 20 53 6f 66   by.the Free Sof
0110: 74 77 61 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e  tware Foundation
0120: 3b 20 65 69 74 68 65 72 20 76 65 72 73 69 6f 6e  ; either version
0130: 20 32 2e 31 20 6f 66 20 74 68 65 20 4c 69 63 65   2.1 of the Lice
0140: 6e 73 65 2c 20 6f 72 20 28 61 74 20 79 6f 75 72  nse, or (at your
0150: 0a 6f 70 74 69 6f 6e 29 20 61 6e 79 20 6c 61 74  .option) any lat
0160: 65 72 20 76 65 72 73 69 6f 6e 2e 0a 0a 54 68 65  er version...The
0170: 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79 20   GNU MP Library 
0180: 69 73 20 64 69 73 74 72 69 62 75 74 65 64 20 69  is distributed i
0190: 6e 20 74 68 65 20 68 6f 70 65 20 74 68 61 74 20  n the hope that 
01a0: 69 74 20 77 69 6c 6c 20 62 65 20 75 73 65 66 75  it will be usefu
01b0: 6c 2c 20 62 75 74 0a 57 49 54 48 4f 55 54 20 41  l, but.WITHOUT A
01c0: 4e 59 20 57 41 52 52 41 4e 54 59 3b 20 77 69 74  NY WARRANTY; wit
01d0: 68 6f 75 74 20 65 76 65 6e 20 74 68 65 20 69 6d  hout even the im
01e0: 70 6c 69 65 64 20 77 61 72 72 61 6e 74 79 20 6f  plied warranty o
01f0: 66 20 4d 45 52 43 48 41 4e 54 41 42 49 4c 49 54  f MERCHANTABILIT
0200: 59 0a 6f 72 20 46 49 54 4e 45 53 53 20 46 4f 52  Y.or FITNESS FOR
0210: 20 41 20 50 41 52 54 49 43 55 4c 41 52 20 50 55   A PARTICULAR PU
0220: 52 50 4f 53 45 2e 20 20 53 65 65 20 74 68 65 20  RPOSE.  See the 
0230: 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e 65 72  GNU Lesser Gener
0240: 61 6c 20 50 75 62 6c 69 63 0a 4c 69 63 65 6e 73  al Public.Licens
0250: 65 20 66 6f 72 20 6d 6f 72 65 20 64 65 74 61 69  e for more detai
0260: 6c 73 2e 0a 0a 59 6f 75 20 73 68 6f 75 6c 64 20  ls...You should 
0270: 68 61 76 65 20 72 65 63 65 69 76 65 64 20 61 20  have received a 
0280: 63 6f 70 79 20 6f 66 20 74 68 65 20 47 4e 55 20  copy of the GNU 
0290: 4c 65 73 73 65 72 20 47 65 6e 65 72 61 6c 20 50  Lesser General P
02a0: 75 62 6c 69 63 20 4c 69 63 65 6e 73 65 0a 61 6c  ublic License.al
02b0: 6f 6e 67 20 77 69 74 68 20 74 68 65 20 47 4e 55  ong with the GNU
02c0: 20 4d 50 20 4c 69 62 72 61 72 79 3b 20 73 65 65   MP Library; see
02d0: 20 74 68 65 20 66 69 6c 65 20 43 4f 50 59 49 4e   the file COPYIN
02e0: 47 2e 4c 49 42 2e 20 20 49 66 20 6e 6f 74 2c 20  G.LIB.  If not, 
02f0: 77 72 69 74 65 20 74 6f 0a 74 68 65 20 46 72 65  write to.the Fre
0300: 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64  e Software Found
0310: 61 74 69 6f 6e 2c 20 49 6e 63 2e 2c 20 35 39 20  ation, Inc., 59 
0320: 54 65 6d 70 6c 65 20 50 6c 61 63 65 20 2d 20 53  Temple Place - S
0330: 75 69 74 65 20 33 33 30 2c 20 42 6f 73 74 6f 6e  uite 330, Boston
0340: 2c 0a 4d 41 20 30 32 31 31 31 2d 31 33 30 37 2c  ,.MA 02111-1307,
0350: 20 55 53 41 2e 0a 2a 2f 0a 0a 2f 2a 20 54 68 69   USA..*/../* Thi
0360: 73 20 66 69 6c 65 20 72 65 71 75 69 72 65 73 20  s file requires 
0370: 74 68 65 20 66 6f 6c 6c 6f 77 69 6e 67 20 68 65  the following he
0380: 61 64 65 72 20 66 69 6c 65 73 3a 20 67 6d 70 2e  ader files: gmp.
0390: 68 20 2a 2f 0a 0a 23 69 66 6e 64 65 66 09 5f 5f  h */..#ifndef.__
03a0: 47 4d 50 53 54 41 54 5f 48 5f 5f 0a 23 64 65 66  GMPSTAT_H__.#def
03b0: 69 6e 65 09 5f 5f 47 4d 50 53 54 41 54 5f 48 5f  ine.__GMPSTAT_H_
03c0: 5f 0a 0a 2f 2a 20 47 6c 6f 62 61 6c 20 64 65 62  _../* Global deb
03d0: 75 67 20 66 6c 61 67 2e 20 20 46 49 58 4d 45 3a  ug flag.  FIXME:
03e0: 20 52 65 6d 6f 76 65 2e 20 2a 2f 0a 65 78 74 65   Remove. */.exte
03f0: 72 6e 20 69 6e 74 20 67 5f 64 65 62 75 67 3b 0a  rn int g_debug;.
0400: 23 64 65 66 69 6e 65 20 44 45 42 55 47 5f 31 20  #define DEBUG_1 
0410: 30 0a 23 64 65 66 69 6e 65 20 44 45 42 55 47 5f  0.#define DEBUG_
0420: 32 20 31 0a 0a 2f 2a 20 4d 61 78 20 6e 75 6d 62  2 1../* Max numb
0430: 65 72 20 6f 66 20 64 69 6d 65 6e 73 69 6f 6e 73  er of dimensions
0440: 20 69 6e 20 73 70 65 63 74 72 61 6c 20 74 65 73   in spectral tes
0450: 74 2e 20 20 46 49 58 4d 45 3a 20 4d 61 6b 77 20  t.  FIXME: Makw 
0460: 64 79 6e 61 6d 69 63 2e 20 2a 2f 0a 23 64 65 66  dynamic. */.#def
0470: 69 6e 65 20 47 4d 50 5f 53 50 45 43 54 5f 4d 41  ine GMP_SPECT_MA
0480: 58 54 20 31 30 0a 0a 76 6f 69 64 0a 6d 70 66 5f  XT 10..void.mpf_
0490: 66 72 65 71 74 20 28 6d 70 66 5f 74 20 4b 70 2c  freqt (mpf_t Kp,
04a0: 0a 09 20 20 20 6d 70 66 5f 74 20 4b 6d 2c 0a 09  ..   mpf_t Km,..
04b0: 20 20 20 6d 70 66 5f 74 20 58 5b 5d 2c 0a 09 20     mpf_t X[],.. 
04c0: 20 20 63 6f 6e 73 74 20 75 6e 73 69 67 6e 65 64    const unsigned
04d0: 20 6c 6f 6e 67 20 69 6e 74 20 6e 29 3b 0a 75 6e   long int n);.un
04e0: 73 69 67 6e 65 64 20 6c 6f 6e 67 20 69 6e 74 0a  signed long int.
04f0: 6d 70 7a 5f 66 72 65 71 74 20 28 6d 70 66 5f 74  mpz_freqt (mpf_t
0500: 20 56 2c 0a 09 20 20 20 6d 70 7a 5f 74 20 58 5b   V,..   mpz_t X[
0510: 5d 2c 0a 09 20 20 20 75 6e 73 69 67 6e 65 64 20  ],..   unsigned 
0520: 69 6e 74 20 69 6d 61 78 2c 0a 09 20 20 20 63 6f  int imax,..   co
0530: 6e 73 74 20 75 6e 73 69 67 6e 65 64 20 6c 6f 6e  nst unsigned lon
0540: 67 20 69 6e 74 20 6e 29 3b 0a 0a 2f 2a 20 4c 6f  g int n);../* Lo
0550: 77 20 6c 65 76 65 6c 20 66 75 6e 63 74 69 6f 6e  w level function
0560: 73 2e 20 2a 2f 0a 76 6f 69 64 0a 6b 73 20 28 6d  s. */.void.ks (m
0570: 70 66 5f 74 20 4b 70 2c 0a 20 20 20 20 6d 70 66  pf_t Kp,.    mpf
0580: 5f 74 20 4b 6d 2c 0a 20 20 20 20 6d 70 66 5f 74  _t Km,.    mpf_t
0590: 20 58 5b 5d 2c 0a 20 20 20 20 76 6f 69 64 20 28   X[],.    void (
05a0: 50 29 20 28 6d 70 66 5f 74 2c 20 6d 70 66 5f 74  P) (mpf_t, mpf_t
05b0: 29 2c 0a 20 20 20 20 63 6f 6e 73 74 20 75 6e 73  ),.    const uns
05c0: 69 67 6e 65 64 20 6c 6f 6e 67 20 69 6e 74 20 6e  igned long int n
05d0: 29 3b 0a 0a 76 6f 69 64 0a 6b 73 5f 74 61 62 6c  );..void.ks_tabl
05e0: 65 20 28 6d 70 66 5f 74 20 70 2c 20 6d 70 66 5f  e (mpf_t p, mpf_
05f0: 74 20 76 61 6c 2c 20 63 6f 6e 73 74 20 75 6e 73  t val, const uns
0600: 69 67 6e 65 64 20 69 6e 74 20 6e 29 3b 0a 0a 76  igned int n);..v
0610: 6f 69 64 0a 78 32 5f 74 61 62 6c 65 20 28 64 6f  oid.x2_table (do
0620: 75 62 6c 65 20 74 5b 5d 2c 0a 09 20 20 75 6e 73  uble t[],..  uns
0630: 69 67 6e 65 64 20 69 6e 74 20 76 29 3b 0a 0a 76  igned int v);..v
0640: 6f 69 64 0a 73 70 65 63 74 72 61 6c 5f 74 65 73  oid.spectral_tes
0650: 74 20 28 6d 70 66 5f 74 20 72 6f 70 5b 5d 2c 20  t (mpf_t rop[], 
0660: 75 6e 73 69 67 6e 65 64 20 69 6e 74 20 54 2c 20  unsigned int T, 
0670: 6d 70 7a 5f 74 20 61 2c 20 6d 70 7a 5f 74 20 6d  mpz_t a, mpz_t m
0680: 29 3b 0a 76 6f 69 64 0a 76 7a 5f 64 6f 74 20 28  );.void.vz_dot (
0690: 6d 70 7a 5f 74 20 72 6f 70 2c 20 6d 70 7a 5f 74  mpz_t rop, mpz_t
06a0: 20 56 31 5b 5d 2c 20 6d 70 7a 5f 74 20 56 32 5b   V1[], mpz_t V2[
06b0: 5d 2c 20 75 6e 73 69 67 6e 65 64 20 69 6e 74 20  ], unsigned int 
06c0: 6e 29 3b 0a 76 6f 69 64 0a 66 5f 66 6c 6f 6f 72  n);.void.f_floor
06d0: 20 28 6d 70 66 5f 74 20 72 6f 70 2c 20 6d 70 66   (mpf_t rop, mpf
06e0: 5f 74 20 6f 70 29 3b 0a 0a 76 6f 69 64 0a 6d 65  _t op);..void.me
06f0: 72 69 74 20 28 6d 70 66 5f 74 20 72 6f 70 2c 20  rit (mpf_t rop, 
0700: 75 6e 73 69 67 6e 65 64 20 69 6e 74 20 74 2c 20  unsigned int t, 
0710: 6d 70 66 5f 74 20 76 2c 20 6d 70 7a 5f 74 20 6d  mpf_t v, mpz_t m
0720: 29 3b 0a 64 6f 75 62 6c 65 0a 6d 65 72 69 74 5f  );.double.merit_
0730: 75 20 28 75 6e 73 69 67 6e 65 64 20 69 6e 74 20  u (unsigned int 
0740: 74 2c 20 6d 70 66 5f 74 20 76 2c 20 6d 70 7a 5f  t, mpf_t v, mpz_
0750: 74 20 6d 29 3b 0a 0a 2f 2a 20 46 72 6f 6d 20 73  t m);../* From s
0760: 65 70 61 72 61 74 65 20 73 6f 75 72 63 65 20 66  eparate source f
0770: 69 6c 65 73 3a 20 2a 2f 0a 76 6f 69 64 20 7a 64  iles: */.void zd
0780: 69 76 5f 72 6f 75 6e 64 20 28 6d 70 7a 5f 74 20  iv_round (mpz_t 
0790: 72 6f 70 2c 20 6d 70 7a 5f 74 20 6e 2c 20 6d 70  rop, mpz_t n, mp
07a0: 7a 5f 74 20 64 29 3b 0a 0a 23 65 6e 64 69 66 20  z_t d);..#endif 
07b0: 2f 2a 20 21 5f 5f 47 4d 50 53 54 41 54 5f 48 5f  /* !__GMPSTAT_H_
07c0: 5f 20 2a 2f 0a                                   _ */.