Annotation of OpenXM_contrib/gmp/mpn/tests/trace.c, Revision 1.1.1.1
1.1 maekawa 1: /* trace.c -- Support for diagnostic traces. */
2:
3: /*
4: Copyright (C) 1999, 2000 Free Software Foundation, Inc.
5:
6: This file is part of the GNU MP Library.
7:
8: The GNU MP Library is free software; you can redistribute it and/or modify
9: it under the terms of the GNU Lesser General Public License as published by
10: the Free Software Foundation; either version 2.1 of the License, or (at your
11: option) any later version.
12:
13: The GNU MP Library is distributed in the hope that it will be useful, but
14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
16: License for more details.
17:
18: You should have received a copy of the GNU Lesser General Public License
19: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
21: MA 02111-1307, USA.
22: */
23:
24: #include <stdio.h>
25: #include "gmp.h"
26: #include "gmp-impl.h"
27: #include "try.h"
28:
29:
30: /* Number base for the various trace printing routines.
31: Set this in main() or with the debugger.
32: If hexadecimal is going to be fed into GNU bc, remember to use -16
33: because bc requires upper case. */
34:
35: int mp_trace_base = 10;
36:
37:
38: /* Print "name=value\n" to stdout for an mpq_t value. */
39: void
40: mpq_trace (const char *name, mpq_srcptr q)
41: {
42: if (name != NULL && name[0] != '\0')
43: printf ("%s=", name);
44:
45: switch (ABS (mp_trace_base)) {
46: case 8: printf ("oct"); break;
47: case 10: break;
48: case 16: printf ("0x"); break;
49: default: printf ("base%d:", ABS (mp_trace_base)); break;
50: }
51:
52: mpq_out_str (stdout, mp_trace_base, q);
53:
54: /* It's not very interesting to know when numbers are unnormalized.
55: mpz's should always be and ought to be checked with ASSERT() if in doubt.
56: mpn's can often be unnormalized without affecting anything. */
57: /* if (!mpq_normalized_p (z)) printf (" (unnorm)"); */
58:
59: printf ("\n");
60: }
61:
62:
63: /* Print "name=value\n" to stdout for an mpz_t value. */
64: void
65: mpz_trace (const char *name, mpz_srcptr z)
66: {
67: mpq_t q;
68: mpz_t one;
69: mp_limb_t one_limb;
70:
71: PTR(one) = &one_limb;
72: SIZ(one) = 1;
73: one_limb = 1;
74:
75: q->_mp_num = *z;
76: q->_mp_den = *one;
77:
78: mpq_trace(name, q);
79: }
80:
81:
82: /* Print "namenum=value\n" to stdout for an mpz_t value.
83: "name" should have a "%d" to get the number. */
84: void
85: mpz_tracen (const char *name, int num, mpz_srcptr z)
86: {
87: if (name != NULL && name[0] != '\0')
88: {
89: printf (name, num);
90: putchar ('=');
91: }
92: mpz_trace (NULL, z);
93: }
94:
95:
96: /* Print "name=value\n" to stdout for an mpn style ptr,size. */
97: void
98: mpn_trace (const char *name, mp_srcptr ptr, mp_size_t size)
99: {
100: mpz_t z;
101: MPN_NORMALIZE (ptr, size);
102: PTR(z) = (mp_ptr) ptr;
103: SIZ(z) = (int) size;
104: mpz_trace (name, z);
105: }
106:
107:
108: /* Print "namenum=value\n" to stdout for an mpn style ptr,size.
109: "name" should have a "%d" to get the number. */
110: void
111: mpn_tracen (const char *name, int num, mp_srcptr ptr, mp_size_t size)
112: {
113: if (name != NULL && name[0] != '\0')
114: {
115: printf (name, num);
116: putchar ('=');
117: }
118: mpn_trace (NULL, ptr, size);
119: }
120:
121:
122: /* Print "namenum=value\n" to stdout for an array of mpn style ptr,size.
123:
124: "a" is an array of pointers, each a[i] is a pointer to "size" many limbs.
125: The formal parameter isn't mp_srcptr because that causes compiler
126: warnings, but the values aren't modified.
127:
128: "name" should have a printf style "%d" to get the array index. */
129:
130: void
131: mpn_tracea (const char *name, const mp_ptr *a, int count, mp_size_t size)
132: {
133: int i;
134: for (i = 0; i < count; i++)
135: mpn_tracen (name, i, a[i], size);
136: }
137:
138:
139: /* Print "value\n" to a file for an mpz_t value. Any previous contents of
140: the file are overwritten, so you need different file names each time this
141: is called.
142:
143: Overwriting the file is a feature, it means you get old data replaced
144: when you run a test program repeatedly. */
145:
146: void
147: mpn_trace_file (const char *filename, mp_srcptr ptr, mp_size_t size)
148: {
149: FILE *fp;
150: mpz_t z;
151:
152: fp = fopen (filename, "w");
153: if (fp == NULL)
154: {
155: perror ("fopen");
156: abort();
157: }
158:
159: MPN_NORMALIZE (ptr, size);
160: PTR(z) = (mp_ptr) ptr;
161: SIZ(z) = (int) size;
162:
163: mpz_out_str (fp, mp_trace_base, z);
164: fprintf (fp, "\n");
165:
166: if (ferror (fp) || fclose (fp) != 0)
167: {
168: printf ("error writing %s\n", filename);
169: abort();
170: }
171: }
172:
173:
174: /* Print "value\n" to a set of files, one file for each element of the given
175: array of mpn style ptr,size. Any previous contents of the files are
176: overwritten, so you need different file names each time this is called.
177: Each file is "filenameN" where N is 0 to count-1.
178:
179: "a" is an array of pointers, each a[i] is a pointer to "size" many limbs.
180: The formal parameter isn't mp_srcptr because that causes compiler
181: warnings, but the values aren't modified.
182:
183: Overwriting the files is a feature, it means you get old data replaced
184: when you run a test program repeatedly. The output style isn't
185: particularly pretty, but at least it gets something out, and you can cat
186: the files into bc, or whatever. */
187:
188: void
189: mpn_tracea_file (const char *filename,
190: const mp_ptr *a, int count, mp_size_t size)
191: {
192: char *s;
193: int i;
194: TMP_DECL (marker);
195:
196: TMP_MARK (marker);
197: s = (char *) TMP_ALLOC (strlen (filename) + 50);
198:
199: for (i = 0; i < count; i++)
200: {
201: sprintf (s, "%s%d", filename, i);
202: mpn_trace_file (s, a[i], size);
203: }
204:
205: TMP_FREE (marker);
206: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to trace.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpn / tests