235 const char *sys_cmd=(
char *)(args->
Data());
240 if (strcmp(sys_cmd,
"nblocks") == 0)
262 r = (ring)
h->Data();
269 if(strcmp(sys_cmd,
"version")==0)
277 if(strcmp(sys_cmd,
"alarm")==0)
284 struct itimerval t,o;
285 memset(&t,0,
sizeof(t));
286 t.it_value.tv_sec =(unsigned)((
unsigned long)
h->Data());
287 setitimer(ITIMER_VIRTUAL,&t,&o);
295 if(strcmp(sys_cmd,
"cpu")==0)
298 #ifdef _SC_NPROCESSORS_ONLN
299 cpu=sysconf(_SC_NPROCESSORS_ONLN);
300 #elif defined(_SC_NPROCESSORS_CONF)
301 cpu=sysconf(_SC_NPROCESSORS_CONF);
303 res->data=(
void *)cpu;
309 if(strcmp(sys_cmd,
"executable")==0)
324 if(strcmp(sys_cmd,
"flatten")==0)
337 if(strcmp(sys_cmd,
"unflatten")==0)
350 if(strcmp(sys_cmd,
"neworder")==0)
364 if(strcmp(sys_cmd,
"nc_hilb") == 0)
372 i = (ideal)
h->Data();
375 WerrorS(
"nc_Hilb:ideal expected");
380 lV = (int)(
long)
h->Data();
383 WerrorS(
"nc_Hilb:int expected");
389 if((
int)(
long)
h->Data() == 1)
391 else if((
int)(
long)
h->Data() == 2)
396 trunDegHs = (int)(
long)
h->Data();
401 WerrorS(
"nc_Hilb:int 1,2, total degree for the truncation, and a string for printing the details are expected");
410 if(strcmp(sys_cmd,
"rcolon") == 0)
415 ideal
i = (ideal)
h->Data();
417 poly
w=(poly)
h->Data();
419 int lV = (int)(
long)
h->Data();
430 if(strcmp(sys_cmd,
"sh")==0)
434 WerrorS(
"shell execution is disallowed in restricted mode");
440 res->data = (
void*)(
long)
system((
char*)(
h->Data()));
447 if(strcmp(sys_cmd,
"reduce_bound")==0)
457 p = (poly)
h->CopyD();
461 pid = (ideal)
h->CopyD();
466 ideal q = (ideal)
h->next->CopyD();
467 int bound = (int)(
long)
h->next->next->Data();
476 if(strcmp(sys_cmd,
"uname")==0)
484 if(strcmp(sys_cmd,
"with")==0)
494 #define TEST_FOR(A) if(strcmp(s,A)==0) res->data=(void *)1; else
495 char *
s=(
char *)
h->Data();
508 #ifdef TEST_MAC_ORDER
513 #ifdef HAVE_DYNAMIC_LOADING
536 if (strcmp(sys_cmd,
"browsers")==0)
546 if (strcmp(sys_cmd,
"pid")==0)
549 res->data=(
void *)(
long) getpid();
554 if (strcmp(sys_cmd,
"getenv")==0)
559 const char *r=
getenv((
char *)
h->Data());
572 if (strcmp(sys_cmd,
"setenv")==0)
579 setenv((
char *)
h->Data(), (
char *)
h->next->Data(), 1);
589 WerrorS(
"setenv not supported on this platform");
595 if (strcmp(sys_cmd,
"Singular") == 0)
604 if (strcmp(sys_cmd,
"SingularLib") == 0)
614 if (strstr(sys_cmd,
"--") == sys_cmd)
616 if (strcmp(sys_cmd,
"--") == 0)
624 Werror(
"Unknown option %s", sys_cmd);
625 WerrorS(
"Use 'system(\"--\");' for listing of available options");
640 const char *r=(
const char*)
feOptSpec[opt].value;
654 WerrorS(
"Need string or int argument to set option value");
657 const char* errormsg;
662 Werror(
"Need string argument to set value of option %s", sys_cmd);
666 if (errormsg !=
NULL)
667 Werror(
"Option '--%s=%d' %s", sys_cmd, (
int) ((
long)
h->Data()), errormsg);
672 if (errormsg !=
NULL)
673 Werror(
"Option '--%s=%s' %s", sys_cmd, (
char*)
h->Data(), errormsg);
680 if (strcmp(sys_cmd,
"HC")==0)
688 if(strcmp(sys_cmd,
"random")==0)
711 if (strcmp(sys_cmd,
"std_syz") == 0)
717 i1=(ideal)
h->CopyD();
723 i2=(int)((
long)
h->Data());
732 if (strcmp(sys_cmd,
"denom_list")==0)
741 if(strcmp(sys_cmd,
"complexNearZero")==0)
748 WerrorS(
"unsupported ground field!");
755 (int)((
long)(
h->next->Data())));
766 if(strcmp(sys_cmd,
"getPrecDigits")==0)
771 WerrorS(
"unsupported ground field!");
782 if(strcmp(sys_cmd,
"lduDecomp")==0)
789 poly
l; poly u; poly prodLU;
790 lduDecomp(aMat, pMat, lMat, dMat, uMat,
l, u, prodLU);
801 res->data = (
char *)L;
811 if(strcmp(sys_cmd,
"lduSolve")==0)
845 poly
l = (poly)
h->next->next->next->next->Data();
846 poly u = (poly)
h->next->next->next->next->next->Data();
847 poly lTimesU = (poly)
h->next->next->next->next->next->next->Data();
848 matrix bVec = (
matrix)
h->next->next->next->next->next->next->next->Data();
852 Werror(
"first matrix (%d x %d) is not quadratic",
858 Werror(
"second matrix (%d x %d) is not quadratic",
864 Werror(
"third matrix (%d x %d) is not quadratic",
870 Werror(
"third matrix (%d x %d) and fourth matrix (%d x %d) %s",
877 Werror(
"fourth matrix (%d x %d) and vector (%d x 1) do not fit",
882 bVec, xVec, homogSolSpace);
905 if (strcmp(sys_cmd,
"shared") == 0)
907 #ifndef SI_COUNTEDREF_AUTOLOAD
914 else if (strcmp(sys_cmd,
"reference") == 0)
916 #ifndef SI_COUNTEDREF_AUTOLOAD
925 #ifdef HAVE_SIMPLEIPC
926 if (strcmp(sys_cmd,
"semaphore")==0)
932 v=(int)(
long)
h->next->next->Data();
939 WerrorS(
"Usage: system(\"semaphore\",<cmd>,int)");
946 if (strcmp(sys_cmd,
"reserve")==0)
954 res->data=(
void*)(
long)
p;
961 if (strcmp(sys_cmd,
"reservedLink")==0)
970 if (strcmp(sys_cmd,
"install")==0)
976 (
int)(
long)
h->next->next->next->Data(),
983 if (strcmp(sys_cmd,
"newstruct")==0)
989 char *n=(
char*)
h->Data();
996 newstruct_desc desc=(newstruct_desc)bb->data;
1000 else Werror(
"'%s' is not a newstruct",n);
1002 else Werror(
"'%s' is not a blackbox object",n);
1008 if (strcmp(sys_cmd,
"blackbox")==0)
1016 if (strcmp(sys_cmd,
"absFact") == 0)
1033 l->m[0].data=(
void *)
f;
1035 l->m[1].data=(
void *)
v;
1037 l->m[2].data=(
void*) mipos;
1039 l->m[3].data=(
void*) (
long) n;
1040 res->data=(
void *)
l;
1049 if (strcmp(sys_cmd,
"LLL") == 0)
1072 #if __FLINT_RELEASE >= 20500
1073 if (strcmp(sys_cmd,
"LLL_Flint") == 0)
1094 WerrorS(
"matrix,int or bigint,int expected");
1099 if(((
int)((
long)(
h->next->Data())) != 0) && (
int)((
long)(
h->next->Data()) != 1))
1101 WerrorS(
"int is different from 0, 1");
1105 if((
long)(
h->next->Data()) == 0)
1120 if((
long)(
h->next->Data()) == 1)
1126 for(
int i = 1;
i<=
m->rows();
i++)
1131 m = singflint_LLL(
m,
T);
1144 for(
int i = 1;
i<=
m->rows();
i++)
1146 m = singflint_LLL(
m,
T);
1167 if(strcmp(sys_cmd,
"pcvLAddL")==0)
1172 if(strcmp(sys_cmd,
"pcvPMulL")==0)
1177 if(strcmp(sys_cmd,
"pcvMinDeg")==0)
1182 if(strcmp(sys_cmd,
"pcvP2CV")==0)
1187 if(strcmp(sys_cmd,
"pcvCV2P")==0)
1192 if(strcmp(sys_cmd,
"pcvDim")==0)
1197 if(strcmp(sys_cmd,
"pcvBasis")==0)
1204 #ifdef HAVE_EIGENVAL
1205 if(strcmp(sys_cmd,
"hessenberg")==0)
1212 #ifdef HAVE_EIGENVAL
1213 if(strcmp(sys_cmd,
"eigenvals")==0)
1220 #ifdef HAVE_EIGENVAL
1221 if(strcmp(sys_cmd,
"rowelim")==0)
1228 #ifdef HAVE_EIGENVAL
1229 if(strcmp(sys_cmd,
"rowcolswap")==0)
1237 if(strcmp(sys_cmd,
"gmsnf")==0)
1244 if(strcmp(sys_cmd,
"contributors") == 0)
1248 "Olaf Bachmann, Michael Brickenstein, Hubert Grassmann, Kai Krueger, Victor Levandovskyy, Wolfgang Neumann, Thomas Nuessler, Wilfred Pohl, Jens Schmidt, Mathias Schulze, Thomas Siebert, Ruediger Stobbe, Moritz Wenk, Tim Wichmann");
1253 #ifdef HAVE_SPECTRUM
1254 if(strcmp(sys_cmd,
"spectrum") == 0)
1268 if(((
long)
h->next->Data())==1L)
1274 if(strcmp(sys_cmd,
"semic") == 0)
1280 if (
h->next->next==
NULL)
1282 else if (
h->next->next->Typ()==
INT_CMD)
1289 if(strcmp(sys_cmd,
"spadd") == 0)
1300 if(strcmp(sys_cmd,
"spmul") == 0)
1312 #define HAVE_SHEAFCOH_TRICKS 1
1314 #ifdef HAVE_SHEAFCOH_TRICKS
1315 if(strcmp(sys_cmd,
"tensorModuleMult")==0)
1321 int m = (int)( (
long)
h->Data() );
1322 ideal
M = (ideal)
h->next->Data();
1333 if (strcmp(sys_cmd,
"twostd") == 0)
1338 I=(ideal)
h->CopyD();
1352 if (strcmp(sys_cmd,
"bracket") == 0)
1357 poly
p=(poly)
h->CopyD();
1359 poly q=(poly)
h->Data();
1370 if (strcmp(sys_cmd,
"env")==0)
1374 ring r = (ring)
h->Data();
1381 WerrorS(
"`system(\"env\",<ring>)` expected");
1389 if (strcmp(sys_cmd,
"opp")==0)
1393 ring r=(ring)
h->Data();
1400 WerrorS(
"`system(\"opp\",<ring>)` expected");
1408 if (strcmp(sys_cmd,
"oppose")==0)
1411 && (
h->next!=
NULL))
1413 ring Rop = (ring)
h->Data();
1426 WerrorS(
"`system(\"oppose\",<ring>,<poly>)` expected");
1436 if (strcmp(sys_cmd,
"walkNextWeight") == 0)
1443 Werror(
"system(\"walkNextWeight\" ...) intvecs not of length %d\n",
1447 res->data = (
void*) walkNextWeight(((
intvec*)
h->Data()),
1449 (ideal)
h->next->next->Data());
1450 if (
res->data ==
NULL ||
res->data == (
void*) 1L)
1466 if (strcmp(sys_cmd,
"walkInitials") == 0)
1470 WerrorS(
"system(\"walkInitials\", ideal) expected");
1473 res->data = (
void*) walkInitials((ideal)
h->Data());
1483 if (strcmp(sys_cmd,
"walkAddIntVec") == 0)
1489 res->data = (
intvec*) walkAddIntVec(arg1, arg2);
1498 #ifdef MwaklNextWeight
1499 if (strcmp(sys_cmd,
"MwalkNextWeight") == 0)
1506 Werror(
"system(\"MwalkNextWeight\" ...) intvecs not of length %d\n",
1512 ideal arg3 = (ideal)
h->next->next->Data();
1523 if(strcmp(sys_cmd,
"Mivdp") == 0)
1527 WerrorS(
"system(\"Mivdp\", int) expected");
1530 if ((
int) ((
long)(
h->Data())) !=
currRing->N)
1532 Werror(
"system(\"Mivdp\" ...) intvecs not of length %d\n",
1536 int arg1 = (int) ((
long)(
h->Data()));
1546 if(strcmp(sys_cmd,
"Mivlp") == 0)
1550 WerrorS(
"system(\"Mivlp\", int) expected");
1553 if ((
int) ((
long)(
h->Data())) !=
currRing->N)
1555 Werror(
"system(\"Mivlp\" ...) intvecs not of length %d\n",
1559 int arg1 = (int) ((
long)(
h->Data()));
1570 if(strcmp(sys_cmd,
"MpDiv") == 0)
1574 poly arg1 = (poly)
h->Data();
1575 poly arg2 = (poly)
h->next->Data();
1576 poly
result = MpDiv(arg1, arg2);
1587 if(strcmp(sys_cmd,
"MpMult") == 0)
1591 poly arg1 = (poly)
h->Data();
1592 poly arg2 = (poly)
h->next->Data();
1593 poly
result = MpMult(arg1, arg2);
1603 if (strcmp(sys_cmd,
"MivSame") == 0)
1624 res->data = (
void*)(
long)
MivSame(arg1, arg2);
1631 if (strcmp(sys_cmd,
"M3ivSame") == 0)
1654 res->data = (
void*)(
long)
M3ivSame(arg1, arg2, arg3);
1661 if(strcmp(sys_cmd,
"MwalkInitialForm") == 0)
1667 Werror(
"system \"MwalkInitialForm\"...) intvec not of length %d\n",
1671 ideal
id = (ideal)
h->Data();
1683 if(strcmp(sys_cmd,
"MivMatrixOrder") == 0)
1687 WerrorS(
"system(\"MivMatrixOrder\",intvec) expected");
1700 if(strcmp(sys_cmd,
"MivMatrixOrderdp") == 0)
1704 WerrorS(
"system(\"MivMatrixOrderdp\",intvec) expected");
1707 int arg1 = (int) ((
long)(
h->Data()));
1717 if(strcmp(sys_cmd,
"MPertVectors") == 0)
1721 ideal arg1 = (ideal)
h->Data();
1723 int arg3 = (int) ((
long)(
h->next->next->Data()));
1733 if(strcmp(sys_cmd,
"MPertVectorslp") == 0)
1737 ideal arg1 = (ideal)
h->Data();
1739 int arg3 = (int) ((
long)(
h->next->next->Data()));
1749 if(strcmp(sys_cmd,
"Mfpertvector") == 0)
1753 ideal arg1 = (ideal)
h->Data();
1764 if(strcmp(sys_cmd,
"MivUnit") == 0)
1768 int arg1 = (int) ((
long)(
h->Data()));
1778 if(strcmp(sys_cmd,
"MivWeightOrderlp") == 0)
1792 if(strcmp(sys_cmd,
"MivWeightOrderdp") == 0)
1796 WerrorS(
"system(\"MivWeightOrderdp\",intvec) expected");
1810 if(strcmp(sys_cmd,
"MivMatrixOrderlp") == 0)
1814 WerrorS(
"system(\"MivMatrixOrderlp\",int) expected");
1817 int arg1 = (int) ((
long)(
h->Data()));
1827 if (strcmp(sys_cmd,
"MkInterRedNextWeight") == 0)
1834 Werror(
"system(\"MkInterRedNextWeight\" ...) intvecs not of length %d\n",
1840 ideal arg3 = (ideal)
h->next->next->Data();
1850 #ifdef MPertNextWeight
1851 if (strcmp(sys_cmd,
"MPertNextWeight") == 0)
1857 Werror(
"system(\"MPertNextWeight\" ...) intvecs not of length %d\n",
1862 ideal arg2 = (ideal)
h->next->Data();
1863 int arg3 = (int)
h->next->next->Data();
1874 #ifdef Mivperttarget
1875 if (strcmp(sys_cmd,
"Mivperttarget") == 0)
1879 ideal arg1 = (ideal)
h->Data();
1880 int arg2 = (int)
h->next->Data();
1891 if (strcmp(sys_cmd,
"Mwalk") == 0)
1898 Werror(
"system(\"Mwalk\" ...) intvecs not of length %d\n",
1902 ideal arg1 = (ideal)
h->CopyD();
1905 ring arg4 = (ring)
h->next->next->next->Data();
1906 int arg5 = (int) (
long)
h->next->next->next->next->Data();
1907 int arg6 = (int) (
long)
h->next->next->next->next->next->Data();
1908 ideal
result = (ideal)
Mwalk(arg1, arg2, arg3, arg4, arg5, arg6);
1918 if (strcmp(sys_cmd,
"Mwalk") == 0)
1927 Werror(
"system(\"Mwalk\" ...) intvecs not of length %d or %d\n",
1931 ideal arg1 = (ideal)
h->Data();
1934 ring arg4 = (ring)
h->next->next->next->Data();
1935 ideal
result = (ideal)
Mwalk(arg1, arg2, arg3,arg4);
1942 if (strcmp(sys_cmd,
"Mpwalk") == 0)
1949 Werror(
"system(\"Mpwalk\" ...) intvecs not of length %d\n",
currRing->N);
1952 ideal arg1 = (ideal)
h->Data();
1953 int arg2 = (int) (
long)
h->next->Data();
1954 int arg3 = (int) (
long)
h->next->next->Data();
1957 int arg6 = (int) (
long)
h->next->next->next->next->next->Data();
1958 int arg7 = (int) (
long)
h->next->next->next->next->next->next->Data();
1959 int arg8 = (int) (
long)
h->next->next->next->next->next->next->next->Data();
1960 ideal
result = (ideal)
Mpwalk(arg1, arg2, arg3, arg4, arg5, arg6, arg7, arg8);
1970 if (strcmp(sys_cmd,
"Mrwalk") == 0)
1979 Werror(
"system(\"Mrwalk\" ...) intvecs not of length %d or %d\n",
1983 ideal arg1 = (ideal)
h->Data();
1986 int arg4 = (int)(
long)
h->next->next->next->Data();
1987 int arg5 = (int)(
long)
h->next->next->next->next->Data();
1988 int arg6 = (int)(
long)
h->next->next->next->next->next->Data();
1989 int arg7 = (int)(
long)
h->next->next->next->next->next->next->Data();
1990 ideal
result = (ideal)
Mrwalk(arg1, arg2, arg3, arg4, arg5, arg6, arg7);
1999 if (strcmp(sys_cmd,
"MAltwalk1") == 0)
2006 Werror(
"system(\"MAltwalk1\" ...) intvecs not of length %d\n",
2010 ideal arg1 = (ideal)
h->Data();
2011 int arg2 = (int) ((
long)(
h->next->Data()));
2012 int arg3 = (int) ((
long)(
h->next->next->Data()));
2025 if (strcmp(sys_cmd,
"Mfwalk_alt") == 0)
2032 Werror(
"system(\"Mfwalk\" ...) intvecs not of length %d\n",
2036 ideal arg1 = (ideal)
h->Data();
2039 int arg4 = (int)
h->next->next->next->Data();
2040 ideal
result = (ideal) Mfwalk_alt(arg1, arg2, arg3, arg4);
2050 if (strcmp(sys_cmd,
"Mfwalk") == 0)
2057 Werror(
"system(\"Mfwalk\" ...) intvecs not of length %d\n",
2061 ideal arg1 = (ideal)
h->Data();
2064 int arg4 = (int)(
long)
h->next->next->next->Data();
2065 int arg5 = (int)(
long)
h->next->next->next->next->Data();
2066 ideal
result = (ideal)
Mfwalk(arg1, arg2, arg3, arg4, arg5);
2075 if (strcmp(sys_cmd,
"Mfrwalk") == 0)
2092 Werror(
"system(\"Mfrwalk\" ...) intvecs not of length %d or %d\n",
2097 ideal arg1 = (ideal)
h->Data();
2100 int arg4 = (int)(
long)
h->next->next->next->Data();
2101 int arg5 = (int)(
long)
h->next->next->next->next->Data();
2102 int arg6 = (int)(
long)
h->next->next->next->next->next->Data();
2103 ideal
result = (ideal)
Mfrwalk(arg1, arg2, arg3, arg4, arg5, arg6);
2110 if (strcmp(sys_cmd,
"Mprwalk") == 0)
2119 Werror(
"system(\"Mrwalk\" ...) intvecs not of length %d or %d\n",
2123 ideal arg1 = (ideal)
h->Data();
2126 int arg4 = (int)(
long)
h->next->next->next->Data();
2127 int arg5 = (int)(
long)
h->next->next->next->next->Data();
2128 int arg6 = (int)(
long)
h->next->next->next->next->next->Data();
2129 int arg7 = (int)(
long)
h->next->next->next->next->next->next->Data();
2130 int arg8 = (int)(
long)
h->next->next->next->next->next->next->next->Data();
2131 int arg9 = (int)(
long)
h->next->next->next->next->next->next->next->next->Data();
2132 ideal
result = (ideal)
Mprwalk(arg1, arg2, arg3, arg4, arg5, arg6, arg7, arg8, arg9);
2142 if (strcmp(sys_cmd,
"TranMImprovwalk") == 0)
2149 Werror(
"system(\"TranMImprovwalk\" ...) intvecs not of length %d\n",
2153 ideal arg1 = (ideal)
h->Data();
2166 if (strcmp(sys_cmd,
"MAltwalk2") == 0)
2173 Werror(
"system(\"MAltwalk2\" ...) intvecs not of length %d\n",
2177 ideal arg1 = (ideal)
h->Data();
2189 if (strcmp(sys_cmd,
"TranMImprovwalk") == 0)
2196 Werror(
"system(\"TranMImprovwalk\" ...) intvecs not of length %d\n",
2200 ideal arg1 = (ideal)
h->Data();
2203 int arg4 = (int) ((
long)(
h->next->next->next->Data()));
2214 if (strcmp(sys_cmd,
"TranMrImprovwalk") == 0)
2219 h->next->next->next ==
NULL ||
h->next->next->next->Typ() !=
INT_CMD ||
2220 h->next->next->next ==
NULL ||
h->next->next->next->next->Typ() !=
INT_CMD ||
2221 h->next->next->next ==
NULL ||
h->next->next->next->next->next->Typ() !=
INT_CMD)
2223 WerrorS(
"system(\"TranMrImprovwalk\", ideal, intvec, intvec) expected");
2229 Werror(
"system(\"TranMrImprovwalk\" ...) intvecs not of length %d\n",
currRing->N);
2232 ideal arg1 = (ideal)
h->Data();
2235 int arg4 = (int)(
long)
h->next->next->next->Data();
2236 int arg5 = (int)(
long)
h->next->next->next->next->Data();
2237 int arg6 = (int)(
long)
h->next->next->next->next->next->Data();
2238 ideal
result = (ideal) TranMrImprovwalk(arg1, arg2, arg3, arg4, arg5, arg6);
2248 #ifndef MAKE_DISTRIBUTION
#define BIMATELEM(M, I, J)
blackbox * getBlackboxStuff(const int t)
return the structure to the type given by t
int blackboxIsCmd(const char *n, int &tok)
used by scanner: returns ROOT_DECL for known types (and the type number in tok)
void printBlackboxTypes()
list all defined type (for debugging)
static CanonicalForm bound(const CFMatrix &M)
void factoryseed(int s)
random seed initializer
char * singclap_neworder(ideal I, const ring r)
matrix singntl_LLL(matrix m, const ring s)
ideal singclap_absFactorize(poly f, ideal &mipos, intvec **exps, int &numFactors, const ring r)
gmp_complex numbers based on
@ n_transExt
used for all transcendental extensions, i.e., the top-most extension in an extension tower is transce...
static FORCE_INLINE n_coeffType getCoeffType(const coeffs r)
Returns the type of coeffs domain.
void countedref_reference_load()
Initialize blackbox types 'reference' and 'shared', or both.
void countedref_shared_load()
matrix evRowElim(matrix M, int i, int j, int k)
matrix evHessenberg(matrix M)
matrix evSwap(matrix M, int i, int j)
lists evEigenvals(matrix M)
const char * feSetOptValue(feOptIndex opt, char *optarg)
feOptIndex feGetOptIndex(const char *name)
static void * feOptValue(feOptIndex opt)
EXTERN_VAR struct fe_option feOptSpec[]
static char * feResource(feResourceConfig config, int warn)
void feStringAppendBrowsers(int warn)
lists gmsNF(ideal p, ideal g, matrix B, int D, int K)
void HilbertSeries_OrbitData(ideal S, int lV, bool IG_CASE, bool mgrad, bool odp, int trunDegHs)
ideal RightColonOperation(ideal S, poly w, int lV)
ideal idXXX(ideal h1, int k)
ideal id_TensorModuleMult(const int m, const ideal M, const ring rRing)
static BOOLEAN length(leftv result, leftv arg)
#define IMATELEM(M, I, J)
BOOLEAN spaddProc(leftv result, leftv first, leftv second)
BOOLEAN semicProc3(leftv res, leftv u, leftv v, leftv w)
BOOLEAN spectrumfProc(leftv result, leftv first)
BOOLEAN spmulProc(leftv result, leftv first, leftv second)
BOOLEAN spectrumProc(leftv result, leftv first)
BOOLEAN semicProc(leftv res, leftv u, leftv v)
poly kNFBound(ideal F, ideal Q, poly p, int bound, int syzComp, int lazyReduce)
poly pOppose(ring Rop_src, poly p, const ring Rop_dst)
opposes a vector p from Rop to currRing (dst!)
poly nc_p_Bracket_qq(poly p, const poly q, const ring r)
returns [p,q], destroys p
bool luSolveViaLDUDecomp(const matrix pMat, const matrix lMat, const matrix dMat, const matrix uMat, const poly l, const poly u, const poly lTimesU, const matrix bVec, matrix &xVec, matrix &H)
Solves the linear system A * x = b, where A is an (m x n)-matrix which is given by its LDU-decomposit...
void lduDecomp(const matrix aMat, matrix &pMat, matrix &lMat, matrix &dMat, matrix &uMat, poly &l, poly &u, poly &lTimesU)
LU-decomposition of a given (m x n)-matrix with performing only those divisions that yield zero remai...
ideal sm_UnFlatten(ideal a, int col, const ring R)
ideal sm_Flatten(ideal a, const ring R)
EXTERN_VAR size_t gmp_output_digits
bool complexNearZero(gmp_complex *c, int digits)
ideal twostd(ideal I)
Compute two-sided GB:
void newstructShow(newstruct_desc d)
BOOLEAN newstruct_set_proc(const char *bbname, const char *func, int args, procinfov pr)
char * omFindExec(const char *name, char *exec)
poly pcvP2CV(poly p, int d0, int d1)
int pcvBasis(lists b, int i, poly m, int d, int n)
int pcvDim(int d0, int d1)
lists pcvPMulL(poly p, lists l1)
poly pcvCV2P(poly cv, int d0, int d1)
lists pcvLAddL(lists l1, lists l2)
void StringSetS(const char *st)
static BOOLEAN rField_is_Ring(const ring r)
static int rBlocks(ring r)
static BOOLEAN rField_is_long_C(const ring r)
static BOOLEAN rField_is_long_R(const ring r)
int simpleipc_cmd(char *cmd, int id, int v)
int ssiReservePort(int clients)
int M3ivSame(intvec *temp, intvec *u, intvec *v)
intvec * MivMatrixOrderlp(int nV)
intvec * MivMatrixOrder(intvec *iv)
intvec * MkInterRedNextWeight(intvec *iva, intvec *ivb, ideal G)
intvec * Mfpertvector(ideal G, intvec *ivtarget)
ideal TranMImprovwalk(ideal G, intvec *curr_weight, intvec *target_tmp, int nP)
ideal Mfwalk(ideal G, intvec *ivstart, intvec *ivtarget, int reduction, int printout)
intvec * MPertVectors(ideal G, intvec *ivtarget, int pdeg)
intvec * MPertVectorslp(ideal G, intvec *ivtarget, int pdeg)
ideal Mprwalk(ideal Go, intvec *orig_M, intvec *target_M, int weight_rad, int op_deg, int tp_deg, int nP, int reduction, int printout)
intvec * MivWeightOrderdp(intvec *ivstart)
intvec * MivMatrixOrderdp(int nV)
intvec * MivWeightOrderlp(intvec *ivstart)
ideal MAltwalk2(ideal Go, intvec *curr_weight, intvec *target_weight)
ideal MAltwalk1(ideal Go, int op_deg, int tp_deg, intvec *curr_weight, intvec *target_weight)
ideal Mrwalk(ideal Go, intvec *orig_M, intvec *target_M, int weight_rad, int pert_deg, int reduction, int printout)
ideal Mfrwalk(ideal G, intvec *ivstart, intvec *ivtarget, int weight_rad, int reduction, int printout)
ideal Mwalk(ideal Go, intvec *orig_M, intvec *target_M, ring baseRing, int reduction, int printout)
ideal Mpwalk(ideal Go, int op_deg, int tp_deg, intvec *curr_weight, intvec *target_weight, int nP, int reduction, int printout)
int MivSame(intvec *u, intvec *v)
ideal MwalkInitialForm(ideal G, intvec *ivw)
intvec * MPertNextWeight(intvec *iva, ideal G, int deg)
intvec * MwalkNextWeight(intvec *curr_weight, intvec *target_weight, ideal G)
intvec * Mivperttarget(ideal G, int ndeg)