My Project
Functions
tesths.cc File Reference
#include "kernel/mod2.h"
#include "omalloc/omalloc.h"
#include "misc/options.h"
#include "factory/factory.h"
#include "kernel/oswrapper/feread.h"
#include "Singular/fevoices.h"
#include "kernel/oswrapper/timer.h"
#include "ipshell.h"
#include "cntrlc.h"
#include "links/silink.h"
#include "ipid.h"
#include "sdb.h"
#include "feOpt.h"
#include "distrib.h"
#include "mmalloc.h"
#include "tok.h"
#include "fegetopt.h"
#include "Singular/countedref.h"
#include "Singular/pyobject_setup.h"
#include <unistd.h>
#include <NTL/config.h>

Go to the source code of this file.

Functions

int siInit (char *)
 
int mmInit (void)
 
int main (int argc, char **argv)
 

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 59 of file tesths.cc.

62 {
63  mmInit();
64  // Don't worry: ifdef OM_NDEBUG, then all these calls are undef'ed
65  omInitRet_2_Info(argv[0]);
67 
68 #ifdef PSINGULAR
69  global_argv0 = argv[0];
70 #endif
71  siInit(argv[0]);
72  init_signals();
73  #ifdef HAVE_NTL
74  #if NTL_MAJOR_VERSION>=10
75  #ifdef NTL_THREAD_BOOST
76  SetNumThreads(feOptValue(FE_OPT_CPUS));
77  #endif
78  #endif
79  #endif
80 
81  // parse command line options
82  int optc, option_index;
83  const char* errormsg;
84  while((optc = fe_getopt_long(argc, argv,
85  SHORT_OPTS_STRING, feOptSpec, &option_index))
86  != EOF)
87  {
88  if (optc == '?' || optc == 0)
89  {
90  fprintf(stderr, "Use '%s --help' for a complete list of options\n", feArgv0);
91  exit(1);
92  }
93 
94  if (optc != LONG_OPTION_RETURN)
95  option_index = feGetOptIndex(optc);
96 
97  assume(option_index >= 0 && option_index < (int) FE_OPT_UNDEF);
98 
99  if (fe_optarg == NULL &&
100  (feOptSpec[option_index].type == feOptBool ||
101  feOptSpec[option_index].has_arg == optional_argument))
102  errormsg = feSetOptValue((feOptIndex) option_index, (int) 1);
103  else
104  errormsg = feSetOptValue((feOptIndex) option_index, fe_optarg);
105 
106  if (errormsg)
107  {
108  if (fe_optarg == NULL)
109  fprintf(stderr, "Error: Option '--%s' %s\n",
110  feOptSpec[option_index].name, errormsg);
111  else
112  fprintf(stderr, "Error: Option '--%s=%s' %s\n",
113  feOptSpec[option_index].name, fe_optarg, errormsg);
114  fprintf(stderr, "Use '%s --help' for a complete list of options\n", feArgv0);
115  exit(1);
116  }
117  if (optc == 'h') exit(0);
118  switch(option_index)
119  {
120  case FE_OPT_DUMP_VERSIONTUPLE:
121  exit(0);
122  break;
123  default: ;
124  }
125  }
126 
127  /* say hello */
128 
129  if (TEST_V_QUIET)
130  {
131  (printf)(
132 " SINGULAR /"
133 #ifndef MAKE_DISTRIBUTION
134 " Development"
135 #endif
136 "\n"
137 " A Computer Algebra System for Polynomial Computations / version %s\n"
138 " 0<\n"
139 " by: W. Decker, G.-M. Greuel, G. Pfister, H. Schoenemann \\ %s\n"
140 "FB Mathematik der Universitaet, D-67653 Kaiserslautern \\\n"
141 , VERSION, VERSION_DATE);
142  if (feOptValue(FE_OPT_NO_SHELL)) WarnS("running in restricted mode:"
143  " shell invocation and links are disallowed");
144  }
145  else
146  {
147  if (feOptValue(FE_OPT_SORT)) On(SW_USE_NTL_SORT);
148  dup2(1,2);
149  /* alternative:
150  * memcpy(stderr,stdout,sizeof(FILE));
151  */
152  }
153 
154 #ifdef SINGULAR_PYOBJECT_SETUP_H
155  pyobject_setup();
156 #endif
157 #ifdef SI_COUNTEDREF_AUTOLOAD
158  countedref_init();
159 #endif
160  errorreported = 0;
161 
162  // -- example for "static" modules ------
163  //load_builtin("huhu.so",FALSE,(SModulFunc_t)huhu_mod_init);
164  //module_help_main("huhu.so","Help for huhu\nhaha\n");
165  //module_help_proc("huhu.so","p","Help for huhu::p\nhaha\n");
166  setjmp(si_start_jmpbuf);
167 
168  // Now, put things on the stack of stuff to do
169  // Last thing to do is to execute given scripts
170  if (fe_optind < argc)
171  {
172  int i = argc - 1;
173  FILE *fd;
174  while (i >= fe_optind)
175  {
176  if ((fd = feFopen(argv[i], "r")) == NULL)
177  {
178  Warn("Can not open %s", argv[i]);
179  }
180  else
181  {
182  fclose(fd);
183  newFile(argv[i]);
184  }
185  i--;
186  }
187  }
188  else
189  {
191  }
192 
193  // before scripts, we execute -c, if it was given
194  if (feOptValue(FE_OPT_EXECUTE) != NULL)
195  newBuffer(omStrDup((char*) feOptValue(FE_OPT_EXECUTE)), BT_execute);
196 
197  // first thing, however, is to load .singularrc from Singularpath
198  // and cwd/$HOME (in that order).
199  if (! feOptValue(FE_OPT_NO_RC))
200  {
201  char buf[MAXPATHLEN];
202  FILE * rc = feFopen("." DIR_SEPP ".singularrc", "r", buf);
203  if (rc == NULL) rc = feFopen("~" DIR_SEPP ".singularrc", "r", buf);
204  if (rc == NULL) rc = feFopen(".singularrc", "r", buf);
205 
206  if (rc != NULL)
207  {
208  if (BVERBOSE(V_LOAD_LIB))
209  Print("// ** executing %s\n", buf);
210  fclose(rc);
211  newFile(buf);
212  }
213  }
214 
215  /* start shell */
217  {
219  char *linkname=(char*) feOptValue(FE_OPT_LINK);
220  if((linkname!=NULL)&&(strcmp(linkname,"ssi")==0))
221  {
222  return ssiBatch((char*) feOptValue(FE_OPT_MPHOST),(char*) feOptValue(FE_OPT_MPPORT));
223  //Print("batch: p:%s, h:%s\n",(char*) feOptValue(FE_OPT_MPPORT),(char*) feOptValue(FE_OPT_MPHOST));
224  //exit(0);
225  }
226  else
227  {
228  Print("** missing arguments: -b requires --link/--MPhost/--MPport\n");
229  return 1;
230  }
231  }
232  setjmp(si_start_jmpbuf);
233  yyparse();
234  m2_end(0);
235  return 0;
236 }
#define TRUE
Definition: auxiliary.h:100
void On(int sw)
switches
int i
Definition: cfEzgcd.cc:132
static const int SW_USE_NTL_SORT
set to 1 to sort factors in a factorization
Definition: cf_defs.h:38
VAR jmp_buf si_start_jmpbuf
Definition: cntrlc.cc:100
VAR BOOLEAN singular_in_batchmode
Definition: cntrlc.cc:70
void init_signals()
init signal handlers and error handling for libraries: NTL, factory
Definition: cntrlc.cc:535
void countedref_init()
Definition: countedref.h:458
#define Print
Definition: emacs.cc:80
#define Warn
Definition: emacs.cc:77
#define WarnS
Definition: emacs.cc:78
char name(const Variable &v)
Definition: factory.h:196
#define VERSION
Definition: factoryconf.h:282
VAR short errorreported
Definition: feFopen.cc:23
FILE * feFopen(const char *path, const char *mode, char *where, short useWerror, short path_only)
Definition: feFopen.cc:47
feOptIndex
Definition: feOptGen.h:15
@ FE_OPT_UNDEF
Definition: feOptGen.h:15
#define LONG_OPTION_RETURN
Definition: feOptTab.h:4
const char * feSetOptValue(feOptIndex opt, char *optarg)
Definition: feOpt.cc:154
const char SHORT_OPTS_STRING[]
Definition: feOpt.cc:29
feOptIndex feGetOptIndex(const char *name)
Definition: feOpt.cc:104
static void * feOptValue(feOptIndex opt)
Definition: feOpt.h:40
EXTERN_VAR struct fe_option feOptSpec[]
Definition: feOpt.h:17
VAR char * feArgv0
Definition: feResource.cc:19
#define DIR_SEPP
Definition: feResource.h:7
int fe_getopt_long(int argc, char *const *argv, const char *options, const struct fe_option *long_options, int *opt_index)
Definition: fegetopt.c:666
VAR int fe_optind
Definition: fegetopt.c:111
VAR char * fe_optarg
Definition: fegetopt.c:96
@ feOptBool
Definition: fegetopt.h:77
#define optional_argument
Definition: fegetopt.h:101
char *(* fe_fgets_stdin)(const char *pr, char *s, int size)
Definition: feread.cc:34
char * fe_fgets_dummy(const char *, char *, int)
Definition: feread.cc:451
void newBuffer(char *s, feBufferTypes t, procinfo *pi, int lineno)
Definition: fevoices.cc:164
VAR Voice * currentVoice
Definition: fevoices.cc:47
BOOLEAN newFile(char *fname)
Definition: fevoices.cc:127
Voice * feInitStdin(Voice *pp)
Definition: fevoices.cc:669
@ BT_execute
Definition: fevoices.h:23
int yyparse(void)
Definition: grammar.cc:2111
#define assume(x)
Definition: mod2.h:387
void m2_end(int i)
Definition: misc_ip.cc:1098
#define omStrDup(s)
Definition: omAllocDecl.h:263
#define omInitGetBackTrace()
#define NULL
Definition: omList.c:12
#define MAXPATHLEN
Definition: omRet2Info.c:22
#define BVERBOSE(a)
Definition: options.h:34
#define TEST_V_QUIET
Definition: options.h:132
#define V_LOAD_LIB
Definition: options.h:46
void pyobject_setup()
initialize blackbox support for pyobject; functionilty os autoloaded on demand
int status int void * buf
Definition: si_signals.h:59
int status int fd
Definition: si_signals.h:59
int siInit(char *)
Definition: misc_ip.cc:1358
int mmInit(void)
Definition: tesths.cc:43
#define omInitRet_2_Info(argv0)
Definition: xalloc.h:322

