My Project
Loading...
Searching...
No Matches
misc_ip.cc
Go to the documentation of this file.
1/*****************************************************************************\
2 * Computer Algebra System SINGULAR
3\*****************************************************************************/
4/** @file misc_ip.cc
5 *
6 * This file provides miscellaneous functionality.
7 *
8 * For more general information, see the documentation in misc_ip.h.
9 *
10 **/
11/*****************************************************************************/
12
13// include header files
14#define PLURAL_INTERNAL_DECLARATIONS 1
15
16#include "kernel/mod2.h"
17#include "misc/sirandom.h"
18#include "omalloc/omalloc.h"
19#include "misc/mylimits.h"
20#include "reporter/si_signals.h"
21#include "factory/factory.h"
22#include "coeffs/si_gmp.h"
23#include "coeffs/coeffs.h"
24#include "coeffs/flintcf_Q.h"
25#include "coeffs/flintcf_Qrat.h"
26#include "coeffs/flintcf_Zn.h"
27#include "coeffs/rmodulon.h"
30#include "polys/nc/gb_hack.h"
31
32#ifdef HAVE_SIMPLEIPC
34#endif
35
36#include "misc_ip.h"
37#include "ipid.h"
38#include "feOpt.h"
39#include "links/silink.h"
40#include "mod_lib.h"
41#include "misc/distrib.h"
42
43#include "misc/options.h"
44#include "misc/intvec.h"
45
48
53
54#include "subexpr.h"
55#include "cntrlc.h"
56#include "ipshell.h"
57
58#include "fehelp.h"
59
60#ifdef HAVE_READLINE
61 #ifdef READLINE_READLINE_H_OK
62 #include <readline/readline.h>
63 #endif
64 #ifndef RL_VERSION_MAJOR
65 #define RL_VERSION_MAJOR 0
66 #endif
67#endif
68
69#ifdef HAVE_NTL
70#include <NTL/version.h>
71#endif
72
73
74void setListEntry(lists L, int index, mpz_t n)
75{ /* assumes n > 0 */
76 /* try to fit nn into an int: */
77 if (mpz_size1(n)<=1)
78 {
79 int ui=(int)mpz_get_si(n);
80 if ((((ui<<3)>>3)==ui)
81 && (mpz_cmp_si(n,(long)ui)==0))
82 {
83 L->m[index].rtyp = INT_CMD; L->m[index].data = (void*)(long)ui;
84 return;
85 }
86 }
87 number nn = mpz2number(n, coeffs_BIGINT);
88 L->m[index].rtyp = BIGINT_CMD; L->m[index].data = (void*)nn;
89}
90
91void setListEntry_ui(lists L, int index, unsigned long ui)
92{ /* assumes n > 0 */
93 /* try to fit nn into an int: */
94 int i=(int)ui;
95 if ((((unsigned long)i)==ui) && (((i<<3)>>3)==i))
96 {
97 L->m[index].rtyp = INT_CMD; L->m[index].data = (void*)(long)i;
98 }
99 else
100 {
101 number nn = n_Init(ui, coeffs_BIGINT);
102 L->m[index].rtyp = BIGINT_CMD; L->m[index].data = (void*)nn;
103 }
104}
105
106/* Factoring with Pollard's rho method. stolen from GMP/demos */
107STATIC_VAR unsigned add[] = {4, 2, 4, 2, 4, 6, 2, 6};
108
109static int factor_using_division (mpz_t t, unsigned int limit,lists primes, int *multiplicities,int &index, unsigned long bound)
110{
111 mpz_t q, r;
112 unsigned long int f;
113 int ai;
114 unsigned *addv = add;
115 unsigned int failures;
116 int bound_not_reached=1;
117
118 mpz_init (q);
119 mpz_init (r);
120
121 f = mpz_scan1 (t, 0);
122 mpz_div_2exp (t, t, f);
123 if (f>0)
124 {
126 multiplicities[index++] = f;
127 }
128
129 f=0;
130 loop
131 {
132 mpz_tdiv_qr_ui (q, r, t, 3);
133 if (mpz_sgn1 (r) != 0)
134 break;
135 mpz_set (t, q);
136 f++;
137 }
138 if (f>0)
139 {
141 multiplicities[index++] = f;
142 }
143 f=0;
144 loop
145 {
146 mpz_tdiv_qr_ui (q, r, t, 5);
147 if (mpz_sgn1 (r) != 0)
148 break;
149 mpz_set (t, q);
150 f++;
151 }
152 if (f>0)
153 {
155 multiplicities[index++] = f;
156 }
157
158 failures = 0;
159 f = 7;
160 ai = 0;
161 unsigned long last_f=0;
162 while (mpz_cmp_ui (t, 1) != 0)
163 {
164 mpz_tdiv_qr_ui (q, r, t, f);
165 if (mpz_sgn1 (r) != 0)
166 {
167 f += addv[ai];
168 if (mpz_cmp_ui (t, f) < 0)
169 break;
170 ai = (ai + 1) & 7;
171 failures++;
172 if (failures > limit)
173 break;
174 if ((bound!=0) && (f>bound))
175 {
176 bound_not_reached=0;
177 break;
178 }
179 }
180 else
181 {
182 mpz_swap (t, q);
183 if (f!=last_f)
184 {
186 multiplicities[index]++;
187 index++;
188 }
189 else
190 {
191 multiplicities[index-1]++;
192 }
193 last_f=f;
194 failures = 0;
195 }
196 }
197
198 mpz_clear (q);
199 mpz_clear (r);
200 //printf("bound=%d,f=%d,failures=%d, reached=%d\n",bound,f,failures,bound_not_reached);
201 return bound_not_reached;
202}
203
204static void factor_using_pollard_rho (mpz_t n, unsigned long a, lists primes, int * multiplicities,int &index)
205{
206 mpz_t x, x1, y, P;
207 mpz_t t1, t2;
208 mpz_t last_f;
209 unsigned long long k, l, i;
210
211 mpz_init (t1);
212 mpz_init (t2);
213 mpz_init_set_ui (last_f, 0);
214 mpz_init_set_ui (y, 2);
215 mpz_init_set_ui (x, 2);
216 mpz_init_set_ui (x1, 2);
217 mpz_init_set_ui (P, 1);
218 k = 1;
219 l = 1;
220
221 while (mpz_cmp_ui (n, 1) != 0)
222 {
223 loop
224 {
225 do
226 {
227 mpz_mul (t1, x, x);
228 mpz_mod (x, t1, n);
229 mpz_add_ui (x, x, a);
230 mpz_sub (t1, x1, x);
231 mpz_mul (t2, P, t1);
232 mpz_mod (P, t2, n);
233
234 if (k % 32 == 1)
235 {
236 mpz_gcd (t1, P, n);
237 if (mpz_cmp_ui (t1, 1) != 0)
238 goto factor_found;
239 mpz_set (y, x);
240 }
241 }
242 while (--k != 0);
243
244 mpz_gcd (t1, P, n);
245 if (mpz_cmp_ui (t1, 1) != 0)
246 goto factor_found;
247
248 mpz_set (x1, x);
249 k = l;
250 l = 2 * l;
251 for (i = 0; i < k; i++)
252 {
253 mpz_mul (t1, x, x);
254 mpz_mod (x, t1, n);
255 mpz_add_ui (x, x, a);
256 }
257 mpz_set (y, x);
258 }
259
260 factor_found:
261 do
262 {
263 mpz_mul (t1, y, y);
264 mpz_mod (y, t1, n);
265 mpz_add_ui (y, y, a);
266 mpz_sub (t1, x1, y);
267 mpz_gcd (t1, t1, n);
268 }
269 while (mpz_cmp_ui (t1, 1) == 0);
270
271 mpz_divexact (n, n, t1); /* divide by t1, before t1 is overwritten */
272
273 if (!mpz_probab_prime_p (t1, 10))
274 {
275 do
276 {
277 mp_limb_t a_limb;
278 mpn_random (&a_limb, (mp_size_t) 1);
279 a = a_limb;
280 }
281 while (a == 0);
282
283 factor_using_pollard_rho (t1, a, primes,multiplicities,index);
284 }
285 else
286 {
287 if (mpz_cmp(t1,last_f)==0)
288 {
289 multiplicities[index-1]++;
290 }
291 else
292 {
293 mpz_set(last_f,t1);
295 multiplicities[index++] = 1;
296 }
297 }
298 mpz_mod (x, x, n);
299 mpz_mod (x1, x1, n);
300 mpz_mod (y, y, n);
301 if (mpz_probab_prime_p (n, 10))
302 {
303 if (mpz_cmp(n,last_f)==0)
304 {
305 multiplicities[index-1]++;
306 }
307 else
308 {
309 mpz_set(last_f,n);
311 multiplicities[index++] = 1;
312 }
313 mpz_set_ui(n,1);
314 break;
315 }
316 }
317
318 mpz_clear (P);
319 mpz_clear (t2);
320 mpz_clear (t1);
321 mpz_clear (x1);
322 mpz_clear (x);
323 mpz_clear (y);
324 mpz_clear (last_f);
325}
326
327static void factor_gmp (mpz_t t,lists primes,int *multiplicities,int &index,unsigned long bound)
328{
329 unsigned int division_limit;
330
331 if (mpz_sgn (t) == 0)
332 return;
333
334 /* Set the trial division limit according the size of t. */
335 division_limit = mpz_sizeinbase (t, 2);
336 if (division_limit > 1000)
337 division_limit = 1000 * 1000;
338 else
339 division_limit = division_limit * division_limit;
340
341 if (factor_using_division (t, division_limit,primes,multiplicities,index,bound))
342 {
343 if (mpz_cmp_ui (t, 1) != 0)
344 {
345 if (mpz_probab_prime_p (t, 10))
346 {
348 multiplicities[index++] = 1;
349 mpz_set_ui(t,1);
350 }
351 else
352 factor_using_pollard_rho (t, 1L, primes,multiplicities,index);
353 }
354 }
355}
356/* n and pBound are assumed to be bigint numbers */
357lists primeFactorisation(const number n, const int pBound)
358{
359 int i;
360 int index=0;
361 mpz_t nn; number2mpz(n, coeffs_BIGINT, nn);
362 lists primes = (lists)omAllocBin(slists_bin); primes->Init(1000);
363 int* multiplicities = (int*)omAlloc0(1000*sizeof(int));
364 int positive=1;
365
366 if (!n_IsZero(n, coeffs_BIGINT))
367 {
369 {
370 positive=-1;
371 mpz_neg(nn,nn);
372 }
373 factor_gmp(nn,primes,multiplicities,index,pBound);
374 }
375
376 lists primesL = (lists)omAllocBin(slists_bin);
377 primesL->Init(index);
378 for (i = 0; i < index; i++)
379 {
380 primesL->m[i].rtyp = primes->m[i].rtyp;
381 primesL->m[i].data = primes->m[i].data;
382 primes->m[i].rtyp=0;
383 primes->m[i].data=NULL;
384 }
385 primes->Clean(NULL);
386
387 lists multiplicitiesL = (lists)omAllocBin(slists_bin);
388 multiplicitiesL->Init(index);
389 for (i = 0; i < index; i++)
390 {
391 multiplicitiesL->m[i].rtyp = INT_CMD;
392 multiplicitiesL->m[i].data = (void*)(long)multiplicities[i];
393 }
394 omFree(multiplicities);
395
397 L->Init(3);
398 if (positive==-1) mpz_neg(nn,nn);
399 L->m[0].rtyp = LIST_CMD; L->m[0].data = (void*)primesL;
400 L->m[1].rtyp = LIST_CMD; L->m[1].data = (void*)multiplicitiesL;
401 setListEntry(L, 2, nn);
402
403 mpz_clear(nn);
404
405 return L;
406}
407
408//#ifdef HAVE_LIBPARSER
409//# include "libparse.h"
410//#endif /* HAVE_LIBPARSER */
411
412
413/*2
414* the renice routine for very large jobs
415* works only on unix machines,
416* testet on : linux, HP 9.0
417*
418*#include <sys/times.h>
419*#include <sys/resource.h>
420*extern "C" int setpriority(int,int,int);
421*void very_nice()
422*{
423*#ifndef NO_SETPRIORITY
424* setpriority(PRIO_PROCESS,0,19);
425*#endif
426* sleep(10);
427*}
428*/
429
430void singular_example(char *str)
431{
432 assume(str!=NULL);
433 char *s=str;
434 while (*s==' ') s++;
435 char *ss=s;
436 while (*ss!='\0') ss++;
437 while (*ss<=' ')
438 {
439 *ss='\0';
440 ss--;
441 }
442 idhdl h=IDROOT->get_level(s,0);
443 if ((h!=NULL) && (IDTYP(h)==PROC_CMD))
444 {
445 char *lib=iiGetLibName(IDPROC(h));
446 if((lib!=NULL)&&(*lib!='\0'))
447 {
448 Print("// proc %s from lib %s\n",s,lib);
450 if (s!=NULL)
451 {
452 if (strlen(s)>5)
453 {
454 iiEStart(s,IDPROC(h));
455 omFree((ADDRESS)s);
456 return;
457 }
458 else omFree((ADDRESS)s);
459 }
460 }
461 }
462 else
463 {
464 char sing_file[MAXPATHLEN];
465 FILE *fd=NULL;
466 char *res_m=feResource('m', 0);
467 if (res_m!=NULL)
468 {
469 sprintf(sing_file, "%s/%s.sing", res_m, s);
470 fd = feFopen(sing_file, "r");
471 }
472 if (fd != NULL)
473 {
474
475 int old_echo = si_echo;
476 int length, got;
477 char* s;
478
479 fseek(fd, 0, SEEK_END);
480 length = ftell(fd);
481 fseek(fd, 0, SEEK_SET);
482 s = (char*) omAlloc((length+20)*sizeof(char));
483 got = fread(s, sizeof(char), length, fd);
484 fclose(fd);
485 if (got != length)
486 {
487 Werror("Error while reading file %s", sing_file);
488 }
489 else
490 {
491 s[length] = '\0';
492 strcat(s, "\n;return();\n\n");
493 si_echo = 2;
494 iiEStart(s, NULL);
495 si_echo = old_echo;
496 }
497 omFree(s);
498 }
499 else
500 {
501 Werror("no example for %s", str);
502 }
503 }
504}
505
506
508{
509 {"prot", Sy_bit(OPT_PROT), ~Sy_bit(OPT_PROT) },
510 {"redSB", Sy_bit(OPT_REDSB), ~Sy_bit(OPT_REDSB) },
511 {"notBuckets", Sy_bit(OPT_NOT_BUCKETS), ~Sy_bit(OPT_NOT_BUCKETS) },
512 {"notSugar", Sy_bit(OPT_NOT_SUGAR), ~Sy_bit(OPT_NOT_SUGAR) },
513 {"interrupt", Sy_bit(OPT_INTERRUPT), ~Sy_bit(OPT_INTERRUPT) },
514 {"sugarCrit", Sy_bit(OPT_SUGARCRIT), ~Sy_bit(OPT_SUGARCRIT) },
515 {"teach", Sy_bit(OPT_DEBUG), ~Sy_bit(OPT_DEBUG) },
516 {"notSyzMinim", Sy_bit(OPT_NO_SYZ_MINIM), ~Sy_bit(OPT_NO_SYZ_MINIM) },
517 /* 9 return SB in syz, quotient, intersect, modulo */
518 {"returnSB", Sy_bit(OPT_RETURN_SB), ~Sy_bit(OPT_RETURN_SB) },
519 {"fastHC", Sy_bit(OPT_FASTHC), ~Sy_bit(OPT_FASTHC) },
520 /* 11-19 sort in L/T */
521 {"staircaseBound",Sy_bit(OPT_STAIRCASEBOUND),~Sy_bit(OPT_STAIRCASEBOUND) },
522 {"multBound", Sy_bit(OPT_MULTBOUND), ~Sy_bit(OPT_MULTBOUND) },
523 {"degBound", Sy_bit(OPT_DEGBOUND), ~Sy_bit(OPT_DEGBOUND) },
524 {"redTailSyz", Sy_bit(OPT_REDTAIL_SYZ), ~Sy_bit(OPT_REDTAIL_SYZ) },
525 /* 25 no redTail(p)/redTail(s) */
526 {"redTail", Sy_bit(OPT_REDTAIL), ~Sy_bit(OPT_REDTAIL) },
527 {"redThrough", Sy_bit(OPT_REDTHROUGH), ~Sy_bit(OPT_REDTHROUGH) },
528 {"lazy", Sy_bit(OPT_OLDSTD), ~Sy_bit(OPT_OLDSTD) },
529 {"intStrategy", Sy_bit(OPT_INTSTRATEGY), ~Sy_bit(OPT_INTSTRATEGY) },
530 {"infRedTail", Sy_bit(OPT_INFREDTAIL), ~Sy_bit(OPT_INFREDTAIL) },
531 /* 30: use not regularity for syz */
532 {"notRegularity",Sy_bit(OPT_NOTREGULARITY), ~Sy_bit(OPT_NOTREGULARITY) },
533 {"weightM", Sy_bit(OPT_WEIGHTM), ~Sy_bit(OPT_WEIGHTM) },
534/*special for "none" and also end marker for showOption:*/
535 {"ne", 0, 0 }
536};
537
539{
540 {"assign_none",Sy_bit(V_ASSIGN_NONE),~Sy_bit(V_ASSIGN_NONE)},
541 {"mem", Sy_bit(V_SHOW_MEM), ~Sy_bit(V_SHOW_MEM) },
542 {"yacc", Sy_bit(V_YACC), ~Sy_bit(V_YACC) },
543 {"redefine", Sy_bit(V_REDEFINE), ~Sy_bit(V_REDEFINE) },
544 {"reading", Sy_bit(V_READING), ~Sy_bit(V_READING) },
545 {"loadLib", Sy_bit(V_LOAD_LIB), ~Sy_bit(V_LOAD_LIB) },
546 {"debugLib", Sy_bit(V_DEBUG_LIB), ~Sy_bit(V_DEBUG_LIB) },
547 {"loadProc", Sy_bit(V_LOAD_PROC), ~Sy_bit(V_LOAD_PROC) },
548 {"defRes", Sy_bit(V_DEF_RES), ~Sy_bit(V_DEF_RES) },
549 {"usage", Sy_bit(V_SHOW_USE), ~Sy_bit(V_SHOW_USE) },
550 {"Imap", Sy_bit(V_IMAP), ~Sy_bit(V_IMAP) },
551 {"prompt", Sy_bit(V_PROMPT), ~Sy_bit(V_PROMPT) },
552 {"length", Sy_bit(V_LENGTH), ~Sy_bit(V_LENGTH) },
553 {"notWarnSB",Sy_bit(V_NSB), ~Sy_bit(V_NSB) },
554 {"contentSB",Sy_bit(V_CONTENTSB), ~Sy_bit(V_CONTENTSB) },
555 {"cancelunit",Sy_bit(V_CANCELUNIT),~Sy_bit(V_CANCELUNIT)},
556 {"modpsolve",Sy_bit(V_MODPSOLVSB),~Sy_bit(V_MODPSOLVSB)},
557 {"geometricSB",Sy_bit(V_UPTORADICAL),~Sy_bit(V_UPTORADICAL)},
558 {"findMonomials",Sy_bit(V_FINDMONOM),~Sy_bit(V_FINDMONOM)},
559 {"coefStrat",Sy_bit(V_COEFSTRAT), ~Sy_bit(V_COEFSTRAT)},
560 {"qringNF", Sy_bit(V_QRING), ~Sy_bit(V_QRING)},
561 {"warn", Sy_bit(V_ALLWARN), ~Sy_bit(V_ALLWARN)},
562 {"intersectSyz",Sy_bit(V_INTERSECT_SYZ), ~Sy_bit(V_INTERSECT_SYZ)},
563 {"intersectElim",Sy_bit(V_INTERSECT_ELIM), ~Sy_bit(V_INTERSECT_ELIM)},
564/*special for "none" and also end marker for showOption:*/
565 {"ne", 0, 0 }
566};
567
569{
570 const char *n;
571 do
572 {
573 if (v->Typ()==STRING_CMD)
574 {
575 n=(const char *)v->CopyD(STRING_CMD);
576 }
577 else
578 {
579 if (v->name==NULL)
580 return TRUE;
581 if (v->rtyp==0)
582 {
583 n=v->name;
584 v->name=NULL;
585 }
586 else
587 {
588 n=omStrDup(v->name);
589 }
590 }
591
592 int i;
593
594 if(strcmp(n,"get")==0)
595 {
596 intvec *w=new intvec(2);
597 (*w)[0]=si_opt_1;
598 (*w)[1]=si_opt_2;
599 res->rtyp=INTVEC_CMD;
600 res->data=(void *)w;
601 goto okay;
602 }
603 if(strcmp(n,"set")==0)
604 {
605 if((v->next!=NULL)
606 &&(v->next->Typ()==INTVEC_CMD))
607 {
608 v=v->next;
609 intvec *w=(intvec*)v->Data();
610 si_opt_1=(*w)[0];
611 si_opt_2=(*w)[1];
612#if 0
616 ) {
617 si_opt_1 &=~Sy_bit(OPT_INTSTRATEGY);
618 }
619#endif
620 goto okay;
621 }
622 }
623 if(strcmp(n,"none")==0)
624 {
625 si_opt_1=0;
626 si_opt_2=0;
627 goto okay;
628 }
629 for (i=0; (i==0) || (optionStruct[i-1].setval!=0); i++)
630 {
631 if (strcmp(n,optionStruct[i].name)==0)
632 {
634 {
636 // optOldStd disables redthrough
639 }
640 else
641 WarnS("cannot set option");
642#if 0
646 ) {
647 test &=~Sy_bit(OPT_INTSTRATEGY);
648 }
649#endif
650 goto okay;
651 }
652 else if ((strncmp(n,"no",2)==0)
653 && (strcmp(n+2,optionStruct[i].name)==0))
654 {
656 {
658 }
659 else
660 WarnS("cannot clear option");
661 goto okay;
662 }
663 }
664 for (i=0; (i==0) || (verboseStruct[i-1].setval!=0); i++)
665 {
666 if (strcmp(n,verboseStruct[i].name)==0)
667 {
669 #ifdef YYDEBUG
670 #if YYDEBUG
671 /*debugging the bison grammar --> grammar.cc*/
673 if (BVERBOSE(V_YACC)) yydebug=1;
674 else yydebug=0;
675 #endif
676 #endif
677 goto okay;
678 }
679 else if ((strncmp(n,"no",2)==0)
680 && (strcmp(n+2,verboseStruct[i].name)==0))
681 {
683 #ifdef YYDEBUG
684 #if YYDEBUG
685 /*debugging the bison grammar --> grammar.cc*/
687 if (BVERBOSE(V_YACC)) yydebug=1;
688 else yydebug=0;
689 #endif
690 #endif
691 goto okay;
692 }
693 }
694 Werror("unknown option `%s`",n);
695 okay:
696 if (currRing != NULL)
699 v=v->next;
700 } while (v!=NULL);
701
702 // set global variable to show memory usage
704 else om_sing_opt_show_mem = 0;
705
706 return FALSE;
707}
708
710{
711 int i;
712 BITSET tmp;
713
714 StringSetS("//options:");
715 if ((si_opt_1!=0)||(si_opt_2!=0))
716 {
717 tmp=si_opt_1;
718 if(tmp)
719 {
720 for (i=0; optionStruct[i].setval!=0; i++)
721 {
722 if (optionStruct[i].setval & tmp)
723 {
725 tmp &=optionStruct[i].resetval;
726 }
727 }
728 for (i=0; i<32; i++)
729 {
730 if (tmp & Sy_bit(i)) StringAppend(" %d",i);
731 }
732 }
733 tmp=si_opt_2;
734 if (tmp)
735 {
736 for (i=0; verboseStruct[i].setval!=0; i++)
737 {
738 if (verboseStruct[i].setval & tmp)
739 {
741 tmp &=verboseStruct[i].resetval;
742 }
743 }
744 for (i=1; i<32; i++)
745 {
746 if (tmp & Sy_bit(i)) StringAppend(" %d",i+32);
747 }
748 }
749 return StringEndS();
750 }
751 StringAppendS(" none");
752 return StringEndS();
753}
754
755/* version strings */
756#ifdef HAVE_FLINT
757extern "C"
758{
759#ifndef __GMP_BITS_PER_MP_LIMB
760#define __GMP_BITS_PER_MP_LIMB GMP_LIMB_BITS
761#endif
762#include <flint/flint.h>
763}
764#endif
765
766#ifndef MAKE_DISTRIBUTION
767const char *singular_date = __DATE__ " " __TIME__;
768#endif
769
770char * versionString(/*const bool bShowDetails = false*/ )
771{
772 StringSetS("");
773 StringAppend("Singular for %s version %s (%d, %d bit) %s",
774 S_UNAME, VERSION, // SINGULAR_VERSION,
775 SINGULAR_VERSION, sizeof(void*)*8,
777 VERSION_DATE);
778#else
779 singular_date);
780#endif
781 StringAppendS("\nwith\n\t");
782
783#if defined(mpir_version)
784 StringAppend("MPIR(%s)~GMP(%s),", mpir_version, gmp_version);
785#elif defined(gmp_version)
786 // #if defined (__GNU_MP_VERSION) && defined (__GNU_MP_VERSION_MINOR)
787 // StringAppend("GMP(%d.%d),",__GNU_MP_VERSION,__GNU_MP_VERSION_MINOR);
788 StringAppend("GMP(%s),", gmp_version);
789#endif
790#ifdef HAVE_NTL
791 StringAppend("NTL(%s),",NTL_VERSION);
792#endif
793
794#ifdef HAVE_FLINT
795 StringAppend("FLINT(%s),",FLINT_VERSION);
796#endif
797// StringAppendS("factory(" FACTORYVERSION "),");
798 StringAppendS("\n\t");
799#ifndef HAVE_OMALLOC
800 StringAppendS("xalloc,");
801#else
802 StringAppendS("omalloc,");
803#endif
804#if defined(HAVE_DYN_RL)
806 StringAppendS("no input,");
807 else if (fe_fgets_stdin==fe_fgets)
808 StringAppendS("fgets,");
810 StringAppend("dynamic readline%d),",RL_VERSION_MAJOR);
811 #ifdef HAVE_FEREAD
813 StringAppendS("emulated readline,");
814 #endif
815 else
816 StringAppendS("unknown fgets method,");
817#else
818 #if defined(HAVE_READLINE) && !defined(FEREAD)
819 StringAppend("static readline(%d),",RL_VERSION_MAJOR);
820 #else
821 #ifdef HAVE_FEREAD
822 StringAppendS("emulated readline,");
823 #else
824 StringAppendS("fgets,");
825 #endif
826 #endif
827#endif
828#ifdef HAVE_PLURAL
829 StringAppendS("Plural,");
830#endif
831#ifdef HAVE_VSPACE
832 #if defined(__GNUC__) && (__GNUC__<9) &&!defined(__clang__)
833 StringAppendS("vspace(1),");
834 #else
835 StringAppendS("vspace(2),");
836 #endif
837#endif
838#ifdef HAVE_DBM
839 StringAppendS("DBM,\n\t");
840#else
841 StringAppendS("\n\t");
842#endif
843#ifdef HAVE_DYNAMIC_LOADING
844 StringAppendS("dynamic modules,");
845#endif
846#ifdef HAVE_DYNANIC_PPROCS
847 StringAppendS("dynamic p_Procs,");
848#endif
849#if YYDEBUG
850 StringAppendS("YYDEBUG=1,");
851#endif
852#ifdef MDEBUG
853 StringAppend("MDEBUG=%d,",MDEBUG);
854#endif
855#ifdef OM_CHECK
856 StringAppend("OM_CHECK=%d,",OM_CHECK);
857#endif
858#ifdef OM_TRACK
859 StringAppend("OM_TRACK=%d,",OM_TRACK);
860#endif
861#ifdef OM_NDEBUG
862 StringAppendS("OM_NDEBUG,");
863#endif
864#ifdef SING_NDEBUG
865 StringAppendS("SING_NDEBUG,");
866#endif
867#ifdef PDEBUG
868 StringAppendS("PDEBUG,");
869#endif
870#ifdef KDEBUG
871 StringAppendS("KDEBUG,");
872#endif
873 StringAppendS("\n\t");
874#ifdef __OPTIMIZE__
875 StringAppendS("CC:OPTIMIZE,");
876#endif
877#ifdef __OPTIMIZE_SIZE__
878 StringAppendS("CC:OPTIMIZE_SIZE,");
879#endif
880#ifdef __NO_INLINE__
881 StringAppendS("CC:NO_INLINE,");
882#endif
883#ifdef HAVE_NTL
884 #ifdef NTL_AVOID_BRANCHING
885 #undef HAVE_GENERIC_ADD
886 #endif
887#endif
888#ifdef HAVE_GENERIC_ADD
889 StringAppendS("GenericAdd,");
890#else
891 StringAppendS("AvoidBranching,");
892#endif
893#ifdef HAVE_GENERIC_MULT
894 StringAppendS("GenericMult,");
895#else
896 StringAppendS("TableMult,");
897#endif
898#ifdef HAVE_INVTABLE
899 StringAppendS("invTable,");
900#else
901 StringAppendS("no invTable,");
902#endif
903 StringAppendS("\n\t");
904#ifdef HAVE_EIGENVAL
905 StringAppendS("eigenvalues,");
906#endif
907#ifdef HAVE_GMS
908 StringAppendS("Gauss-Manin system,");
909#endif
910#ifdef HAVE_RATGRING
911 StringAppendS("ratGB,");
912#endif
913 StringAppend("random=%d\n",siRandomStart);
914
915#define SI_SHOW_BUILTIN_MODULE(name) StringAppend(" %s", #name);
916 StringAppendS("built-in modules: {");
918 StringAppendS("}\n");
919#undef SI_SHOW_BUILTIN_MODULE
920
921 StringAppend("AC_CONFIGURE_ARGS = %s,\n"
922 "CC = %s,FLAGS : %s,\n"
923 "CXX = %s,FLAGS : %s,\n"
924 "DEFS : %s,CPPFLAGS : %s,\n"
925 "LDFLAGS : %s,LIBS : %s "
926#ifdef __GNUC__
927 "(ver: " __VERSION__ ")"
928#endif
929 "\n",AC_CONFIGURE_ARGS, CC,CFLAGS " " PTHREAD_CFLAGS,
930 CXX,CXXFLAGS " " PTHREAD_CFLAGS, DEFS,CPPFLAGS, LDFLAGS,
931 LIBS " " PTHREAD_LIBS);
934 StringAppendS("\n");
935 return StringEndS();
936}
937
938#ifdef PDEBUG
939#if (OM_TRACK > 2) && defined(OM_TRACK_CUSTOM)
940void p_SetRingOfLeftv(leftv l, ring r)
941{
942 switch(l->rtyp)
943 {
944 case INT_CMD:
945 case BIGINT_CMD:
946 case IDHDL:
947 case DEF_CMD:
948 break;
949 case POLY_CMD:
950 case VECTOR_CMD:
951 {
952 poly p=(poly)l->data;
953 while(p!=NULL) { p_SetRingOfLm(p,r); pIter(p); }
954 break;
955 }
956 case IDEAL_CMD:
957 case MODUL_CMD:
958 case MATRIX_CMD:
959 {
960 ideal I=(ideal)l->data;
961 int i;
962 for(i=IDELEMS(I)-1;i>=0;i--)
963 {
964 poly p=I->m[i];
965 while(p!=NULL) { p_SetRingOfLm(p,r); pIter(p); }
966 }
967 break;
968 }
969 case COMMAND:
970 {
971 command d=(command)l->data;
972 p_SetRingOfLeftv(&d->arg1, r);
973 if (d->argc>1) p_SetRingOfLeftv(&d->arg2, r);
974 if (d->argc>2) p_SetRingOfLeftv(&d->arg3, r);
975 break;
976 }
977 default:
978 printf("type %d not yet implementd in p_SetRingOfLeftv\n",l->rtyp);
979 break;
980 }
981}
982#endif
983#endif
984
985#if 0 /* debug only */
986void listall(int showproc)
987{
988 idhdl hh=basePack->idroot;
989 PrintS("====== Top ==============\n");
990 while (hh!=NULL)
991 {
992 if (showproc || (IDTYP(hh)!=PROC_CMD))
993 {
994 if (IDDATA(hh)==(void *)currRing) PrintS("(R)");
995 else if (IDDATA(hh)==(void *)currPack) PrintS("(P)");
996 else PrintS(" ");
997 Print("::%s, typ %s level %d data %lx",
998 IDID(hh),Tok2Cmdname(IDTYP(hh)),IDLEV(hh),(long)IDDATA(hh));
999 if (IDTYP(hh)==RING_CMD)
1000 Print(" ref: %d\n",IDRING(hh)->ref);
1001 else
1002 PrintLn();
1003 }
1004 hh=IDNEXT(hh);
1005 }
1006 hh=basePack->idroot;
1007 while (hh!=NULL)
1008 {
1009 if (IDDATA(hh)==(void *)basePack)
1010 Print("(T)::%s, typ %s level %d data %lx\n",
1011 IDID(hh),Tok2Cmdname(IDTYP(hh)),IDLEV(hh),(long)IDDATA(hh));
1012 else
1013 if ((IDTYP(hh)==RING_CMD)
1014 || (IDTYP(hh)==PACKAGE_CMD))
1015 {
1016 Print("====== %s ==============\n",IDID(hh));
1017 idhdl h2=IDRING(hh)->idroot;
1018 while (h2!=NULL)
1019 {
1020 if (showproc || (IDTYP(h2)!=PROC_CMD))
1021 {
1022 if ((IDDATA(h2)==(void *)currRing)
1023 && (IDTYP(h2)==RING_CMD))
1024 PrintS("(R)");
1025 else if (IDDATA(h2)==(void *)currPack) PrintS("(P)");
1026 else PrintS(" ");
1027 Print("%s::%s, typ %s level %d data %lx\n",
1028 IDID(hh),IDID(h2),Tok2Cmdname(IDTYP(h2)),IDLEV(h2),(long)IDDATA(h2));
1029 }
1030 h2=IDNEXT(h2);
1031 }
1032 }
1033 hh=IDNEXT(hh);
1034 }
1035 Print("currRing:%lx, currPack:%lx,basePack:%lx\n",(long)currRing,(long)currPack,(long)basePack);
1037}
1038#endif
1039
1040#ifndef SING_NDEBUG
1041void checkall()
1042{
1043 idhdl hh=basePack->idroot;
1044 while (hh!=NULL)
1045 {
1046 omCheckAddr(hh);
1047 omCheckAddr((ADDRESS)IDID(hh));
1048 if (RingDependend(IDTYP(hh)))
1049 {
1050 Print("%s typ %d in Top (should be in ring)\n",IDID(hh),IDTYP(hh));
1051 }
1052 hh=IDNEXT(hh);
1053 }
1054 hh=basePack->idroot;
1055 while (hh!=NULL)
1056 {
1057 if (IDTYP(hh)==PACKAGE_CMD)
1058 {
1059 idhdl h2=NULL;
1060 if (IDPACKAGE(hh)!=NULL)
1061 h2=IDPACKAGE(hh)->idroot;
1062 if (IDPACKAGE(hh)!=basePack)
1063 {
1064 while (h2!=NULL)
1065 {
1066 omCheckAddr(h2);
1067 omCheckAddr((ADDRESS)IDID(h2));
1068 if (RingDependend(IDTYP(h2)))
1069 {
1070 Print("%s typ %d in %s (should be in ring)\n",IDID(h2),IDTYP(h2),IDID(hh));
1071 }
1072 h2=IDNEXT(h2);
1073 }
1074 }
1075 }
1076 hh=IDNEXT(hh);
1077 }
1078}
1079#endif
1080
1081extern "C"
1082int singular_fstat(int fd, struct stat *buf)
1083{
1084 return si_fstat(fd,buf);
1085}
1086
1087/*2
1088* the global exit routine of Singular
1089*/
1090extern "C" {
1091/* Note: We cannot use a mutex here because mutexes are not async-safe, but
1092 * m2_end is called by sig_term_hdl(). Anyway, the race condition in the first
1093 * few lines of m2_end() should not matter.
1094 */
1096
1097void m2_end(int i)
1098{
1099 if (!m2_end_called)
1100 {
1102 EXTERN_VAR FILE* File_Log;
1105 if (File_Log!=NULL)
1106 {
1107 fclose(File_Log);
1108 File_Log=NULL;
1109 if (File_Log_written==FALSE) // remove empty logs
1110 {
1111 int pid=getpid();
1112 char buf[20];
1113 sprintf(buf,"/tmp/sing_log.%d",pid);
1114 remove(buf);
1115 }
1116 }
1118#ifdef HAVE_SIMPLEIPC
1119 for (int j = SIPC_MAX_SEMAPHORES-1; j >= 0; j--)
1120 {
1121 if (semaphore[j] != NULL)
1122 {
1123 while (sem_acquired[j] > 0)
1124 {
1125#if PORTABLE_SEMAPHORES
1126 sem_post(semaphore[j]->sig);
1127#else
1128 sem_post(semaphore[j]);
1129#endif
1130 sem_acquired[j]--;
1131 }
1132 }
1133 }
1134#endif // HAVE_SIMPLEIPC
1135 monitor(NULL,0);
1136#ifdef PAGE_TEST
1137 mmEndStat();
1138#endif
1141 {
1143 while(hh!=NULL)
1144 {
1145 //Print("close %s\n",hh->l->name);
1146 slPrepClose(hh->l);
1147 hh=(link_list)hh->next;
1148 }
1150
1151 idhdl h = currPack->idroot;
1152 while(h != NULL)
1153 {
1154 if(IDTYP(h) == LINK_CMD)
1155 {
1156 idhdl hh=h->next;
1157 //Print("kill %s\n",IDID(h));
1158 killhdl(h, currPack);
1159 h = hh;
1160 }
1161 else
1162 {
1163 h = h->next;
1164 }
1165 }
1166 hh=ssiToBeClosed;
1167 while(hh!=NULL)
1168 {
1169 //Print("close %s\n",hh->l->name);
1170 slClose(hh->l);
1171 hh=ssiToBeClosed;
1172 }
1173 }
1175 {
1176 if (i<=0)
1177 {
1178 //extern long all_farey;
1179 //extern long farey_cnt;
1180 //if (all_farey!=0L) printf("farey:%ld, cnt=%ld\n",all_farey,farey_cnt);
1181 if (TEST_V_QUIET)
1182 {
1183 if (i==0)
1184 printf("Auf Wiedersehen.\n");
1185 else
1186 printf("\n$Bye.\n");
1187 }
1188 //#ifdef sun
1189 // #ifndef __svr4__
1190 // _cleanup();
1191 // _exit(0);
1192 // #endif
1193 //#endif
1194 i=0;
1195 }
1196 else
1197 {
1198 printf("\nhalt %d\n",i);
1199 }
1200 }
1201 exit(i);
1202 }
1203}
1204}
1205
1206extern "C"
1207{
1209 {
1210 fprintf(stderr, "\nSingular error: no more memory\n");
1211 omPrintStats(stderr);
1212 m2_end(14);
1213 /* should never get here */
1214 exit(1);
1215 }
1216}
1217
1218#ifdef HAVE_FLINT
1221//STATIC_VAR n_coeffType n_FlintQrat=n_unknown;
1223{
1224 const short t[]={2,INT_CMD,STRING_CMD};
1225 if (iiCheckTypes(a,t,1))
1226 {
1228 p.ch=(int)(long)a->Data();
1229 p.name=(char*)a->next->Data();
1230 res->rtyp=CRING_CMD;
1231 res->data=(void*)nInitChar(n_FlintZn,(void*)&p);
1232 return FALSE;
1233 }
1234 return TRUE;
1235}
1237{
1238 const short t[]={1,STRING_CMD};
1239 if (iiCheckTypes(a,t,1))
1240 {
1241 char* p;
1242 p=(char*)a->Data();
1243 res->rtyp=CRING_CMD;
1244 res->data=(void*)nInitChar(n_FlintQ,(void*)p);
1245 return FALSE;
1246 }
1247 return TRUE;
1248}
1249#if __FLINT_RELEASE >= 20503
1250static BOOLEAN ii_FlintQrat_init(leftv res,leftv a)
1251{
1252 if (a==NULL)
1253 {
1254 WerrorS("at least one name required");
1255 return TRUE;
1256 }
1257 QaInfo par;
1258 #ifdef QA_DEBUG
1259 par.C=r->cf;
1260 a=a->next;
1261 #endif
1262 par.N=a->listLength();
1263 par.names=(char**)omAlloc(par.N*sizeof(char*));
1264 int i=0;
1265 while(a!=NULL)
1266 {
1267 par.names[i]=omStrDup(a->Name());
1268 i++;
1269 a=a->next;
1270 }
1271 res->rtyp=CRING_CMD;
1272 res->data=(void*)nInitChar(n_FlintQrat,&par);
1273 for(i=par.N-1;i>=0;i--)
1274 {
1275 omFree(par.names[i]);
1276 }
1277 omFreeSize(par.names,par.N*sizeof(char*));
1278 return FALSE;
1279}
1280#endif
1281extern "C" int flint_mod_init(SModulFunctions* psModulFunctions)
1282{
1283 package save=currPack;
1286 if (n_FlintQ!=n_unknown)
1287 {
1288 iiAddCproc("kernel","flintQp",FALSE,ii_FlintQ_init);
1290 }
1291#if __FLINT_RELEASE >= 20503
1292 iiAddCproc("kernel","flintQ",FALSE,ii_FlintQrat_init);
1294#endif
1296 if (n_FlintZn!=n_unknown)
1297 {
1298 iiAddCproc("kernel","flintZn",FALSE,ii_FlintZn_init);
1300 }
1301 currPack=save;
1302 return MAX_TOK;
1303}
1304#endif
1305
1307{
1308 short float_len=3;
1309 short float_len2=SHORT_REAL_LENGTH;
1310 coeffs cf=NULL;
1311 if ((pnn!=NULL) && (pnn->Typ()==INT_CMD))
1312 {
1313 float_len=(int)(long)pnn->Data();
1314 float_len2=float_len;
1315 pnn=pnn->next;
1316 if ((pnn!=NULL) && (pnn->Typ()==INT_CMD))
1317 {
1318 float_len2=(int)(long)pnn->Data();
1319 pnn=pnn->next;
1320 }
1321 }
1322 if (float_len2 <= (short)SHORT_REAL_LENGTH)
1323 cf=nInitChar(n_R, NULL);
1324 else // longR or longC?
1325 {
1326 LongComplexInfo param;
1327 param.float_len = si_min (float_len, 32767);
1328 param.float_len2 = si_min (float_len2, 32767);
1329 cf = nInitChar(n_long_R, (void*)&param);
1330 }
1331 res->rtyp=CRING_CMD;
1332 res->data=cf;
1333 return cf==NULL;
1334}
1336{
1337 leftv h=args;
1338 coeffs *c=NULL;
1339 coeffs cf=NULL;
1340 int i=0;
1341 if (h==NULL) goto crossprod_error;
1342 while (h!=NULL)
1343 {
1344 if (h->Typ()!=CRING_CMD) goto crossprod_error;
1345 i++;
1346 h=h->next;
1347 }
1348 c=(coeffs*)omAlloc0((i+1)*sizeof(coeffs));
1349 h=args;
1350 i=0;
1351 while (h!=NULL)
1352 {
1353 c[i]=(coeffs)h->CopyD();
1354 i++;
1355 h=h->next;
1356 }
1358 res->data=cf;
1359 res->rtyp=CRING_CMD;
1360 return FALSE;
1361
1362 crossprod_error:
1363 WerrorS("expected `crossprod(coeffs, ...)`");
1364 return TRUE;
1365}
1366/*2
1367* initialize components of Singular
1368*/
1369static void callWerrorS(const char *s) { WerrorS(s); }
1370void siInit(char *name)
1371{
1372// memory initialization: -----------------------------------------------
1373 om_Opts.OutOfMemoryFunc = omSingOutOfMemoryFunc;
1374#ifndef OM_NDEBUG
1375#ifndef __OPTIMIZE__
1376 om_Opts.ErrorHook = dErrorBreak;
1377#else
1378 om_Opts.Keep = 0; /* !OM_NDEBUG, __OPTIMIZE__*/
1379#endif
1380#else
1381 om_Opts.Keep = 0; /* OM_NDEBUG */
1382#endif
1383 omInitInfo();
1384// factory
1385#ifndef HAVE_NTL
1386 extern void initPT();
1387 initPT();
1388#endif
1389// options ---------------------------------------------------------------
1390 si_opt_1=0;
1391// interpreter tables etc.: -----------------------------------------------
1392 memset(&sLastPrinted,0,sizeof(sleftv));
1394
1395 extern int iiInitArithmetic(); iiInitArithmetic(); // iparith.cc
1396
1397 basePack=(package)omAlloc0(sizeof(*basePack));
1399 idhdl h;
1400 h=enterid("Top", 0, PACKAGE_CMD, &IDROOT, FALSE);
1402 IDPACKAGE(h)->language = LANG_TOP;
1403 currPackHdl=h;
1404 basePackHdl=h;
1405
1406 coeffs_BIGINT = nInitChar(n_Q,(void*)1);
1407
1408#if 1
1409 // def HAVE_POLYEXTENSIONS
1410 if(TRUE)
1411 {
1412 n_coeffType type;
1413 #ifdef SINGULAR_4_2
1415 assume(type == n_polyExt);
1416 #endif
1417
1418 type = nRegister(n_algExt, naInitChar);
1419 assume(type == n_algExt);
1420
1422 assume(type == n_transExt);
1423
1424 (void)type;
1425 }
1426#endif
1427
1428// random generator: -----------------------------------------------
1429 int t=initTimer();
1430 if (t==0) t=1;
1431 initRTimer();
1432 siSeed=t;
1433 factoryseed(t);
1434 siRandomStart=t;
1435 feOptSpec[FE_OPT_RANDOM].value = (void*) ((long)siRandomStart);
1436
1437// ressource table: ----------------------------------------------------
1438 // Don't worry: ifdef OM_NDEBUG, then all these calls are undef'ed
1439 // hack such that all shared' libs in the bindir are loaded correctly
1441
1442// singular links: --------------------------------------------------
1444 myynest=0;
1445// how many processes ? -----------------------------------------------------
1446 int cpus=2;
1447 int cpu_n;
1448 #ifdef _SC_NPROCESSORS_ONLN
1449 if ((cpu_n=sysconf(_SC_NPROCESSORS_ONLN))>cpus) cpus=cpu_n;
1450 #elif defined(_SC_NPROCESSORS_CONF)
1451 if ((cpu_n=sysconf(_SC_NPROCESSORS_CONF))>cpus) cpus=cpu_n;
1452 #endif
1453 feSetOptValue(FE_OPT_CPUS, cpus);
1454// how many threads ? -----------------------------------------------------
1455 feSetOptValue(FE_OPT_THREADS, cpus);
1456
1457// default coeffs
1458 {
1459 idhdl h;
1460 h=enterid("QQ",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1461 IDDATA(h)=(char*)nInitChar(n_Q,NULL);
1462 h=enterid("ZZ",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1463 IDDATA(h)=(char*)nInitChar(n_Z,NULL);
1465 iiAddCproc("kernel","crossprod",FALSE,iiCrossProd);
1466 iiAddCproc("kernel","Float",FALSE,iiFloat);
1467 //h=enterid("RR",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1468 //IDDATA(h)=(char*)nInitChar(n_R,NULL);
1469 //h=enterid("CC",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1470 //IDDATA(h)=(char*)nInitChar(n_long_C,NULL);
1471 }
1472// setting routines for PLURAL QRINGS:
1473// allowing to use libpolys without libSingular(kStd)
1474#ifdef HAVE_PLURAL
1475 nc_NF=k_NF;
1481#endif
1482// loading standard.lib -----------------------------------------------
1483 if (! feOptValue(FE_OPT_NO_STDLIB))
1484 {
1485 BITSET save1,save2;
1486 SI_SAVE_OPT(save1,save2);
1487 si_opt_2 &= ~Sy_bit(V_LOAD_LIB);
1488 iiLibCmd("standard.lib", TRUE,TRUE,TRUE);
1489 SI_RESTORE_OPT(save1,save2);
1490 }
1491 // interpreter error handling
1492 #ifndef __CYGWIN__
1493 factoryError=callWerrorS; // to honour later changes of variable WerrorS
1494 #endif
1495 errorreported = 0;
1496}
BOOLEAN naInitChar(coeffs cf, void *infoStruct)
Initialize the coeffs object.
Definition: algext.cc:1388
BOOLEAN n2pInitChar(coeffs cf, void *infoStruct)
Definition: algext.cc:1642
int BOOLEAN
Definition: auxiliary.h:87
#define TRUE
Definition: auxiliary.h:100
#define FALSE
Definition: auxiliary.h:96
void * ADDRESS
Definition: auxiliary.h:119
static int si_min(const int a, const int b)
Definition: auxiliary.h:125
int l
Definition: cfEzgcd.cc:100
int i
Definition: cfEzgcd.cc:132
int k
Definition: cfEzgcd.cc:99
Variable x
Definition: cfModGcd.cc:4082
int p
Definition: cfModGcd.cc:4078
CanonicalForm cf
Definition: cfModGcd.cc:4083
CanonicalForm test
Definition: cfModGcd.cc:4096
void initPT()
static CanonicalForm bound(const CFMatrix &M)
Definition: cf_linsys.cc:460
void factoryseed(int s)
random seed initializer
Definition: cf_random.cc:189
VAR void(* factoryError)(const char *s)
Definition: cf_util.cc:80
FILE * f
Definition: checklibs.c:9
char name() const
Definition: variable.cc:122
Variable next() const
Definition: variable.h:52
Definition: idrec.h:35
Definition: intvec.h:23
Class used for (list of) interpreter objects.
Definition: subexpr.h:83
int Typ()
Definition: subexpr.cc:1019
int rtyp
Definition: subexpr.h:91
void * Data()
Definition: subexpr.cc:1162
leftv next
Definition: subexpr.h:86
const char * Name()
Definition: subexpr.h:120
int listLength()
Definition: subexpr.cc:51
void * data
Definition: subexpr.h:88
Definition: lists.h:24
sleftv * m
Definition: lists.h:46
INLINE_THIS void Init(int l=0)
VAR BOOLEAN singular_in_batchmode
Definition: cntrlc.cc:62
VAR int siRandomStart
Definition: cntrlc.cc:93
Coefficient rings, fields and other domains suitable for Singular polynomials.
static FORCE_INLINE void number2mpz(number n, coeffs c, mpz_t m)
Definition: coeffs.h:984
static FORCE_INLINE number mpz2number(mpz_t m, coeffs c)
Definition: coeffs.h:985
n_coeffType
Definition: coeffs.h:27
@ n_R
single prescision (6,6) real numbers
Definition: coeffs.h:31
@ n_FlintQrat
rational funtion field over Q
Definition: coeffs.h:47
@ n_polyExt
used to represent polys as coeffcients
Definition: coeffs.h:34
@ n_Q
rational (GMP) numbers
Definition: coeffs.h:30
@ n_algExt
used for all algebraic extensions, i.e., the top-most extension in an extension tower is algebraic
Definition: coeffs.h:35
@ n_Zn
only used if HAVE_RINGS is defined
Definition: coeffs.h:44
@ n_long_R
real floating point (GMP) numbers
Definition: coeffs.h:33
@ n_transExt
used for all transcendental extensions, i.e., the top-most extension in an extension tower is transce...
Definition: coeffs.h:38
@ n_unknown
Definition: coeffs.h:28
@ n_Z
only used if HAVE_RINGS is defined
Definition: coeffs.h:43
@ n_nTupel
n-tupel of cf: ZZ/p1,... ZZ/pn, R, long_R
Definition: coeffs.h:42
short float_len2
additional char-flags, rInit
Definition: coeffs.h:102
static FORCE_INLINE BOOLEAN n_GreaterZero(number n, const coeffs r)
ordered fields: TRUE iff 'n' is positive; in Z/pZ: TRUE iff 0 < m <= roundedBelow(p/2),...
Definition: coeffs.h:491
coeffs nInitChar(n_coeffType t, void *parameter)
one-time initialisations for new coeffs in case of an error return NULL
Definition: numbers.cc:413
static FORCE_INLINE BOOLEAN n_IsZero(number n, const coeffs r)
TRUE iff 'n' represents the zero element.
Definition: coeffs.h:461
static FORCE_INLINE number n_Init(long i, const coeffs r)
a number representing i in the given coeff field/ring r
Definition: coeffs.h:535
short float_len
additional char-flags, rInit
Definition: coeffs.h:101
#define MAKE_DISTRIBUTION
Definition: distrib.h:1
#define Print
Definition: emacs.cc:80
#define WarnS
Definition: emacs.cc:78
#define StringAppend
Definition: emacs.cc:79
const CanonicalForm int s
Definition: facAbsFact.cc:51
const CanonicalForm int const CFList const Variable & y
Definition: facAbsFact.cc:53
CanonicalForm res
Definition: facAbsFact.cc:60
const CanonicalForm & w
Definition: facAbsFact.cc:51
const Variable & v
< [in] a sqrfree bivariate poly
Definition: facBivar.h:39
int j
Definition: facHensel.cc:110
FILE * feFopen(const char *path, const char *mode, char *where, short useWerror, short path_only)
Definition: feFopen.cc:47
VAR short errorreported
Definition: feFopen.cc:23
void WerrorS(const char *s)
Definition: feFopen.cc:24
const char * feSetOptValue(feOptIndex opt, char *optarg)
Definition: feOpt.cc:154
static void * feOptValue(feOptIndex opt)
Definition: feOpt.h:40
EXTERN_VAR struct fe_option feOptSpec[]
Definition: feOpt.h:17
static char * feResource(feResourceConfig config, int warn)
Definition: feResource.cc:236
void feInitResources(const char *argv0)
Definition: feResource.cc:163
VAR int si_echo
Definition: febase.cc:35
VAR int myynest
Definition: febase.cc:41
void monitor(void *F, int mode)
Definition: febase.cc:68
void * value
Definition: fegetopt.h:93
void feStringAppendBrowsers(int warn)
Definition: fehelp.cc:341
char * fe_fgets_dummy(const char *, char *, int)
Definition: feread.cc:455
char *(* fe_fgets_stdin)(const char *pr, char *s, int size)
Definition: feread.cc:32
char * fe_fgets(const char *pr, char *s, int size)
Definition: feread.cc:309
char * fe_fgets_stdin_drl(const char *pr, char *s, int size)
Definition: feread.cc:269
char * fe_fgets_stdin_emu(const char *pr, char *s, int size)
Definition: feread.cc:253
void fe_reset_input_mode()
Definition: fereadl.c:831
VAR FILE * File_Profiling
Definition: fevoices.cc:32
VAR BOOLEAN File_Log_written
Definition: fevoices.cc:34
VAR FILE * File_Log
Definition: fevoices.cc:33
BOOLEAN flintQ_InitChar(coeffs cf, void *infoStruct)
Definition: flintcf_Q.cc:561
coeffs flintQInitCfByName(char *s, n_coeffType n)
Definition: flintcf_Q.cc:533
BOOLEAN flintZn_InitChar(coeffs cf, void *infoStruct)
Definition: flintcf_Zn.cc:486
coeffs flintZnInitCfByName(char *s, n_coeffType n)
Definition: flintcf_Zn.cc:432
EXTERN_VAR BBA_Proc gnc_gr_bba
Definition: gb_hack.h:10
EXTERN_VAR BBA_Proc gnc_gr_mora
Definition: gb_hack.h:10
EXTERN_VAR BBA_Proc sca_gr_bba
Definition: gb_hack.h:10
EXTERN_VAR NF_Proc nc_NF
Definition: gb_hack.h:9
EXTERN_VAR BBA_Proc sca_mora
Definition: gb_hack.h:10
EXTERN_VAR BBA_Proc sca_bba
Definition: gb_hack.h:10
STATIC_VAR unsigned short primes[]
primes, primes_len: used to step through possible extensions
const char * Tok2Cmdname(int tok)
Definition: gentable.cc:140
static int RingDependend(int t)
Definition: gentable.cc:28
#define STATIC_VAR
Definition: globaldefs.h:7
#define EXTERN_VAR
Definition: globaldefs.h:6
ideal k_gnc_gr_mora(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Definition: gr_kstd2.cc:1283
ideal k_gnc_gr_bba(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Definition: gr_kstd2.cc:1030
VAR int yydebug
Definition: grammar.cc:1805
@ IDEAL_CMD
Definition: grammar.cc:284
@ MATRIX_CMD
Definition: grammar.cc:286
@ PROC_CMD
Definition: grammar.cc:280
@ MODUL_CMD
Definition: grammar.cc:287
@ VECTOR_CMD
Definition: grammar.cc:292
@ POLY_CMD
Definition: grammar.cc:289
@ RING_CMD
Definition: grammar.cc:281
static BOOLEAN length(leftv result, leftv arg)
Definition: interval.cc:257
int iiInitArithmetic()
initialisation of arithmetic structured data
Definition: iparith.cc:9738
idhdl enterid(const char *s, int lev, int t, idhdl *root, BOOLEAN init, BOOLEAN search)
Definition: ipid.cc:279
VAR package basePack
Definition: ipid.cc:58
VAR package currPack
Definition: ipid.cc:57
void killhdl(idhdl h, package proot)
Definition: ipid.cc:414
VAR idhdl currPackHdl
Definition: ipid.cc:55
VAR idhdl basePackHdl
Definition: ipid.cc:56
VAR coeffs coeffs_BIGINT
Definition: ipid.cc:50
#define IDNEXT(a)
Definition: ipid.h:118
ip_command * command
Definition: ipid.h:23
#define IDDATA(a)
Definition: ipid.h:126
#define IDPROC(a)
Definition: ipid.h:140
#define IDID(a)
Definition: ipid.h:122
#define IDROOT
Definition: ipid.h:19
unsigned resetval
Definition: ipid.h:154
#define IDPACKAGE(a)
Definition: ipid.h:139
#define IDLEV(a)
Definition: ipid.h:121
unsigned setval
Definition: ipid.h:153
#define IDRING(a)
Definition: ipid.h:127
#define IDTYP(a)
Definition: ipid.h:119
int iiAddCproc(const char *libname, const char *procname, BOOLEAN pstatic, BOOLEAN(*func)(leftv res, leftv v))
Definition: iplib.cc:1063
BOOLEAN iiEStart(char *example, procinfo *pi)
Definition: iplib.cc:754
BOOLEAN iiLibCmd(const char *newlib, BOOLEAN autoexport, BOOLEAN tellerror, BOOLEAN force)
Definition: iplib.cc:884
char * iiGetLibProcBuffer(procinfo *pi, int part)
Definition: iplib.cc:197
void iiCheckPack(package &p)
Definition: ipshell.cc:1630
BOOLEAN iiCheckTypes(leftv args, const short *type_list, int report)
check a list of arguemys against a given field of types return TRUE if the types match return FALSE (...
Definition: ipshell.cc:6572
static char * iiGetLibName(const procinfov pi)
find the library of an proc
Definition: ipshell.h:66
STATIC_VAR Poly * h
Definition: janet.cc:971
ideal k_sca_bba(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Modified modern Sinuglar Buchberger's algorithm.
Definition: sca.cc:368
ideal k_sca_mora(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Modified modern Sinuglar Mora's algorithm.
Definition: sca.cc:885
ideal k_sca_gr_bba(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Modified Plural's Buchberger's algorithmus.
Definition: sca.cc:95
poly k_NF(ideal F, ideal Q, poly p, int syzComp, int lazyReduce, const ring _currRing)
NOTE: this is just a wrapper which sets currRing for the actual kNF call.
Definition: kstd1.cc:3396
VAR BITSET validOpts
Definition: kstd1.cc:60
VAR omBin slists_bin
Definition: lists.cc:23
static int factor_using_division(mpz_t t, unsigned int limit, lists primes, int *multiplicities, int &index, unsigned long bound)
Definition: misc_ip.cc:109
int singular_fstat(int fd, struct stat *buf)
Definition: misc_ip.cc:1082
#define SI_SHOW_BUILTIN_MODULE(name)
const struct soptionStruct verboseStruct[]
Definition: misc_ip.cc:538
static void factor_using_pollard_rho(mpz_t n, unsigned long a, lists primes, int *multiplicities, int &index)
Definition: misc_ip.cc:204
static BOOLEAN iiCrossProd(leftv res, leftv args)
Definition: misc_ip.cc:1335
void siInit(char *name)
Definition: misc_ip.cc:1370
char * versionString()
Definition: misc_ip.cc:770
static BOOLEAN ii_FlintQ_init(leftv res, leftv a)
Definition: misc_ip.cc:1236
int flint_mod_init(SModulFunctions *psModulFunctions)
Definition: misc_ip.cc:1281
STATIC_VAR n_coeffType n_FlintQ
Definition: misc_ip.cc:1220
lists primeFactorisation(const number n, const int pBound)
Factorises a given bigint number n into its prime factors less than or equal to a given bound,...
Definition: misc_ip.cc:357
void omSingOutOfMemoryFunc()
Definition: misc_ip.cc:1208
BOOLEAN setOption(leftv res, leftv v)
Definition: misc_ip.cc:568
static void factor_gmp(mpz_t t, lists primes, int *multiplicities, int &index, unsigned long bound)
Definition: misc_ip.cc:327
static BOOLEAN iiFloat(leftv res, leftv pnn)
Definition: misc_ip.cc:1306
void setListEntry(lists L, int index, mpz_t n)
Definition: misc_ip.cc:74
void singular_example(char *str)
Definition: misc_ip.cc:430
static BOOLEAN ii_FlintZn_init(leftv res, leftv a)
Definition: misc_ip.cc:1222
const struct soptionStruct optionStruct[]
Definition: misc_ip.cc:507
void setListEntry_ui(lists L, int index, unsigned long ui)
Definition: misc_ip.cc:91
STATIC_VAR n_coeffType n_FlintZn
Definition: misc_ip.cc:1219
char * showOption()
Definition: misc_ip.cc:709
static void callWerrorS(const char *s)
Definition: misc_ip.cc:1369
volatile BOOLEAN m2_end_called
Definition: misc_ip.cc:1095
STATIC_VAR unsigned add[]
Definition: misc_ip.cc:107
void m2_end(int i)
Definition: misc_ip.cc:1097
This file provides miscellaneous functionality.
void dErrorBreak(void)
#define SEEK_SET
Definition: mod2.h:115
#define assume(x)
Definition: mod2.h:389
#define SINGULAR_VERSION
Definition: mod2.h:87
#define SEEK_END
Definition: mod2.h:111
#define MDEBUG
Definition: mod2.h:180
#define SI_FOREACH_BUILTIN(add)
Data for type_of_LIB to determine built-in modules, use add(name) to add built-in library to macro.
Definition: mod_lib.h:17
#define pIter(p)
Definition: monomials.h:37
#define p_SetRingOfLm(p, r)
Definition: monomials.h:144
slists * lists
Definition: mpr_numeric.h:146
The main handler for Singular numbers which are suitable for Singular polynomials.
void nRegisterCfByName(cfInitCfByNameProc p, n_coeffType n)
Definition: numbers.cc:636
n_coeffType nRegister(n_coeffType n, cfInitCharProc p)
Definition: numbers.cc:595
#define SHORT_REAL_LENGTH
Definition: numbers.h:57
#define omStrDup(s)
Definition: omAllocDecl.h:263
#define omFreeSize(addr, size)
Definition: omAllocDecl.h:260
#define omCheckAddr(addr)
Definition: omAllocDecl.h:328
#define omAlloc(size)
Definition: omAllocDecl.h:210
#define omAllocBin(bin)
Definition: omAllocDecl.h:205
#define omFree(addr)
Definition: omAllocDecl.h:261
#define omAlloc0(size)
Definition: omAllocDecl.h:211
#define omFreeBinAddr(addr)
Definition: omAllocDecl.h:258
#define NULL
Definition: omList.c:12
omOpts_t om_Opts
Definition: omOpts.c:13
#define MAXPATHLEN
Definition: omRet2Info.c:22
int om_sing_opt_show_mem
#define OM_TRACK
Definition: omalloc_debug.c:10
#define OM_CHECK
Definition: omalloc_debug.c:15
VAR unsigned si_opt_2
Definition: options.c:6
VAR unsigned si_opt_1
Definition: options.c:5
#define OPT_SUGARCRIT
Definition: options.h:81
#define OPT_PROT
Definition: options.h:76
#define OPT_INFREDTAIL
Definition: options.h:95
#define V_QRING
Definition: options.h:42
#define SI_SAVE_OPT(A, B)
Definition: options.h:20
#define V_DEF_RES
Definition: options.h:50
#define OPT_INTSTRATEGY
Definition: options.h:93
#define V_INTERSECT_SYZ
Definition: options.h:69
#define TEST_OPT_INTSTRATEGY
Definition: options.h:111
#define BVERBOSE(a)
Definition: options.h:35
#define OPT_WEIGHTM
Definition: options.h:98
#define V_SHOW_MEM
Definition: options.h:43
#define OPT_REDTAIL_SYZ
Definition: options.h:88
#define V_ALLWARN
Definition: options.h:67
#define TEST_V_QUIET
Definition: options.h:135
#define V_READING
Definition: options.h:46
#define V_COEFSTRAT
Definition: options.h:61
#define V_ASSIGN_NONE
Definition: options.h:70
#define V_CONTENTSB
Definition: options.h:56
#define OPT_REDTAIL
Definition: options.h:92
#define V_DEBUG_LIB
Definition: options.h:48
#define V_MODPSOLVSB
Definition: options.h:58
#define OPT_RETURN_SB
Definition: options.h:85
#define V_LOAD_PROC
Definition: options.h:49
#define OPT_NOT_SUGAR
Definition: options.h:79
#define V_UPTORADICAL
Definition: options.h:59
#define V_YACC
Definition: options.h:44
#define OPT_REDTHROUGH
Definition: options.h:83
#define OPT_REDSB
Definition: options.h:77
#define Sy_bit(x)
Definition: options.h:31
#define OPT_NOTREGULARITY
Definition: options.h:97
#define V_FINDMONOM
Definition: options.h:60
#define V_CANCELUNIT
Definition: options.h:57
#define V_REDEFINE
Definition: options.h:45
#define V_INTERSECT_ELIM
Definition: options.h:68
#define OPT_DEBUG
Definition: options.h:82
#define V_LOAD_LIB
Definition: options.h:47
#define OPT_MULTBOUND
Definition: options.h:90
#define V_NSB
Definition: options.h:55
#define V_LENGTH
Definition: options.h:64
#define OPT_STAIRCASEBOUND
Definition: options.h:89
#define OPT_INTERRUPT
Definition: options.h:80
#define OPT_DEGBOUND
Definition: options.h:91
#define OPT_NOT_BUCKETS
Definition: options.h:78
#define V_IMAP
Definition: options.h:53
#define OPT_FASTHC
Definition: options.h:86
#define TEST_RINGDEP_OPTS
Definition: options.h:101
#define V_PROMPT
Definition: options.h:54
#define OPT_OLDSTD
Definition: options.h:87
#define V_SHOW_USE
Definition: options.h:52
#define SI_RESTORE_OPT(A, B)
Definition: options.h:23
#define OPT_NO_SYZ_MINIM
Definition: options.h:84
static int index(p_Length length, p_Ord ord)
Definition: p_Procs_Impl.h:592
VAR ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
Definition: polys.cc:13
void StringSetS(const char *st)
Definition: reporter.cc:128
void StringAppendS(const char *st)
Definition: reporter.cc:107
void PrintS(const char *s)
Definition: reporter.cc:284
void feStringAppendResources(int warn)
Definition: reporter.cc:398
char * StringEndS()
Definition: reporter.cc:151
void PrintLn()
Definition: reporter.cc:310
void Werror(const char *fmt,...)
Definition: reporter.cc:189
static BOOLEAN rField_has_simple_inverse(const ring r)
Definition: ring.h:548
#define rField_is_Ring(R)
Definition: ring.h:485
coeffs nrnInitCfByName(char *s, n_coeffType)
Definition: rmodulon.cc:35
VAR sipc_sem_t * semaphore[SIPC_MAX_SEMAPHORES]
Definition: semaphore.c:24
VAR int sem_acquired[SIPC_MAX_SEMAPHORES]
Definition: semaphore.c:25
#define mpz_size1(A)
Definition: si_gmp.h:17
#define mpz_sgn1(A)
Definition: si_gmp.h:18
int status int void * buf
Definition: si_signals.h:59
int status int fd
Definition: si_signals.h:59
#define IDELEMS(i)
Definition: simpleideals.h:23
#define SIPC_MAX_SEMAPHORES
Definition: simpleipc.h:10
VAR int siSeed
Definition: sirandom.c:30
ip_package * package
Definition: structs.h:43
#define BITSET
Definition: structs.h:16
#define loop
Definition: structs.h:75
INST_VAR sleftv sLastPrinted
Definition: subexpr.cc:46
@ LANG_TOP
Definition: subexpr.h:22
int name
New type name for int.
Definition: templateForC.h:21
int initTimer()
Definition: timer.cc:67
void initRTimer()
Definition: timer.cc:156
#define IDHDL
Definition: tok.h:31
@ BIGINT_CMD
Definition: tok.h:38
@ CRING_CMD
Definition: tok.h:56
@ LIST_CMD
Definition: tok.h:118
@ INTVEC_CMD
Definition: tok.h:101
@ PACKAGE_CMD
Definition: tok.h:149
@ DEF_CMD
Definition: tok.h:58
@ LINK_CMD
Definition: tok.h:117
@ STRING_CMD
Definition: tok.h:185
@ INT_CMD
Definition: tok.h:96
@ MAX_TOK
Definition: tok.h:218
#define NONE
Definition: tok.h:221
#define COMMAND
Definition: tok.h:29
BOOLEAN ntInitChar(coeffs cf, void *infoStruct)
Initialize the coeffs object.
Definition: transext.cc:2636
#define omPrintStats(F)
Definition: xalloc.h:231
#define omInitInfo()
Definition: xalloc.h:228