◆ mmInit()

int mmInit ( void  )

Definition at line 43 of file tesths.cc.

44 {
45 #ifndef X_OMALLOC
46 #if defined(OMALLOC_USES_MALLOC)
47  /* in mmstd.c, for some architectures freeSize() unconditionally uses the *system* free() */
48  /* sage ticket 5344: http://trac.sagemath.org/sage_trac/ticket/5344 */
49  /* do not rely on the default in Singular as libsingular may be different */
50  mp_set_memory_functions(omMallocFunc,omReallocSizeFunc,omFreeSizeFunc);
51 #else
52  mp_set_memory_functions(malloc,reallocSize,freeSize);
53 #endif
54 #endif
55  return 1;
56 }
#define freeSize
Definition: omAllocFunc.c:15
#define reallocSize
Definition: omAllocFunc.c:17
void * omMallocFunc(size_t size)
void * omReallocSizeFunc(void *old_addr, size_t old_size, size_t new_size)
void omFreeSizeFunc(void *addr, size_t size)
void * malloc(size_t size)
Definition: omalloc.c:92

◆ siInit()

int siInit ( char *  name)

Definition at line 1358 of file misc_ip.cc.

1359 {
1360 // memory initialization: -----------------------------------------------
1361  om_Opts.OutOfMemoryFunc = omSingOutOfMemoryFunc;
1362 #ifndef OM_NDEBUG
1363 #ifndef __OPTIMIZE__
1364  om_Opts.ErrorHook = dErrorBreak;
1365 #else
1366  om_Opts.Keep = 0; /* !OM_NDEBUG, __OPTIMIZE__*/
1367 #endif
1368 #else
1369  om_Opts.Keep = 0; /* OM_NDEBUG */
1370 #endif
1371  omInitInfo();
1372 // factory
1373 #ifndef HAVE_NTL
1374  extern void initPT();
1375  initPT();
1376 #endif
1377 // options ---------------------------------------------------------------
1378  si_opt_1=0;
1379 // interpreter tables etc.: -----------------------------------------------
1380  memset(&sLastPrinted,0,sizeof(sleftv));
1382 
1383  extern int iiInitArithmetic(); iiInitArithmetic(); // iparith.cc
1384 
1385  basePack=(package)omAlloc0(sizeof(*basePack));
1387  idhdl h;
1388  h=enterid("Top", 0, PACKAGE_CMD, &IDROOT, FALSE);
1389  IDPACKAGE(h)=basePack;
1390  IDPACKAGE(h)->language = LANG_TOP;
1391  currPackHdl=h;
1392  basePackHdl=h;
1393 
1394  coeffs_BIGINT = nInitChar(n_Q,(void*)1);
1395 
1396 #if 1
1397  // def HAVE_POLYEXTENSIONS
1398  if(TRUE)
1399  {
1400  n_coeffType type;
1401  #ifdef SINGULAR_4_2
1402  type = nRegister(n_polyExt, n2pInitChar);
1403  assume(type == n_polyExt);
1404  #endif
1405 
1406  type = nRegister(n_algExt, naInitChar);
1407  assume(type == n_algExt);
1408 
1409  type = nRegister(n_transExt, ntInitChar);
1410  assume(type == n_transExt);
1411 
1412  (void)type;
1413  }
1414 #endif
1415 
1416 // random generator: -----------------------------------------------
1417  int t=initTimer();
1418  if (t==0) t=1;
1419  initRTimer();
1420  siSeed=t;
1421  factoryseed(t);
1422  siRandomStart=t;
1423  feOptSpec[FE_OPT_RANDOM].value = (void*) ((long)siRandomStart);
1424 
1425 // ressource table: ----------------------------------------------------
1426  // Don't worry: ifdef OM_NDEBUG, then all these calls are undef'ed
1427  // hack such that all shared' libs in the bindir are loaded correctly
1429 
1430 // singular links: --------------------------------------------------
1431  slStandardInit();
1432  myynest=0;
1433 // how many processes ? -----------------------------------------------------
1434  int cpus=2;
1435  int cpu_n;
1436  #ifdef _SC_NPROCESSORS_ONLN
1437  if ((cpu_n=sysconf(_SC_NPROCESSORS_ONLN))>cpus) cpus=cpu_n;
1438  #elif defined(_SC_NPROCESSORS_CONF)
1439  if ((cpu_n=sysconf(_SC_NPROCESSORS_CONF))>cpus) cpus=cpu_n;
1440  #endif
1441  feSetOptValue(FE_OPT_CPUS, cpus);
1442 // how many threads ? -----------------------------------------------------
1443  feSetOptValue(FE_OPT_THREADS, cpus);
1444 
1445 // default coeffs
1446  {
1447  idhdl h;
1448  h=enterid("QQ",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1449  IDDATA(h)=(char*)nInitChar(n_Q,NULL);
1450  h=enterid("ZZ",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1451  IDDATA(h)=(char*)nInitChar(n_Z,NULL);
1452  nRegisterCfByName(nrnInitCfByName,n_Zn); // and n_Znm
1453  iiAddCproc("kernel","crossprod",FALSE,iiCrossProd);
1454  iiAddCproc("kernel","Float",FALSE,iiFloat);
1455  //h=enterid("RR",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1456  //IDDATA(h)=(char*)nInitChar(n_R,NULL);
1457  //h=enterid("CC",0/*level*/, CRING_CMD,&(basePack->idroot),FALSE /*init*/,FALSE /*search*/);
1458  //IDDATA(h)=(char*)nInitChar(n_long_C,NULL);
1459  }
1460 // setting routines for PLURAL QRINGS:
1461 // allowing to use libpolys without libSingular(kStd)
1462 #ifdef HAVE_PLURAL
1463  nc_NF=k_NF;
1469 #endif
1470 // loading standard.lib -----------------------------------------------
1471  if (! feOptValue(FE_OPT_NO_STDLIB))
1472  {
1473  BITSET save1,save2;
1474  SI_SAVE_OPT(save1,save2);
1475  si_opt_2 &= ~Sy_bit(V_LOAD_LIB);
1476  iiLibCmd("standard.lib", TRUE,TRUE,TRUE);
1477  SI_RESTORE_OPT(save1,save2);
1478  }
1479  // interpreter error handling
1480  #ifndef __CYGWIN__
1481  factoryError=callWerrorS; // to honour later changes of variable WerrorS
1482  #endif
1483  errorreported = 0;
1484 }
BOOLEAN naInitChar(coeffs cf, void *infoStruct)
Initialize the coeffs object.
Definition: algext.cc:1373
BOOLEAN n2pInitChar(coeffs cf, void *infoStruct)
Definition: algext.cc:1627
#define FALSE
Definition: auxiliary.h:96
void initPT()
void factoryseed(int s)
random seed initializer
Definition: cf_random.cc:189
VAR void(* factoryError)(const char *s)
Definition: cf_util.cc:80
Definition: idrec.h:35
Class used for (list of) interpreter objects.
Definition: subexpr.h:83
int rtyp
Definition: subexpr.h:91
VAR int siRandomStart
Definition: cntrlc.cc:101
n_coeffType
Definition: coeffs.h:28
@ n_polyExt
used to represent polys as coeffcients
Definition: coeffs.h:35
@ n_Q
rational (GMP) numbers
Definition: coeffs.h:31
@ n_algExt
used for all algebraic extensions, i.e., the top-most extension in an extension tower is algebraic
Definition: coeffs.h:36
@ n_Zn
only used if HAVE_RINGS is defined
Definition: coeffs.h:45
@ n_transExt
used for all transcendental extensions, i.e., the top-most extension in an extension tower is transce...
Definition: coeffs.h:39
@ n_Z
only used if HAVE_RINGS is defined
Definition: coeffs.h:44
coeffs nInitChar(n_coeffType t, void *parameter)
one-time initialisations for new coeffs in case of an error return NULL
Definition: numbers.cc:353
void feInitResources(const char *argv0)
Definition: feResource.cc:170
VAR int myynest
Definition: febase.cc:41
void * value
Definition: fegetopt.h:93
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
ideal k_gnc_gr_mora(const ideal F, const ideal Q, const intvec *, const intvec *, kStrategy strat, const ring _currRing)
Definition: gr_kstd2.cc:1289
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
int iiInitArithmetic()
initialisation of arithmetic structured data
Definition: iparith.cc:9662
idhdl enterid(const char *s, int lev, int t, idhdl *root, BOOLEAN init, BOOLEAN search)
Definition: ipid.cc:278
VAR package basePack
Definition: ipid.cc:58
VAR package currPack
Definition: ipid.cc:57
VAR idhdl currPackHdl
Definition: ipid.cc:55
VAR idhdl basePackHdl
Definition: ipid.cc:56
VAR coeffs coeffs_BIGINT
Definition: ipid.cc:50
#define IDDATA(a)
Definition: ipid.h:126
#define IDROOT
Definition: ipid.h:19
#define IDPACKAGE(a)
Definition: ipid.h:139
int iiAddCproc(const char *libname, const char *procname, BOOLEAN pstatic, BOOLEAN(*func)(leftv res, leftv v))
Definition: iplib.cc:1059
BOOLEAN iiLibCmd(const char *newlib, BOOLEAN autoexport, BOOLEAN tellerror, BOOLEAN force)
Definition: iplib.cc:880
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:3369
static BOOLEAN iiCrossProd(leftv res, leftv args)
Definition: misc_ip.cc:1323
void omSingOutOfMemoryFunc()
Definition: misc_ip.cc:1196
static BOOLEAN iiFloat(leftv res, leftv pnn)
Definition: misc_ip.cc:1294
static void callWerrorS(const char *s)
Definition: misc_ip.cc:1357
void dErrorBreak()
Definition: dError.cc:139
void nRegisterCfByName(cfInitCfByNameProc p, n_coeffType n)
Definition: numbers.cc:582
n_coeffType nRegister(n_coeffType n, cfInitCharProc p)
Definition: numbers.cc:541
#define omAlloc0(size)
Definition: omAllocDecl.h:211
omOpts_t om_Opts
Definition: omOpts.c:13
VAR unsigned si_opt_2
Definition: options.c:6
VAR unsigned si_opt_1
Definition: options.c:5
#define SI_SAVE_OPT(A, B)
Definition: options.h:20
#define Sy_bit(x)
Definition: options.h:31
#define SI_RESTORE_OPT(A, B)
Definition: options.h:23
coeffs nrnInitCfByName(char *s, n_coeffType)
Definition: rmodulon.cc:35
VAR int siSeed
Definition: sirandom.c:30
ip_package * package
Definition: structs.h:48
#define BITSET
Definition: structs.h:20
INST_VAR sleftv sLastPrinted
Definition: subexpr.cc:46
@ LANG_TOP
Definition: subexpr.h:22
int initTimer()
Definition: timer.cc:67
void initRTimer()
Definition: timer.cc:156
@ CRING_CMD
Definition: tok.h:56
@ PACKAGE_CMD
Definition: tok.h:149
#define NONE
Definition: tok.h:221
BOOLEAN ntInitChar(coeffs cf, void *infoStruct)
Initialize the coeffs object.
Definition: transext.cc:2544
#define omInitInfo()
Definition: xalloc.h:272