source: git/Singular/iparith.cc @ f7920ef

spielwiese
Last change on this file since f7920ef was f7920ef, checked in by Hans Schönemann <hannes@…>, 25 years ago
*hannes: fixed sign bug in HIGHCORNER_M git-svn-id: file:///usr/local/Singular/svn/trunk@2913 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 164.8 KB
Line 
1/****************************************
2*  Computer Algebra System SINGULAR     *
3****************************************/
4/* $Id: iparith.cc,v 1.136 1999-03-15 09:20:59 Singular Exp $ */
5
6/*
7* ABSTRACT: table driven kernel interface, used by interpreter
8*/
9
10#include <stdlib.h>
11#include <string.h>
12#include <ctype.h>
13#include <stdio.h>
14#include <time.h>
15#include <unistd.h>
16
17#include "mod2.h"
18#include "tok.h"
19#include "ipid.h"
20#include "intvec.h"
21#include "mmemory.h"
22#include "febase.h"
23#include "longalg.h"
24#include "polys.h"
25#include "ideals.h"
26#include "matpol.h"
27#include "kstd1.h"
28#include "timer.h"
29#include "ring.h"
30#include "subexpr.h"
31#include "lists.h"
32#include "numbers.h"
33#include "stairc.h"
34#include "maps.h"
35#include "syz.h"
36#include "weight.h"
37#include "ipconv.h"
38#include "ipprint.h"
39#include "sing_dld.h"
40#include "attrib.h"
41#include "silink.h"
42#include "sparsmat.h"
43#ifdef HAVE_FACTORY
44#include "clapsing.h"
45#include "kstdfac.h"
46#endif
47#ifdef HAVE_FGLM
48#include "fglm.h"
49#endif
50#include "ipshell.h"
51
52/*=============== types =====================*/
53struct sValCmdTab
54{
55  short cmd;
56  short start;
57};
58
59typedef sValCmdTab jjValCmdTab[];
60
61/* ifdef GENTABLE: definitions are in ipshell.h */
62#ifndef GENTABLE
63typedef char * (*Proc1)(char *);
64struct sValCmd1
65{
66  proc1 p;
67  short cmd;
68  short res;
69  short arg;
70};
71
72typedef BOOLEAN (*proc2)(leftv,leftv,leftv);
73struct sValCmd2
74{
75  proc2 p;
76  short cmd;
77  short res;
78  short arg1;
79  short arg2;
80};
81
82typedef BOOLEAN (*proc3)(leftv,leftv,leftv,leftv);
83struct sValCmd3
84{
85  proc3 p;
86  short cmd;
87  short res;
88  short arg1;
89  short arg2;
90  short arg3;
91};
92struct sValCmdM
93{
94  proc1 p;
95  short cmd;
96  short res;
97  short number_of_args; /* -1: any, -2: any >0, .. */
98};
99#endif
100
101/*============= proc =======================*/
102static BOOLEAN jjLOAD(leftv res, leftv v, BOOLEAN autoexport = FALSE);
103static int iiTabIndex(const jjValCmdTab dArithTab, const int len, const int op);
104#ifdef MDEBUG
105#define jjMakeSub(A) jjDBMakeSub(A,__FILE__,__LINE__)
106static Subexpr jjDBMakeSub(leftv e,char *f, int l);
107#else
108static Subexpr jjMakeSub(leftv e);
109#endif
110
111/*============= vars ======================*/
112extern int cmdtok;
113extern BOOLEAN expected_parms;
114
115int iiOp; /* the current operation*/
116
117#ifdef GENTABLE
118cmdnames cmds[] =
119{  // name-string alias tokval          toktype
120  { "$INVALID$",   0, -1,                 0},
121  { "and",         0, '&' ,               '&'},
122  { "attrib",      0, ATTRIB_CMD ,        CMD_123},
123  { "bareiss",     0, BAREISS_CMD ,       CMD_1},
124  { "betti",       0, BETTI_CMD ,         CMD_1},
125  { "break",       0, BREAK_CMD ,         BREAK_CMD},
126  { "char",        0, CHARACTERISTIC_CMD ,CMD_1},
127  { "char_series", 0, CHAR_SERIES_CMD ,   CMD_1},
128  { "charstr",     0, CHARSTR_CMD ,       CMD_1},
129  { "cleardenom",  0, CONTENT_CMD ,       CMD_1},
130  { "close",       0, CLOSE_CMD ,         CMD_1},
131  { "coef",        0, COEF_CMD ,          CMD_M},
132  { "coeffs",      0, COEFFS_CMD ,        CMD_23},
133  { "continue",    0, CONTINUE_CMD ,      CONTINUE_CMD},
134  { "contract",    0, CONTRACT_CMD ,      CMD_2},
135  { "dbprint",     0, DBPRINT_CMD ,       CMD_M},
136  { "def",         0, DEF_CMD ,           ROOT_DECL},
137  { "defined",     0, DEFINED_CMD ,       CMD_1},
138  { "deg",         0, DEG_CMD ,           CMD_12},
139  { "degree",      0, DEGREE_CMD ,        CMD_1},
140  { "delete",      0, DELETE_CMD ,        CMD_2},
141  { "det",         0, DET_CMD ,           CMD_1},
142  { "diff",        0, DIFF_CMD ,          CMD_2},
143  { "dim",         0, DIM_CMD ,           CMD_1},
144  { "div",         0, INTDIV ,            INTDIV},
145#ifdef DRING
146  { "dring",       0, DRING_CMD ,         DRING_CMD},
147#endif
148  { "dump",        0, DUMP_CMD,           CMD_1},
149  { "extgcd",      0, EXTGCD_CMD ,        CMD_2},
150  { "EXTGCD",      2, EXTGCD_CMD ,        CMD_2},
151  { "eliminate",   0, ELIMINATION_CMD,    CMD_23},
152  { "else",        0, ELSE_CMD ,          ELSE_CMD},
153  { "eval",        0, EVAL ,              EVAL},
154  { "example",     0, EXAMPLE_CMD ,       EXAMPLE_CMD},
155  { "execute",     0, EXECUTE_CMD ,       EXECUTE_CMD},
156  { "export",      0, EXPORT_CMD ,        EXPORT_CMD},
157#ifdef HAVE_NAMESPACES
158  { "exportto",    0, EXPORTTO_CMD ,      CMD_M},
159#endif
160  { "factorize",   0, FAC_CMD ,           CMD_12},
161  { "fetch",       0, FETCH_CMD ,         CMD_2},
162  { "fglm",        0, FGLM_CMD ,          CMD_2},
163  { "find",        0, FIND_CMD ,          CMD_23},
164  { "finduni",     0, FINDUNI_CMD,        CMD_1},
165  { "forif",       0, IF_CMD ,            IF_CMD},
166  { "freemodule",  0, FREEMODULE_CMD ,    CMD_1},
167  { "facstd",      0, FACSTD_CMD ,        CMD_12},
168  { "gen",         0, E_CMD ,             CMD_1},
169  { "getdump",     0, GETDUMP_CMD,        CMD_1},
170  { "gcd",         0, GCD_CMD ,           CMD_2},
171  { "GCD",         2, GCD_CMD ,           CMD_2},
172  { "hilb",        0, HILBERT_CMD ,       CMD_12},
173  { "highcorner",  0, HIGHCORNER_CMD,     CMD_1},
174  { "homog",       0, HOMOG_CMD ,         CMD_12},
175  { "ideal",       0, IDEAL_CMD ,         IDEAL_CMD},
176  { "if",          0, IF_CMD ,            IF_CMD},
177  { "imap",        0, IMAP_CMD ,          CMD_2},
178#ifdef HAVE_NAMESPACES
179  { "importfrom",  0, IMPORTFROM_CMD ,    CMD_M},
180#endif
181  { "indepSet",    0, INDEPSET_CMD ,      CMD_12},
182  { "insert",      0, INSERT_CMD ,        CMD_23},
183  { "int",         0, INT_CMD ,           ROOT_DECL},
184  { "interred",    0, INTERRED_CMD ,      CMD_1},
185  { "intersect",   0, INTERSECT_CMD ,     CMD_M},
186  { "intmat",      0, INTMAT_CMD ,        INTMAT_CMD},
187  { "intvec",      0, INTVEC_CMD ,        ROOT_DECL_LIST},
188  { "jacob",       0, JACOB_CMD ,         CMD_1},
189  { "jet",         0, JET_CMD ,           CMD_23},
190  { "kbase",       0, KBASE_CMD ,         CMD_12},
191  { "keepring",    0, KEEPRING_CMD ,      KEEPRING_CMD},
192  { "kill",        0, KILL_CMD ,          KILL_CMD},
193  { "killattrib",  0, KILLATTR_CMD ,      CMD_12},
194  { "koszul",      0, KOSZUL_CMD ,        CMD_23},
195  { "lead",        0, LEAD_CMD ,          CMD_1},
196  { "leadcoef",    0, LEADCOEF_CMD ,      CMD_1},
197  { "leadexp",     0, LEADEXP_CMD ,       CMD_1},
198  { "LIB",         0, LIB_CMD ,           SYSVAR},
199  { "lift",        0, LIFT_CMD ,          CMD_23},
200  { "liftstd",     0, LIFTSTD_CMD ,       CMD_2},
201  { "link",        0, LINK_CMD ,          ROOT_DECL},
202  { "listvar",     0, LISTVAR_CMD ,       LISTVAR_CMD},
203  { "list",        0, LIST_CMD ,          ROOT_DECL_LIST},
204  { "load",        0, LOAD_CMD ,          CMD_12},
205  { "lres",        0, LRES_CMD ,          CMD_2},
206  { "map",         0, MAP_CMD ,           RING_DECL},
207  { "matrix",      0, MATRIX_CMD ,        MATRIX_CMD},
208  { "maxideal",    0, MAXID_CMD ,         CMD_1},
209  { "memory",      0, MEMORY_CMD ,        CMD_1},
210  { "minbase",     0, MINBASE_CMD ,       CMD_1},
211  { "minor",       0, MINOR_CMD ,         CMD_2},
212  { "minres",      0, MINRES_CMD ,        CMD_1},
213  { "mod",         0, '%',                '%'},
214  { "module",      0, MODUL_CMD ,         MODUL_CMD},
215  { "modulo",      0, MODULO_CMD ,        CMD_2},
216  { "monitor",     0, MONITOR_CMD ,       CMD_12},
217  { "mult",        0, MULTIPLICITY_CMD ,  CMD_1},
218  #ifdef OLD_RES
219  { "mres",        0, MRES_CMD ,          CMD_23},
220  #else
221  { "mres",        0, MRES_CMD ,          CMD_2},
222  #endif
223  { "mstd",        0, MSTD_CMD ,          CMD_1},
224  { "nameof",      0, NAMEOF_CMD ,        CMD_1},
225  { "names",       0, NAMES_CMD ,         CMD_M},
226  { "ncols",       0, COLS_CMD ,          CMD_1},
227  { "not",         0, NOT ,               NOT},
228  { "npars",       0, NPARS_CMD ,         CMD_1},
229  #ifdef OLD_RES
230  { "nres",        0, RES_CMD ,           CMD_23},
231  #else
232  { "nres",        0, RES_CMD ,           CMD_2},
233  #endif
234  { "nrows",       0, ROWS_CMD ,          CMD_1},
235  { "number",      0, NUMBER_CMD ,        RING_DECL},
236  { "nvars",       0, NVARS_CMD ,         CMD_1},
237  { "open",        0, OPEN_CMD ,          CMD_1},
238  { "option",      0, OPTION_CMD ,        CMD_M},
239  { "or",          0, '|' ,               '|'},
240  { "ord",         0, ORD_CMD ,           CMD_1},
241  { "ordstr",      0, ORDSTR_CMD ,        CMD_1},
242  { "package",     0, PACKAGE_CMD ,       ROOT_DECL},
243  { "par",         0, PAR_CMD ,           CMD_1},
244  { "parameter",   0, PARAMETER ,         PARAMETER},
245  { "pardeg",      0, PARDEG_CMD ,        CMD_1},
246  { "parstr",      0, PARSTR_CMD ,        CMD_12},
247  { "poly",        0, POLY_CMD ,          RING_DECL},
248  { "preimage",    0, PREIMAGE_CMD ,      CMD_3},
249  { "prime",       0, PRIME_CMD ,         CMD_1},
250  { "print",       0, PRINT_CMD ,         CMD_12},
251  { "prune",       0, PRUNE_CMD ,         CMD_1},
252  { "proc",        0, PROC_CMD ,          PROC_CMD},
253  { "qhweight",    0, QHWEIGHT_CMD ,      CMD_1},
254  { "qring",       0, QRING_CMD ,         ROOT_DECL},
255  { "quote",       0, QUOTE ,             QUOTE},
256  { "quotient",    0, QUOTIENT_CMD ,      CMD_2},
257  { "random",      0, RANDOM_CMD ,        CMD_23},
258  { "read",        0, READ_CMD ,          CMD_12},
259  { "reduce",      0, REDUCE_CMD ,        CMD_23},
260  { "regularity",  0, REGULARITY_CMD ,    CMD_1},
261  { "reservedName",0, RESERVEDNAME_CMD ,  CMD_M},
262  { "resolution",  0, RESOLUTION_CMD ,    RING_DECL},
263  { "resultant",   0, RESULTANT_CMD,      CMD_3},
264  { "return",      0, RETURN ,            RETURN},
265  { "RETURN",      0, END_GRAMMAR ,       RETURN},
266  { "ring",        0, RING_CMD ,          RING_CMD},
267  { "rvar",        0, IS_RINGVAR ,        CMD_1},
268  { "setring",     0, SETRING_CMD ,       SETRING_CMD},
269  { "simplify",    0, SIMPLIFY_CMD ,      CMD_2},
270  { "size",        0, COUNT_CMD ,         CMD_1},
271  { "sortvec",     0, SORTVEC_CMD ,       CMD_1},
272  #ifdef OLD_RES
273  { "sres",        0, SRES_CMD ,          CMD_23},
274  #else
275  { "sres",        0, SRES_CMD ,          CMD_2},
276  #endif
277  { "status",      0, STATUS_CMD,         CMD_M},
278  { "std",         0, STD_CMD ,           CMD_123},
279  { "string",      0, STRING_CMD ,        ROOT_DECL_LIST},
280  { "subst",       0, SUBST_CMD ,         CMD_3},
281  { "system",      0, SYSTEM_CMD,         CMD_M},
282  { "syz",         0, SYZYGY_CMD ,        CMD_1},
283  { "test",        0, TEST_CMD ,          CMD_M},
284  { "trace",       0, TRACE_CMD ,         CMD_1},
285  { "transpose",   0, TRANSPOSE_CMD ,     CMD_1},
286  { "type",        0, TYPE_CMD ,          TYPE_CMD},
287  { "typeof",      0, TYPEOF_CMD ,        CMD_1},
288#ifdef HAVE_NAMESPACES
289  { "unload",      0, UNLOAD_CMD ,        CMD_M},
290#endif
291  { "var",         0, VAR_CMD ,           CMD_1},
292  { "varstr",      0, VARSTR_CMD ,        CMD_12},
293  { "vdim",        0, VDIM_CMD ,          CMD_1},
294  { "vector",      0, VECTOR_CMD ,        RING_DECL},
295  { "wedge",       0, WEDGE_CMD ,         CMD_2},
296  { "weight",      0, WEIGHT_CMD ,        CMD_1},
297  { "whileif",     0, IF_CMD ,            IF_CMD},
298  { "write",       0, WRITE_CMD ,         CMD_M},
299/* delete for next version:*/
300  { "IN",          1, LEAD_CMD ,          CMD_1},
301  { "NF",          1, REDUCE_CMD ,        CMD_23},
302  { "multiplicity",1, MULTIPLICITY_CMD ,  CMD_1},
303  { "verbose",     2, OPTION_CMD ,        CMD_M},
304//  { "rank",        1, ROWS_CMD ,          CMD_1},
305//  { "Current",     0, -1 ,                SYSVAR},
306//  { "Top",         0, -1 ,                SYSVAR},
307//  { "Up",          0, -1 ,                SYSVAR},
308
309/* set sys vars*/
310//#ifdef SRING
311  { "alternating", 0, VALTVARS ,          SYSVAR},
312//#endif
313  { "degBound",    0, VMAXDEG ,           SYSVAR},
314  { "echo",        0, VECHO ,             SYSVAR},
315  { "minpoly",     0, VMINPOLY ,          SYSVAR},
316  { "multBound",   0, VMAXMULT ,          SYSVAR},
317  { "noether",     0, VNOETHER ,          SYSVAR},
318  { "pagelength",  0, VPAGELENGTH ,       SYSVAR},
319  { "pagewidth",   0, VCOLMAX ,           SYSVAR},
320  { "printlevel",  0, VPRINTLEVEL ,       SYSVAR},
321  { "short",       0, VSHORTOUT ,         SYSVAR},
322  { "timer",       0, VTIMER ,            SYSVAR},
323  { "rtimer",      0, VRTIMER,            SYSVAR},
324  { "TRACE",       0, TRACE ,             SYSVAR},
325  { "voice",       0, VOICE ,             SYSVAR},
326
327/* other reserved words:scanner.l */
328  { "pause",       2, -1 ,             0},
329  { "while",       0, -1 ,             0},
330  { "for",         0, -1 ,             0},
331  { "help",        0, -1 ,             0},
332  { "newline",     0, -1 ,             0},
333  { "exit",        0, -1 ,             0},
334  { "quit",        0, -1 ,             0},
335/* end of list marker */
336  { NULL, 0, 0, 0}
337};
338#endif
339
340/*=================== operations with 2 args.: static proc =================*/
341static BOOLEAN jjOP_IV_I(leftv res, leftv u, leftv v)
342{
343  intvec* aa= (intvec *)u->CopyD(INTVEC_CMD);
344  int bb = (int)(v->Data());
345  if (errorreported) return TRUE;
346  switch (iiOp)
347  {
348    case '+': (*aa) += bb; break;
349    case '-': (*aa) -= bb; break;
350    case '*': (*aa) *= bb; break;
351    case '/':
352    case INTDIV: (*aa) /= bb; break;
353    case '%': (*aa) %= bb; break;
354  }
355  res->data=(char *)aa;
356  return FALSE;
357}
358static BOOLEAN jjOP_I_IV(leftv res, leftv u, leftv v)
359{
360  return jjOP_IV_I(res,v,u);
361}
362static BOOLEAN jjOP_IM_I(leftv res, leftv u, leftv v)
363{
364  intvec* aa= (intvec *)u->CopyD(INTVEC_CMD);
365  int bb = (int)(v->Data());
366  int i=min(aa->rows(),aa->cols());
367  switch (iiOp)
368  {
369    case '+': for (;i>0;i--) IMATELEM(*aa,i,i) += bb;
370              break;
371    case '-': for (;i>0;i--) IMATELEM(*aa,i,i) -= bb;
372              break;
373  }
374  res->data=(char *)aa;
375  return FALSE;
376}
377static BOOLEAN jjOP_I_IM(leftv res, leftv u, leftv v)
378{
379  return jjOP_IM_I(res,v,u);
380}
381static BOOLEAN jjDOTDOT(leftv res, leftv u, leftv v)
382{
383  res->data=(char *)new intvec((int)u->Data(),(int)v->Data());
384  return FALSE;
385}
386static void jjEQUAL_REST(leftv res,leftv u,leftv v)
387{
388  if ((res->data) && (u->next!=NULL) && (v->next!=NULL))
389  {
390    int save_iiOp=iiOp;
391    if (iiOp==NOTEQUAL)
392      iiExprArith2(res,u->next,EQUAL_EQUAL,v->next);
393    else
394      iiExprArith2(res,u->next,iiOp,v->next);
395    iiOp=save_iiOp;
396  }
397  if (iiOp==NOTEQUAL) res->data=(char *)(!(int)res->data);
398}
399static BOOLEAN jjCOMPARE_IV(leftv res, leftv u, leftv v)
400{
401  intvec*    a = (intvec * )(u->Data());
402  intvec*    b = (intvec * )(v->Data());
403  int r=a->compare(b);
404  switch  (iiOp)
405  {
406    case '<':
407      res->data  = (char *) (r<0);
408      break;
409    case '>':
410      res->data  = (char *) (r>0);
411      break;
412    case LE:
413      res->data  = (char *) (r<=0);
414      break;
415    case GE:
416      res->data  = (char *) (r>=0);
417      break;
418    case EQUAL_EQUAL:
419    case NOTEQUAL: /* negation handled by jjEQUAL_REST */
420      res->data  = (char *) (r==0);
421      break;
422  }
423  jjEQUAL_REST(res,u,v);
424  return r==-2;
425}
426static BOOLEAN jjCOMPARE_IV_I(leftv res, leftv u, leftv v)
427{
428  intvec* a = (intvec * )(u->Data());
429  int     b = (int )     (v->Data());
430  int r=a->compare(b);
431  switch  (iiOp)
432  {
433    case '<':
434      res->data  = (char *) (r<0);
435      break;
436    case '>':
437      res->data  = (char *) (r>0);
438      break;
439    case LE:
440      res->data  = (char *) (r<=0);
441      break;
442    case GE:
443      res->data  = (char *) (r>=0);
444      break;
445    case EQUAL_EQUAL:
446    case NOTEQUAL: /* negation handled by jjEQUAL_REST */
447      res->data  = (char *) (r==0);
448      break;
449  }
450  jjEQUAL_REST(res,u,v);
451  return FALSE;
452}
453static BOOLEAN jjCOMPARE_P(leftv res, leftv u, leftv v)
454{
455  poly p=(poly)u->Data();
456  poly q=(poly)v->Data();
457  int r=pComp(p,q);
458  if (r==0)
459  {
460    /* compare lead coeffs */
461    number h=nSub(pGetCoeff(p),pGetCoeff(q));
462    r = -1+nIsZero(h)+2*nGreaterZero(h); /* -1: <, 0:==, 1: > */
463    nDelete(&h);
464  }
465  else if (p==NULL)
466  {
467    if (q==NULL)
468    {
469      /* compare 0, 0 */
470      r=0;
471    }
472    else if(pIsConstant(q))
473    {
474      /* compare 0, const */
475      r = 1-2*nGreaterZero(pGetCoeff(q)); /* -1: <, 1: > */
476    }
477  }
478  else if (q==NULL)
479  {
480    if (pIsConstant(p))
481    {
482      /* compare const, 0 */
483      r = -1+2*nGreaterZero(pGetCoeff(p)); /* -1: <, 1: > */
484    }
485  }
486  switch  (iiOp)
487  {
488    case '<':
489      res->data  = (char *) (r < 0);
490      break;
491    case '>':
492      res->data  = (char *) (r > 0);
493      break;
494    case LE:
495      res->data  = (char *) (r <= 0);
496      break;
497    case GE:
498      res->data  = (char *) (r >= 0);
499      break;
500    //case EQUAL_EQUAL:
501    //case NOTEQUAL: /* negation handled by jjEQUAL_REST */
502    //  res->data  = (char *) (r == 0);
503    //  break;
504  }
505  jjEQUAL_REST(res,u,v);
506  return FALSE;
507}
508static BOOLEAN jjCOMPARE_S(leftv res, leftv u, leftv v)
509{
510  char*    a = (char * )(u->Data());
511  char*    b = (char * )(v->Data());
512  int result = strcmp(a,b);
513  switch  (iiOp)
514  {
515    case '<':
516      res->data  = (char *) (result  < 0);
517      break;
518    case '>':
519      res->data  = (char *) (result  > 0);
520      break;
521    case LE:
522      res->data  = (char *) (result  <= 0);
523      break;
524    case GE:
525      res->data  = (char *) (result  >= 0);
526      break;
527    case EQUAL_EQUAL:
528    case NOTEQUAL: /* negation handled by jjEQUAL_REST */
529      res->data  = (char *) (result  == 0);
530      break;
531  }
532  jjEQUAL_REST(res,u,v);
533  return FALSE;
534}
535static BOOLEAN jjOP_REST(leftv res, leftv u, leftv v)
536{
537  if (u->Next()!=NULL)
538  {
539    u=u->next;
540    res->next = (leftv)Alloc(sizeof(sleftv));
541    return iiExprArith2(res->next,u,iiOp,v);
542  }
543  else if (v->Next()!=NULL)
544  {
545    v=v->next;
546    res->next = (leftv)Alloc(sizeof(sleftv));
547    return iiExprArith2(res->next,u,iiOp,v);
548  }
549  return FALSE;
550}
551static BOOLEAN jjPOWER_I(leftv res, leftv u, leftv v)
552{
553  int b=(int)u->Data();
554  int e=(int)v->Data();
555  int rc = 1;
556  BOOLEAN overflow=FALSE;
557  if ((e >= 0)&&(b!=0))
558  {
559    int oldrc;
560    while ((e--)!=0)
561    {
562      oldrc=rc;
563      rc *= b;
564      if (!overflow)
565      {
566        if(rc/b!=oldrc) overflow=TRUE;
567      }
568    }
569    if (overflow)
570      Warn("int overflow(^), result may be wrong");
571    res->data = (char *)rc;
572    if (u!=NULL) return jjOP_REST(res,u,v);
573    return FALSE;
574  }
575  else
576  {
577    WerrorS("exponent must be non-negative");
578    return TRUE;
579  }
580}
581static BOOLEAN jjPOWER_N(leftv res, leftv u, leftv v)
582{
583  int e=(int)v->Data();
584  number n=(number)u->CopyD();
585  if (e<0)
586  {
587    number m=nInvers(n);
588    nDelete(&n);
589    n=m;
590    e=-e;
591  }
592  nPower(n,e,(number*)&res->data);
593  nDelete(&n);
594  if (u!=NULL) return jjOP_REST(res,u,v);
595  return FALSE;
596}
597static BOOLEAN jjPOWER_P(leftv res, leftv u, leftv v)
598{
599  res->data = (char *)pPower((poly)u->CopyD(POLY_CMD),(int)v->Data());
600  if (u!=NULL) return jjOP_REST(res,u,v);
601  return FALSE;
602}
603static BOOLEAN jjPOWER_ID(leftv res, leftv u, leftv v)
604{
605  res->data = (char *)idPower((ideal)(u->Data()),(int)(v->Data()));
606  if (u!=NULL) return jjOP_REST(res,u,v);
607  return FALSE;
608}
609static BOOLEAN jjPLUSMINUS_Gen(leftv res, leftv u, leftv v)
610{
611  u=u->next;
612  v=v->next;
613  if (u==NULL)
614  {
615    if (v==NULL) return FALSE;      /* u==NULL, v==NULL */
616    if (iiOp=='-')                  /* u==NULL, v<>NULL, iiOp=='-'*/
617    {
618      if (res->next==NULL)
619        res->next = (leftv)Alloc(sizeof(sleftv));
620      return iiExprArith1(res->next,v,'-');
621    }
622    loop                            /* u==NULL, v<>NULL, iiOp=='+' */
623    {
624      res->next = (leftv)Alloc0(sizeof(sleftv));
625      res=res->next;
626      res->data = v->CopyD();
627      res->rtyp = v->Typ();
628      v=v->next;
629      if (v==NULL) return FALSE;
630    }
631  }
632  if (v!=NULL)                     /* u<>NULL, v<>NULL */
633  {
634    res->next = (leftv)Alloc(sizeof(sleftv));
635    return iiExprArith2(res->next,u,iiOp,v);
636  }
637  loop                             /* u<>NULL, v==NULL */
638  {
639    res->next = (leftv)Alloc0(sizeof(sleftv));
640    res=res->next;
641    res->data = u->CopyD();
642    res->rtyp = u->Typ();
643    u=u->next;
644    if (u==NULL) return FALSE;
645  }
646}
647static BOOLEAN jjCOLCOL(leftv res, leftv u, leftv v)
648{
649#ifdef HAVE_NAMESPACES
650  BOOLEAN iiReLoadLib(idhdl packhdl);
651  BOOLEAN iiTryLoadLib(leftv v, char *id);
652  idhdl packhdl;
653
654  switch(u->Typ())
655  {
656      case 0:
657        //Print("%s of type 'ANY'. Trying load.\n", v->name);
658        if(iiTryLoadLib(u, u->name))
659        {
660          Werror("'%s' no such package", u->name);
661          return TRUE;
662        }
663        // else: use next case !!! no break !!!
664
665      case PACKAGE_CMD:
666        packhdl = (idhdl)u->data;
667        if(!IDPACKAGE(packhdl)->loaded)
668        {
669          //if(iiReLoadLib(packhdl))
670          //  Werror("unable to reload package '%s'", IDID(packhdl));
671        }
672        if(v->rtyp == IDHDL)
673        {
674          v->name = mstrdup(v->name);
675        }
676        namespaceroot->push( IDPACKAGE(packhdl), IDID(packhdl));
677        syMake(v, v->name, packhdl);
678        memcpy(res, v, sizeof(sleftv));
679        memset(v, 0, sizeof(sleftv));
680        namespaceroot->pop();
681        break;
682
683      case DEF_CMD:
684        break;
685
686      default:
687        WerrorS("<package>::<id> expected");
688        return TRUE;
689  }
690#else /* HAVE_NAMESPACES */
691  WerrorS("package is not supported in this version");
692#endif /* HAVE_NAMESPACES */
693  return FALSE;
694}
695static BOOLEAN jjPLUS_I(leftv res, leftv u, leftv v)
696{
697  unsigned int a=(unsigned int)u->Data();
698  unsigned int b=(unsigned int)v->Data();
699  unsigned int c=a+b;
700  res->data = (char *)c;
701  if (((Sy_bit(31)&a)==(Sy_bit(31)&b))&&((Sy_bit(31)&a)!=(Sy_bit(31)&c)))
702  {
703    Warn("int overflow(+), result may be wrong");
704  }
705  return jjPLUSMINUS_Gen(res,u,v);
706}
707static BOOLEAN jjPLUS_N(leftv res, leftv u, leftv v)
708{
709  res->data = (char *)(nAdd((number)u->Data(), (number)v->Data()));
710  return jjPLUSMINUS_Gen(res,u,v);
711}
712static BOOLEAN jjPLUS_P(leftv res, leftv u, leftv v)
713{
714  res->data = (char *)(pAdd((poly)u->CopyD(POLY_CMD) , (poly)v->CopyD(POLY_CMD)));
715  return jjPLUSMINUS_Gen(res,u,v);
716}
717static BOOLEAN jjPLUS_IV(leftv res, leftv u, leftv v)
718{
719  res->data = (char *)ivAdd((intvec*)(u->Data()), (intvec*)(v->Data()));
720  if (res->data==NULL)
721  {
722     WerrorS("intmat size not compatible");
723     return TRUE;
724  }
725  return jjPLUSMINUS_Gen(res,u,v);
726}
727static BOOLEAN jjPLUS_MA(leftv res, leftv u, leftv v)
728{
729  res->data = (char *)(mpAdd((matrix)u->Data() , (matrix)v->Data()));
730  if (res->data==NULL)
731  {
732     WerrorS("matrix size not compatible");
733     return TRUE;
734  }
735  return jjPLUSMINUS_Gen(res,u,v);
736}
737static BOOLEAN jjPLUS_MA_P(leftv res, leftv u, leftv v)
738{
739  matrix m=(matrix)u->Data();
740  matrix p= mpInitP(m->nrows,m->ncols,(poly)(v->CopyD()));
741  if (iiOp=='+')
742    res->data = (char *)mpAdd(m , p);
743  else
744    res->data = (char *)mpSub(m , p);
745  idDelete((ideal *)&p);
746  return jjPLUSMINUS_Gen(res,u,v);
747}
748static BOOLEAN jjPLUS_P_MA(leftv res, leftv u, leftv v)
749{
750  return jjPLUS_MA_P(res,v,u);
751}
752static BOOLEAN jjPLUS_S(leftv res, leftv u, leftv v)
753{
754  char*    a = (char * )(u->Data());
755  char*    b = (char * )(v->Data());
756  char*    r = (char * )AllocL(strlen(a) + strlen(b) + 1);
757  strcpy(r,a);
758  strcat(r,b);
759  res->data=r;
760  return jjPLUSMINUS_Gen(res,u,v);
761}
762static BOOLEAN jjPLUS_ID(leftv res, leftv u, leftv v)
763{
764  res->data = (char *)idAdd((ideal)u->Data(),(ideal)v->Data());
765  return jjPLUSMINUS_Gen(res,u,v);
766}
767static BOOLEAN jjMINUS_I(leftv res, leftv u, leftv v)
768{
769  unsigned int a=(unsigned int)u->Data();
770  unsigned int b=(unsigned int)v->Data();
771  unsigned int c=a-b;
772  if (((Sy_bit(31)&a)!=(Sy_bit(31)&b))&&((Sy_bit(31)&a)!=(Sy_bit(31)&c)))
773  {
774    Warn("int overflow(-), result may be wrong");
775  }
776  res->data = (char *)c;
777  return jjPLUSMINUS_Gen(res,u,v);
778}
779static BOOLEAN jjMINUS_N(leftv res, leftv u, leftv v)
780{
781  res->data = (char *)(nSub((number)u->Data(), (number)v->Data()));
782  return jjPLUSMINUS_Gen(res,u,v);
783}
784static BOOLEAN jjMINUS_P(leftv res, leftv u, leftv v)
785{
786  res->data = (char *)(pSub((poly)u->CopyD(POLY_CMD) , (poly)v->CopyD(POLY_CMD)));
787  return jjPLUSMINUS_Gen(res,u,v);
788}
789static BOOLEAN jjMINUS_IV(leftv res, leftv u, leftv v)
790{
791  res->data = (char *)ivSub((intvec*)(u->Data()), (intvec*)(v->Data()));
792  if (res->data==NULL)
793  {
794     WerrorS("intmat size not compatible");
795     return TRUE;
796  }
797  return jjPLUSMINUS_Gen(res,u,v);
798}
799static BOOLEAN jjMINUS_MA(leftv res, leftv u, leftv v)
800{
801  res->data = (char *)(mpSub((matrix)u->Data() , (matrix)v->Data()));
802  if (res->data==NULL)
803  {
804     WerrorS("matrix size not compatible");
805     return TRUE;
806  }
807  return jjPLUSMINUS_Gen(res,u,v);
808}
809static BOOLEAN jjTIMES_I(leftv res, leftv u, leftv v)
810{
811  int a=(int)u->Data();
812  int b=(int)v->Data();
813  int c=a * b;
814  if ((b!=0) && (c/b !=a))
815    Warn("int overflow(*), result may be wrong");
816  res->data = (char *)c;
817  if ((u->Next()!=NULL) || (v->Next()!=NULL))
818    return jjOP_REST(res,u,v);
819  return FALSE;
820}
821static BOOLEAN jjTIMES_N(leftv res, leftv u, leftv v)
822{
823  res->data = (char *)(nMult( (number)u->Data(), (number)v->Data()));
824  if ((v->next!=NULL) || (u->next!=NULL))
825    return jjOP_REST(res,u,v);
826  return FALSE;
827}
828static BOOLEAN jjTIMES_P(leftv res, leftv u, leftv v)
829{
830  poly a;
831  poly b;
832
833  if (v->next==NULL)
834  {
835    a=(poly)u->CopyD();
836    if (u->next==NULL)
837    {
838      b=(poly)v->CopyD();
839      res->data = (char *)(pMult( a, b));
840      return FALSE;
841    }
842    // u->next exists: copy v
843    b=pCopy((poly)v->Data());
844    res->data = (char *)(pMult( a, b));
845    return jjOP_REST(res,u,v);
846  }
847  // v->next exists: copy u
848  a=pCopy((poly)u->Data());
849  b=(poly)v->CopyD();
850  res->data = (char *)(pMult( a, b));
851  return jjOP_REST(res,u,v);
852}
853static BOOLEAN jjTIMES_ID(leftv res, leftv u, leftv v)
854{
855  res->data = (char *)idMult((ideal)u->Data(),(ideal)v->Data());
856  if ((v->next!=NULL) || (u->next!=NULL))
857    return jjOP_REST(res,u,v);
858  return FALSE;
859}
860static BOOLEAN jjTIMES_IV(leftv res, leftv u, leftv v)
861{
862  res->data = (char *)ivMult((intvec*)(u->Data()), (intvec*)(v->Data()));
863  if (res->data==NULL)
864  {
865     WerrorS("intmat size not compatible");
866     return TRUE;
867  }
868  if ((v->next!=NULL) || (u->next!=NULL))
869    return jjOP_REST(res,u,v);
870  return FALSE;
871}
872static BOOLEAN jjTIMES_MA_P1(leftv res, leftv u, leftv v)
873{
874  poly p=(poly)v->CopyD(POLY_CMD);
875  int r=pMaxComp(p);/* recompute the rank for the case ideal*vector*/
876  ideal I= (ideal)mpMultP((matrix)u->CopyD(MATRIX_CMD),p);
877  if (r>0) I->rank=r;
878  res->data = (char *)I;
879  return FALSE;
880}
881static BOOLEAN jjTIMES_MA_P2(leftv res, leftv u, leftv v)
882{
883  return jjTIMES_MA_P1(res,v,u);
884}
885static BOOLEAN jjTIMES_MA_N1(leftv res, leftv u, leftv v)
886{
887  number n=(number)v->CopyD(NUMBER_CMD);
888  poly p=pOne();
889  pSetCoeff(p,n);
890  res->data = (char *)mpMultP((matrix)u->CopyD(MATRIX_CMD),p);
891  return FALSE;
892}
893static BOOLEAN jjTIMES_MA_N2(leftv res, leftv u, leftv v)
894{
895  return jjTIMES_MA_N1(res,v,u);
896}
897static BOOLEAN jjTIMES_MA_I1(leftv res, leftv u, leftv v)
898{
899  res->data = (char *)mpMultI((matrix)u->CopyD(MATRIX_CMD),(int)v->Data());
900  return FALSE;
901}
902static BOOLEAN jjTIMES_MA_I2(leftv res, leftv u, leftv v)
903{
904  return jjTIMES_MA_I1(res,v,u);
905}
906static BOOLEAN jjTIMES_MA(leftv res, leftv u, leftv v)
907{
908  res->data = (char *)mpMult((matrix)u->Data(),(matrix)v->Data());
909  if (res->data==NULL)
910  {
911     WerrorS("matrix size not compatible");
912     return TRUE;
913  }
914  if ((v->next!=NULL) || (u->next!=NULL))
915    return jjOP_REST(res,u,v);
916  return FALSE;
917}
918static BOOLEAN jjGE_I(leftv res, leftv u, leftv v)
919{
920  res->data = (char *)((int)u->Data() >= (int)v->Data());
921  return FALSE;
922}
923static BOOLEAN jjGE_N(leftv res, leftv u, leftv v)
924{
925  number h=nSub((number)u->Data(),(number)v->Data());
926  res->data = (char *) (nGreaterZero(h)||(nIsZero(h)));
927  nDelete(&h);
928  return FALSE;
929}
930static BOOLEAN jjGT_I(leftv res, leftv u, leftv v)
931{
932  res->data = (char *)((int)u->Data()>(int)v->Data());
933  return FALSE;
934}
935static BOOLEAN jjGT_N(leftv res, leftv u, leftv v)
936{
937  number h=nSub((number)u->Data(),(number)v->Data());
938  res->data = (char *) (nGreaterZero(h)&&(!nIsZero(h)));
939  nDelete(&h);
940  return FALSE;
941}
942static BOOLEAN jjLE_I(leftv res, leftv u, leftv v)
943{
944  res->data = (char *)((int)u->Data() <= (int)v->Data());
945  return FALSE;
946}
947static BOOLEAN jjLE_N(leftv res, leftv u, leftv v)
948{
949  number h=nSub((number)v->Data(),(number)u->Data());
950  res->data = (char *) (nGreaterZero(h)||nIsZero(h));
951  nDelete(&h);
952  return FALSE;
953}
954static BOOLEAN jjLT_I(leftv res, leftv u, leftv v)
955{
956  res->data = (char *)((int)u->Data() < (int)v->Data());
957  return FALSE;
958}
959static BOOLEAN jjLT_N(leftv res, leftv u, leftv v)
960{
961  number h=nSub((number)v->Data(),(number)u->Data());
962  res->data = (char *) (nGreaterZero(h)&&(!nIsZero(h)));
963  nDelete(&h);
964  return FALSE;
965}
966static BOOLEAN jjDIV_I(leftv res, leftv u, leftv v)
967{
968  int a= (int) u->Data();
969  int b= (int) v->Data();
970  if (b==0)
971  {
972    WerrorS("div. by 0");
973    return TRUE;
974  }
975  //int c=a%ABS(b);
976  //if(c<0) c+=ABS(b);
977  //res->data = (char *)((a-c) / b);
978  res->data = (char *)(a / b);
979  return FALSE;
980}
981static BOOLEAN jjDIV_N(leftv res, leftv u, leftv v)
982{
983  number q=(number)v->Data();
984  if (nIsZero(q))
985  {
986    WerrorS("div. by 0");
987    return TRUE;
988  }
989  q = nDiv((number)u->Data(),q);
990  nNormalize(q);
991  res->data = (char *)q;
992  return FALSE;
993}
994static BOOLEAN jjDIV_P(leftv res, leftv u, leftv v)
995{
996  poly q=(poly)v->Data();
997  if (q==NULL)
998  {
999    WerrorS("div. by 0");
1000    return TRUE;
1001  }
1002  if (pNext(q)!=NULL)
1003  {
1004#ifdef HAVE_FACTORY
1005    res->data=(void*)(singclap_pdivide((poly)(u->Data()),
1006                                       (poly)(v->Data())));
1007#else
1008    WerrorS("division only by a monomial");
1009    return TRUE;
1010#endif
1011  }
1012  else
1013  {
1014    res->data = (char *)pDivideM((poly)u->CopyD(POLY_CMD),pHead(q));
1015  }
1016  return FALSE;
1017}
1018static BOOLEAN jjMOD_I(leftv res, leftv u, leftv v)
1019{
1020  int a=(int)u->Data();
1021  int b=ABS((int)v->Data());
1022  if (errorreported) return TRUE;
1023  if (b==0)
1024  {
1025    WerrorS("div. by 0");
1026    return TRUE;
1027  }
1028  int c=a%b;
1029  if(c<0) c+=b;
1030  res->data = (char *)c;
1031  return FALSE;
1032}
1033static BOOLEAN jjEQUAL_I(leftv res, leftv u, leftv v)
1034{
1035  res->data = (char *)((int)u->Data() == (int)v->Data());
1036  jjEQUAL_REST(res,u,v);
1037  return FALSE;
1038}
1039static BOOLEAN jjEQUAL_Ma(leftv res, leftv u, leftv v)
1040{
1041  res->data = (char *)mpEqual((matrix)u->Data(),(matrix)v->Data());
1042  jjEQUAL_REST(res,u,v);
1043  return FALSE;
1044}
1045static BOOLEAN jjEQUAL_N(leftv res, leftv u, leftv v)
1046{
1047  res->data = (char *)nEqual((number)u->Data(),(number)v->Data());
1048  jjEQUAL_REST(res,u,v);
1049  return FALSE;
1050}
1051static BOOLEAN jjEQUAL_P(leftv res, leftv u, leftv v)
1052{
1053  poly p=(poly)u->Data();
1054  poly q=(poly)v->Data();
1055  if (p==NULL)
1056  {
1057    res->data=(char *)(q==NULL);
1058  }
1059  else
1060  {
1061    if (q==NULL) res->data=(char *)FALSE;
1062    else
1063    {
1064      int r=pComp(p,q);
1065      if (r==0)
1066      {
1067        p=pSub(pCopy(p),pCopy(q));
1068        res->data = (char *) (p==NULL);
1069        pDelete(&p);
1070      }
1071      else
1072        res->data = (char *) FALSE;
1073    }
1074  }
1075  jjEQUAL_REST(res,u,v);
1076  return FALSE;
1077}
1078static BOOLEAN jjAND_I(leftv res, leftv u, leftv v)
1079{
1080  res->data = (char *)((int)u->Data() && (int)v->Data());
1081  return FALSE;
1082}
1083static BOOLEAN jjOR_I(leftv res, leftv u, leftv v)
1084{
1085  res->data = (char *)((int)u->Data() || (int)v->Data());
1086  return FALSE;
1087}
1088static BOOLEAN jjINDEX_I(leftv res, leftv u, leftv v)
1089{
1090  res->rtyp=u->rtyp; u->rtyp=0;
1091  res->data=u->data; u->data=NULL;
1092  res->name=u->name; u->name=NULL;
1093#ifdef HAVE_NAMESPACES
1094  res->packhdl=u->packhdl; u->packhdl=NULL;
1095  res->req_packhdl=u->req_packhdl; u->req_packhdl=NULL;
1096#endif /* HAVE_NAMESPACES */
1097  res->e=u->e;       u->e=NULL;
1098  if (res->e==NULL) res->e=jjMakeSub(v);
1099  else
1100  {
1101    Subexpr sh=res->e;
1102    while (sh->next != NULL) sh=sh->next;
1103    sh->next=jjMakeSub(v);
1104  }
1105  return FALSE;
1106}
1107static BOOLEAN jjINDEX_IV(leftv res, leftv u, leftv v)
1108{
1109  if ((u->rtyp!=IDHDL)||(u->e!=NULL))
1110  {
1111    WerrorS("indexed object must have a name");
1112    return TRUE;
1113  }
1114  intvec * iv=(intvec *)v->Data();
1115  leftv p=NULL;
1116  int i;
1117  sleftv t;
1118
1119  memset(&t,0,sizeof(t));
1120  t.rtyp=INT_CMD;
1121  for (i=0;i<iv->length(); i++)
1122  {
1123    t.data=(char *)(*iv)[i];
1124    if (p==NULL)
1125    {
1126      p=res;
1127    }
1128    else
1129    {
1130      p->next=(leftv)Alloc0(sizeof(sleftv));
1131      p=p->next;
1132    }
1133    p->rtyp=IDHDL;
1134    p->data=u->data;
1135    p->name=u->name;
1136#ifdef HAVE_NAMESPACES
1137    p->packhdl=u->packhdl;
1138#endif /* HAVE_NAMESPACES */
1139    p->flag|=u->flag;
1140    p->e=jjMakeSub(&t);
1141  }
1142  u->rtyp=0;
1143  u->data=NULL;
1144  u->name=NULL;
1145  return FALSE;
1146}
1147static BOOLEAN jjINDEX_P(leftv res, leftv u, leftv v)
1148{
1149  poly p=(poly)u->Data();
1150  int i=(int)v->Data();
1151  int j=0;
1152  while (p!=NULL)
1153  {
1154    j++;
1155    if (j==i)
1156    {
1157      res->data=(char *)pHead(p);
1158      return FALSE;
1159    }
1160    pIter(p);
1161  }
1162  return FALSE;
1163}
1164static BOOLEAN jjINDEX_P_IV(leftv res, leftv u, leftv v)
1165{
1166  poly p=(poly)u->Data();
1167  poly r=NULL;
1168  intvec *iv=(intvec *)v->CopyD();
1169  int i;
1170  int sum=0;
1171  for(i=iv->length()-1;i>=0;i--)
1172    sum+=(*iv)[i];
1173  int j=0;
1174  while ((p!=NULL) && (sum>0))
1175  {
1176    j++;
1177    for(i=iv->length()-1;i>=0;i--)
1178    {
1179      if (j==(*iv)[i])
1180      {
1181        r=pAdd(r,pHead(p));
1182        sum-=j;
1183        (*iv)[i]=0;
1184        break;
1185      }
1186    }
1187    pIter(p);
1188  }
1189  delete iv;
1190  res->data=(char *)r;
1191  return FALSE;
1192}
1193static BOOLEAN jjINDEX_V(leftv res, leftv u, leftv v)
1194{
1195  poly p=(poly)u->CopyD();
1196  poly r=p; // pointer to the beginning of component i
1197  poly o=NULL;
1198  int i=(int)v->Data();
1199  while (p!=NULL)
1200  {
1201    if (pGetComp(p)!=i)
1202    {
1203      if (r==p) r=pNext(p);
1204      if (o!=NULL)
1205      {
1206        pDelete1(&pNext(o));
1207        p=pNext(o);
1208      }
1209      else
1210        pDelete1(&p);
1211    }
1212    else
1213    {
1214      pSetComp(p, 0);
1215      o=p;
1216      p=pNext(o);
1217    }
1218  }
1219  res->data=(char *)r;
1220  return FALSE;
1221}
1222static BOOLEAN jjINDEX_V_IV(leftv res, leftv u, leftv v)
1223{
1224  poly p=(poly)u->CopyD();
1225  if (p!=NULL)
1226  {
1227    poly r=pOne();
1228    poly hp=r;
1229    intvec *iv=(intvec *)v->Data();
1230    int i;
1231    loop
1232    {
1233      for(i=0;i<iv->length();i++)
1234      {
1235        if (pGetComp(p)==(*iv)[i])
1236        {
1237          poly h;
1238          pSplit(p,&h);
1239          pNext(hp)=p;
1240          p=h;
1241          pIter(hp);
1242          break;
1243        }
1244      }
1245      if (p==NULL) break;
1246      if (i==iv->length())
1247      {
1248        pDelete1(&p);
1249        if (p==NULL) break;
1250      }
1251    }
1252    pDelete1(&r);
1253    res->data=(char *)r;
1254  }
1255  return FALSE;
1256}
1257static BOOLEAN jjKLAMMER(leftv res, leftv u, leftv v)
1258{
1259  if(u->name==NULL) return TRUE;
1260  char * nn = (char *)AllocL(strlen(u->name) + 13);
1261  sprintf(nn,"%s(%d)",u->name,(int)v->Data());
1262  FreeL((ADDRESS)u->name);
1263  u->name=NULL;
1264  char *n=mstrdup(nn);
1265  FreeL((ADDRESS)nn);
1266#ifdef HAVE_NAMESPACES
1267  if(u->req_packhdl != NULL)
1268  {
1269    namespaceroot->push( IDPACKAGE(u->req_packhdl), IDID(u->req_packhdl));
1270    syMake(res,n);
1271    namespaceroot->pop();
1272  }
1273  else
1274#endif /* HAVE_NAMESPACES */
1275    syMake(res,n);
1276  return FALSE;
1277}
1278static BOOLEAN jjKLAMMER_IV(leftv res, leftv u, leftv v)
1279{
1280  intvec * iv=(intvec *)v->Data();
1281  leftv p=NULL;
1282  int i;
1283  char *n;
1284#ifdef HAVE_NAMESPACES
1285  BOOLEAN needpop=FALSE;
1286
1287  if(u->req_packhdl != NULL)
1288  {
1289    namespaceroot->push( IDPACKAGE(u->req_packhdl), IDID(u->req_packhdl));
1290    needpop = TRUE;
1291  }
1292#endif /* HAVE_NAMESPACES */
1293
1294  for (i=0;i<iv->length(); i++)
1295  {
1296    if (p==NULL)
1297    {
1298      p=res;
1299    }
1300    else
1301    {
1302      p->next=(leftv)Alloc0(sizeof(sleftv));
1303      p=p->next;
1304    }
1305    n = (char *)AllocL(strlen(u->name) + 6);
1306    sprintf(n,"%s(%d)",u->name,(*iv)[i]);
1307    syMake(p,n);
1308  }
1309  FreeL((ADDRESS)u->name);
1310  u->name=NULL;
1311#ifdef HAVE_NAMESPACES
1312  if(needpop) namespaceroot->pop();
1313#endif /* HAVE_NAMESPACES */
1314  return FALSE;
1315}
1316static BOOLEAN jjPROC(leftv res, leftv u, leftv v)
1317{
1318#ifdef HAVE_NAMESPACES
1319  leftv sl = iiMake_proc((idhdl)u->data,u,v);
1320#else /* HAVE_NAMESPACES */
1321  leftv sl = iiMake_proc((idhdl)u->data,v);
1322#endif /* HAVE_NAMESPACES */
1323  if (sl==NULL)
1324  {
1325    return TRUE;
1326  }
1327  else
1328  {
1329    memcpy(res,sl,sizeof(sleftv));
1330  }
1331  return FALSE;
1332}
1333static BOOLEAN jjMAP(leftv res, leftv u, leftv v)
1334{
1335  //Print("try to map %s with %s\n",$3.Name(),$1.Name());
1336  leftv sl=NULL;
1337  if ((v->e==NULL)&&(v->name!=NULL))
1338  {
1339    map m=(map)u->Data();
1340    sl=iiMap(m,v->name);
1341  }
1342  else
1343  {
1344    Werror("%s(<name>) expected",u->Name());
1345  }
1346  if (sl==NULL) return TRUE;
1347  memcpy(res,sl,sizeof(sleftv));
1348  Free((ADDRESS)sl,sizeof(*sl));
1349  return FALSE;
1350}
1351static BOOLEAN jjCALL2MANY(leftv res, leftv u, leftv v)
1352{
1353  u->next=(leftv)Alloc(sizeof(sleftv));
1354  memcpy(u->next,v,sizeof(sleftv));
1355  BOOLEAN r=iiExprArithM(res,u,iiOp);
1356  v->rtyp=0; v->data=NULL; // iiExprArithM did the CleanUp
1357  return r;
1358}
1359static BOOLEAN jjCOEF(leftv res, leftv u, leftv v)
1360{
1361  res->data=(char *)mpCoeffProc((poly)u->Data(),(poly)v->Data());
1362  return FALSE;
1363}
1364static BOOLEAN jjCOEFFS_Id(leftv res, leftv u, leftv v)
1365{
1366  int i=pVar((poly)v->Data());
1367  if (i==0)
1368  {
1369    WerrorS("ringvar expected");
1370    return TRUE;
1371  }
1372  res->data=(char *)mpCoeffs((ideal)u->CopyD(),i);
1373  return FALSE;
1374}
1375static BOOLEAN jjCONTRACT(leftv res, leftv u, leftv v)
1376{
1377  res->data=(char *)idDiffOp((ideal)u->Data(),(ideal)v->Data(),FALSE);
1378  return FALSE;
1379}
1380static BOOLEAN jjDEG_IV(leftv res, leftv u, leftv v)
1381{
1382  short *iv=iv2array((intvec *)v->Data());
1383  res->data = (char *)pDegW((poly)u->Data(),iv);
1384  Free((ADDRESS)iv,(pVariables+1)*sizeof(short));
1385  return FALSE;
1386}
1387static BOOLEAN jjDIFF_P(leftv res, leftv u, leftv v)
1388{
1389  int i=pVar((poly)v->Data());
1390  if (i==0)
1391  {
1392    WerrorS("ringvar expected");
1393    return TRUE;
1394  }
1395  res->data=(char *)pDiff((poly)(u->Data()),i);
1396  return FALSE;
1397}
1398static BOOLEAN jjDIFF_ID(leftv res, leftv u, leftv v)
1399{
1400  int i=pVar((poly)v->Data());
1401  if (i==0)
1402  {
1403    WerrorS("ringvar expected");
1404    return TRUE;
1405  }
1406  res->data=(char *)idDiff((matrix)(u->Data()),i);
1407  return FALSE;
1408}
1409static BOOLEAN jjDIFF_ID_ID(leftv res, leftv u, leftv v)
1410{
1411  res->data=(char *)idDiffOp((ideal)u->Data(),(ideal)v->Data());
1412  return FALSE;
1413}
1414
1415static BOOLEAN jjELIMIN(leftv res, leftv u, leftv v)
1416{
1417  res->data=(char *)idElimination((ideal)u->Data(),(poly)v->Data());
1418  setFlag(res,FLAG_STD);
1419  return FALSE;
1420}
1421static BOOLEAN jjDIM2(leftv res, leftv v, leftv w)
1422{
1423  assumeStdFlag(v);
1424  if(currQuotient==NULL)
1425    res->data = (char *)scDimInt((ideal)(v->Data()),(ideal)w->Data());
1426  else
1427  {
1428    ideal q=idSimpleAdd(currQuotient,(ideal)w->Data());
1429    res->data = (char *)scDimInt((ideal)(v->Data()),q);
1430    idDelete(&q);
1431  }
1432  return FALSE;
1433}
1434static BOOLEAN jjEXTGCD_I(leftv res, leftv u, leftv v)
1435{
1436  int p0=ABS((int)u->Data()),p1=ABS((int)v->Data());
1437  int f0 = 1, f1 = 0, g0 = 0, g1 = 1, q, r;
1438
1439  while ( p1!=0 )
1440  {
1441    q=p0 / p1;
1442    r=p0 % p1;
1443    p0 = p1; p1 = r;
1444    r = g0 - g1 * q;
1445    g0 = g1; g1 = r;
1446    r = f0 - f1 * q;
1447    f0 = f1; f1 = r;
1448  }
1449  int a = f0;
1450  int b = g0;
1451  if ( (int)u->Data() < 0 ) a=-a;
1452  if ( (int)v->Data() < 0 ) b=-b;
1453  lists L=(lists)Alloc(sizeof(slists));
1454  L->Init(3);
1455  L->m[0].rtyp=INT_CMD;
1456  L->m[0].data=(void *)p0;
1457  L->m[1].rtyp=INT_CMD;
1458  L->m[1].data=(void *)a;
1459  L->m[2].rtyp=INT_CMD;
1460  L->m[2].data=(void *)b;
1461  res->rtyp=LIST_CMD;
1462  res->data=(char *)L;
1463  return FALSE;
1464}
1465#ifdef HAVE_FACTORY
1466static BOOLEAN jjFACSTD2(leftv res, leftv v, leftv w)
1467{
1468  res->data=(void *)kStdfac((ideal)v->Data(),NULL,testHomog,NULL,
1469           (ideal)w->Data());
1470  setFlag(res,FLAG_STD);
1471  return FALSE;
1472}
1473#endif
1474static BOOLEAN jjFETCH(leftv res, leftv u, leftv v)
1475{
1476  ring r=(ring)u->Data();
1477  idhdl w;
1478
1479  if ((iiOp!=IMAP_CMD)
1480  && ((rInternalChar(currRing) != rInternalChar(r))
1481    || ((currRing->N != r->N)&& (iiOp==FETCH_CMD)))
1482  )
1483    goto err_fetch;
1484  if ((w=r->idroot->get(v->Name(),myynest))!=NULL)
1485  {
1486    int *perm=NULL;
1487    int *par_perm=NULL;
1488    int par_perm_size=0;
1489    BOOLEAN bo;
1490    if (iiOp==IMAP_CMD)
1491    {
1492      if (!nSetMap(rInternalChar(r),r->parameter,rPar(r),r->minpoly))
1493      {
1494        if (iiOp != IMAP_CMD) goto err_fetch;
1495        // Allow imap to be make an exception only for:
1496        if ( (rField_is_Q_a(r) &&  // Q(a..) -> Q(a..) || Q || Zp || Zp(a)
1497              (rField_is_Q() || rField_is_Q_a() ||
1498               (rField_is_Zp() || rField_is_Zp_a())))
1499             ||
1500             (rField_is_Zp_a(r) &&  // Zp(a..) -> Zp(a..) || Zp
1501              (rField_is_Zp(currRing, rInternalChar(r)) ||
1502               rField_is_Zp_a(currRing, rInternalChar(r)))) )
1503        {
1504          par_perm_size=rPar(r);
1505          BITSET save_test=test;
1506          naSetChar(rInternalChar(r),TRUE,r->parameter,rPar(r));
1507          nSetChar(rInternalChar(currRing),TRUE,currRing->parameter,
1508                   rPar(currRing));
1509          test=save_test;
1510        }
1511        else
1512        {
1513          goto err_fetch;
1514        }
1515      }
1516      perm=(int *)Alloc0((r->N+1)*sizeof(int));
1517      if (par_perm_size!=0)
1518        par_perm=(int *)Alloc0(par_perm_size*sizeof(int));
1519      maFindPerm(r->names,       r->N,       r->parameter,        r->P,
1520                 currRing->names,currRing->N,currRing->parameter, currRing->P,
1521                 perm,par_perm, currRing->ch);
1522    }
1523    sleftv tmpW;
1524    memset(&tmpW,0,sizeof(sleftv));
1525    tmpW.rtyp=IDTYP(w);
1526    tmpW.data=IDDATA(w);
1527    if ((bo=maApplyFetch(iiOp,NULL,res,&tmpW, r, perm,par_perm,par_perm_size)))
1528    {
1529      Werror("cannot map %s of type %s(%d)",v->name, Tok2Cmdname(w->typ),w->typ);
1530    }
1531    if (perm!=NULL)
1532      Free((ADDRESS)perm,(r->N+1)*sizeof(int));
1533    if (par_perm!=NULL)
1534      Free((ADDRESS)par_perm,par_perm_size*sizeof(int));
1535    return bo;
1536  }
1537  else
1538  {
1539    Werror("%s undefined in %s",v->Fullname(),u->Fullname());
1540  }
1541  return TRUE;
1542err_fetch:
1543  Werror("no identity map from %s",u->Fullname());
1544  return TRUE;
1545}
1546static BOOLEAN jjFIND2(leftv res, leftv u, leftv v)
1547{
1548  /*4
1549  * look for the substring what in the string where
1550  * return the position of the first char of what in where
1551  * or 0
1552  */
1553  char *where=(char *)u->Data();
1554  char *what=(char *)v->Data();
1555  char *found = strstr(where,what);
1556  if (found != NULL)
1557  {
1558    res->data=(char *)((found-where)+1);
1559  }
1560  /*else res->data=NULL;*/
1561  return FALSE;
1562}
1563static BOOLEAN jjGCD_I(leftv res, leftv u, leftv v)
1564{
1565  int p0=ABS((int)u->Data()),p1=ABS((int)v->Data());
1566  int r;
1567
1568  while ( p1!=0 )
1569  {
1570    r=p0 % p1;
1571    p0 = p1; p1 = r;
1572  }
1573  res->rtyp=INT_CMD;
1574  res->data=(char *)p0;
1575  return FALSE;
1576}
1577static BOOLEAN jjHILBERT2(leftv res, leftv u, leftv v)
1578{
1579  assumeStdFlag(u);
1580  intvec *module_w=(intvec*)atGet(u,"isHomog");
1581  intvec *iv=hFirstSeries((ideal)u->Data(),module_w,currQuotient);
1582  switch((int)v->Data())
1583  {
1584    case 1:
1585      res->data=(void *)iv;
1586      return FALSE;
1587    case 2:
1588      res->data=(void *)hSecondSeries(iv);
1589      delete iv;
1590      return FALSE;
1591  }
1592  WerrorS(feNotImplemented);
1593  delete iv;
1594  return TRUE;
1595}
1596static BOOLEAN jjHOMOG_P(leftv res, leftv u, leftv v)
1597{
1598  int i=pVar((poly)v->Data());
1599  if (i==0)
1600  {
1601    WerrorS("ringvar expected");
1602    return TRUE;
1603  }
1604  res->data = (char *)pHomogen((poly)u->Data(),i);
1605  return FALSE;
1606}
1607static BOOLEAN jjHOMOG_ID(leftv res, leftv u, leftv v)
1608{
1609  int i=pVar((poly)v->Data());
1610  if (i==0)
1611  {
1612    WerrorS("ringvar expected");
1613    return TRUE;
1614  }
1615  res->data = (char *)idHomogen((ideal)u->Data(),i);
1616  return FALSE;
1617}
1618static BOOLEAN jjINDEPSET2(leftv res, leftv u, leftv v)
1619{
1620  assumeStdFlag(u);
1621  res->data=(void *)scIndIndset((ideal)(u->Data()),(int)(v->Data()),
1622                    currQuotient);
1623  return FALSE;
1624}
1625static BOOLEAN jjINTERSECT(leftv res, leftv u, leftv v)
1626{
1627  res->data=(char *)idSect((ideal)u->Data(),(ideal)v->Data());
1628  setFlag(res,FLAG_STD);
1629  return FALSE;
1630}
1631static BOOLEAN jjJET_P(leftv res, leftv u, leftv v)
1632{
1633  res->data = (char *)pJet((poly)u->Data(),(int)v->Data());
1634  return FALSE;
1635}
1636static BOOLEAN jjJET_ID(leftv res, leftv u, leftv v)
1637{
1638  res->data = (char *)idJet((ideal)u->Data(),(int)v->Data());
1639  return FALSE;
1640}
1641static BOOLEAN jjKBASE2(leftv res, leftv u, leftv v)
1642{
1643  assumeStdFlag(u);
1644  res->data = (char *)scKBase((int)v->Data(),
1645                              (ideal)(u->Data()),currQuotient);
1646  return FALSE;
1647}
1648static BOOLEAN jjKoszul(leftv res, leftv u, leftv v)
1649{
1650  return mpKoszul(res, u,v);
1651}
1652static BOOLEAN jjKoszul_Id(leftv res, leftv u, leftv v)
1653{
1654  sleftv h;
1655  memset(&h,0,sizeof(sleftv));
1656  h.rtyp=INT_CMD;
1657  h.data=(void *)IDELEMS((ideal)v->Data());
1658  return mpKoszul(res, &h, u, v);
1659}
1660static BOOLEAN jjLIFT(leftv res, leftv u, leftv v)
1661{
1662  ideal m;
1663  int ul= IDELEMS((ideal)u->Data());
1664  int vl= IDELEMS((ideal)v->Data());
1665  if (hasFlag(u,FLAG_STD))
1666  {
1667    m = idLift((ideal)u->Data(),(ideal)v->Data());
1668    res->data = (char *)idModule2formatedMatrix(m,ul,vl);
1669  }
1670  else
1671  {
1672    // matrix ma=mpNew(1,1);
1673    // ideal mo=idLiftStd((ideal)u->CopyD(),currQuotient,&ma,testHomog);
1674    // m = idLift(mo,(ideal)v->Data());
1675    // matrix r=idModule2formatedMatrix(m, IDELEMS(mo),vl);
1676    // idDelete(&mo);
1677    // // idDelete(&m); already done by idModule2formatedMatrix
1678    // res->data=(char *)mpMult(ma,r);
1679    // idDelete((ideal *)&ma);
1680    // idDelete((ideal *)&r);
1681    m = idLiftNonStB((ideal)u->Data(),(ideal)v->Data());
1682    res->data = (char *)idModule2formatedMatrix(m,ul,vl);
1683  }
1684  return FALSE;
1685}
1686static BOOLEAN jjLIFTSTD(leftv res, leftv u, leftv v)
1687{
1688  if ((v->rtyp!=IDHDL)||(v->e!=NULL)) return TRUE;
1689  idhdl h=(idhdl)v->data;
1690  res->data = (char *)idLiftStd((ideal)u->CopyD(), currQuotient,
1691              &(h->data.umatrix),testHomog);
1692  setFlag(res,FLAG_STD);
1693  return FALSE;
1694}
1695static BOOLEAN jjMINOR(leftv res, leftv u, leftv v)
1696{
1697  res->data = (char *)idMinors((matrix)u->Data(),(int)v->Data());
1698  return FALSE;
1699}
1700static BOOLEAN jjMODULO(leftv res, leftv u, leftv v)
1701{
1702  res->data = (char *)idModulo((ideal)u->Data(),(ideal)v->Data());
1703  return FALSE;
1704}
1705static BOOLEAN jjMONITOR2(leftv res, leftv u,leftv v)
1706{
1707  char *opt=(char *)v->Data();
1708  int mode=0;
1709  while(*opt!='\0')
1710  {
1711    if (*opt=='i') mode |= PROT_I;
1712    else if (*opt=='o') mode |= PROT_O;
1713    opt++;
1714  }
1715  monitor((char *)(u->Data()),mode);
1716  return FALSE;
1717}
1718static BOOLEAN jjPARSTR2(leftv res, leftv u, leftv v)
1719{
1720  idhdl h=(idhdl)u->data;
1721  int i=(int)v->Data();
1722  int p=0;
1723  if ((0<i)
1724  && (IDRING(h)->parameter!=NULL)
1725  && (i<=(p=rPar(IDRING(h)))))
1726    res->data=mstrdup(IDRING(h)->parameter[i-1]);
1727  else
1728  {
1729    Werror("par number %d out of range 1..%d",i,p);
1730    return TRUE;
1731  }
1732  return FALSE;
1733}
1734static BOOLEAN jjQUOT(leftv res, leftv u, leftv v)
1735{
1736  res->data = (char *)idQuot((ideal)u->Data(),(ideal)v->Data(),
1737    hasFlag(u,FLAG_STD),u->Typ()==v->Typ());
1738  idDelMultiples((ideal)(res->data));
1739  return FALSE;
1740}
1741static BOOLEAN jjRANDOM(leftv res, leftv u, leftv v)
1742{
1743  int i=(int)u->Data();
1744  int j=(int)v->Data();
1745#ifdef buildin_rand
1746  res->data =(char *)((i > j) ? i : (siRand() % (j-i+1)) + i);
1747#else
1748  res->data =(char *)((i > j) ? i : (rand() % (j-i+1)) + i);
1749#endif
1750  return FALSE;
1751}
1752static BOOLEAN jjREAD2(leftv res, leftv u, leftv v)
1753{
1754  si_link l=(si_link)u->Data();
1755  leftv r=slRead(l,v);
1756  if (r==NULL)
1757  {
1758    const char *s;
1759    if ((l!=NULL)&&(l->name!=NULL)) s=l->name;
1760    else                            s=sNoName;
1761    Werror("cannot read from `%s`",s);
1762    return TRUE;
1763  }
1764  memcpy(res,r,sizeof(sleftv));
1765  Free((ADDRESS)r,sizeof(sleftv));
1766  return FALSE;
1767}
1768static BOOLEAN jjREDUCE_P(leftv res, leftv u, leftv v)
1769{
1770  assumeStdFlag(v);
1771  res->data = (char *)kNF((ideal)v->Data(),currQuotient,(poly)u->Data());
1772  return FALSE;
1773}
1774static BOOLEAN jjREDUCE_ID(leftv res, leftv u, leftv v)
1775{
1776  assumeStdFlag(v);
1777  res->data = (char *)kNF((ideal)v->Data(),currQuotient,(ideal)u->Data());
1778  return FALSE;
1779}
1780static BOOLEAN jjRES(leftv res, leftv u, leftv v)
1781{
1782  int maxl=(int)v->Data();
1783  if (maxl<0)
1784  {
1785    WerrorS("length for res must not be negative");
1786    return TRUE;
1787  }
1788  int l=0;
1789  //resolvente r;
1790  syStrategy r;
1791  intvec **weights=NULL;
1792  int wmaxl=maxl;
1793  ideal u_id=(ideal)u->Data();
1794  maxl--;
1795  if ((maxl==-1) && (iiOp!=MRES_CMD))
1796    maxl = pVariables-1;
1797  if ((iiOp == RES_CMD) || (iiOp == MRES_CMD))
1798  {
1799    //if (BTEST1(28))
1800    //{
1801    //  r=syMinRes(u_id,maxl,&l, iiOp==MRES_CMD);
1802    //}
1803    //else
1804    {
1805      intvec * iv=(intvec*)atGet(u,"isHomog");
1806      //if (iv!=NULL)
1807      //{
1808      //  weights = (intvec**)Alloc0(sizeof(intvec*));
1809      //  weights[0] = ivCopy(iv);
1810      //  l=1;
1811      //}
1812      //r=syResolvente(u_id,maxl,&l, &weights, iiOp==MRES_CMD);
1813      r=syResolution(u_id,maxl, iv, iiOp==MRES_CMD);
1814    }
1815  }
1816  else if (iiOp==SRES_CMD)
1817  //  r=sySchreyerResolvente(u_id,maxl+1,&l);
1818    r=sySchreyer(u_id,maxl+1);
1819  else /* LRES */
1820  {
1821    int dummy;
1822    if((currQuotient!=NULL)||
1823    (!idHomIdeal (u_id,NULL)))
1824    {
1825       WerrorS
1826       ("`lres` can not be called in q ring or with inhomogeneous input");
1827       return TRUE;
1828    }
1829    r=syLaScala3(u_id,&dummy);
1830  }
1831  if (r==NULL) return TRUE;
1832  //res->data=(void *)liMakeResolv(r,l,wmaxl,u->Typ(),weights);
1833  r->list_length=wmaxl;
1834  res->data=(void *)r;
1835  return FALSE;
1836}
1837static BOOLEAN jjRSUM(leftv res, leftv u, leftv v)
1838{
1839  ring r;
1840  int i=rSum((ring)u->Data(),(ring)v->Data(),r);
1841  res->data = (char *)r;
1842  return (i==-1);
1843}
1844#define SIMPL_LMDIV 32
1845#define SIMPL_LMEQ  16
1846#define SIMPL_MULT 8
1847#define SIMPL_EQU  4
1848#define SIMPL_NULL 2
1849#define SIMPL_NORM 1
1850static BOOLEAN jjSIMPL_ID(leftv res, leftv u, leftv v)
1851{
1852  int sw = (int)v->Data();
1853  ideal id = (ideal)u->CopyD();
1854  if (sw & SIMPL_LMDIV)
1855  {
1856    idDelDiv(id);
1857  }
1858  if (sw & SIMPL_LMEQ)
1859  {
1860    idDelLmEquals(id);
1861  }
1862  if (sw & SIMPL_NULL)
1863  {
1864    idSkipZeroes(id);
1865  }
1866  if (sw & SIMPL_NORM)
1867  {
1868    idNorm(id);
1869  }
1870  if (sw & SIMPL_MULT)
1871  {
1872    idDelMultiples(id);
1873  }
1874  else if(sw & SIMPL_EQU)
1875  {
1876    idDelEquals(id);
1877  }
1878  res->data = (char * )id;
1879  return FALSE;
1880}
1881static BOOLEAN jjSTATUS2(leftv res, leftv u, leftv v)
1882{
1883  res->data = mstrdup(slStatus((si_link) u->Data(), (char *) v->Data()));
1884  return FALSE;
1885}
1886static BOOLEAN jjSIMPL_P(leftv res, leftv u, leftv v)
1887{
1888  int sw = (int)v->Data();
1889  poly p = (poly)u->CopyD();
1890  if (sw & SIMPL_NORM)
1891  {
1892    pNorm(p);
1893  }
1894  res->data = (char * )p;
1895  return FALSE;
1896}
1897static BOOLEAN jjSTD_HILB(leftv res, leftv u, leftv v)
1898{
1899  ideal result;
1900  intvec *w=(intvec *)atGet(u,"isHomog");
1901  tHomog hom=testHomog;
1902  if (w!=NULL)
1903  {
1904    w=ivCopy(w);
1905    hom=isHomog;
1906  }
1907  result=kStd((ideal)(u->Data()),currQuotient,hom,&w,(intvec *)v->Data());
1908  idSkipZeroes(result);
1909  res->data = (char *)result;
1910  setFlag(res,FLAG_STD);
1911  if (w!=NULL) atSet(res,mstrdup("isHomog"),w,INTVEC_CMD);
1912  return FALSE;
1913}
1914static BOOLEAN jjSTD_1(leftv res, leftv u, leftv v)
1915{
1916  ideal result;
1917  intvec *w=(intvec *)atGet(u,"isHomog");
1918  tHomog hom=testHomog;
1919  if (w!=NULL)
1920  {
1921    w=ivCopy(w);
1922    hom=isHomog;
1923  }
1924  ideal i1=(ideal)u->Data();
1925  ideal i0=idInit(1,i1->rank);
1926  i0->m[0]=(poly)v->Data();
1927  i1=idSimpleAdd(i1,i0);
1928  i0->m[0]=NULL;
1929  idDelete(&i0);
1930  BITSET save_test=test;
1931  test|=Sy_bit(OPT_SB_1);
1932  result=kStd(i1,currQuotient,hom,&w);
1933  test=save_test;
1934  idDelete(&i1);
1935  idSkipZeroes(result);
1936  res->data = (char *)result;
1937  setFlag(res,FLAG_STD);
1938  if (w!=NULL) atSet(res,mstrdup("isHomog"),w,INTVEC_CMD);
1939  return FALSE;
1940}
1941static BOOLEAN jjVARSTR2(leftv res, leftv u, leftv v)
1942{
1943  idhdl h=(idhdl)u->data;
1944  int i=(int)v->Data();
1945  if ((0<i) && (i<=IDRING(h)->N))
1946    res->data=mstrdup(IDRING(h)->names[i-1]);
1947  else
1948  {
1949    Werror("var number %d out of range 1..%d",i,IDRING(h)->N);
1950    return TRUE;
1951  }
1952  return FALSE;
1953}
1954static BOOLEAN jjWEDGE(leftv res, leftv u, leftv v)
1955{
1956  res->data = (char *)mpWedge((matrix)u->Data(),(int)v->Data());
1957  return FALSE;
1958}
1959#define jjWRONG2 (proc2)jjWRONG
1960#define jjWRONG3 (proc3)jjWRONG
1961static BOOLEAN jjWRONG(leftv res, leftv u)
1962{
1963  return TRUE;
1964}
1965
1966static BOOLEAN jjLOAD_E(leftv res, leftv v, leftv u)
1967{
1968  char * s=(char *)u->Data();
1969  if(strcmp(s, "with")==0)
1970    return jjLOAD(res, v, TRUE);
1971  WerrorS("invalid second argument");
1972  WerrorS("load(\"libname\" [,\"with\"]);");
1973  return TRUE;
1974}
1975
1976/*=================== operations with 2 args.: table =================*/
1977
1978struct sValCmd2 dArith2[]=
1979{
1980// operations:
1981// proc        cmd              res             arg1        arg2
1982 {jjCOLCOL,    COLONCOLON,     ANY_TYPE,       DEF_CMD,    DEF_CMD      PROFILER}
1983,{jjPLUS_I,    '+',            INT_CMD,        INT_CMD,    INT_CMD        PROFILER}
1984,{jjPLUS_N,    '+',            NUMBER_CMD,     NUMBER_CMD, NUMBER_CMD     PROFILER}
1985,{jjPLUS_P,    '+',            POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
1986,{jjPLUS_P,    '+',            VECTOR_CMD,     VECTOR_CMD, VECTOR_CMD PROFILER}
1987,{jjPLUS_ID,   '+',            IDEAL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
1988,{jjPLUS_ID,   '+',            MODUL_CMD,      MODUL_CMD,  MODUL_CMD PROFILER}
1989,{jjPLUS_P_MA, '+',            MATRIX_CMD,     POLY_CMD,   MATRIX_CMD PROFILER}
1990,{jjPLUS_MA_P, '+',            MATRIX_CMD,     MATRIX_CMD, POLY_CMD PROFILER}
1991,{jjPLUS_MA,   '+',            MATRIX_CMD,     MATRIX_CMD, MATRIX_CMD PROFILER}
1992,{jjPLUS_S,    '+',            STRING_CMD,     STRING_CMD, STRING_CMD PROFILER}
1993,{jjOP_IV_I,   '+',            INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
1994,{jjOP_I_IV,   '+',            INTVEC_CMD,     INT_CMD,    INTVEC_CMD PROFILER}
1995,{jjOP_IM_I,   '+',            INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
1996,{jjOP_I_IM,   '+',            INTMAT_CMD,     INT_CMD,    INTMAT_CMD PROFILER}
1997,{jjPLUS_IV,   '+',            INTVEC_CMD,     INTVEC_CMD, INTVEC_CMD PROFILER}
1998,{jjPLUS_IV,   '+',            INTMAT_CMD,     INTMAT_CMD, INTMAT_CMD PROFILER}
1999,{lAdd,        '+',            LIST_CMD,       LIST_CMD,   LIST_CMD PROFILER}
2000,{jjRSUM,      '+',            RING_CMD,       RING_CMD,   RING_CMD PROFILER}
2001,{jjMINUS_I,   '-',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2002,{jjMINUS_N,   '-',            NUMBER_CMD,     NUMBER_CMD, NUMBER_CMD PROFILER}
2003,{jjMINUS_P,   '-',            POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2004,{jjMINUS_P,   '-',            VECTOR_CMD,     VECTOR_CMD, VECTOR_CMD PROFILER}
2005,{jjPLUS_MA_P, '-',            MATRIX_CMD,     MATRIX_CMD, POLY_CMD PROFILER}
2006,{jjMINUS_MA,  '-',            MATRIX_CMD,     MATRIX_CMD, MATRIX_CMD PROFILER}
2007,{jjOP_IV_I,   '-',            INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
2008,{jjOP_IM_I,   '-',            INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
2009,{jjMINUS_IV,  '-',            INTVEC_CMD,     INTVEC_CMD, INTVEC_CMD PROFILER}
2010,{jjMINUS_IV,  '-',            INTMAT_CMD,     INTMAT_CMD, INTMAT_CMD PROFILER}
2011,{jjTIMES_I,   '*',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2012,{jjTIMES_N,   '*',            NUMBER_CMD,     NUMBER_CMD, NUMBER_CMD PROFILER}
2013,{jjTIMES_P,   '*',            POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2014,{jjTIMES_P,   '*',            VECTOR_CMD,     POLY_CMD,   VECTOR_CMD PROFILER}
2015,{jjTIMES_P,   '*',            VECTOR_CMD,     VECTOR_CMD, POLY_CMD PROFILER}
2016,{jjTIMES_MA_P1,'*',           IDEAL_CMD,      IDEAL_CMD,  POLY_CMD PROFILER}
2017,{jjTIMES_MA_P2,'*',           IDEAL_CMD,      POLY_CMD,   IDEAL_CMD PROFILER}
2018,{jjTIMES_ID,  '*',            IDEAL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
2019,{jjTIMES_MA_P1,'*',           MODUL_CMD,      IDEAL_CMD,  VECTOR_CMD PROFILER}
2020,{jjTIMES_MA_P2,'*',           MODUL_CMD,      VECTOR_CMD, IDEAL_CMD PROFILER}
2021,{jjTIMES_ID,  '*',            MODUL_CMD,      IDEAL_CMD,  MODUL_CMD PROFILER}
2022,{jjTIMES_ID,  '*',            MODUL_CMD,      MODUL_CMD,  IDEAL_CMD PROFILER}
2023,{jjTIMES_MA_P1,'*',           MATRIX_CMD,     MATRIX_CMD, POLY_CMD PROFILER}
2024,{jjTIMES_MA_P2,'*',           MATRIX_CMD,     POLY_CMD,   MATRIX_CMD PROFILER}
2025,{jjTIMES_MA_N1,'*',           MATRIX_CMD,     MATRIX_CMD, NUMBER_CMD PROFILER}
2026,{jjTIMES_MA_N2,'*',           MATRIX_CMD,     NUMBER_CMD, MATRIX_CMD PROFILER}
2027,{jjTIMES_MA_I1,'*',           MATRIX_CMD,     MATRIX_CMD, INT_CMD PROFILER}
2028,{jjTIMES_MA_I2,'*',           MATRIX_CMD,     INT_CMD,    MATRIX_CMD PROFILER}
2029,{jjTIMES_MA,  '*',            MATRIX_CMD,     MATRIX_CMD, MATRIX_CMD PROFILER}
2030,{jjOP_IV_I,   '*',            INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
2031,{jjOP_I_IV,   '*',            INTVEC_CMD,     INT_CMD,    INTVEC_CMD PROFILER}
2032,{jjOP_IV_I,   '*',            INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
2033,{jjOP_I_IV,   '*',            INTMAT_CMD,     INT_CMD,    INTMAT_CMD PROFILER}
2034,{jjTIMES_IV,  '*',            INTVEC_CMD,     INTMAT_CMD, INTVEC_CMD PROFILER}
2035,{jjTIMES_IV,  '*',            INTMAT_CMD,     INTMAT_CMD, INTMAT_CMD PROFILER}
2036,{jjTIMES_IV,  '*',            INTMAT_CMD,     INTVEC_CMD, INTMAT_CMD PROFILER}
2037,{jjDIV_N,     '/',            NUMBER_CMD,     NUMBER_CMD, NUMBER_CMD PROFILER}
2038,{jjDIV_P,     '/',            POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2039,{jjDIV_P,     '/',            VECTOR_CMD,     VECTOR_CMD, POLY_CMD PROFILER}
2040,{jjDIV_I,     '/',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2041,{jjOP_IV_I,   '/',            INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
2042,{jjOP_IV_I,   '/',            INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
2043,{jjDIV_I,     INTDIV,         INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2044,{jjOP_IV_I,   INTDIV,         INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
2045,{jjOP_IV_I,   INTDIV,         INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
2046,{jjMOD_I,     '%',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2047,{jjOP_IV_I,   '%',            INTVEC_CMD,     INTVEC_CMD, INT_CMD PROFILER}
2048,{jjOP_IV_I,   '%',            INTMAT_CMD,     INTMAT_CMD, INT_CMD PROFILER}
2049,{jjPOWER_I,   '^',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2050,{jjPOWER_N,   '^',            NUMBER_CMD,     NUMBER_CMD, INT_CMD PROFILER}
2051,{jjPOWER_P,   '^',            POLY_CMD,       POLY_CMD,   INT_CMD PROFILER}
2052,{jjPOWER_ID,  '^',            IDEAL_CMD,      IDEAL_CMD,  INT_CMD PROFILER}
2053,{jjLE_I,      LE,             INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2054,{jjLE_N,      LE,             INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2055,{jjCOMPARE_S, LE,             INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2056,{jjCOMPARE_IV_I,LE,           INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2057,{jjCOMPARE_IV,LE,             INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2058,{jjCOMPARE_P, LE,             INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2059,{jjCOMPARE_P, LE,             INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2060,{jjLT_I,      '<',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2061,{jjLT_N,      '<',            INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2062,{jjCOMPARE_IV_I,'<',           INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2063,{jjCOMPARE_IV,'<',            INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2064,{jjCOMPARE_S, '<',            INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2065,{jjCOMPARE_P, '<',            INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2066,{jjCOMPARE_P, '<',            INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2067,{jjGE_I,      GE,             INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2068,{jjGE_N,      GE,             INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2069,{jjCOMPARE_S, GE,             INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2070,{jjCOMPARE_IV_I,GE,           INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2071,{jjCOMPARE_IV,GE,             INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2072,{jjCOMPARE_P, GE,             INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2073,{jjCOMPARE_P, GE,             INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2074,{jjGT_I,      '>',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2075,{jjGT_N,      '>',            INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2076,{jjCOMPARE_S, '>',            INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2077,{jjCOMPARE_IV_I,'>',          INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2078,{jjCOMPARE_IV,'>',            INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2079,{jjCOMPARE_P, '>',            INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2080,{jjCOMPARE_P, '>',            INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2081,{jjAND_I,     '&',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2082,{jjOR_I,      '|',            INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2083,{jjEQUAL_I,   EQUAL_EQUAL,    INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2084,{jjEQUAL_N,   EQUAL_EQUAL,    INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2085,{jjCOMPARE_S, EQUAL_EQUAL,    INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2086,{jjEQUAL_P,   EQUAL_EQUAL,    INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2087,{jjEQUAL_P,   EQUAL_EQUAL,    INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2088,{jjCOMPARE_IV_I,EQUAL_EQUAL,  INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2089,{jjCOMPARE_IV,EQUAL_EQUAL,    INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2090,{jjCOMPARE_IV,EQUAL_EQUAL,    INT_CMD,        INTMAT_CMD, INTMAT_CMD PROFILER}
2091,{jjEQUAL_Ma,  EQUAL_EQUAL,    INT_CMD,        MATRIX_CMD, MATRIX_CMD PROFILER}
2092,{jjWRONG2,    EQUAL_EQUAL,    0,              IDEAL_CMD,  IDEAL_CMD PROFILER}
2093,{jjWRONG2,    EQUAL_EQUAL,    0,              MODUL_CMD,  MODUL_CMD PROFILER}
2094,{jjEQUAL_I,   NOTEQUAL,       INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2095,{jjEQUAL_N,   NOTEQUAL,       INT_CMD,        NUMBER_CMD, NUMBER_CMD PROFILER}
2096,{jjCOMPARE_S, NOTEQUAL,       INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2097,{jjEQUAL_P,   NOTEQUAL,       INT_CMD,        POLY_CMD,   POLY_CMD PROFILER}
2098,{jjEQUAL_P,   NOTEQUAL,       INT_CMD,        VECTOR_CMD, VECTOR_CMD PROFILER}
2099,{jjCOMPARE_IV,NOTEQUAL,       INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2100,{jjCOMPARE_IV,NOTEQUAL,       INT_CMD,        INTMAT_CMD, INTMAT_CMD PROFILER}
2101,{jjEQUAL_Ma,  NOTEQUAL,       INT_CMD,        MATRIX_CMD, MATRIX_CMD PROFILER}
2102,{jjWRONG2,    NOTEQUAL,       0,              IDEAL_CMD,  IDEAL_CMD PROFILER}
2103,{jjWRONG2,    NOTEQUAL,       0,              MODUL_CMD,  MODUL_CMD PROFILER}
2104,{jjDOTDOT,    DOTDOT,         INTVEC_CMD,     INT_CMD,    INT_CMD PROFILER}
2105,{jjINDEX_I,   '[',            INT_CMD,        INTVEC_CMD, INT_CMD PROFILER}
2106,{jjINDEX_IV,  '[',            INT_CMD,        INTVEC_CMD, INTVEC_CMD PROFILER}
2107,{jjINDEX_I,   '[',            POLY_CMD,       IDEAL_CMD,  INT_CMD PROFILER}
2108,{jjINDEX_I,   '[',            POLY_CMD,       MAP_CMD,    INT_CMD PROFILER}
2109,{jjINDEX_IV,  '[',            POLY_CMD,       IDEAL_CMD,  INTVEC_CMD PROFILER}
2110,{jjINDEX_I,   '[',            VECTOR_CMD,     MODUL_CMD,  INT_CMD PROFILER}
2111,{jjINDEX_IV,  '[',            VECTOR_CMD,     MODUL_CMD,  INTVEC_CMD PROFILER}
2112,{jjINDEX_I,   '[',            STRING_CMD,     STRING_CMD, INT_CMD PROFILER}
2113,{jjINDEX_IV,  '[',            STRING_CMD,     STRING_CMD, INTVEC_CMD PROFILER}
2114,{jjINDEX_I,   '[',            ANY_TYPE/*set by p*/,LIST_CMD, INT_CMD PROFILER}
2115,{jjINDEX_IV,  '[',            ANY_TYPE/*set by p*/,LIST_CMD, INTVEC_CMD PROFILER}
2116,{jjINDEX_P,   '[',            POLY_CMD,       POLY_CMD,   INT_CMD PROFILER}
2117,{jjINDEX_P_IV,'[',            POLY_CMD,       POLY_CMD,   INTVEC_CMD PROFILER}
2118,{jjINDEX_V,   '[',            POLY_CMD,       VECTOR_CMD, INT_CMD PROFILER}
2119,{jjINDEX_V_IV,'[',            VECTOR_CMD,     VECTOR_CMD, INTVEC_CMD PROFILER}
2120,{jjPROC,      '(',            ANY_TYPE/*set by p*/,PROC_CMD, DEF_CMD PROFILER}
2121,{jjMAP,       '(',            ANY_TYPE/*set by p*/,MAP_CMD, DEF_CMD PROFILER}
2122,{jjKLAMMER,   '(',            ANY_TYPE/*set by p*/,ANY_TYPE, INT_CMD PROFILER}
2123,{jjKLAMMER_IV,'(',            ANY_TYPE/*set by p*/,ANY_TYPE, INTVEC_CMD PROFILER}
2124// and the procedures with 2 arguments:
2125,{atATTRIB2,   ATTRIB_CMD,     NONE/*set by p*/,DEF_CMD,   STRING_CMD PROFILER}
2126,{jjCOEF,      COEF_CMD,       MATRIX_CMD,     POLY_CMD,   POLY_CMD PROFILER}
2127,{jjCOEFFS_Id, COEFFS_CMD,     MATRIX_CMD,     IDEAL_CMD,  POLY_CMD PROFILER}
2128,{jjCOEFFS_Id, COEFFS_CMD,     MATRIX_CMD,     MODUL_CMD,  POLY_CMD PROFILER}
2129,{jjCONTRACT,  CONTRACT_CMD,   MATRIX_CMD,     IDEAL_CMD,  IDEAL_CMD PROFILER}
2130,{jjDEG_IV,    DEG_CMD,        INT_CMD,        POLY_CMD,   INTVEC_CMD PROFILER}
2131,{lDelete,     DELETE_CMD,     LIST_CMD,       LIST_CMD,   INT_CMD PROFILER}
2132,{jjDIFF_P,    DIFF_CMD,       POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2133,{jjDIFF_P,    DIFF_CMD,       VECTOR_CMD,     VECTOR_CMD, POLY_CMD PROFILER}
2134,{jjDIFF_ID,   DIFF_CMD,       IDEAL_CMD,      IDEAL_CMD,  POLY_CMD PROFILER}
2135,{jjDIFF_ID_ID,DIFF_CMD,       MATRIX_CMD,     IDEAL_CMD,  IDEAL_CMD PROFILER}
2136,{jjDIFF_ID,   DIFF_CMD,       MODUL_CMD,      MODUL_CMD,  POLY_CMD PROFILER}
2137,{jjDIFF_ID,   DIFF_CMD,       MATRIX_CMD,     MATRIX_CMD, POLY_CMD PROFILER}
2138,{jjDIM2,      DIM_CMD,        INT_CMD,        IDEAL_CMD,  IDEAL_CMD PROFILER}
2139,{jjDIM2,      DIM_CMD,        INT_CMD,        MODUL_CMD,  IDEAL_CMD PROFILER}
2140,{jjELIMIN,    ELIMINATION_CMD,IDEAL_CMD,      IDEAL_CMD,  POLY_CMD PROFILER}
2141,{jjELIMIN,    ELIMINATION_CMD,MODUL_CMD,      MODUL_CMD,  POLY_CMD PROFILER}
2142,{jjEXTGCD_I,  EXTGCD_CMD,     LIST_CMD,       INT_CMD,    INT_CMD PROFILER}
2143#ifdef HAVE_FACTORY
2144,{jjEXTGCD_P,  EXTGCD_CMD,     LIST_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2145,{jjSQR_FREE_DEC,FAC_CMD,      IDEAL_CMD,      POLY_CMD,   INT_CMD PROFILER}
2146,{jjFACSTD2,    FACSTD_CMD,    LIST_CMD,       IDEAL_CMD,  IDEAL_CMD PROFILER}
2147#else
2148,{jjWRONG2,    EXTGCD_CMD,     LIST_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2149,{jjWRONG2,    FAC_CMD,        IDEAL_CMD,      POLY_CMD,   INT_CMD PROFILER}
2150,{jjWRONG2,    FACSTD_CMD,     LIST_CMD,       IDEAL_CMD,  IDEAL_CMD PROFILER}
2151#endif
2152,{jjFETCH,     FETCH_CMD,      ANY_TYPE/*set by p*/,RING_CMD,  ANY_TYPE PROFILER}
2153,{jjFETCH,     FETCH_CMD,      ANY_TYPE/*set by p*/,QRING_CMD, ANY_TYPE PROFILER}
2154#ifdef HAVE_FGLM
2155,{fglmProc,    FGLM_CMD,       IDEAL_CMD,      RING_CMD,   DEF_CMD PROFILER}
2156,{fglmProc,    FGLM_CMD,       IDEAL_CMD,      QRING_CMD,  DEF_CMD PROFILER}
2157#else
2158,{jjWRONG2,    FGLM_CMD,       IDEAL_CMD,      RING_CMD,   DEF_CMD PROFILER}
2159,{jjWRONG2,    FGLM_CMD,       IDEAL_CMD,      QRING_CMD,  DEF_CMD PROFILER}
2160#endif
2161,{jjFIND2,     FIND_CMD,       INT_CMD,        STRING_CMD, STRING_CMD PROFILER}
2162,{jjGCD_I,     GCD_CMD,        INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2163#ifdef HAVE_FACTORY
2164,{jjGCD_P,     GCD_CMD,        POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2165#else
2166,{jjWRONG2,    GCD_CMD,        POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2167#endif
2168,{jjHILBERT2,  HILBERT_CMD,    INTVEC_CMD,     IDEAL_CMD,  INT_CMD PROFILER}
2169,{jjHILBERT2,  HILBERT_CMD,    INTVEC_CMD,     MODUL_CMD,  INT_CMD PROFILER}
2170,{jjHOMOG_P,   HOMOG_CMD,      POLY_CMD,       POLY_CMD,   POLY_CMD PROFILER}
2171,{jjHOMOG_P,   HOMOG_CMD,      VECTOR_CMD,     VECTOR_CMD, POLY_CMD PROFILER}
2172,{jjHOMOG_ID,  HOMOG_CMD,      IDEAL_CMD,      IDEAL_CMD,  POLY_CMD PROFILER}
2173,{jjHOMOG_ID,  HOMOG_CMD,      MODUL_CMD,      MODUL_CMD,  POLY_CMD PROFILER}
2174,{jjCALL2MANY, IDEAL_CMD,      IDEAL_CMD,      DEF_CMD,    DEF_CMD PROFILER}
2175,{jjFETCH,     IMAP_CMD,       ANY_TYPE/*set by p*/,RING_CMD,  ANY_TYPE PROFILER}
2176,{jjFETCH,     IMAP_CMD,       ANY_TYPE/*set by p*/,QRING_CMD, ANY_TYPE PROFILER}
2177,{jjINDEPSET2, INDEPSET_CMD,   LIST_CMD,       IDEAL_CMD,  INT_CMD PROFILER}
2178,{lInsert,     INSERT_CMD,     LIST_CMD,       LIST_CMD,   DEF_CMD PROFILER}
2179,{jjINTERSECT, INTERSECT_CMD,  IDEAL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
2180,{jjINTERSECT, INTERSECT_CMD,  MODUL_CMD,      MODUL_CMD,  MODUL_CMD PROFILER}
2181,{jjJET_P,     JET_CMD,        POLY_CMD,       POLY_CMD,   INT_CMD PROFILER}
2182,{jjJET_ID,    JET_CMD,        IDEAL_CMD,      IDEAL_CMD,  INT_CMD PROFILER}
2183,{jjJET_P,     JET_CMD,        VECTOR_CMD,     VECTOR_CMD, INT_CMD PROFILER}
2184,{jjJET_ID,    JET_CMD,        MODUL_CMD,      MODUL_CMD,  INT_CMD PROFILER}
2185,{jjKBASE2,    KBASE_CMD,      IDEAL_CMD,      IDEAL_CMD,  INT_CMD PROFILER}
2186,{jjKBASE2,    KBASE_CMD,      MODUL_CMD,      MODUL_CMD,  INT_CMD PROFILER}
2187,{atKILLATTR2, KILLATTR_CMD,   NONE,           IDHDL,      STRING_CMD PROFILER}
2188,{jjKoszul,    KOSZUL_CMD,     MATRIX_CMD,     INT_CMD,    INT_CMD PROFILER}
2189,{jjKoszul_Id, KOSZUL_CMD,     MATRIX_CMD,     INT_CMD,    IDEAL_CMD PROFILER}
2190,{jjLIFT,      LIFT_CMD,       MATRIX_CMD,     IDEAL_CMD,  IDEAL_CMD PROFILER}
2191,{jjLIFT,      LIFT_CMD,       MATRIX_CMD,     MODUL_CMD,  MODUL_CMD PROFILER}
2192,{jjLIFTSTD,   LIFTSTD_CMD,    IDEAL_CMD,      IDEAL_CMD,  MATRIX_CMD PROFILER}
2193,{jjLIFTSTD,   LIFTSTD_CMD,    MODUL_CMD,      MODUL_CMD,  MATRIX_CMD PROFILER}
2194,{jjCALL2MANY, LIST_CMD,       LIST_CMD,       DEF_CMD,    DEF_CMD PROFILER}
2195,{jjRES,       LRES_CMD,       RESOLUTION_CMD, IDEAL_CMD,  INT_CMD PROFILER}
2196,{jjRES,       LRES_CMD,       RESOLUTION_CMD, MODUL_CMD,  INT_CMD PROFILER}
2197,{jjMINOR,     MINOR_CMD,      IDEAL_CMD,      MATRIX_CMD, INT_CMD PROFILER}
2198,{jjCALL2MANY, MODUL_CMD,      MODUL_CMD,      DEF_CMD,    DEF_CMD PROFILER}
2199,{jjMODULO,    MODULO_CMD,     MODUL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
2200,{jjMODULO,    MODULO_CMD,     MODUL_CMD,      MODUL_CMD,  MODUL_CMD PROFILER}
2201,{jjMONITOR2,  MONITOR_CMD,    NONE,           STRING_CMD, STRING_CMD PROFILER}
2202//,{jjRES,       MRES_CMD,       LIST_CMD,       IDEAL_CMD,  INT_CMD PROFILER}
2203//,{jjRES,       MRES_CMD,       LIST_CMD,       MODUL_CMD,  INT_CMD PROFILER}
2204,{jjRES,       MRES_CMD,       RESOLUTION_CMD, IDEAL_CMD,  INT_CMD PROFILER}
2205,{jjRES,       MRES_CMD,       RESOLUTION_CMD, MODUL_CMD,  INT_CMD PROFILER}
2206,{jjPARSTR2,   PARSTR_CMD,     STRING_CMD,     RING_CMD,   INT_CMD PROFILER}
2207,{jjPARSTR2,   PARSTR_CMD,     STRING_CMD,     QRING_CMD,  INT_CMD PROFILER}
2208,{jjPRINT_FORMAT, PRINT_CMD,   NONE,           DEF_CMD,    STRING_CMD PROFILER}
2209,{jjQUOT,      QUOTIENT_CMD,   IDEAL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
2210,{jjQUOT,      QUOTIENT_CMD,   MODUL_CMD,      MODUL_CMD,  IDEAL_CMD PROFILER}
2211,{jjQUOT,      QUOTIENT_CMD,   IDEAL_CMD,      MODUL_CMD,  MODUL_CMD PROFILER}
2212,{jjRANDOM,    RANDOM_CMD,     INT_CMD,        INT_CMD,    INT_CMD PROFILER}
2213,{jjREAD2,     READ_CMD,       STRING_CMD,     LINK_CMD,   STRING_CMD PROFILER}
2214,{jjREDUCE_P,  REDUCE_CMD,     POLY_CMD,       POLY_CMD,   IDEAL_CMD PROFILER}
2215,{jjREDUCE_P,  REDUCE_CMD,     VECTOR_CMD,     VECTOR_CMD, IDEAL_CMD PROFILER}
2216,{jjREDUCE_P,  REDUCE_CMD,     VECTOR_CMD,     VECTOR_CMD, MODUL_CMD PROFILER}
2217,{jjREDUCE_ID, REDUCE_CMD,     IDEAL_CMD,      IDEAL_CMD,  IDEAL_CMD PROFILER}
2218,{jjREDUCE_ID, REDUCE_CMD,     MODUL_CMD,      MODUL_CMD,  MODUL_CMD PROFILER}
2219,{jjREDUCE_ID, REDUCE_CMD,     MODUL_CMD,      MODUL_CMD,  IDEAL_CMD PROFILER}
2220//,{jjRES,       RES_CMD,        LIST_CMD,       IDEAL_CMD,  INT_CMD PROFILER}
2221//,{jjRES,       RES_CMD,        LIST_CMD,       MODUL_CMD,  INT_CMD PROFILER}
2222,{jjRES,       RES_CMD,        RESOLUTION_CMD, IDEAL_CMD,  INT_CMD PROFILER}
2223,{jjRES,       RES_CMD,        RESOLUTION_CMD, MODUL_CMD,  INT_CMD PROFILER}
2224,{jjSTATUS2,   STATUS_CMD,     STRING_CMD,     LINK_CMD,   STRING_CMD PROFILER}
2225,{jjSIMPL_P,   SIMPLIFY_CMD,   POLY_CMD,       POLY_CMD,   INT_CMD PROFILER}
2226,{jjSIMPL_P,   SIMPLIFY_CMD,   VECTOR_CMD,     VECTOR_CMD, INT_CMD PROFILER}
2227,{jjSIMPL_ID,  SIMPLIFY_CMD,   IDEAL_CMD,      IDEAL_CMD,  INT_CMD PROFILER}
2228,{jjSIMPL_ID,  SIMPLIFY_CMD,   MODUL_CMD,      MODUL_CMD,  INT_CMD PROFILER}
2229//,{jjRES,       SRES_CMD,       LIST_CMD,       IDEAL_CMD,  INT_CMD PROFILER}
2230//,{jjRES,       SRES_CMD,       LIST_CMD,       MODUL_CMD,  INT_CMD PROFILER}
2231,{jjRES,       SRES_CMD,       RESOLUTION_CMD, IDEAL_CMD,  INT_CMD PROFILER}
2232,{jjRES,       SRES_CMD,       RESOLUTION_CMD, MODUL_CMD,  INT_CMD PROFILER}
2233,{jjCALL2MANY, SYSTEM_CMD,     ANY_TYPE/*set by p*/,STRING_CMD, DEF_CMD PROFILER}
2234,{jjSTD_1,     STD_CMD,        IDEAL_CMD,      IDEAL_CMD,  POLY_CMD PROFILER}
2235,{jjSTD_1,     STD_CMD,        MODUL_CMD,      MODUL_CMD,  VECTOR_CMD PROFILER}
2236,{jjSTD_HILB,  STD_CMD,        IDEAL_CMD,      IDEAL_CMD,  INTVEC_CMD PROFILER}
2237,{jjSTD_HILB,  STD_CMD,        MODUL_CMD,      MODUL_CMD,  INTVEC_CMD PROFILER}
2238,{jjVARSTR2,   VARSTR_CMD,     STRING_CMD,     RING_CMD,   INT_CMD PROFILER}
2239,{jjVARSTR2,   VARSTR_CMD,     STRING_CMD,     QRING_CMD,  INT_CMD PROFILER}
2240,{jjWEDGE,     WEDGE_CMD,      MATRIX_CMD,     MATRIX_CMD, INT_CMD PROFILER}
2241,{jjLOAD_E,    LOAD_CMD,       NONE,           STRING_CMD, STRING_CMD PROFILER}
2242,{NULL,        0,              0,              0,          0 PROFILER}
2243};
2244/*=================== operations with 1 arg.: static proc =================*/
2245static BOOLEAN jjDUMMY(leftv res, leftv u)
2246{
2247  res->data = (char *)u->CopyD();
2248  return FALSE;
2249}
2250static BOOLEAN jjNULL(leftv res, leftv u)
2251{
2252  return FALSE;
2253}
2254//static BOOLEAN jjPLUSPLUS(leftv res, leftv u)
2255//{
2256//  res->data = (char *)((int)u->Data()+1);
2257//  return FALSE;
2258//}
2259//static BOOLEAN jjMINUSMINUS(leftv res, leftv u)
2260//{
2261//  res->data = (char *)((int)u->Data()-1);
2262//  return FALSE;
2263//}
2264static BOOLEAN jjPLUSPLUS(leftv res, leftv u)
2265{
2266  if (IDTYP((idhdl)u->data)==INT_CMD)
2267  {
2268    if (iiOp==PLUSPLUS) IDINT((idhdl)u->data)++;
2269    else                IDINT((idhdl)u->data)--;
2270    return FALSE;
2271  }
2272  return TRUE;
2273}
2274static BOOLEAN jjUMINUS_I(leftv res, leftv u)
2275{
2276  res->data = (char *)(-(int)u->Data());
2277  return FALSE;
2278}
2279static BOOLEAN jjUMINUS_N(leftv res, leftv u)
2280{
2281  number n=(number)u->CopyD();
2282  n=nNeg(n);
2283  res->data = (char *)n;
2284  return FALSE;
2285}
2286static BOOLEAN jjUMINUS_P(leftv res, leftv u)
2287{
2288  res->data = (char *)pNeg((poly)u->CopyD(POLY_CMD));
2289  return FALSE;
2290}
2291static BOOLEAN jjUMINUS_MA(leftv res, leftv u)
2292{
2293  poly m1=pISet(-1);
2294  res->data = (char *)mpMultP((matrix)u->CopyD(MATRIX_CMD),m1);
2295  return FALSE;
2296}
2297static BOOLEAN jjUMINUS_IV(leftv res, leftv u)
2298{
2299  intvec *iv=(intvec *)u->CopyD(INTVEC_CMD);
2300  (*iv)*=(-1);
2301  res->data = (char *)iv;
2302  return FALSE;
2303}
2304static BOOLEAN jjPROC1(leftv res, leftv u)
2305{
2306  if ((u->rtyp!=IDHDL) || (u->e!=NULL))
2307    return TRUE;
2308#ifdef HAVE_NAMESPACES
2309  leftv sl = iiMake_proc((idhdl) u->data,u,NULL);
2310#else /* HAVE_NAMESPACES */
2311  leftv sl = iiMake_proc((idhdl) u->data,NULL);
2312#endif /* HAVE_NAMESPACES */
2313  if (sl==NULL)
2314  {
2315    return TRUE;
2316  }
2317  else
2318  {
2319    memcpy(res,sl,sizeof(sleftv));
2320  }
2321  return FALSE;
2322}
2323static BOOLEAN jjBAREISS(leftv res, leftv v)
2324{
2325  //matrix m=(matrix)v->Data();
2326  //lists l=mpBareiss(m,FALSE);
2327  lists l=smCallNewBareiss((ideal)v->Data(),0,0);
2328  res->data = (char *)l;
2329  return FALSE;
2330}
2331static BOOLEAN jjBAREISS_IM(leftv res, leftv v)
2332{
2333  intvec *m=(intvec *)v->CopyD();
2334  ivTriangMat(m,1,1);
2335  res->data = (char *)m;
2336  return FALSE;
2337}
2338static BOOLEAN jjCALL1MANY(leftv res, leftv u)
2339{
2340  return iiExprArithM(res,u,iiOp);
2341}
2342static BOOLEAN jjCHAR(leftv res, leftv v)
2343{
2344  res->data = (char *)rChar((ring)v->Data());
2345  return FALSE;
2346}
2347static BOOLEAN jjCOLS(leftv res, leftv v)
2348{
2349  res->data = (char *)MATCOLS((matrix)(v->Data()));
2350  return FALSE;
2351}
2352static BOOLEAN jjCOLS_IV(leftv res, leftv v)
2353{
2354  res->data = (char *)((intvec*)(v->Data()))->cols();
2355  return FALSE;
2356}
2357static BOOLEAN jjCONTENT(leftv res, leftv v)
2358{
2359  poly p=(poly)v->CopyD();
2360  if (p!=NULL) pCleardenom(p);
2361  res->data = (char *)p;
2362  return FALSE;
2363}
2364static BOOLEAN jjCOUNT_N(leftv res, leftv v)
2365{
2366  res->data = (char *)nSize((number)v->Data());
2367  return FALSE;
2368}
2369static BOOLEAN jjCOUNT_L(leftv res, leftv v)
2370{
2371  lists l=(lists)v->Data();
2372  res->data = (char *)(l->nr+1);
2373  return FALSE;
2374}
2375static BOOLEAN jjCOUNT_M(leftv res, leftv v)
2376{
2377  matrix m=(matrix)v->Data();
2378  res->data = (char *)(MATROWS(m)*MATCOLS(m));
2379  return FALSE;
2380}
2381static BOOLEAN jjCOUNT_IV(leftv res, leftv v)
2382{
2383  res->data = (char *)((intvec*)(v->Data()))->length();
2384  return FALSE;
2385}
2386static BOOLEAN jjDEG(leftv res, leftv v)
2387{
2388  int dummy;
2389  poly p=(poly)v->Data();
2390  if (p!=NULL) res->data = (char *)pLDeg(p,&dummy);
2391  else res->data=(char *)-1;
2392  return FALSE;
2393}
2394static BOOLEAN jjDEGREE(leftv res, leftv v)
2395{
2396  assumeStdFlag(v);
2397  scDegree((ideal)v->Data(),currQuotient);
2398  return FALSE;
2399}
2400static BOOLEAN jjDEFINED(leftv res, leftv v)
2401{
2402  if ((v->rtyp==IDHDL)
2403  && ((myynest==IDLEV((idhdl)v->data))||(0==IDLEV((idhdl)v->data))))
2404  {
2405    res->data=(void *)(IDLEV((idhdl)v->data)+1);
2406  }
2407  else if (v->rtyp!=0) res->data=(void *)(-1);
2408  return FALSE;
2409}
2410#ifdef HAVE_FACTORY
2411static BOOLEAN jjDET(leftv res, leftv v)
2412{
2413  int i,j;
2414  matrix m=(matrix)v->Data();
2415  if (currRing->parameter==NULL)
2416  {
2417    for(i=m->rows();i>0;i--)
2418    {
2419      for(j=m->cols();j>0;j--)
2420      {
2421        if((MATELEM(m,i,j)!=NULL)
2422        && (!pIsConstant(MATELEM(m,i,j))))
2423        {
2424          goto nonconst;
2425        }
2426      }
2427    }
2428    res->data = (char *)singclap_det(m);
2429    return FALSE;
2430  }
2431nonconst:
2432//  res->data = (char *)mpDet(m);
2433  poly p = mpDetBareiss(m);
2434#ifdef PDEBUG
2435  poly old = mpDet(m);
2436  if (!pEqual(p,old))
2437  {
2438    WerrorS("error in mpDetBareiss");
2439  }
2440#endif
2441  res ->data = (char *)p;
2442  return FALSE;
2443}
2444static BOOLEAN jjDET_I(leftv res, leftv v)
2445{
2446  intvec * m=(intvec*)v->Data();
2447  int i,j;
2448  i=m->rows();j=m->cols();
2449  if(i==j)
2450    res->data = (char *)singclap_det_i(m);
2451  else
2452  {
2453    Werror("det of %d x %d intmat",i,j);
2454    return TRUE;
2455  }
2456  return FALSE;
2457}
2458#endif
2459static BOOLEAN jjDIM(leftv res, leftv v)
2460{
2461  assumeStdFlag(v);
2462  res->data = (char *)scDimInt((ideal)(v->Data()),currQuotient);
2463  return FALSE;
2464}
2465static BOOLEAN jjDUMP(leftv res, leftv v)
2466{
2467  si_link l = (si_link)v->Data();
2468  if (slDump(l))
2469  {
2470    const char *s;
2471    if ((l!=NULL)&&(l->name!=NULL)) s=l->name;
2472    else                            s=sNoName;
2473    Werror("cannot dump to `%s`",s);
2474    return TRUE;
2475  }
2476  else
2477    return FALSE;
2478}
2479static BOOLEAN jjE(leftv res, leftv v)
2480{
2481  res->data = (char *)pOne();
2482  pSetComp((poly)res->data,(int)v->Data());
2483  return FALSE;
2484}
2485#ifdef HAVE_FACTORY
2486static BOOLEAN jjFACSTD(leftv res, leftv v)
2487{
2488  res->data=(void *)kStdfac((ideal)v->Data(),NULL,testHomog,NULL);
2489  setFlag(res,FLAG_STD);
2490  return FALSE;
2491}
2492#endif
2493static BOOLEAN jjGETDUMP(leftv res, leftv v)
2494{
2495  si_link l = (si_link)v->Data();
2496  if (slGetDump(l))
2497  {
2498    const char *s;
2499    if ((l!=NULL)&&(l->name!=NULL)) s=l->name;
2500    else                            s=sNoName;
2501    Werror("cannot get dump from `%s`",s);
2502    return TRUE;
2503  }
2504  else
2505    return FALSE;
2506}
2507static BOOLEAN jjHIGHCORNER(leftv res, leftv v)
2508{
2509  assumeStdFlag(v);
2510  ideal I=(ideal)v->Data();
2511  res->data=(void *)iiHighCorner(I,0);
2512  return FALSE;
2513}
2514static BOOLEAN jjHIGHCORNER_M(leftv res, leftv v)
2515{
2516  assumeStdFlag(v);
2517  intvec *w=(intvec*)atGet(v,"isHomog");
2518  intvec *module_w=NULL;
2519  ideal I=(ideal)v->Data();
2520  int i;
2521  poly p=NULL,po=NULL;
2522  int rk=idRankFreeModule(I);
2523  if (w!=NULL)
2524    module_w=new intvec(*w);
2525  else
2526    module_w = new intvec(rk);
2527  for(i=rk;i>0;i--)
2528  {
2529    p=iiHighCorner(I,i);
2530    if (p==NULL)
2531    {
2532      Werror("module must be zero-dimensional");
2533      delete module_w;
2534      return TRUE;
2535    }
2536    if (po==NULL)
2537    {
2538      po=p;
2539    }
2540    else
2541    {
2542      // now po!=NULL, p!=NULL
2543      int d=(pFDeg(po)+(*module_w)[pGetComp(po)-1] - pFDeg(p)-(*module_w)[i-1]);
2544      if (d==0)
2545        d=pComp0(p,po);
2546      if (d > 0)
2547      {
2548        pDelete(&po); po=p;
2549      }
2550      else // (d < 0)
2551      {
2552        pDelete(&p);
2553      }
2554    }
2555  }
2556  delete module_w;
2557  res->data=(void *)po;
2558  return FALSE;
2559}
2560static BOOLEAN jjHILBERT(leftv res, leftv v)
2561{
2562  assumeStdFlag(v);
2563  intvec *module_w=(intvec*)atGet(v,"isHomog");
2564  //scHilbertPoly((ideal)v->Data(),currQuotient);
2565  hLookSeries((ideal)v->Data(),module_w,currQuotient);
2566  return FALSE;
2567}
2568static BOOLEAN jjHILBERT_IV(leftv res, leftv v)
2569{
2570  res->data=(void *)hSecondSeries((intvec *)v->Data());
2571  return FALSE;
2572}
2573static BOOLEAN jjHOMOG1(leftv res, leftv v)
2574{
2575  intvec *w;
2576  res->data=(void *)idHomModule((ideal)v->Data(),currQuotient,&w);
2577  if ((v->rtyp==IDHDL)&&(w!=NULL))
2578    atSet((idhdl)v->data,mstrdup("isHomog"),w,INTVEC_CMD);
2579  return FALSE;
2580}
2581static BOOLEAN jjIDEAL_Ma(leftv res, leftv v)
2582{
2583  matrix mat=(matrix)v->CopyD();
2584  IDELEMS((ideal)mat)=MATCOLS(mat)*MATROWS(mat);
2585  MATROWS(mat)=1;
2586  mat->rank=1;
2587  res->data=(char *)mat;
2588  return FALSE;
2589}
2590static BOOLEAN jjIDEAL_Map(leftv res, leftv v)
2591{
2592  map m=(map)v->CopyD();
2593  FreeL((ADDRESS)m->preimage);
2594  m->preimage=NULL;
2595  ideal I=(ideal)m;
2596  I->rank=1;
2597  res->data=(char *)I;
2598  return FALSE;
2599}
2600static BOOLEAN jjIDEAL_R(leftv res, leftv v)
2601{
2602  if (currRing!=NULL)
2603  {
2604    ring q=(ring)v->Data();
2605    if ((rInternalChar(q)==rInternalChar(currRing))
2606    && (q->N==currRing->N)) /* && varnames equ. */
2607    {
2608      if (q->qideal==NULL)
2609        res->data=(char *)idInit(1,1);
2610      else
2611        res->data=(char *)idCopy(q->qideal);
2612      return FALSE;
2613    }
2614  }
2615  return TRUE;
2616}
2617static BOOLEAN jjIm2Iv(leftv res, leftv v)
2618{
2619  intvec *iv = (intvec *)v->CopyD();
2620  iv->makeVector();
2621  res->data = iv;
2622  return FALSE;
2623}
2624static BOOLEAN jjINDEPSET(leftv res, leftv v)
2625{
2626  assumeStdFlag(v);
2627  res->data=(void *)scIndIntvec((ideal)(v->Data()),currQuotient);
2628  return FALSE;
2629}
2630static BOOLEAN jjINTERRED(leftv res, leftv v)
2631{
2632  ideal result=kInterRed((ideal)(v->Data()),currQuotient);
2633  //if (result==NULL) return TRUE;
2634  idSkipZeroes(result);
2635  res->data = result;
2636  return FALSE;
2637}
2638static BOOLEAN jjIS_RINGVAR_P(leftv res, leftv v)
2639{
2640  res->data = (char *)pVar((poly)v->Data());
2641  return FALSE;
2642}
2643static BOOLEAN jjIS_RINGVAR_S(leftv res, leftv v)
2644{
2645  res->data = (char *)(rIsRingVar((char *)v->Data())+1);
2646  return FALSE;
2647}
2648static BOOLEAN jjIS_RINGVAR0(leftv res, leftv v)
2649{
2650  res->data = (char *)0;
2651  return FALSE;
2652}
2653static BOOLEAN jjJACOB_P(leftv res, leftv v)
2654{
2655  ideal i=idInit(pVariables,1);
2656  int k;
2657  for (k=pVariables;k>0;k--)
2658  {
2659    i->m[k-1]=pDiff(pCopy((poly)(v->Data())),k);
2660  }
2661  res->data = (char *)i;
2662  return FALSE;
2663}
2664static BOOLEAN jjKBASE(leftv res, leftv v)
2665{
2666  assumeStdFlag(v);
2667  res->data = (char *)scKBase(-1,(ideal)(v->Data()),currQuotient);
2668  return FALSE;
2669}
2670static BOOLEAN jjKLAMMER_LIB(leftv res, leftv u)
2671{
2672  char * s=(char *)u->Data();
2673  return iiLibCmd(s);
2674}
2675#ifdef MDEBUG
2676static BOOLEAN jjpHead(leftv res, leftv v)
2677{
2678  res->data=(char *)pHead((poly)v->Data());
2679  return FALSE;
2680}
2681#endif
2682static BOOLEAN jjL2R(leftv res, leftv v)
2683{
2684  res->data=(char *)syConvList((lists)v->Data());
2685  return FALSE;
2686}
2687static BOOLEAN jjLEADCOEF(leftv res, leftv v)
2688{
2689  poly p=(poly)v->Data();
2690  if (p==NULL)
2691  {
2692    res->data=(char *)nInit(0);
2693  }
2694  else
2695  {
2696    res->data=(char *)nCopy(pGetCoeff(p));
2697  }
2698  return FALSE;
2699}
2700static BOOLEAN jjLEADEXP(leftv res, leftv v)
2701{
2702  poly p=(poly)v->Data();
2703  int s=pVariables;
2704  if (v->Typ()==VECTOR_CMD) s++;
2705  intvec *iv=new intvec(s);
2706  if (p!=NULL)
2707  {
2708    for(int i = pVariables;i;i--)
2709    {
2710      (*iv)[i-1]=pGetExp(p,i);
2711    }
2712    if (s!=pVariables)
2713      (*iv)[pVariables]=pGetComp(p);
2714  }
2715  res->data=(char *)iv;
2716  return FALSE;
2717}
2718static BOOLEAN jjLIB(leftv res, leftv v)
2719{
2720  return iiLibCmd((char *)v->CopyD());
2721}
2722static BOOLEAN jjMEMORY(leftv res, leftv v)
2723{
2724#ifdef MM_COUNT
2725  switch(((int)v->Data()))
2726  {
2727  case 0:
2728    res->data = (char *)mmMemUsed();
2729    break;
2730  case 1:
2731    res->data = (char *)mmMemAlloc();
2732    break;
2733#ifdef HAVE_SBRK
2734  case 2:
2735    res->data = (char *)mmMemPhysical();
2736    break;
2737#endif
2738  default:
2739#ifndef MAKE_DISRIBUTION
2740    mmPrintStat();
2741#endif
2742    res->data = (char *)0;
2743  }
2744  return FALSE;
2745#else
2746  res->data = (char *)0;
2747  return FALSE;
2748#endif
2749}
2750static BOOLEAN jjMONITOR1(leftv res, leftv v)
2751{
2752  monitor((char *)(v->Data()),PROT_I);
2753  return FALSE;
2754}
2755static BOOLEAN jjMSTD(leftv res, leftv v)
2756{
2757  res->data=(void *)min_std((ideal)v->Data(),currQuotient,testHomog,NULL);
2758  return FALSE;
2759}
2760static BOOLEAN jjMULT(leftv res, leftv v)
2761{
2762  assumeStdFlag(v);
2763  res->data = (char *)scMultInt((ideal)(v->Data()),currQuotient);
2764  return FALSE;
2765}
2766static BOOLEAN jjNAMEOF(leftv res, leftv v)
2767{
2768  res->data = (char *)v->name;
2769  if (res->data==NULL) res->data=mstrdup("");
2770  v->name=NULL;
2771  return FALSE;
2772}
2773static BOOLEAN jjNAMES(leftv res, leftv v)
2774{
2775  res->data=ipNameList(((ring)v->Data())->idroot);
2776  return FALSE;
2777}
2778static BOOLEAN jjNVARS(leftv res, leftv v)
2779{
2780  res->data = (char *)(((ring)(v->Data()))->N);
2781  return FALSE;
2782}
2783static BOOLEAN jjOpenClose(leftv res, leftv v)
2784{
2785  si_link l=(si_link)v->Data();
2786  if (iiOp==OPEN_CMD) return slOpen(l, SI_LINK_OPEN);
2787  else                return slClose(l);
2788}
2789static BOOLEAN jjORD(leftv res, leftv v)
2790{
2791  poly p=(poly)v->Data();
2792  res->data=(char *)( p==NULL ? -1 : pFDeg(p) );
2793  return FALSE;
2794}
2795static BOOLEAN jjPAR1(leftv res, leftv v)
2796{
2797  int i=(int)v->Data();
2798  int p=0;
2799  p=rPar(currRing);
2800  if ((0<i) && (i<=p))
2801  {
2802    res->data=(char *)nPar(i);
2803  }
2804  else
2805  {
2806    Werror("par number %d out of range 1..%d",i,p);
2807    return TRUE;
2808  }
2809  return FALSE;
2810}
2811static BOOLEAN jjPARDEG(leftv res, leftv v)
2812{
2813  res->data = (char *)nParDeg((number)v->Data());
2814  return FALSE;
2815}
2816static BOOLEAN jjPARSTR1(leftv res, leftv v)
2817{
2818  if (currRing==NULL)
2819  {
2820    WerrorS("no ring active");
2821    return TRUE;
2822  }
2823  int i=(int)v->Data();
2824  int p=0;
2825  if ((0<i) && (currRing->parameter!=NULL) && (i<=(p=rPar(currRing))))
2826    res->data=mstrdup(currRing->parameter[i-1]);
2827  else
2828  {
2829    Werror("par number %d out of range 1..%d",i,p);
2830    return TRUE;
2831  }
2832  return FALSE;
2833}
2834static BOOLEAN jjP2I(leftv res, leftv v)
2835{
2836  poly p=(poly)v->Data();
2837  if (p==NULL) { /*res->data=(char *)0;*/ return FALSE; }
2838  if ((pNext(p)!=NULL)|| (!pIsConstant(p)))
2839  {
2840    WerrorS("poly must be constant");
2841    return TRUE;
2842  }
2843  res->data = (char *)nInt(pGetCoeff(p));
2844  return FALSE;
2845}
2846static BOOLEAN jjPRIME(leftv res, leftv v)
2847{
2848  int i = IsPrime((int)(v->Data()));
2849  res->data = (char *)(i > 1 ? i : 2);
2850  return FALSE;
2851}
2852static BOOLEAN jjP2N(leftv res, leftv v)
2853{
2854  number n;
2855  poly p;
2856  if (((p=(poly)v->Data())!=NULL)
2857  && (pIsConstant(p)))
2858  {
2859    n=nCopy(pGetCoeff(p));
2860  }
2861  else
2862  {
2863    n=nInit(0);
2864  }
2865  res->data = (char *)n;
2866  return FALSE;
2867}
2868static BOOLEAN jjRESERVEDNAME(leftv res, leftv v)
2869{
2870  char *s= (char *)v->Data();
2871  int i = 1;
2872  int l = strlen(s);
2873  while (cmds[i].tokval!=0)
2874  {
2875    if (strcmp(s, cmds[i].name) == 0)
2876    {
2877      res->data = (char *)1;
2878      return FALSE;
2879    }
2880    i++;
2881  }
2882  //res->data = (char *)0;
2883  return FALSE;
2884}
2885static BOOLEAN jjREAD(leftv res, leftv v)
2886{
2887  return jjREAD2(res,v,NULL);
2888}
2889static BOOLEAN jjREGULARITY(leftv res, leftv v)
2890{
2891  res->data = (char *)iiRegularity((lists)v->Data());
2892  return FALSE;
2893}
2894static BOOLEAN jjROWS(leftv res, leftv v)
2895{
2896  ideal i = (ideal)v->Data();
2897  res->data = (char *)i->rank;
2898  return FALSE;
2899}
2900static BOOLEAN jjROWS_IV(leftv res, leftv v)
2901{
2902  res->data = (char *)((intvec*)(v->Data()))->rows();
2903  return FALSE;
2904}
2905static BOOLEAN jjRPAR(leftv res, leftv v)
2906{
2907  res->data = (char *)rPar(((ring)v->Data()));
2908  return FALSE;
2909}
2910static BOOLEAN jjSTD(leftv res, leftv v)
2911{
2912  ideal result;
2913  intvec *w=(intvec *)atGet(v,"isHomog");
2914  tHomog hom=testHomog;
2915  if (w!=NULL)
2916  {
2917    w=ivCopy(w);
2918    hom=isHomog;
2919  }
2920  //if (hasFlag(v,FLAG_STD))
2921  //  result=stdred((ideal)v->Data(),currQuotient,hom,&w);
2922  //else
2923    result=kStd((ideal)(v->Data()),currQuotient,hom,&w);
2924  idSkipZeroes(result);
2925  res->data = (char *)result;
2926  setFlag(res,FLAG_STD);
2927  if (w!=NULL) atSet(res,mstrdup("isHomog"),w,INTVEC_CMD);
2928  return FALSE;
2929}
2930static BOOLEAN jjSTRING(leftv res, leftv v)
2931{
2932  res->data = (char *)v->String();
2933  if (res->data==NULL)
2934  {
2935    Werror("cannot convert %s to string",Tok2Cmdname(v->Typ()));
2936    return TRUE;
2937  }
2938  return FALSE;
2939}
2940static BOOLEAN jjSort_Id(leftv res, leftv v)
2941{
2942  res->data = (char *)idSort((ideal)v->Data());
2943  return FALSE;
2944}
2945static BOOLEAN jjSTRING_PROC(leftv res, leftv v)
2946{
2947  procinfov pi = IDPROC((idhdl)v->data);
2948  if((pi->language == LANG_SINGULAR) && (pi->data.s.body!=NULL))
2949  //if(pi->language == LANG_SINGULAR)
2950  {
2951    //if(pi->data.s.body==NULL)
2952    //  iiGetLibProcBuffer(IDPROC((idhdl)v->data));
2953    res->data=mstrdup(pi->data.s.body);
2954  }
2955  else
2956    res->data=mstrdup("");
2957  return FALSE;
2958}
2959static BOOLEAN jjSYZYGY(leftv res, leftv v)
2960{
2961  intvec *w=NULL;
2962  res->data = (char *)idSyzygies((ideal)v->Data(),currQuotient,testHomog,&w);
2963  if (w!=NULL) delete w;
2964  return FALSE;
2965}
2966static BOOLEAN jjTRACE_IV(leftv res, leftv v)
2967{
2968  res->data = (char *)ivTrace((intvec*)(v->Data()));
2969  return FALSE;
2970}
2971static BOOLEAN jjTRANSP_IV(leftv res, leftv v)
2972{
2973  res->data = (char *)ivTranp((intvec*)(v->Data()));
2974  return FALSE;
2975}
2976static BOOLEAN jjTYPEOF(leftv res, leftv v)
2977{
2978  switch ((int)v->data)
2979  {
2980    case INT_CMD:     res->data=mstrdup("int"); break;
2981    case POLY_CMD:   res->data=mstrdup("poly"); break;
2982    case VECTOR_CMD:  res->data=mstrdup("vector"); break;
2983    case STRING_CMD:  res->data=mstrdup("string"); break;
2984    case INTVEC_CMD:  res->data=mstrdup("intvec"); break;
2985    case IDEAL_CMD:   res->data=mstrdup("ideal"); break;
2986    case MATRIX_CMD: res->data=mstrdup("matrix"); break;
2987    case MODUL_CMD:   res->data=mstrdup("module"); break;
2988    case MAP_CMD:     res->data=mstrdup("map"); break;
2989    case PROC_CMD:    res->data=mstrdup("proc"); break;
2990    case RING_CMD:    res->data=mstrdup("ring"); break;
2991    case QRING_CMD:   res->data=mstrdup("qring"); break;
2992    case INTMAT_CMD:  res->data=mstrdup("intmat"); break;
2993    case NUMBER_CMD:  res->data=mstrdup("number"); break;
2994    case LIST_CMD:   res->data=mstrdup("list"); break;
2995    case PACKAGE_CMD: res->data=mstrdup("package"); break;
2996    case LINK_CMD:   res->data=mstrdup("link"); break;
2997    case RESOLUTION_CMD:res->data=mstrdup("resolution");break;
2998    case DEF_CMD:
2999    case NONE:    res->data=mstrdup("none"); break;
3000    default:       res->data=mstrdup("?unknown type?");
3001  }
3002  return FALSE;
3003}
3004static BOOLEAN jjVAR1(leftv res, leftv v)
3005{
3006  int i=(int)v->Data();
3007  if ((0<i) && (i<=currRing->N))
3008  {
3009    poly p=pOne();
3010    pSetExp(p,i,1);
3011    pSetm(p);
3012    res->data=(char *)p;
3013  }
3014  else
3015  {
3016    Werror("var number %d out of range 1..%d",i,currRing->N);
3017    return TRUE;
3018  }
3019  return FALSE;
3020}
3021static BOOLEAN jjVARSTR1(leftv res, leftv v)
3022{
3023  if (currRing==NULL)
3024  {
3025    WerrorS("no ring active");
3026    return TRUE;
3027  }
3028  int i=(int)v->Data();
3029  if ((0<i) && (i<=currRing->N))
3030    res->data=mstrdup(currRing->names[i-1]);
3031  else
3032  {
3033    Werror("var number %d out of range 1..%d",i,currRing->N);
3034    return TRUE;
3035  }
3036  return FALSE;
3037}
3038static BOOLEAN jjVDIM(leftv res, leftv v)
3039{
3040  assumeStdFlag(v);
3041  res->data = (char *)scMult0Int((ideal)v->Data(),currQuotient);
3042  return FALSE;
3043}
3044
3045static BOOLEAN jjLOAD1(leftv res, leftv v)
3046{
3047  return jjLOAD(res, v);
3048}
3049
3050static BOOLEAN jjLOAD(leftv res, leftv v, BOOLEAN autoexport)
3051{
3052  char * s=mstrdup((char *)v->Data());
3053  char libnamebuf[256];
3054  lib_types LT = type_of_LIB(s, libnamebuf);
3055  BOOLEAN result = TRUE;
3056#ifdef HAVE_DYNAMIC_LOADING
3057  extern BOOLEAN load_modules(char *newlib, char *fullpath, BOOLEAN tellerror);
3058#endif /* HAVE_DYNAMIC_LOADING */
3059
3060  switch(LT)
3061  {
3062      default:
3063      case LT_NONE:
3064        Werror("%s: unknown type", s);
3065        break;
3066
3067      case LT_SINGULAR:
3068#ifdef HAVE_NAMESPACES
3069        result = iiLibCmd(s, autoexport);
3070#else
3071        result = iiLibCmd(s);
3072#endif
3073        break;
3074
3075      case LT_ELF:
3076#ifdef HAVE_DYNAMIC_LOADING
3077        result = load_modules(s, libnamebuf, autoexport);
3078#else /* HAVE_DYNAMIC_LOADING */
3079        WerrorS("Dynamic modules are not supported by this version of Singular");
3080#endif /* HAVE_DYNAMIC_LOADING */
3081        break;
3082
3083      case LT_HPUX:
3084#ifdef HAVE_DYNAMIC_LOADING
3085        result = load_modules(s, libnamebuf, autoexport);
3086#else /* HAVE_DYNAMIC_LOADING */
3087        WerrorS("Dynamic modules are not supported by this version of Singular");
3088#endif /* HAVE_DYNAMIC_LOADING */
3089        break;
3090  }
3091  return result;
3092}
3093
3094/*=================== operations with 1 arg.: table =================*/
3095
3096#ifdef INIT_BUG
3097#define XS(A) -((short)A)
3098#define jjstrlen       (proc1)1
3099#define jjpLength      (proc1)2
3100#define jjidElem       (proc1)3
3101#define jjmpDetBareiss (proc1)4
3102#define jjidFreeModule (proc1)5
3103#define jjidVec2Ideal  (proc1)6
3104#define jjrCharStr     (proc1)7
3105#ifndef MDEBUG
3106#define jjpHead        (proc1)8
3107#endif
3108#define jjidHead       (proc1)9
3109#define jjidMaxIdeal   (proc1)10
3110#define jjidMinBase    (proc1)11
3111#define jjsyMinBase    (proc1)12
3112#define jjpMaxComp     (proc1)13
3113#define jjmpTrace      (proc1)14
3114#define jjmpTransp     (proc1)15
3115#define jjrOrdStr      (proc1)16
3116#define jjrVarStr      (proc1)18
3117#define jjrParStr      (proc1)19
3118#define jjidMinEmbedding (proc1)20
3119#define jjBETTI_R        (proc1)21
3120#define jjCOUNT_R        (proc1)22
3121#define jjDIM_R          (proc1)23
3122#define jjMINRES_R       (proc1)24
3123#define jjidTransp       (proc1)25
3124
3125extern struct sValCmd1 dArith1[];
3126void jjInitTab1()
3127{
3128  int i=0;
3129  for (;dArith1[i].cmd!=0;i++)
3130  {
3131    if (dArith1[i].res<0)
3132    {
3133      switch ((int)dArith1[i].p)
3134      {
3135        case (int)jjstrlen:       dArith1[i].p=(proc1)strlen; break;
3136        case (int)jjpLength:      dArith1[i].p=(proc1)pLength; break;
3137        case (int)jjidElem:       dArith1[i].p=(proc1)idElem; break;
3138        case (int)jjidVec2Ideal:  dArith1[i].p=(proc1)idVec2Ideal; break;
3139#ifndef HAVE_FACTORY
3140        case (int)jjmpDetBareiss: dArith1[i].p=(proc1)mpDetBareiss; break;
3141#endif
3142        case (int)jjidFreeModule: dArith1[i].p=(proc1)idFreeModule; break;
3143        case (int)jjrCharStr:     dArith1[i].p=(proc1)rCharStr; break;
3144#ifndef MDEBUG
3145        case (int)jjpHead:        dArith1[i].p=(proc1)pHeadProc; break;
3146#endif
3147        case (int)jjidHead:       dArith1[i].p=(proc1)idHead; break;
3148        case (int)jjidMaxIdeal:   dArith1[i].p=(proc1)idMaxIdeal; break;
3149        case (int)jjidMinBase:    dArith1[i].p=(proc1)idMinBase; break;
3150        case (int)jjsyMinBase:    dArith1[i].p=(proc1)syMinBase; break;
3151        case (int)jjpMaxComp:     dArith1[i].p=(proc1)pMaxComp; break;
3152        case (int)jjmpTrace:      dArith1[i].p=(proc1)mpTrace; break;
3153        case (int)jjmpTransp:     dArith1[i].p=(proc1)mpTransp; break;
3154        case (int)jjrOrdStr:      dArith1[i].p=(proc1)rOrdStr; break;
3155        case (int)jjrVarStr:      dArith1[i].p=(proc1)rVarStr; break;
3156        case (int)jjrParStr:      dArith1[i].p=(proc1)rParStr; break;
3157        case (int)jjidMinEmbedding: dArith1[i].p=(proc1)idMinEmbedding; break;
3158        case (int)jjBETTI_R:      dArith1[i].p=(proc1)syBettiOfComputation; break;
3159        case (int)jjCOUNT_R:      dArith1[i].p=(proc1)syLength; break;
3160        case (int)jjDIM_R:        dArith1[i].p=(proc1)syDim; break;
3161        case (int)jjMINRES_R:     dArith1[i].p=(proc1)syMinimize; break;
3162        case (int)jjidTransp:     dArith1[i].p=(proc1)idTransp; break;
3163#ifdef GENTABLE
3164        default: Werror("missing proc1-definition for %d",(int)dArith1[i].p);
3165#endif
3166      }
3167    }
3168  }
3169}
3170#else
3171#if defined(PROC_BUG)
3172#define XS(A) A
3173static BOOLEAN jjstrlen(leftv res, leftv v)
3174{
3175  res->data = (char *)strlen((char *)v->Data());
3176  return FALSE;
3177}
3178static BOOLEAN jjpLength(leftv res, leftv v)
3179{
3180  res->data = (char *)pLength((poly)v->Data());
3181  return FALSE;
3182}
3183static BOOLEAN jjidElem(leftv res, leftv v)
3184{
3185  res->data = (char *)idElem((ideal)v->Data());
3186  return FALSE;
3187}
3188static BOOLEAN jjmpDetBareiss(leftv res, leftv v)
3189{
3190  res->data = (char *)mpDetBareiss((matrix)v->Data());
3191  return FALSE;
3192}
3193static BOOLEAN jjidFreeModule(leftv res, leftv v)
3194{
3195  res->data = (char *)idFreeModule((int)v->Data());
3196  return FALSE;
3197}
3198static BOOLEAN jjidVec2Ideal(leftv res, leftv v)
3199{
3200  res->data = (char *)idVec2Ideal((poly)v->Data());
3201  return FALSE;
3202}
3203static BOOLEAN jjrCharStr(leftv res, leftv v)
3204{
3205  res->data = rCharStr((ring)v->Data());
3206  return FALSE;
3207}
3208#ifndef MDEBUG
3209static BOOLEAN jjpHead(leftv res, leftv v)
3210{
3211  res->data = (char *)pHead((poly)v->Data());
3212  return FALSE;
3213}
3214#endif
3215static BOOLEAN jjidHead(leftv res, leftv v)
3216{
3217  res->data = (char *)idHead((ideal)v->Data());
3218  return FALSE;
3219}
3220static BOOLEAN jjidMaxIdeal(leftv res, leftv v)
3221{
3222  res->data = (char *)idMaxIdeal((int)v->Data());
3223  return FALSE;
3224}
3225static BOOLEAN jjidMinBase(leftv res, leftv v)
3226{
3227  res->data = (char *)idMinBase((ideal)v->Data());
3228  return FALSE;
3229}
3230static BOOLEAN jjsyMinBase(leftv res, leftv v)
3231{
3232  res->data = (char *)syMinBase((ideal)v->Data());
3233  return FALSE;
3234}
3235static BOOLEAN jjpMaxComp(leftv res, leftv v)
3236{
3237  res->data = (char *)pMaxComp((poly)v->Data());
3238  return FALSE;
3239}
3240static BOOLEAN jjmpTrace(leftv res, leftv v)
3241{
3242  res->data = (char *)mpTrace((matrix)v->Data());
3243  return FALSE;
3244}
3245static BOOLEAN jjmpTransp(leftv res, leftv v)
3246{
3247  res->data = (char *)mpTransp((matrix)v->Data());
3248  return FALSE;
3249}
3250static BOOLEAN jjrOrdStr(leftv res, leftv v)
3251{
3252  res->data = rOrdStr((ring)v->Data());
3253  return FALSE;
3254}
3255static BOOLEAN jjrVarStr(leftv res, leftv v)
3256{
3257  res->data = rVarStr((ring)v->Data());
3258  return FALSE;
3259}
3260static BOOLEAN jjrParStr(leftv res, leftv v)
3261{
3262  res->data = rParStr((ring)v->Data());
3263  return FALSE;
3264}
3265static BOOLEAN jjidMinEmbedding(leftv res, leftv v)
3266{
3267  res->data = (char *)idMinEmbedding((ideal)v->Data());
3268  return FALSE;
3269}
3270static BOOLEAN jjBETTI_R(leftv res, leftv v)
3271{
3272  res->data=(char *)syBettiOfComputation((syStrategy)v->Data());
3273  return FALSE;
3274}
3275static BOOLEAN jjCOUNT_R(leftv res, leftv v)
3276{
3277  res->data=(char *)syLength((syStrategy)v->Data());
3278  return FALSE;
3279}
3280static BOOLEAN jjDIM_R(leftv res, leftv v)
3281{
3282  res->data = (char *)syDim((syStrategy)v->Data());
3283  return FALSE;
3284}
3285static BOOLEAN jjMINRES_R(leftv res, leftv v)
3286{
3287  res->data=(char *)syMinimize((syStrategy)v->Data());
3288  return FALSE;
3289}
3290static BOOLEAN jjidTransp(leftv res, leftv v)
3291{
3292  res->data = (char *)idTransp((ideal)v->Data());
3293  return FALSE;
3294}
3295#else
3296#define XS(A)          -((short)A)
3297#define jjstrlen       (proc1)strlen
3298#define jjpLength      (proc1)pLength
3299#define jjidElem       (proc1)idElem
3300#define jjmpDetBareiss (proc1)mpDetBareiss
3301#define jjidFreeModule (proc1)idFreeModule
3302#define jjidVec2Ideal  (proc1)idVec2Ideal
3303#define jjrCharStr     (proc1)rCharStr
3304#ifndef MDEBUG
3305#define jjpHead        (proc1)pHead
3306#endif
3307#define jjidHead       (proc1)idHead
3308#define jjidMaxIdeal   (proc1)idMaxIdeal
3309#define jjidMinBase    (proc1)idMinBase
3310#define jjsyMinBase    (proc1)syMinBase
3311#define jjpMaxComp     (proc1)pMaxComp
3312#define jjmpTrace      (proc1)mpTrace
3313#define jjmpTransp     (proc1)mpTransp
3314#define jjrOrdStr      (proc1)rOrdStr
3315#define jjrVarStr      (proc1)rVarStr
3316#define jjrParStr      (proc1)rParStr
3317#define jjidMinEmbedding (proc1)idMinEmbedding
3318#define jjBETTI_R        (proc1)syBettiOfComputation
3319#define jjCOUNT_R        (proc1)syLength
3320#define jjDIM_R        (proc1)syDim
3321#define jjMINRES_R     (proc1)syMinimize
3322#define jjidTransp     (proc1)idTransp
3323#endif
3324#endif
3325BOOLEAN jjnInt(leftv res, leftv u)
3326{
3327  number n=(number)u->CopyD();
3328  res->data=(char *)nInt(n);
3329  nDelete(&n);
3330  return FALSE;
3331}
3332#define s short
3333struct sValCmd1 dArith1[]=
3334{
3335// operations:
3336// proc         cmd               res             arg
3337// {jjPLUSPLUS,   PLUSPLUS,        INT_CMD,        INT_CMD }
3338//,{jjMINUSMINUS, MINUSMINUS,      INT_CMD,        INT_CMD }
3339 {jjPLUSPLUS,   PLUSPLUS,        NONE,           IDHDL }
3340,{jjPLUSPLUS,   MINUSMINUS,      NONE,           IDHDL }
3341,{jjUMINUS_I,   '-',             INT_CMD,        INT_CMD }
3342,{jjUMINUS_N,   '-',             NUMBER_CMD,     NUMBER_CMD }
3343,{jjUMINUS_P,   '-',             POLY_CMD,       POLY_CMD }
3344,{jjUMINUS_P,   '-',             VECTOR_CMD,     VECTOR_CMD }
3345,{jjUMINUS_MA,  '-',             MATRIX_CMD,     MATRIX_CMD }
3346,{jjUMINUS_IV,  '-',             INTVEC_CMD,     INTVEC_CMD }
3347,{jjUMINUS_IV,  '-',             INTMAT_CMD,     INTMAT_CMD }
3348,{jjPROC1,      '(',             ANY_TYPE/*set by p*/,PROC_CMD }
3349,{jjKLAMMER_LIB,'(',             NONE,           STRING_CMD }
3350// and the procedures with 1 argument:
3351,{atATTRIB1,    ATTRIB_CMD,      NONE,           DEF_CMD }
3352,{jjBAREISS_IM, BAREISS_CMD,     INTMAT_CMD,     INTMAT_CMD }
3353,{jjBAREISS,    BAREISS_CMD,     LIST_CMD,       MODUL_CMD }
3354,{jjBETTI,      BETTI_CMD,       INTMAT_CMD,     LIST_CMD }
3355,{jjBETTI_R,    BETTI_CMD,       XS(INTMAT_CMD), RESOLUTION_CMD }
3356,{jjCHAR,       CHARACTERISTIC_CMD, INT_CMD,     RING_CMD }
3357,{jjCHAR,       CHARACTERISTIC_CMD, INT_CMD,     QRING_CMD }
3358#ifdef HAVE_FACTORY
3359,{jjCHARSERIES, CHAR_SERIES_CMD, MATRIX_CMD,     IDEAL_CMD }
3360#else
3361,{jjWRONG,      CHAR_SERIES_CMD, MATRIX_CMD,     IDEAL_CMD }
3362#endif
3363,{jjrCharStr,   CHARSTR_CMD,     XS(STRING_CMD), RING_CMD }
3364,{jjrCharStr,   CHARSTR_CMD,     XS(STRING_CMD), QRING_CMD }
3365,{jjOpenClose,  CLOSE_CMD,       NONE,           LINK_CMD }
3366,{jjWRONG,      COLS_CMD,        0,              VECTOR_CMD }
3367,{jjCOLS,       COLS_CMD,        INT_CMD,        MATRIX_CMD }
3368,{jjCOLS,       COLS_CMD,        INT_CMD,        IDEAL_CMD }
3369,{jjCOLS,       COLS_CMD,        INT_CMD,        MODUL_CMD }
3370,{jjCOLS_IV,    COLS_CMD,        INT_CMD,        INTMAT_CMD }
3371,{jjWRONG,      COLS_CMD,        0,              INTVEC_CMD }
3372,{jjCONTENT,    CONTENT_CMD,     POLY_CMD,       POLY_CMD }
3373,{jjCONTENT,    CONTENT_CMD,     VECTOR_CMD,     VECTOR_CMD }
3374,{jjCOUNT_N,    COUNT_CMD,       INT_CMD,        NUMBER_CMD }
3375,{jjCOUNT_R,    COUNT_CMD,       XS(INT_CMD),    RESOLUTION_CMD }
3376,{jjstrlen,     COUNT_CMD,       XS(INT_CMD),    STRING_CMD }
3377,{jjpLength,    COUNT_CMD,       XS(INT_CMD),    POLY_CMD }
3378,{jjpLength,    COUNT_CMD,       XS(INT_CMD),    VECTOR_CMD }
3379,{jjidElem,     COUNT_CMD,       XS(INT_CMD),    IDEAL_CMD }
3380,{jjidElem,     COUNT_CMD,       XS(INT_CMD),    MODUL_CMD }
3381,{jjCOUNT_M,    COUNT_CMD,       INT_CMD,        MATRIX_CMD }
3382,{jjCOUNT_IV,   COUNT_CMD,       INT_CMD,        INTVEC_CMD }
3383,{jjCOUNT_IV,   COUNT_CMD,       INT_CMD,        INTMAT_CMD }
3384,{jjCOUNT_L,    COUNT_CMD,       INT_CMD,        LIST_CMD }
3385,{jjNULL,       DEF_CMD,         DEF_CMD,        INT_CMD }
3386,{jjWRONG,      DEF_CMD,         0,              ANY_TYPE }
3387,{jjDEG,        DEG_CMD,         INT_CMD,        POLY_CMD }
3388,{jjDEG,        DEG_CMD,         INT_CMD,        VECTOR_CMD }
3389,{jjDEGREE,     DEGREE_CMD,      NONE,           IDEAL_CMD }
3390,{jjDEGREE,     DEGREE_CMD,      NONE,           MODUL_CMD }
3391,{jjDEFINED,    DEFINED_CMD,     INT_CMD,        DEF_CMD }
3392#ifdef HAVE_FACTORY
3393,{jjDET_I,      DET_CMD,         INT_CMD,        INTMAT_CMD }
3394,{jjDET,        DET_CMD,         POLY_CMD,       MATRIX_CMD }
3395#else
3396,{jjWRONG,      DET_CMD,         INT_CMD,        INTMAT_CMD }
3397,{jjmpDetBareiss,DET_CMD,        XS(POLY_CMD),   MATRIX_CMD }
3398#endif
3399,{jjDIM,        DIM_CMD,         INT_CMD,        IDEAL_CMD }
3400,{jjDIM,        DIM_CMD,         INT_CMD,        MODUL_CMD }
3401,{jjDIM_R,      DIM_CMD,         XS(INT_CMD),    RESOLUTION_CMD }
3402,{jjDUMP,       DUMP_CMD,        NONE,           LINK_CMD }
3403,{jjE,          E_CMD,           VECTOR_CMD,     INT_CMD }
3404#ifdef HAVE_FACTORY
3405,{jjFAC_P,      FAC_CMD,         LIST_CMD,       POLY_CMD }
3406#else
3407,{jjWRONG,      FAC_CMD,         LIST_CMD,       POLY_CMD }
3408#endif
3409#ifdef HAVE_FGLM
3410,{findUniProc,  FINDUNI_CMD,     IDEAL_CMD,      IDEAL_CMD}
3411#else
3412,{jjWRONG,      FINDUNI_CMD,     IDEAL_CMD,      IDEAL_CMD}
3413#endif
3414,{jjidFreeModule,FREEMODULE_CMD, XS(MODUL_CMD),  INT_CMD }
3415#ifdef HAVE_FACTORY
3416,{jjFACSTD,     FACSTD_CMD,      LIST_CMD,       IDEAL_CMD }
3417#else
3418,{jjWRONG,      FACSTD_CMD,      LIST_CMD,       IDEAL_CMD }
3419#endif
3420,{jjGETDUMP,    GETDUMP_CMD,     NONE,           LINK_CMD }
3421,{jjHIGHCORNER, HIGHCORNER_CMD,  POLY_CMD,       IDEAL_CMD }
3422,{jjHIGHCORNER_M, HIGHCORNER_CMD,VECTOR_CMD,     MODUL_CMD }
3423,{jjHILBERT,    HILBERT_CMD,     NONE,           IDEAL_CMD }
3424,{jjHILBERT,    HILBERT_CMD,     NONE,           MODUL_CMD }
3425,{jjHILBERT_IV, HILBERT_CMD,     INTVEC_CMD,     INTVEC_CMD }
3426,{jjHOMOG1,     HOMOG_CMD,       INT_CMD,        IDEAL_CMD }
3427,{jjHOMOG1,     HOMOG_CMD,       INT_CMD,        MODUL_CMD }
3428,{jjidVec2Ideal,IDEAL_CMD,       XS(IDEAL_CMD),  VECTOR_CMD }
3429,{jjIDEAL_Ma,   IDEAL_CMD,       IDEAL_CMD,      MATRIX_CMD }
3430,{jjIDEAL_R,    IDEAL_CMD,       IDEAL_CMD,      QRING_CMD }
3431,{jjIDEAL_R,    IDEAL_CMD,       IDEAL_CMD,      RING_CMD }
3432,{jjIDEAL_Map,  IDEAL_CMD,       IDEAL_CMD,      MAP_CMD }
3433,{jjDUMMY,      IDEAL_CMD,       IDEAL_CMD,      IDEAL_CMD }
3434,{jjINDEPSET,   INDEPSET_CMD,    INTVEC_CMD,     IDEAL_CMD }
3435,{jjDUMMY,      INT_CMD,         INT_CMD,        INT_CMD }
3436,{jjnInt,       INT_CMD,         INT_CMD,        NUMBER_CMD }
3437,{jjP2I,        INT_CMD,         INT_CMD,        POLY_CMD }
3438,{jjINTERRED,   INTERRED_CMD,    IDEAL_CMD,      IDEAL_CMD }
3439,{jjINTERRED,   INTERRED_CMD,    MODUL_CMD,      MODUL_CMD }
3440,{jjDUMMY,      INTMAT_CMD,      INTMAT_CMD,     INTMAT_CMD }
3441,{jjIm2Iv,      INTVEC_CMD,      INTVEC_CMD,     INTMAT_CMD }
3442,{jjDUMMY,      INTVEC_CMD,      INTVEC_CMD,     INTVEC_CMD }
3443,{jjIS_RINGVAR_P, IS_RINGVAR,    INT_CMD,        POLY_CMD }
3444,{jjIS_RINGVAR_S, IS_RINGVAR,    INT_CMD,        STRING_CMD }
3445,{jjIS_RINGVAR0,IS_RINGVAR,      INT_CMD,        ANY_TYPE }
3446,{jjJACOB_P,    JACOB_CMD,       IDEAL_CMD,      POLY_CMD }
3447,{mpJacobi,     JACOB_CMD,       MATRIX_CMD,     IDEAL_CMD }
3448,{jjKBASE,      KBASE_CMD,       IDEAL_CMD,      IDEAL_CMD }
3449,{jjKBASE,      KBASE_CMD,       MODUL_CMD,      MODUL_CMD }
3450,{atKILLATTR1,  KILLATTR_CMD,    NONE,           IDHDL }
3451#ifdef MDEBUG
3452,{jjpHead,      LEAD_CMD,        POLY_CMD,       POLY_CMD }
3453#else
3454,{jjpHead,      LEAD_CMD,        XS(POLY_CMD),   POLY_CMD }
3455#endif
3456,{jjidHead,     LEAD_CMD,        XS(IDEAL_CMD),  IDEAL_CMD }
3457#ifdef MDEBUG
3458,{jjpHead,      LEAD_CMD,        VECTOR_CMD,     VECTOR_CMD }
3459#else
3460,{jjpHead,      LEAD_CMD,        XS(VECTOR_CMD), VECTOR_CMD }
3461#endif
3462,{jjidHead,     LEAD_CMD,        XS(MODUL_CMD),  MODUL_CMD }
3463,{jjLEADCOEF,   LEADCOEF_CMD,    NUMBER_CMD,     POLY_CMD }
3464,{jjLEADCOEF,   LEADCOEF_CMD,    NUMBER_CMD,     VECTOR_CMD }
3465,{jjLEADEXP,    LEADEXP_CMD,     INTVEC_CMD,     POLY_CMD }
3466,{jjLEADEXP,    LEADEXP_CMD,     INTVEC_CMD,     VECTOR_CMD }
3467,{jjLIB,        LIB_CMD,         NONE,           STRING_CMD }
3468,{jjCALL1MANY,  LIST_CMD,        LIST_CMD,       DEF_CMD }
3469,{jjWRONG,      MAP_CMD,         0,              ANY_TYPE}
3470,{jjDUMMY,      MATRIX_CMD,      MATRIX_CMD,     MATRIX_CMD }
3471,{jjidMaxIdeal, MAXID_CMD,       XS(IDEAL_CMD),  INT_CMD }
3472,{jjMEMORY,     MEMORY_CMD,      INT_CMD,        INT_CMD }
3473,{jjidMinBase,  MINBASE_CMD,     XS(IDEAL_CMD),  IDEAL_CMD }
3474,{jjidMinBase,  MINBASE_CMD,     XS(MODUL_CMD),  MODUL_CMD }
3475,{jjMINRES,     MINRES_CMD,      LIST_CMD,       LIST_CMD }
3476,{jjMINRES_R,   MINRES_CMD,      XS(RESOLUTION_CMD), RESOLUTION_CMD }
3477,{jjDUMMY,      MODUL_CMD,       MODUL_CMD,      MODUL_CMD }
3478,{jjMONITOR1,   MONITOR_CMD,     NONE,           STRING_CMD }
3479,{jjMULT,       MULTIPLICITY_CMD,  INT_CMD,      IDEAL_CMD }
3480,{jjMULT,       MULTIPLICITY_CMD,  INT_CMD,      MODUL_CMD }
3481,{jjMSTD,       MSTD_CMD,        LIST_CMD,       IDEAL_CMD }
3482,{jjNAMEOF,     NAMEOF_CMD,      STRING_CMD,     ANY_TYPE }
3483,{jjNAMES,      NAMES_CMD,       LIST_CMD,       PACKAGE_CMD }
3484,{jjNAMES,      NAMES_CMD,       LIST_CMD,       RING_CMD }
3485,{jjNAMES,      NAMES_CMD,       LIST_CMD,       QRING_CMD }
3486,{jjDUMMY,      NUMBER_CMD,      NUMBER_CMD,     NUMBER_CMD }
3487,{jjP2N,        NUMBER_CMD,      NUMBER_CMD,     POLY_CMD }
3488,{jjRPAR,       NPARS_CMD,       INT_CMD,        RING_CMD }
3489,{jjRPAR,       NPARS_CMD,       INT_CMD,        QRING_CMD }
3490,{jjNVARS,      NVARS_CMD,       INT_CMD,        RING_CMD }
3491,{jjNVARS,      NVARS_CMD,       INT_CMD,        QRING_CMD }
3492,{jjOpenClose,  OPEN_CMD,        NONE,           LINK_CMD }
3493,{jjORD,        ORD_CMD,         INT_CMD,        POLY_CMD }
3494,{jjORD,        ORD_CMD,         INT_CMD,        VECTOR_CMD }
3495,{jjrOrdStr,    ORDSTR_CMD,      XS(STRING_CMD), RING_CMD }
3496,{jjrOrdStr,    ORDSTR_CMD,      XS(STRING_CMD), QRING_CMD }
3497,{jjPAR1,       PAR_CMD,         NUMBER_CMD,     INT_CMD }
3498,{jjPARDEG,     PARDEG_CMD,      INT_CMD,        NUMBER_CMD }
3499,{jjPARSTR1,    PARSTR_CMD,      STRING_CMD,     INT_CMD }
3500,{jjrParStr,    PARSTR_CMD,      XS(STRING_CMD), RING_CMD }
3501,{jjrParStr,    PARSTR_CMD,      XS(STRING_CMD), QRING_CMD }
3502,{jjDUMMY,      POLY_CMD,        POLY_CMD,       POLY_CMD }
3503,{jjPRIME,      PRIME_CMD,       INT_CMD,        INT_CMD }
3504,{jjPRINT_GEN,  PRINT_CMD,       NONE,           INT_CMD }
3505,{jjPRINT_INTVEC,PRINT_CMD,      NONE,           INTVEC_CMD }
3506,{jjPRINT_INTMAT,PRINT_CMD,      NONE,           INTMAT_CMD }
3507,{jjPRINT_GEN,  PRINT_CMD,       NONE,           STRING_CMD }
3508,{jjPRINT_GEN,  PRINT_CMD,       NONE,           POLY_CMD }
3509,{jjPRINT_V,    PRINT_CMD,       NONE,           VECTOR_CMD }
3510,{jjPRINT_GEN,  PRINT_CMD,       NONE,           IDEAL_CMD }
3511,{jjPRINT_MA,   PRINT_CMD,       NONE,           MATRIX_CMD }
3512,{jjPRINT_LIST, PRINT_CMD,       NONE,           LIST_CMD }
3513,{jjPRINT_GEN,  PRINT_CMD,       NONE,           RESOLUTION_CMD }
3514,{jjidMinEmbedding, PRUNE_CMD,   XS(MODUL_CMD),  MODUL_CMD }
3515,{kQHWeight,    QHWEIGHT_CMD,    INTVEC_CMD,     IDEAL_CMD }
3516,{kQHWeight,    QHWEIGHT_CMD,    INTVEC_CMD,     MODUL_CMD }
3517,{jjWRONG,      QRING_CMD,       0,              ANY_TYPE}
3518,{jjREAD,       READ_CMD,        STRING_CMD,     LINK_CMD }
3519,{jjREGULARITY, REGULARITY_CMD,  INT_CMD,        LIST_CMD }
3520,{jjRESERVEDNAME,RESERVEDNAME_CMD, INT_CMD,      STRING_CMD }
3521,{jjL2R,        RESOLUTION_CMD,  RESOLUTION_CMD, LIST_CMD }
3522,{jjDUMMY,      RESOLUTION_CMD,  RESOLUTION_CMD, RESOLUTION_CMD }
3523,{jjWRONG,      ROWS_CMD,        0,              POLY_CMD }
3524,{jjpMaxComp,   ROWS_CMD,        XS(INT_CMD),    VECTOR_CMD }
3525,{jjROWS,       ROWS_CMD,        INT_CMD,        MODUL_CMD }
3526,{jjROWS,       ROWS_CMD,        INT_CMD,        MATRIX_CMD }
3527,{jjROWS_IV,    ROWS_CMD,        INT_CMD,        INTMAT_CMD }
3528,{jjCOUNT_IV,   ROWS_CMD,        INT_CMD,        INTVEC_CMD }
3529,{jjSort_Id,    SORTVEC_CMD,     INTVEC_CMD,     IDEAL_CMD }
3530,{jjSort_Id,    SORTVEC_CMD,     INTVEC_CMD,     MODUL_CMD }
3531,{jjSTD,        STD_CMD,         IDEAL_CMD,      IDEAL_CMD }
3532,{jjSTD,        STD_CMD,         MODUL_CMD,      MODUL_CMD }
3533,{jjDUMMY,      STRING_CMD,      STRING_CMD,     STRING_CMD }
3534,{jjSTRING,     STRING_CMD,      STRING_CMD,     INT_CMD }
3535,{jjSTRING,     STRING_CMD,      STRING_CMD,     POLY_CMD }
3536,{jjSTRING,     STRING_CMD,      STRING_CMD,     VECTOR_CMD }
3537,{jjSTRING,     STRING_CMD,      STRING_CMD,     NUMBER_CMD }
3538,{jjSTRING,     STRING_CMD,      STRING_CMD,     MATRIX_CMD }
3539,{jjSTRING,     STRING_CMD,      STRING_CMD,     MODUL_CMD }
3540,{jjSTRING,     STRING_CMD,      STRING_CMD,     IDEAL_CMD }
3541,{jjSTRING,     STRING_CMD,      STRING_CMD,     MAP_CMD }
3542,{jjSTRING,     STRING_CMD,      STRING_CMD,     INTVEC_CMD }
3543,{jjSTRING,     STRING_CMD,      STRING_CMD,     INTMAT_CMD }
3544,{jjSTRING_PROC,STRING_CMD,      STRING_CMD,     PROC_CMD }
3545,{jjSTRING,     STRING_CMD,      STRING_CMD,     RING_CMD }
3546,{jjSTRING,     STRING_CMD,      STRING_CMD,     QRING_CMD }
3547,{jjSYSTEM,     SYSTEM_CMD,      NONE,           STRING_CMD }
3548,{jjSYZYGY,     SYZYGY_CMD,      MODUL_CMD,      IDEAL_CMD }
3549,{jjSYZYGY,     SYZYGY_CMD,      MODUL_CMD,      MODUL_CMD }
3550,{jjWRONG,      TRACE_CMD,       0,              INTVEC_CMD }
3551,{jjWRONG,      TRACE_CMD,       0,              IDEAL_CMD }
3552,{jjTRACE_IV,   TRACE_CMD,       INT_CMD,        INTMAT_CMD }
3553,{jjmpTrace,    TRACE_CMD,       XS(POLY_CMD),   MATRIX_CMD }
3554,{jjTRANSP_IV,  TRANSPOSE_CMD,   INTMAT_CMD,     INTVEC_CMD }
3555,{jjTRANSP_IV,  TRANSPOSE_CMD,   INTMAT_CMD,     INTMAT_CMD }
3556,{jjmpTransp,   TRANSPOSE_CMD,   XS(MATRIX_CMD), MATRIX_CMD }
3557,{jjidTransp,   TRANSPOSE_CMD,   XS(MODUL_CMD),  MODUL_CMD }
3558,{jjTYPEOF,     TYPEOF_CMD,      STRING_CMD,     ANY_TYPE }
3559,{jjDUMMY,      VECTOR_CMD,      VECTOR_CMD,     VECTOR_CMD }
3560,{jjVDIM,       VDIM_CMD,        INT_CMD,        IDEAL_CMD }
3561,{jjVDIM,       VDIM_CMD,        INT_CMD,        MODUL_CMD }
3562,{jjVAR1,       VAR_CMD,         POLY_CMD,       INT_CMD }
3563,{jjVARSTR1,    VARSTR_CMD,      STRING_CMD,     INT_CMD }
3564,{jjrVarStr,    VARSTR_CMD,      XS(STRING_CMD), RING_CMD }
3565,{jjrVarStr,    VARSTR_CMD,      XS(STRING_CMD), QRING_CMD }
3566,{kWeight,      WEIGHT_CMD,      INTVEC_CMD,     IDEAL_CMD }
3567,{kWeight,      WEIGHT_CMD,      INTVEC_CMD,     MODUL_CMD }
3568,{jjLOAD1,      LOAD_CMD,        NONE,           STRING_CMD }
3569,{NULL,         0,               0,              0}
3570};
3571#undef s
3572/*=================== operations with 3 args.: static proc =================*/
3573static BOOLEAN jjBRACK_S(leftv res, leftv u, leftv v,leftv w)
3574{
3575  char *s= (char *)u->Data();
3576  int   r = (int)v->Data();
3577  int   c = (int)w->Data();
3578  int l = strlen(s);
3579
3580  if ( (r<1) || (r>l) || (c<0) )
3581  {
3582    Werror("wrong range[%d,%d] in string %s",r,c,u->Fullname());
3583    return TRUE;
3584  }
3585  res->data = (char *)AllocL(c+1);
3586  sprintf((char *)res->data,"%-*.*s",c,c,s+r-1);
3587  return FALSE;
3588}
3589static BOOLEAN jjBRACK_Im(leftv res, leftv u, leftv v,leftv w)
3590{
3591  intvec *iv = (intvec *)u->Data();
3592  int   r = (int)v->Data();
3593  int   c = (int)w->Data();
3594  if ((r<1)||(r>iv->rows())||(c<1)||(c>iv->cols()))
3595  {
3596    Werror("wrong range[%d,%d] in intmat %s(%d x %d)",
3597           r,c,u->Fullname(),iv->rows(),iv->cols());
3598    return TRUE;
3599  }
3600  res->data=u->data;
3601  u->data=NULL;
3602  res->rtyp=u->rtyp;
3603  u->rtyp=0;
3604  res->name=u->name;
3605  u->name=NULL;
3606  Subexpr e=jjMakeSub(v);
3607          e->next=jjMakeSub(w);
3608  if (u->e==NULL) res->e=e;
3609  else
3610  {
3611    Subexpr h=u->e;
3612    while (h->next!=NULL) h=h->next;
3613    h->next=e;
3614    res->e=u->e;
3615    u->e=NULL;
3616  }
3617  return FALSE;
3618}
3619static BOOLEAN jjBRACK_Ma(leftv res, leftv u, leftv v,leftv w)
3620{
3621  matrix m= (matrix)u->Data();
3622  int   r = (int)v->Data();
3623  int   c = (int)w->Data();
3624  //Print("gen. elem %d, %d\n",r,c);
3625  if ((r<1)||(r>MATROWS(m))||(c<1)||(c>MATCOLS(m)))
3626  {
3627    Werror("wrong range[%d,%d] in matrix %s(%d x %d)",r,c,u->Fullname(),
3628      MATROWS(m),MATCOLS(m));
3629    return TRUE;
3630  }
3631  res->data=u->data;
3632  u->data=NULL;
3633  res->rtyp=u->rtyp;
3634  u->rtyp=0;
3635  res->name=u->name;
3636  u->name=NULL;
3637  Subexpr e=jjMakeSub(v);
3638          e->next=jjMakeSub(w);
3639  if (u->e==NULL)
3640    res->e=e;
3641  else
3642  {
3643    Subexpr h=u->e;
3644    while (h->next!=NULL) h=h->next;
3645    h->next=e;
3646    res->e=u->e;
3647    u->e=NULL;
3648  }
3649  return FALSE;
3650}
3651static BOOLEAN jjBRACK_Ma_I_IV(leftv res, leftv u, leftv v,leftv w)
3652{
3653  sleftv t;
3654  sleftv ut;
3655  leftv p=NULL;
3656  intvec *iv=(intvec *)w->Data();
3657  int l;
3658  BOOLEAN nok;
3659
3660  if ((u->rtyp!=IDHDL)||(u->e!=NULL))
3661  {
3662    WerrorS("cannot build expression lists from unnamed objects");
3663    return TRUE;
3664  }
3665  memcpy(&ut,u,sizeof(ut));
3666  memset(&t,0,sizeof(t));
3667  t.rtyp=INT_CMD;
3668  for (l=0;l< iv->length(); l++)
3669  {
3670    t.data=(char *)(*iv)[l];
3671    if (p==NULL)
3672    {
3673      p=res;
3674    }
3675    else
3676    {
3677      p->next=(leftv)Alloc0(sizeof(sleftv));
3678      p=p->next;
3679    }
3680    memcpy(u,&ut,sizeof(ut));
3681    if (u->Typ() == MATRIX_CMD)
3682      nok=jjBRACK_Ma(p,u,v,&t);
3683    else /* INTMAT_CMD */
3684      nok=jjBRACK_Im(p,u,v,&t);
3685    if (nok)
3686    {
3687      while (res->next!=NULL)
3688      {
3689        p=res->next->next;
3690        Free((ADDRESS)res->next,sizeof(sleftv));
3691        // res->e aufraeumen !!!!
3692        res->next=p;
3693      }
3694      return TRUE;
3695    }
3696  }
3697  return FALSE;
3698}
3699static BOOLEAN jjBRACK_Ma_IV_I(leftv res, leftv u, leftv v,leftv w)
3700{
3701  sleftv t;
3702  sleftv ut;
3703  leftv p=NULL;
3704  intvec *iv=(intvec *)v->Data();
3705  int l;
3706  BOOLEAN nok;
3707
3708  if ((u->rtyp!=IDHDL)||(u->e!=NULL))
3709  {
3710    WerrorS("cannot build expression lists from unnamed objects");
3711    return TRUE;
3712  }
3713  memcpy(&ut,u,sizeof(ut));
3714  memset(&t,0,sizeof(t));
3715  t.rtyp=INT_CMD;
3716  for (l=0;l< iv->length(); l++)
3717  {
3718    t.data=(char *)((*iv)[l]);
3719    if (p==NULL)
3720    {
3721      p=res;
3722    }
3723    else
3724    {
3725      p->next=(leftv)Alloc0(sizeof(sleftv));
3726      p=p->next;
3727    }
3728    memcpy(u,&ut,sizeof(ut));
3729    if (u->Typ() == MATRIX_CMD)
3730      nok=jjBRACK_Ma(p,u,&t,w);
3731    else /* INTMAT_CMD */
3732      nok=jjBRACK_Im(p,u,&t,w);
3733    if (nok)
3734    {
3735      while (res->next!=NULL)
3736      {
3737        p=res->next->next;
3738        Free((ADDRESS)res->next,sizeof(sleftv));
3739        // res->e aufraeumen !!
3740        res->next=p;
3741      }
3742      return TRUE;
3743    }
3744  }
3745  return FALSE;
3746}
3747static BOOLEAN jjBRACK_Ma_IV_IV(leftv res, leftv u, leftv v,leftv w)
3748{
3749  sleftv t1,t2,ut;
3750  leftv p=NULL;
3751  intvec *vv=(intvec *)v->Data();
3752  intvec *wv=(intvec *)w->Data();
3753  int vl;
3754  int wl;
3755  BOOLEAN nok;
3756
3757  if ((u->rtyp!=IDHDL)||(u->e!=NULL))
3758  {
3759    WerrorS("cannot build expression lists from unnamed objects");
3760    return TRUE;
3761  }
3762  memcpy(&ut,u,sizeof(ut));
3763  memset(&t1,0,sizeof(sleftv));
3764  memset(&t2,0,sizeof(sleftv));
3765  t1.rtyp=INT_CMD;
3766  t2.rtyp=INT_CMD;
3767  for (vl=0;vl< vv->length(); vl++)
3768  {
3769    t1.data=(char *)((*vv)[vl]);
3770    for (wl=0;wl< wv->length(); wl++)
3771    {
3772      t2.data=(char *)((*wv)[wl]);
3773      if (p==NULL)
3774      {
3775        p=res;
3776      }
3777      else
3778      {
3779        p->next=(leftv)Alloc0(sizeof(sleftv));
3780        p=p->next;
3781      }
3782      memcpy(u,&ut,sizeof(ut));
3783      if (u->Typ() == MATRIX_CMD)
3784        nok=jjBRACK_Ma(p,u,&t1,&t2);
3785      else /* INTMAT_CMD */
3786        nok=jjBRACK_Im(p,u,&t1,&t2);
3787      if (nok)
3788      {
3789        res->CleanUp();
3790        return TRUE;
3791      }
3792    }
3793  }
3794  return FALSE;
3795}
3796static BOOLEAN jjPROC3(leftv res, leftv u, leftv v, leftv w)
3797{
3798  v->next=(leftv)Alloc(sizeof(sleftv));
3799  memcpy(v->next,w,sizeof(sleftv));
3800  BOOLEAN r=iiExprArith2(res,u,'(',v);
3801  v->rtyp=0; v->data=NULL;
3802  w->rtyp=0; w->data=NULL;
3803  return r;
3804}
3805static BOOLEAN jjCALL3MANY(leftv res, leftv u, leftv v, leftv w)
3806{
3807  u->next=(leftv)Alloc(sizeof(sleftv));
3808  memcpy(u->next,v,sizeof(sleftv));
3809  u->next->next=(leftv)Alloc(sizeof(sleftv));
3810  memcpy(u->next->next,w,sizeof(sleftv));
3811  BOOLEAN r=iiExprArithM(res,u,iiOp);
3812  v->rtyp=0; v->data=NULL; // iiExprArithM did the CleanUp
3813  w->rtyp=0; w->data=NULL;
3814  return r;
3815}
3816static BOOLEAN jjCOEFFS3_Id(leftv res, leftv u, leftv v, leftv w)
3817{
3818  if ((w->rtyp!=IDHDL)||(w->e!=NULL))
3819  {
3820    WerrorS("3rd argument must be a name of a matrix");
3821    return TRUE;
3822  }
3823  ideal i=(ideal)u->Data();
3824  int rank=(int)i->rank;
3825  BOOLEAN r=jjCOEFFS_Id(res,u,v);
3826  if (r) return TRUE;
3827  mpMonomials((matrix)res->data, rank, pVar((poly)v->Data()),(matrix)w->Data());
3828  return FALSE;
3829}
3830static BOOLEAN jjCOEFFS3_KB(leftv res, leftv u, leftv v, leftv w)
3831{
3832  res->data=(void*)idCoeffOfKBase((ideal)(u->Data()),
3833           (ideal)(v->Data()),(poly)(w->Data()));
3834  return FALSE;
3835}
3836static BOOLEAN jjCOEFFS3_P(leftv res, leftv u, leftv v, leftv w)
3837{
3838  if ((w->rtyp!=IDHDL)||(w->e!=NULL))
3839  {
3840    WerrorS("3rd argument must be a name of a matrix");
3841    return TRUE;
3842  }
3843  poly p=(poly)u->CopyD();
3844  ideal i=idInit(1,1);
3845  i->m[0]=p;
3846  sleftv t;
3847  memset(&t,0,sizeof(t));
3848  t.data=(char *)i;
3849  t.rtyp=IDEAL_CMD;
3850  int rank=1;
3851  if (u->Typ()==VECTOR_CMD)
3852  {
3853    i->rank=rank=pMaxComp(p);
3854    t.rtyp=MODUL_CMD;
3855  }
3856  BOOLEAN r=jjCOEFFS_Id(res,&t,v);
3857  t.CleanUp();
3858  if (r) return TRUE;
3859  mpMonomials((matrix)res->data, rank, pVar((poly)v->Data()),(matrix)w->Data());
3860  return FALSE;
3861}
3862static BOOLEAN jjELIMIN_HILB(leftv res, leftv u, leftv v, leftv w)
3863{
3864  res->data=(char *)idElimination((ideal)u->Data(),(poly)v->Data(),
3865    (intvec *)w->Data());
3866  setFlag(res,FLAG_STD);
3867  return FALSE;
3868}
3869static BOOLEAN jjFIND3(leftv res, leftv u, leftv v, leftv w)
3870{
3871  /*4
3872  * look for the substring what in the string where
3873  * starting at position n
3874  * return the position of the first char of what in where
3875  * or 0
3876  */
3877  int n=(int)w->Data();
3878  char *where=(char *)u->Data();
3879  char *what=(char *)v->Data();
3880  char *found;
3881  if ((1>n)||(n>(int)strlen(where)))
3882  {
3883    Werror("start position %d out of range",n);
3884    return TRUE;
3885  }
3886  found = strchr(where+n-1,*what);
3887  if (*(what+1)!='\0')
3888  {
3889    while((found !=NULL) && (strncmp(found+1,what+1,strlen(what+1))!=0))
3890    {
3891      found=strchr(found+1,*what);
3892    }
3893  }
3894  if (found != NULL)
3895  {
3896    res->data=(char *)((found-where)+1);
3897  }
3898  return FALSE;
3899}
3900static BOOLEAN jjJET_P_IV(leftv res, leftv u, leftv v, leftv w)
3901{
3902  short *iw=iv2array((intvec *)w->Data());
3903  res->data = (char *)pJetW((poly)u->Data(),(int)v->Data(),iw);
3904  Free((ADDRESS)iw,(pVariables+1)*sizeof(short));
3905  return FALSE;
3906}
3907static BOOLEAN jjJET_ID_IV(leftv res, leftv u, leftv v, leftv w)
3908{
3909  res->data = (char *)idJetW((ideal)u->Data(),(int)v->Data(),
3910                             (intvec *)w->Data());
3911  return FALSE;
3912}
3913static BOOLEAN jjPREIMAGE(leftv res, leftv u, leftv v, leftv w)
3914{
3915  idhdl h;
3916  ring rr;
3917  map mapping;
3918
3919  if ((v->name==NULL) || (w->name==NULL))
3920  {
3921    WerrorS("2nd/3rd arguments must have names");
3922    return TRUE;
3923  }
3924  rr=(ring)u->Data();
3925  const char *ring_name=u->Name();
3926  if ((h=rr->idroot->get(v->name,myynest))!=NULL)
3927  {
3928    if (h->typ==MAP_CMD)
3929    {
3930      mapping=IDMAP(h);
3931    }
3932    else if (h->typ==IDEAL_CMD)
3933    {
3934      mapping=IDMAP(h);
3935    }
3936    else
3937    {
3938      Werror("`%s` is no map nor ideal",IDID(h));
3939      return TRUE;
3940    }
3941  }
3942  else
3943  {
3944    Werror("`%s` is not defined in `%s`",v->name,ring_name);
3945    return TRUE;
3946  }
3947  if ((h=rr->idroot->get(w->name,myynest))!=NULL)
3948  {
3949    if (h->typ==IDEAL_CMD)
3950    {
3951      res->data=(char *)maGetPreimage(rr,mapping,IDIDEAL(h));
3952    }
3953    else
3954    {
3955      Werror("`%s` is no ideal",IDID(h));
3956      return TRUE;
3957    }
3958  }
3959  else
3960  {
3961    Werror("`%s` is not defined in `%s`",w->name,ring_name);
3962    return TRUE;
3963  }
3964  return FALSE;
3965}
3966static BOOLEAN jjRANDOM_Im(leftv res, leftv u, leftv v, leftv w)
3967{
3968  int di, k;
3969  int i=(int)u->Data();
3970  int r=(int)v->Data();
3971  int c=(int)w->Data();
3972  intvec *iv = new intvec(r, c, 0);
3973  if (iv->rows()==0)
3974  {
3975    delete iv;
3976    return TRUE;
3977  }
3978  if (i!=0)
3979  {
3980    if (i<0) i = -i;
3981    di = 2 * i + 1;
3982    for (k=0; k<iv->length(); k++)
3983    {
3984#ifdef buildin_rand
3985      (*iv)[k] = ((siRand() % di) - i);
3986#else
3987      (*iv)[k] = ((rand() % di) - i);
3988#endif
3989    }
3990  }
3991  res->data = (char *)iv;
3992  return FALSE;
3993}
3994static BOOLEAN jjSUBST_Test(leftv v,leftv w,
3995  int &ringvar, poly &monomexpr)
3996{
3997  monomexpr=(poly)w->Data();
3998  if (pLength(monomexpr)>1)
3999  {
4000    Werror("`%s` substitutes a ringvar only by a term",
4001      Tok2Cmdname(SUBST_CMD));
4002    return TRUE;
4003  }
4004  if (!(ringvar=pVar((poly)v->Data())))
4005  {
4006    WerrorS("ringvar expected");
4007    return TRUE;
4008  }
4009  return FALSE;
4010}
4011static BOOLEAN jjSUBST_P(leftv res, leftv u, leftv v,leftv w)
4012{
4013  int ringvar;
4014  poly monomexpr;
4015  BOOLEAN nok=jjSUBST_Test(v,w,ringvar,monomexpr);
4016  if (nok) return TRUE;
4017  res->data = pSubst((poly)u->CopyD(POLY_CMD),ringvar,monomexpr);
4018  return FALSE;
4019}
4020static BOOLEAN jjSUBST_Id(leftv res, leftv u, leftv v,leftv w)
4021{
4022  int ringvar;
4023  poly monomexpr;
4024  BOOLEAN nok=jjSUBST_Test(v,w,ringvar,monomexpr);
4025  if (nok) return TRUE;
4026  res->data = idSubst((ideal)u->CopyD(IDEAL_CMD),ringvar,monomexpr);
4027  return FALSE;
4028}
4029static BOOLEAN jjMATRIX_Id(leftv res, leftv u, leftv v,leftv w)
4030{
4031  matrix m=mpNew((int)v->Data(),(int)w->Data());
4032  ideal I=(ideal)u->CopyD(IDEAL_CMD);
4033  int i=min(IDELEMS(I),(int)v->Data()*(int)w->Data());
4034  //for(i=i-1;i>=0;i--)
4035  //{
4036  //  m->m[i]=I->m[i];
4037  //  I->m[i]=NULL;
4038  //}
4039  memcpy4(m->m,I->m,i*sizeof(poly));
4040  memset(I->m,0,i*sizeof(poly));
4041  idDelete(&I);
4042  res->data = (char *)m;
4043  return FALSE;
4044}
4045static BOOLEAN jjMATRIX_Mo(leftv res, leftv u, leftv v,leftv w)
4046{
4047  res->data = (char *)idModule2formatedMatrix((ideal)u->CopyD(MODUL_CMD),
4048           (int)v->Data(),(int)w->Data());
4049  return FALSE;
4050}
4051static BOOLEAN jjMATRIX_Ma(leftv res, leftv u, leftv v,leftv w)
4052{
4053  matrix m=mpNew((int)v->Data(),(int)w->Data());
4054  matrix I=(matrix)u->CopyD(MATRIX_CMD);
4055  int r=min(MATROWS(I),(int)v->Data());
4056  int c=min(MATCOLS(I),(int)w->Data());
4057  int i,j;
4058  for(i=r;i>0;i--)
4059  {
4060    for(j=c;j>0;j--)
4061    {
4062      MATELEM(m,i,j)=MATELEM(I,i,j);
4063      MATELEM(I,i,j)=NULL;
4064    }
4065  }
4066  idDelete((ideal *)&I);
4067  res->data = (char *)m;
4068  return FALSE;
4069}
4070static BOOLEAN jjLIFT3(leftv res, leftv u, leftv v, leftv w)
4071{
4072  int sw=(int)w->Data();
4073  if (sw==0)
4074  {
4075    return jjLIFT(res,u,v);
4076  }
4077  else
4078  {
4079    ideal m;
4080    ideal ui=(ideal)u->Data();
4081    ideal vi=(ideal)v->Data();
4082    int ul= IDELEMS(ui);
4083    int vl= IDELEMS(vi);
4084    m = idLiftNonStB(ui,vi, sw);
4085    res->data = (char *)idModule2formatedMatrix(m,ul,vl);
4086  }
4087  return FALSE;
4088}
4089static BOOLEAN jjREDUCE3_P(leftv res, leftv u, leftv v, leftv w)
4090{
4091  assumeStdFlag(v);
4092  res->data = (char *)kNF((ideal)v->Data(),currQuotient,(poly)u->Data(),
4093    0,(int)w->Data());
4094  return FALSE;
4095}
4096static BOOLEAN jjREDUCE3_ID(leftv res, leftv u, leftv v, leftv w)
4097{
4098  assumeStdFlag(v);
4099  res->data = (char *)kNF((ideal)v->Data(),currQuotient,(ideal)u->Data(),
4100    0,(int)w->Data());
4101  return FALSE;
4102}
4103#ifdef OLD_RES
4104static BOOLEAN jjRES3(leftv res, leftv u, leftv v, leftv w)
4105{
4106  int maxl=(int)v->Data();
4107  int l=0;
4108  resolvente r;
4109  intvec **weights=NULL;
4110  int wmaxl=maxl;
4111  maxl--;
4112  if ((maxl==-1) && (iiOp!=MRES_CMD))
4113    maxl = pVariables-1;
4114  if ((iiOp == RES_CMD) || (iiOp == MRES_CMD))
4115  {
4116    if (BTEST1(28))
4117    {
4118      r=syMinRes((ideal)u->Data(),maxl,&l,iiOp==MRES_CMD);
4119    }
4120    else
4121    {
4122      intvec * iv=(intvec*)atGet(u,"isHomog");
4123      if (iv!=NULL)
4124      {
4125        weights = (intvec**)Alloc0(sizeof(intvec*));
4126        weights[0] = ivCopy(iv);
4127        l=1;
4128      }
4129      r=syResolvente((ideal)u->Data(),maxl,&l, &weights, iiOp==MRES_CMD);
4130    }
4131  }
4132  else
4133    r=sySchreyerResolvente((ideal)u->Data(),maxl+1,&l);
4134  if (r==NULL) return TRUE;
4135  int t3=u->Typ();
4136  iiMakeResolv(r,l,wmaxl,w->name,t3,weights);
4137  return FALSE;
4138  return TRUE;
4139}
4140#endif
4141static BOOLEAN jjSTATUS3(leftv res, leftv u, leftv v, leftv w)
4142{
4143  int yes;
4144  jjSTATUS2(res, u, v);
4145  yes = (strcmp((char *) res->data, (char *) w->Data()) == 0);
4146  FreeL((ADDRESS) res->data);
4147  res->data = (void *) yes;
4148  return FALSE;
4149}
4150static BOOLEAN jjSTD_HILB_W(leftv res, leftv u, leftv v, leftv w)
4151{
4152  ideal result;
4153  intvec *ww=(intvec *)atGet(u,"isHomog");
4154  tHomog hom=testHomog;
4155  if (ww!=NULL)
4156  {
4157    ww=ivCopy(ww);
4158    hom=isHomog;
4159  }
4160  result=kStd((ideal)(u->Data()),
4161              currQuotient,
4162              hom,
4163              &ww,                  // module weights
4164              (intvec *)v->Data(),  // hilbert series
4165              0,0,                  // syzComp, newIdeal
4166              (intvec *)w->Data()); // weights of vars
4167  idSkipZeroes(result);
4168  res->data = (char *)result;
4169  setFlag(res,FLAG_STD);
4170  if (ww!=NULL) atSet(res,mstrdup("isHomog"),ww,INTVEC_CMD);
4171  return FALSE;
4172}
4173
4174/*=================== operations with 3 args.: table =================*/
4175struct sValCmd3 dArith3[]=
4176{
4177// operations:
4178// proc             cmd          res         arg1        arg2        arg3
4179 {jjBRACK_S,        '[',        STRING_CMD, STRING_CMD, INT_CMD,    INT_CMD }
4180,{jjBRACK_Im,       '[',        INT_CMD,    INTMAT_CMD, INT_CMD,    INT_CMD }
4181,{jjBRACK_Ma_I_IV,  '[',        INT_CMD,    INTMAT_CMD, INT_CMD,    INTVEC_CMD }
4182,{jjBRACK_Ma_IV_I,  '[',        INT_CMD,    INTMAT_CMD, INTVEC_CMD, INT_CMD }
4183,{jjBRACK_Ma_IV_IV, '[',        INT_CMD,    INTMAT_CMD, INTVEC_CMD, INTVEC_CMD }
4184,{jjBRACK_Ma,       '[',        POLY_CMD,   MATRIX_CMD, INT_CMD,    INT_CMD }
4185,{jjBRACK_Ma_I_IV,  '[',        POLY_CMD,   MATRIX_CMD, INT_CMD,    INTVEC_CMD }
4186,{jjBRACK_Ma_IV_I,  '[',        POLY_CMD,   MATRIX_CMD, INTVEC_CMD, INT_CMD }
4187,{jjBRACK_Ma_IV_IV, '[',        POLY_CMD,   MATRIX_CMD, INTVEC_CMD, INTVEC_CMD }
4188,{jjPROC3,          '(',        ANY_TYPE,   PROC_CMD,   DEF_CMD,    DEF_CMD }
4189,{atATTRIB3,        ATTRIB_CMD, NONE,       IDHDL,      STRING_CMD, DEF_CMD }
4190,{jjCOEFFS3_P,      COEFFS_CMD, MATRIX_CMD, POLY_CMD,   POLY_CMD,   MATRIX_CMD }
4191,{jjCOEFFS3_P,      COEFFS_CMD, MATRIX_CMD, VECTOR_CMD, POLY_CMD,   MATRIX_CMD }
4192,{jjCOEFFS3_Id,     COEFFS_CMD, MATRIX_CMD, IDEAL_CMD,  POLY_CMD,   MATRIX_CMD }
4193,{jjCOEFFS3_Id,     COEFFS_CMD, MATRIX_CMD, MODUL_CMD,  POLY_CMD,   MATRIX_CMD }
4194,{jjCOEFFS3_KB,     COEFFS_CMD, MATRIX_CMD, IDEAL_CMD,  IDEAL_CMD,  POLY_CMD }
4195,{jjCOEFFS3_KB,     COEFFS_CMD, MATRIX_CMD, MODUL_CMD,  MODUL_CMD,  POLY_CMD }
4196,{jjELIMIN_HILB,    ELIMINATION_CMD,IDEAL_CMD, IDEAL_CMD, POLY_CMD, INTVEC_CMD }
4197,{jjFIND3,          FIND_CMD,   INT_CMD,    STRING_CMD, STRING_CMD, INT_CMD }
4198,{jjCALL3MANY,      IDEAL_CMD,  IDEAL_CMD,  DEF_CMD,    DEF_CMD,    DEF_CMD }
4199//,{jjCALL3MANY,      INTERSECT_CMD,  NONE,   DEF_CMD,    DEF_CMD,    DEF_CMD }
4200,{lInsert3,         INSERT_CMD, LIST_CMD,   LIST_CMD,   DEF_CMD,    INT_CMD }
4201,{jjCALL3MANY,      INTVEC_CMD, INTVEC_CMD, DEF_CMD,    DEF_CMD,    DEF_CMD }
4202,{jjJET_P_IV,       JET_CMD,    POLY_CMD,   POLY_CMD,   INT_CMD,    INTVEC_CMD }
4203,{jjJET_ID_IV,      JET_CMD,    IDEAL_CMD,  IDEAL_CMD,  INT_CMD,    INTVEC_CMD }
4204,{jjJET_P_IV,       JET_CMD,    VECTOR_CMD, VECTOR_CMD, INT_CMD,    INTVEC_CMD }
4205,{jjJET_ID_IV,      JET_CMD,    MODUL_CMD,  MODUL_CMD,  INT_CMD,    INTVEC_CMD }
4206,{mpKoszul,         KOSZUL_CMD, MATRIX_CMD, INT_CMD,    INT_CMD,    IDEAL_CMD }
4207,{jjCALL3MANY,      LIST_CMD,   LIST_CMD,   DEF_CMD,    DEF_CMD,    DEF_CMD }
4208,{jjMATRIX_Id,      MATRIX_CMD, MATRIX_CMD, IDEAL_CMD,  INT_CMD,    INT_CMD }
4209,{jjMATRIX_Mo,      MATRIX_CMD, MATRIX_CMD, MODUL_CMD,  INT_CMD,    INT_CMD }
4210,{jjMATRIX_Ma,      MATRIX_CMD, MATRIX_CMD, MATRIX_CMD, INT_CMD,    INT_CMD }
4211,{jjCALL3MANY,      MODUL_CMD,  MODUL_CMD,  DEF_CMD,    DEF_CMD,    DEF_CMD }
4212#ifdef OLD_RES
4213,{jjRES3,           MRES_CMD,   NONE,       IDEAL_CMD,  INT_CMD,    ANY_TYPE }
4214,{jjRES3,           MRES_CMD,   NONE,       MODUL_CMD,  INT_CMD,    ANY_TYPE }
4215#endif
4216,{jjLIFT3,          LIFT_CMD,   MATRIX_CMD, IDEAL_CMD,  IDEAL_CMD,  INT_CMD }
4217,{jjPREIMAGE,       PREIMAGE_CMD, IDEAL_CMD, RING_CMD,  ANY_TYPE,   ANY_TYPE }
4218,{jjPREIMAGE,       PREIMAGE_CMD, IDEAL_CMD, QRING_CMD, ANY_TYPE,   ANY_TYPE }
4219,{jjRANDOM_Im,      RANDOM_CMD, INTMAT_CMD, INT_CMD,    INT_CMD,    INT_CMD }
4220,{jjREDUCE3_P,      REDUCE_CMD, POLY_CMD,   POLY_CMD,   IDEAL_CMD,  INT_CMD }
4221,{jjREDUCE3_P,      REDUCE_CMD, VECTOR_CMD, VECTOR_CMD, IDEAL_CMD,  INT_CMD }
4222,{jjREDUCE3_P,      REDUCE_CMD, VECTOR_CMD, VECTOR_CMD, MODUL_CMD,  INT_CMD }
4223,{jjREDUCE3_ID,     REDUCE_CMD, IDEAL_CMD,  IDEAL_CMD,  IDEAL_CMD,  INT_CMD }
4224,{jjREDUCE3_ID,     REDUCE_CMD, MODUL_CMD,  MODUL_CMD,  MODUL_CMD,  INT_CMD }
4225,{jjREDUCE3_ID,     REDUCE_CMD, MODUL_CMD,  MODUL_CMD,  IDEAL_CMD,  INT_CMD }
4226#ifdef OLD_RES
4227,{jjRES3,           RES_CMD,    NONE,       IDEAL_CMD,  INT_CMD,    ANY_TYPE }
4228,{jjRES3,           RES_CMD,    NONE,       MODUL_CMD,  INT_CMD,    ANY_TYPE }
4229#endif
4230#ifdef HAVE_FACTORY
4231,{jjRESULTANT,      RESULTANT_CMD, POLY_CMD,POLY_CMD,   POLY_CMD,   POLY_CMD }
4232#else
4233,{jjWRONG3,         RESULTANT_CMD, POLY_CMD,POLY_CMD,   POLY_CMD,   POLY_CMD }
4234#endif
4235#ifdef OLD_RES
4236,{jjRES3,           SRES_CMD,   NONE,       IDEAL_CMD,  INT_CMD,    ANY_TYPE }
4237,{jjRES3,           SRES_CMD,   NONE,       MODUL_CMD,  INT_CMD,    ANY_TYPE }
4238#endif
4239,{jjSTATUS3,        STATUS_CMD, INT_CMD,    LINK_CMD,   STRING_CMD, STRING_CMD}
4240,{jjSTD_HILB_W,     STD_CMD,    IDEAL_CMD,  IDEAL_CMD,  INTVEC_CMD, INTVEC_CMD}
4241,{jjSTD_HILB_W,     STD_CMD,    MODUL_CMD,  MODUL_CMD,  INTVEC_CMD, INTVEC_CMD}
4242,{jjSUBST_P,        SUBST_CMD,  POLY_CMD,   POLY_CMD,   POLY_CMD,   POLY_CMD }
4243,{jjSUBST_P,        SUBST_CMD,  VECTOR_CMD, VECTOR_CMD, POLY_CMD,   POLY_CMD }
4244,{jjSUBST_Id,       SUBST_CMD,  IDEAL_CMD,  IDEAL_CMD,  POLY_CMD,   POLY_CMD }
4245,{jjSUBST_Id,       SUBST_CMD,  MODUL_CMD,  MODUL_CMD,  POLY_CMD,   POLY_CMD }
4246,{jjSUBST_Id,       SUBST_CMD,  MATRIX_CMD, MATRIX_CMD, POLY_CMD,   POLY_CMD }
4247,{jjCALL3MANY,      SYSTEM_CMD, NONE,       STRING_CMD, DEF_CMD,    DEF_CMD }
4248,{NULL,             0,          0,          0,          0,          0 }
4249};
4250/*=================== operations with many arg.: static proc =================*/
4251static BOOLEAN jjCALL1ARG(leftv res, leftv v)
4252{
4253  return iiExprArith1(res,v,iiOp);
4254}
4255static BOOLEAN jjCALL2ARG(leftv res, leftv u)
4256{
4257  leftv v=u->next;
4258  u->next=NULL;
4259  BOOLEAN b=iiExprArith2(res,u,iiOp,v);
4260  u->next=v;
4261  return b;
4262}
4263static BOOLEAN jjCALL3ARG(leftv res, leftv u)
4264{
4265  leftv v = u->next;
4266  leftv w = v->next;
4267  u->next = NULL;
4268  v->next = NULL;
4269  BOOLEAN b = iiExprArith3(res, iiOp, u, v, w);
4270  u->next = v;
4271  v->next = w;
4272  return b;
4273}
4274
4275static BOOLEAN jjCOEF_M(leftv res, leftv v)
4276{
4277  if((v->Typ() != VECTOR_CMD)
4278  || (v->next->Typ() != POLY_CMD)
4279  || (v->next->next->Typ() != MATRIX_CMD)
4280  || (v->next->next->next->Typ() != MATRIX_CMD))
4281     return TRUE;
4282  if (v->next->next->rtyp!=IDHDL) return TRUE;
4283  idhdl c=(idhdl)v->next->next->data;
4284  if (v->next->next->next->rtyp!=IDHDL) return TRUE;
4285  idhdl m=(idhdl)v->next->next->next->data;
4286  idDelete((ideal *)&(c->data.uideal));
4287  idDelete((ideal *)&(m->data.uideal));
4288  mpCoef2((poly)v->Data(),(poly)v->next->Data(),
4289    (matrix *)&(c->data.umatrix),(matrix *)&(m->data.umatrix));
4290  return FALSE;
4291}
4292static BOOLEAN jjIDEAL_PL(leftv res, leftv v)
4293{
4294  int s=1;
4295  leftv h=v;
4296  if (h!=NULL) s=exprlist_length(h);
4297  ideal id=idInit(s,1);
4298  int rank=1;
4299  int i=0;
4300  poly p=NULL;
4301  while (h!=NULL)
4302  {
4303    switch(h->Typ())
4304    {
4305      case POLY_CMD:
4306      {
4307        p=(poly)h->CopyD();
4308        break;
4309      }
4310      case INT_CMD:
4311      {
4312        number n=nInit((int)h->Data());
4313        if (!nIsZero(n))
4314        {
4315          p=pOne();
4316          pSetCoeff(p,n);
4317        }
4318        else
4319        {
4320          //p=NULL;
4321          nDelete(&n);
4322        }
4323        break;
4324      }
4325      case NUMBER_CMD:
4326      {
4327        number n=(number)h->CopyD();
4328        if (!nIsZero(n))
4329        {
4330          p=pOne();
4331          pSetCoeff(p,n);
4332        }
4333        else
4334        {
4335          //p=NULL;
4336          nDelete(&n);
4337        }
4338        break;
4339      }
4340      case VECTOR_CMD:
4341      {
4342        p=(poly)h->CopyD();
4343        if (iiOp!=MODUL_CMD)
4344        {
4345          idDelete(&id);
4346          pDelete(&p);
4347          return TRUE;
4348        }
4349        rank=max(rank,pMaxComp(p));
4350        break;
4351      }
4352      default:
4353      {
4354        idDelete(&id);
4355        return TRUE;
4356      }
4357    }
4358    if ((iiOp==MODUL_CMD)&&(p!=NULL)&&(pGetComp(p)==0))
4359    {
4360      pSetCompP(p,1);
4361    }
4362    id->m[i]=p;
4363    i++;
4364    h=h->next;
4365  }
4366  id->rank=rank;
4367  res->data=(char *)id;
4368  return FALSE;
4369}
4370static BOOLEAN jjINTERSECT_PL(leftv res, leftv v)
4371{
4372  leftv h=v;
4373  int l=v->listLength();
4374  resolvente r=(resolvente)Alloc0(l*sizeof(ideal));
4375  BOOLEAN *copied=(BOOLEAN *)Alloc0(l*sizeof(BOOLEAN));
4376  int t=0;
4377  // try to convert to IDEAL_CMD
4378  while (h!=NULL)
4379  {
4380    if (iiTestConvert(h->Typ(),IDEAL_CMD)!=0)
4381    {
4382      t=IDEAL_CMD;
4383    }
4384    else break;
4385    h=h->next;
4386  }
4387  // if failure, try MODUL_CMD
4388  if (t==0)
4389  {
4390    h=v;
4391    while (h!=NULL)
4392    {
4393      if (iiTestConvert(h->Typ(),MODUL_CMD)!=0)
4394      {
4395        t=MODUL_CMD;
4396      }
4397      else break;
4398      h=h->next;
4399    }
4400  }
4401  // check for success  in converting
4402  if (t==0)
4403  {
4404    WerrorS("cannot convert to ideal or module");
4405    return TRUE;
4406  }
4407  // call idMultSect
4408  h=v;
4409  int i=0;
4410  sleftv tmp;
4411  while (h!=NULL)
4412  {
4413    if (h->Typ()==t)
4414    {
4415      r[i]=(ideal)h->Data(); /*no copy*/
4416      h=h->next;
4417    }
4418    else if(iiConvert(h->Typ(),t,iiTestConvert(h->Typ(),t),h,&tmp))
4419    {
4420      Free((ADDRESS)copied,l*sizeof(BOOLEAN));
4421      Free((ADDRESS)r,l*sizeof(ideal));
4422      Werror("cannot convert arg. %d to %s",i+1,Tok2Cmdname(t));
4423      return TRUE;
4424    }
4425    else
4426    {
4427      r[i]=(ideal)tmp.Data(); /*now it's a copy*/
4428      copied[i]=TRUE;
4429      h=tmp.next;
4430    }
4431    i++;
4432  }
4433  res->rtyp=t;
4434  res->data=(char *)idMultSect(r,i);
4435  while(i>0)
4436  {
4437    i--;
4438    if (copied[i]) idDelete(&(r[i]));
4439  }
4440  Free((ADDRESS)copied,l*sizeof(BOOLEAN));
4441  Free((ADDRESS)r,l*sizeof(ideal));
4442  return FALSE;
4443}
4444static BOOLEAN jjINTVEC_PL(leftv res, leftv v)
4445{
4446  int i=0;
4447  leftv h=v;
4448  if (h!=NULL) i=exprlist_length(h);
4449  intvec *iv=new intvec(i);
4450  i=0;
4451  while (h!=NULL)
4452  {
4453    if(h->Typ()==INT_CMD)
4454    {
4455      (*iv)[i]=(int)h->Data();
4456    }
4457    else
4458    {
4459      delete iv;
4460      return TRUE;
4461    }
4462    i++;
4463    h=h->next;
4464  }
4465  res->data=(char *)iv;
4466  return FALSE;
4467}
4468static BOOLEAN jjKLAMMER_PL(leftv res, leftv u)
4469{
4470  leftv v=u->next;
4471  BOOLEAN b;
4472  if(v==NULL)
4473    b=iiExprArith1(res,u,iiOp);
4474  else
4475  {
4476    u->next=NULL;
4477    b=iiExprArith2(res,u,iiOp,v);
4478    u->next=v;
4479  }
4480  return b;
4481}
4482static BOOLEAN jjLIST_PL(leftv res, leftv v)
4483{
4484  int sl=0;
4485  if (v!=NULL) sl = v->listLength();
4486  lists L;
4487  if((sl==1)&&(v->Typ()==RESOLUTION_CMD))
4488  {
4489    L=syConvRes((syStrategy)v->Data());
4490  }
4491  else
4492  {
4493    L=(lists)Alloc(sizeof(slists));
4494    leftv h=NULL;
4495    int i;
4496    int rt;
4497
4498    L->Init(sl);
4499    for (i=0;i<sl;i++)
4500    {
4501      if (h!=NULL) { /* e.g. not in the first step:
4502                     * h is the pointer to the old sleftv,
4503                     * v is the pointer to the next sleftv
4504                     * (in this moment) */
4505                     h->next=v;
4506                   }
4507      h=v;
4508      v=v->next;
4509      h->next=NULL;
4510      rt=h->Typ();
4511      if (rt==0)
4512      {
4513        L->Clean();
4514        Werror("`%s` is undefined",h->Fullname());
4515        return TRUE;
4516      }
4517      if ((rt==RING_CMD)||(rt==QRING_CMD))
4518      {
4519        L->m[i].rtyp=rt;
4520        L->m[i].data=h->Data();
4521        ((ring)L->m[i].data)->ref++;
4522      }
4523      else
4524        L->m[i].Copy(h);
4525    }
4526  }
4527  res->data=(char *)L;
4528  return FALSE;
4529}
4530static BOOLEAN jjNAMES0(leftv res, leftv v)
4531{
4532  res->data=(void *)ipNameList(IDROOT);
4533  return FALSE;
4534}
4535static BOOLEAN jjOPTION_PL(leftv res, leftv v)
4536{
4537  if(v==NULL)
4538  {
4539    res->data=(char *)showOption();
4540    return FALSE;
4541  }
4542  res->rtyp=NONE;
4543  return setOption(res,v);
4544}
4545static BOOLEAN jjRESERVED0(leftv res, leftv v)
4546{
4547  int i=1;
4548  loop
4549  {
4550    Print("%-20s",cmds[i].name);
4551    i++;
4552    if(cmds[i].name==NULL)
4553      break;
4554    if ((i%3)==1) PrintLn();
4555  }
4556  PrintLn();
4557  return FALSE;
4558}
4559static BOOLEAN jjSTRING_PL(leftv res, leftv v)
4560{
4561  if ((v!=NULL)&&(v->next==NULL))
4562    return iiExprArith1(res,v,iiOp);
4563  int sl=0;
4564  leftv h=v,nachher;
4565  sleftv tmpR;
4566  BOOLEAN bo;
4567  char *sum=NULL;
4568  while(h!=NULL)
4569  {
4570    /* convert to string =================================*/
4571    nachher=h->next;
4572    h->next=NULL;
4573    bo=iiExprArith1(&tmpR,h,iiOp);
4574    if(bo)
4575    {
4576      h->next=nachher;
4577      FreeL(sum);
4578      return TRUE;
4579    }
4580    /* append to 'sum' =================================*/
4581    if(sum==NULL)
4582    {
4583      sum=(char *)tmpR.data;
4584      sl=strlen(sum);
4585    }
4586    else
4587    {
4588      int nl=strlen((char *)tmpR.data);
4589      char *s=(char *)AllocL(sl+nl+1);
4590      memcpy(s,sum,sl);
4591      memcpy(s+sl,(char *)tmpR.data,nl);
4592      s[sl+nl]='\0';
4593      sl+=nl;
4594      FreeL(sum);
4595      sum=s;
4596    }
4597    /* next sleftv =================================*/
4598    h->next=nachher;
4599    h=nachher;
4600  }
4601  if (sum==NULL) sum=mstrdup("");
4602  res->data=(char *)sum;
4603  return FALSE;
4604}
4605static BOOLEAN jjTEST(leftv res, leftv v)
4606{
4607  do
4608  {
4609    if (v->Typ()!=INT_CMD)
4610      return TRUE;
4611    test_cmd((int)v->Data());
4612    v=v->next;
4613  }
4614  while (v!=NULL);
4615  return FALSE;
4616}
4617
4618#ifndef __MWERKS__
4619static BOOLEAN jjSTATUS_M(leftv res, leftv v)
4620{
4621  if ((v->Typ() != LINK_CMD) ||
4622      (v->next->Typ() != STRING_CMD) ||
4623      (v->next->next->Typ() != STRING_CMD) ||
4624      (v->next->next->next->Typ() != INT_CMD))
4625    return TRUE;
4626  jjSTATUS3(res, v, v->next, v->next->next);
4627#if defined(HAVE_USLEEP)
4628  if (((int) res->data) == 0)
4629  {
4630    int i_s = (int) v->next->next->next->Data();
4631    if (i_s > 0)
4632    {
4633      usleep((int) v->next->next->next->Data());
4634      jjSTATUS3(res, v, v->next, v->next->next);
4635    }
4636  }
4637#elif defined(HAVE_SLEEP)
4638  if (((int) res->data) == 0)
4639  {
4640    int i_s = (int) v->next->next->next->Data();
4641    if (i_s > 0)
4642    {
4643      sleep((is - 1)/1000000 + 1);
4644      jjSTATUS3(res, v, v->next, v->next->next);
4645    }
4646  }
4647#endif
4648  return FALSE;
4649}
4650#endif
4651
4652#ifdef HAVE_NAMESPACES
4653static BOOLEAN jjIMPORTFROM(leftv res, leftv v);
4654
4655static BOOLEAN jjEXPORTTO(leftv res, leftv v)
4656{
4657  BOOLEAN nok=TRUE;
4658  leftv u=v;
4659  if(u->rtyp==NSHDL)
4660  {
4661    namehdl ns = (namehdl)(u->data);
4662    idhdl h = namespaceroot->root->get(ns->name, 0, TRUE);
4663    //Print("Export to '%s', lev %d\n", ns->name, ns->myynest);
4664    while(v->next!=NULL)
4665    {
4666      nok = iiInternalExport(v->next, ns->myynest, h);
4667      if(nok) { return nok; }
4668
4669      v = v->next;
4670    }
4671    return FALSE;
4672  }
4673  if(u->Typ()==PACKAGE_CMD)
4674  {
4675    //PrintS("export to package\n");
4676    while(v->next!=NULL)
4677    {
4678      nok = iiInternalExport(v->next, 0, u->data);
4679      if(nok) return nok;
4680      v = v->next;
4681    }
4682    return FALSE;
4683  }
4684  return TRUE;
4685}
4686#endif /* HAVE_NAMESPACES */
4687#ifdef HAVE_NAMESPACES
4688static BOOLEAN jjIMPORTFROM(leftv res, leftv v)
4689{
4690  BOOLEAN nok=FALSE;
4691  PrintS("jjIMPORT_FROM()\n");
4692  if(v->rtyp==NSHDL)
4693  {
4694    PrintS("Import from toplevel\n");
4695//While-schleife!!!
4696    return FALSE;
4697  }
4698  if(v->Typ()==PACKAGE_CMD)
4699  {
4700    Print("Import from package %s\n", v->name);
4701    while(v->next!=NULL)
4702    {
4703      //nok = iiInternalImport(v->next, 0, v->data);
4704      if(nok) return nok;
4705      v = v->next;
4706    }
4707    return FALSE;
4708  }
4709  return TRUE;
4710}
4711#endif
4712#ifdef HAVE_NAMESPACES
4713static BOOLEAN jjUNLOAD(leftv res, leftv v)
4714{
4715  if(v->Typ()==PACKAGE_CMD)
4716  {
4717    char *typ;
4718    idhdl h = (idhdl)v->data;
4719    package d=(package)v->Data();
4720    switch (d->language)
4721    {
4722      case LANG_C:
4723        typ="object";
4724        break;
4725      case LANG_SINGULAR:
4726      case LANG_NONE:
4727      default:
4728        killhdl(h);
4729    }
4730    return FALSE;
4731  }
4732  return TRUE;
4733}
4734#endif
4735/*=================== operations with many arg.: table =================*/
4736/* number_of_args:  -1: any, -2: any >0, .. */
4737struct sValCmdM dArithM[]=
4738{
4739// operations:
4740// proc         cmd               res            number_of_args
4741 {jjKLAMMER_PL,  '(',           ANY_TYPE,           -2  }
4742,{jjCALL2ARG,  COEF_CMD,        MATRIX_CMD,         2  }
4743,{jjCOEF_M,    COEF_CMD,        NONE,               4  }
4744,{jjDBPRINT,   DBPRINT_CMD,     NONE,               -2 }
4745,{jjCALL1ARG,  IDEAL_CMD,       IDEAL_CMD,          1  }
4746,{jjIDEAL_PL,  IDEAL_CMD,       IDEAL_CMD,          -1 }
4747,{jjCALL2ARG,  INTERSECT_CMD,   IDEAL_CMD,          2  }
4748,{jjINTERSECT_PL,INTERSECT_CMD, IDEAL_CMD,          -2 }
4749,{jjCALL1ARG,  INTVEC_CMD,      INTVEC_CMD,         1  }
4750,{jjINTVEC_PL, INTVEC_CMD,      INTVEC_CMD,         -1 }
4751,{jjLIST_PL,   LIST_CMD,        LIST_CMD,           -1 }
4752,{jjCALL1ARG,  MODUL_CMD,       MODUL_CMD,          1  }
4753,{jjIDEAL_PL,  MODUL_CMD,       MODUL_CMD,          -1 }
4754,{jjCALL1ARG,  NAMES_CMD,       LIST_CMD,            1 }
4755,{jjNAMES0,    NAMES_CMD,       LIST_CMD,            0 }
4756,{jjOPTION_PL, OPTION_CMD,      STRING_CMD/*or set by p*/,-1 }
4757,{jjCALL1ARG,  RESERVEDNAME_CMD, INT_CMD,            1 }
4758,{jjRESERVED0, RESERVEDNAME_CMD, NONE,               0 }
4759,{jjSTRING_PL, STRING_CMD,      STRING_CMD,         -1 }
4760,{jjSYSTEM,    SYSTEM_CMD,      NONE/*or set by p*/,-2 }
4761,{jjTEST,      TEST_CMD,        NONE,               -2 }
4762,{iiWRITE,     WRITE_CMD,       NONE,               -2 }
4763,{jjCALL2ARG,  STATUS_CMD,      STRING_CMD,          2 }
4764,{jjCALL3ARG,  STATUS_CMD,      INT_CMD,             3 }
4765#ifndef __MWERKS__
4766,{jjSTATUS_M,  STATUS_CMD,      INT_CMD,             4 }
4767#endif
4768#ifdef HAVE_NAMESPACES
4769,{jjIMPORTFROM,IMPORTFROM_CMD,  ANY_TYPE,           -2 }
4770,{jjEXPORTTO,  EXPORTTO_CMD,    NONE,               -2 }
4771,{jjUNLOAD,    UNLOAD_CMD,      NONE,               -2 }
4772#endif /* HAVE_NAMESPACES */
4773,{NULL,        0,               0,                  0  }
4774};
4775#ifdef MDEBUG
4776static Subexpr jjDBMakeSub(leftv e,char *f, int l)
4777#else
4778static Subexpr jjMakeSub(leftv e)
4779#endif
4780{
4781  // assume: e->Typ()==INT_CMD
4782#ifdef MDEBUG
4783  Subexpr r=(Subexpr)mmDBAllocBlock0(sizeof(*r),f,l);
4784#else
4785  Subexpr r=(Subexpr)Alloc0(sizeof(*r));
4786#endif
4787  r->start =(int)e->Data();
4788  return r;
4789}
4790
4791/*================ generating tables ============================*/
4792#ifdef GENTABLE
4793extern struct sValAssign dAssign[];
4794struct sValCmdTab dArithTab1[]={ {0,0}};
4795#define JJTAB1LEN 0
4796struct sValCmdTab dArithTab2[]={ {0,0}};
4797#define JJTAB2LEN 0
4798void ttGen1()
4799{
4800  FILE *outfile = myfopen("iparith.inc","w");
4801  int i,j,l1=0,l2=0;
4802  currRing=(ring)Alloc(sizeof(*currRing));
4803  fprintf(outfile,
4804  "/****************************************\n"
4805  "*  Computer Algebra System SINGULAR     *\n"
4806  "****************************************/\n\n");
4807/*-------------------------------------------------------------------*/
4808  fprintf(outfile,"// syntax table for Singular\n//\n");
4809  fprintf(outfile,"// - search for an exact match of the argument types\n");
4810  fprintf(outfile,"// - otherwise search for the first possibility\n");
4811  fprintf(outfile,"//   with converted types of the arguments\n");
4812  fprintf(outfile,"// - otherwise report an error\n//\n");
4813
4814  int op;
4815  i=0;
4816  while ((op=dArith1[i].cmd)!=0)
4817  {
4818    if (dArith1[i].p==jjWRONG)
4819      fprintf(outfile,"// DUMMY ");
4820    char *s = iiTwoOps(op);
4821    fprintf(outfile,"// operation: %s (%s)  ->  %s\n",
4822          s,
4823          Tok2Cmdname(dArith1[i].arg),
4824          Tok2Cmdname(ABS(dArith1[i].res)));
4825    i++;
4826  }
4827  fprintf(outfile,"/*---------------------------------------------*/\n");
4828  i=0;
4829  while ((op=dArith2[i].cmd)!=0)
4830  {
4831    if (dArith2[i].p==jjWRONG2)
4832      fprintf(outfile,"// DUMMY ");
4833    char *s = iiTwoOps(op);
4834    fprintf(outfile,"// operation: %s (%s, %s)  ->  %s\n",
4835          s,
4836          Tok2Cmdname(dArith2[i].arg1),
4837          Tok2Cmdname(dArith2[i].arg2),
4838          Tok2Cmdname(dArith2[i].res));
4839    i++;
4840  }
4841  fprintf(outfile,"/*---------------------------------------------*/\n");
4842  i=0;
4843  while ((op=dArith3[i].cmd)!=0)
4844  {
4845    char *s = iiTwoOps(op);
4846    if (dArith3[i].p==jjWRONG3)
4847      fprintf(outfile,"// DUMMY ");
4848    fprintf(outfile,"// operation: %s (%s, %s, %s)  ->  %s\n",
4849          s,
4850          Tok2Cmdname(dArith3[i].arg1),
4851          Tok2Cmdname(dArith3[i].arg2),
4852          Tok2Cmdname(dArith3[i].arg3),
4853          Tok2Cmdname(dArith3[i].res));
4854    i++;
4855  }
4856  fprintf(outfile,"/*---------------------------------------------*/\n");
4857  i=0;
4858  while ((op=dArithM[i].cmd)!=0)
4859  {
4860    char *s = iiTwoOps(op);
4861    fprintf(outfile,"// operation: %s (...)  ->  %s",
4862          s,
4863          Tok2Cmdname(dArithM[i].res));
4864    switch(dArithM[i].number_of_args)
4865    {
4866      case -2:
4867         fprintf(outfile," ( number of arguments >0 )\n");
4868         break;
4869      case -1:
4870         fprintf(outfile," ( any number of arguments )\n");
4871         break;
4872      default:
4873         fprintf(outfile," ( %d arguments )\n",dArithM[i].number_of_args);
4874         break;
4875    }
4876    i++;
4877  }
4878  fprintf(outfile,"/*---------------------------------------------*/\n");
4879  i=0;
4880  while ((op=dAssign[i].res)!=0)
4881  {
4882    fprintf(outfile,"// assign: %s =  %s\n",
4883          Tok2Cmdname(op/*dAssign[i].res*/),
4884          Tok2Cmdname(dAssign[i].arg));
4885    i++;
4886  }
4887/*-------------------------------------------------------------------*/
4888  fprintf(outfile,"/*---------------------------------------------*/\n");
4889  for (j=257;j<=MAX_TOK+1;j++)
4890  {
4891    for(i=257;i<=MAX_TOK+1;i++)
4892    {
4893      if ((i!=j) && (j!=IDHDL) && (j!=DEF_CMD) && (j!=ANY_TYPE)
4894      && iiTestConvert(i,j))
4895      {
4896        fprintf(outfile,"// convert %s -> %s\n",
4897          Tok2Cmdname(i), Tok2Cmdname(j));
4898        if (j==ANY_TYPE) break;
4899      }
4900    }
4901  }
4902  fprintf(outfile,"/*---------------------------------------------*/\n");
4903  char ops[]="=><+*/[.^,%(;";
4904  for(i=0;ops[i]!='\0';i++)
4905    fprintf(outfile,"// token %d : %c\n", (int)ops[i], ops[i]);
4906  for (i=257;i<=MAX_TOK;i++)
4907  {
4908    char *s=iiTwoOps(i);
4909    if (s[0]!='$')
4910    {
4911      fprintf(outfile,"// token %d : %s\n", i, s);
4912    }
4913  }
4914/*-------------------------------------------------------------------*/
4915  fprintf(outfile,"/*--max. token: %d, gr: %d --*/\n",MAX_TOK,UMINUS);
4916/*-------------------------------------------------------------------*/
4917  fprintf(outfile,"/*---------------------------------------------*/\n");
4918  fprintf(outfile,
4919  "struct sValCmdTab dArithTab1[]=\n"
4920  "{\n");
4921  for (j=1;j<=MAX_TOK+1;j++)
4922  {
4923    for(i=0;dArith1[i].cmd!=0;i++)
4924    {
4925      if (dArith1[i].cmd==j)
4926      {
4927        fprintf(outfile," { %d,%d },\n",j,i);
4928        l1++;
4929        break;
4930      }
4931    }
4932  }
4933  fprintf(outfile," { 10000,0 }\n};\n");
4934  fprintf(outfile,"#define JJTAB1LEN %d\n",l1);
4935/*-------------------------------------------------------------------*/
4936  fprintf(outfile,
4937  "struct sValCmdTab dArithTab2[]=\n"
4938  "{\n");
4939  for (j=1;j<=MAX_TOK+1;j++)
4940  {
4941    for(i=0;dArith2[i].cmd!=0;i++)
4942    {
4943      if (dArith2[i].cmd==j)
4944      {
4945        fprintf(outfile," { %d,%d },\n",j,i);
4946        l2++;
4947        break;
4948      }
4949    }
4950  }
4951  fprintf(outfile," { 10000,0 }\n};\n");
4952  fprintf(outfile,"#define JJTAB2LEN %d\n",l2);
4953  fclose(outfile);
4954}
4955/*-------------------------------------------------------------------*/
4956void ttGen2()
4957{
4958  FILE *outfile = myfopen("iparith.inc","a");
4959  fprintf(outfile,
4960  "/****************************************\n"
4961  "*  Computer Algebra System SINGULAR     *\n"
4962  "****************************************/\n\n");
4963/*-------------------------------------------------------------------*/
4964  fprintf(outfile,"// identifier table for Singular\n//\n");
4965
4966  fprintf(outfile,
4967  "cmdnames cmds[] =\n"
4968  "{  // name-string     alias  tokval toktype\n"
4969  "{ \"$INVALID$\",            0,  -1, 0},\n");
4970  int i=1;
4971  int m=-1;
4972  int id_nr=0;
4973  BOOLEAN f=FALSE;
4974  loop
4975  {
4976    while (cmds[i].tokval!=0)
4977    {
4978      if ((cmds[i].tokval!=-1)&&(cmds[i].name!=NULL))
4979      {
4980        if(m==-1)
4981        {
4982          m=i;
4983          f=TRUE;
4984        }
4985        else if(strcmp(cmds[m].name,cmds[i].name)>0)
4986        {
4987          m=i;
4988          f=TRUE;
4989        }
4990      }
4991      i++;
4992    }
4993    if(f)
4994    {
4995      id_nr++;
4996      if(cmds[m].tokval==VRTIMER) fprintf(outfile,"#ifdef HAVE_RTIMER\n");
4997      fprintf(outfile,"  {\"%s\", %*d, %3d, ",cmds[m].name,
4998                                             20-strlen(cmds[m].name),
4999                                             cmds[m].alias,
5000                                             cmds[m].tokval);
5001      switch(cmds[m].toktype)
5002      {
5003        case CMD_1:            fprintf(outfile,"CMD_1 },\n"); break;
5004        case CMD_2:            fprintf(outfile,"CMD_2 },\n"); break;
5005        case CMD_3:            fprintf(outfile,"CMD_3 },\n"); break;
5006        case CMD_12:           fprintf(outfile,"CMD_12 },\n"); break;
5007        case CMD_123 :         fprintf(outfile,"CMD_123 },\n"); break;
5008        case CMD_23:           fprintf(outfile,"CMD_23 },\n"); break;
5009        case CMD_M:            fprintf(outfile,"CMD_M },\n"); break;
5010        case SYSVAR:           fprintf(outfile,"SYSVAR },\n"); break;
5011        case ROOT_DECL:        fprintf(outfile,"ROOT_DECL },\n"); break;
5012        case ROOT_DECL_LIST:   fprintf(outfile,"ROOT_DECL_LIST },\n"); break;
5013        case RING_DECL:        fprintf(outfile,"RING_DECL },\n"); break;
5014        case NONE:             fprintf(outfile,"NONE },\n"); break;
5015        default:               if((cmds[m].toktype>' ')
5016                               &&(cmds[m].toktype<127))
5017                               {
5018                                 fprintf(outfile,"'%c' },\n",cmds[m].toktype);
5019                               }
5020                               else
5021                               {
5022                                 fprintf(outfile,"%d },\n",cmds[m].toktype);
5023                               }
5024                               break;
5025      }
5026      if(cmds[m].tokval==VRTIMER) fprintf(outfile,"#endif\n");
5027      cmds[m].name=NULL;
5028      m=-1;
5029      i=1;
5030      f=FALSE;
5031    }
5032    else break;
5033  }
5034  fprintf(outfile,
5035"/* list of scanner identifiers/only for feread/reservedName */\n");
5036  f=FALSE;
5037  i=1;m=-1;
5038  loop
5039  {
5040    while (cmds[i].tokval!=0)
5041    {
5042      if (cmds[i].name!=NULL)
5043      {
5044        if(m==-1)
5045        {
5046          m=i;
5047          f=TRUE;
5048        }
5049        else if(strcmp(cmds[m].name,cmds[i].name)>0)
5050        {
5051          m=i;
5052          f=TRUE;
5053        }
5054      }
5055      i++;
5056    }
5057    if(f)
5058    {
5059      fprintf(outfile,"  {\"%s\", %*d,  -1, 0 },\n",cmds[m].name,
5060                                             20-strlen(cmds[m].name),
5061                                             0/*cmds[m].alias*/
5062                                             /*-1 cmds[m].tokval*/
5063                                             /*0 cmds[m].toktype*/);
5064      cmds[m].name=NULL;
5065      m=-1;
5066      i=1;
5067      f=FALSE;
5068    }
5069    else break;
5070  }
5071  fprintf(outfile,
5072"/* end of list marker */\n"
5073"  { NULL, 0, 0, 0}\n"
5074"};\n"
5075"#ifdef HAVE_RTIMER\n"
5076"#define LAST_IDENTIFIER %d\n"
5077"#else\n"
5078"#define LAST_IDENTIFIER %d\n"
5079"#endif\n",id_nr,id_nr-1);
5080  fclose(outfile);
5081}
5082/*-------------------------------------------------------------------*/
5083#if 0
5084void ttGen3()
5085{
5086  FILE *outfile = myfopen("mpsr_tok.inc","w");
5087  fprintf(outfile,
5088  "/****************************************\n"
5089  "*  Computer Algebra System SINGULAR     *\n"
5090  "****************************************/\n\n");
5091/*-------------------------------------------------------------------*/
5092  fprintf(outfile,"// token table for Singular\n//\n");
5093
5094  fprintf(outfile,
5095  "short vtok[] =\n"
5096  "{\n");
5097  // operations with 1 arg: ===========================================
5098  int i=0;
5099  while (dArith1[i].cmd!=0)
5100  {
5101    if ((dArith1[i].p!=jjWRONG)
5102    &&((i==0)||(dArith1[i].cmd!=dArith1[i-1].cmd)))
5103    {
5104      fprintf(outfile,"  %d,\n",dArith1[i].cmd);
5105    }
5106    i++;
5107  }
5108  // operations with 2 args: ===========================================
5109  i=0;
5110  while (dArith2[i].cmd!=0)
5111  {
5112    if ((dArith2[i].p!=jjWRONG2)
5113    &&((i==0)||(dArith2[i].cmd!=dArith2[i-1].cmd)))
5114    {
5115      fprintf(outfile,"  %d,\n",dArith2[i].cmd);
5116    }
5117    i++;
5118  }
5119  // operations with 3 args: ===========================================
5120  i=0;
5121  while (dArith3[i].cmd!=0)
5122  {
5123    if (
5124    ((i==0)||(dArith3[i].cmd!=dArith3[i-1].cmd)))
5125    {
5126      fprintf(outfile,"  %d,\n",dArith3[i].cmd);
5127    }
5128    i++;
5129  }
5130  // operations with many args: ===========================================
5131  i=0;
5132  while (dArithM[i].cmd!=0)
5133  {
5134    if (
5135    ((i==0)||(dArithM[i].cmd!=dArithM[i-1].cmd)))
5136    {
5137      fprintf(outfile,"  %d,\n",dArithM[i].cmd);
5138    }
5139    i++;
5140  }
5141  // ====================================================================
5142  fprintf(outfile,
5143  "/* end of list marker */\n"
5144  " %d };\n",MAX_TOK);
5145  fclose(outfile);
5146}
5147#endif
5148/*-------------------------------------------------------------------*/
5149#else
5150#include "iparith.inc"
5151#endif
5152
5153/*=================== operations with 2 args. ============================*/
5154
5155BOOLEAN iiExprArith2(leftv res, leftv a, int op, leftv b, BOOLEAN proccall)
5156{
5157#ifndef GENTABLE
5158  memset(res,0,sizeof(sleftv));
5159
5160  if (!errorreported)
5161  {
5162#ifdef SIQ
5163    if (siq>0)
5164    {
5165      //Print("siq:%d\n",siq);
5166      command d=(command)Alloc0(sizeof(ip_command));
5167      memcpy(&d->arg1,a,sizeof(sleftv));
5168      memcpy(&d->arg2,b,sizeof(sleftv));
5169      d->argc=2;
5170      d->op=op;
5171      res->data=(char *)d;
5172      res->rtyp=COMMAND;
5173      return FALSE;
5174    }
5175#endif
5176    int at=a->Typ();
5177    int bt=b->Typ();
5178    int i=iiTabIndex(dArithTab2,JJTAB2LEN,op);
5179
5180    iiOp=op;
5181    while (dArith2[i].cmd==op)
5182    {
5183      if ((at==dArith2[i].arg1)
5184      && (bt==dArith2[i].arg2))
5185      {
5186        res->rtyp=dArith2[i].res;
5187        if (dArith2[i].p(res,a,b))
5188        {
5189          break;// leave loop, goto error handling
5190        }
5191        a->CleanUp();
5192        b->CleanUp();
5193        //Print("op: %d,result typ:%d\n",op,res->rtyp);
5194        return FALSE;
5195      }
5196      i++;
5197    }
5198    // implicite type conversion ----------------------------------------------
5199    if (dArith2[i].cmd!=op)
5200    {
5201      int ai,bi;
5202      leftv an = (leftv)Alloc0(sizeof(sleftv));
5203      leftv bn = (leftv)Alloc0(sizeof(sleftv));
5204      BOOLEAN failed=FALSE;
5205      i=iiTabIndex(dArithTab2,JJTAB2LEN,op);
5206      //Print("op: %c, type: %s %s\n",op,Tok2Cmdname(at),Tok2Cmdname(bt));
5207      while (dArith2[i].cmd==op)
5208      {
5209        //Print("test %s %s\n",Tok2Cmdname(dArith2[i].arg1),Tok2Cmdname(dArith2[i].arg2));
5210        if ((ai=iiTestConvert(at,dArith2[i].arg1))!=0)
5211        {
5212          if ((bi=iiTestConvert(bt,dArith2[i].arg2))!=0)
5213          {
5214            res->rtyp=dArith2[i].res;
5215            failed= ((iiConvert(at,dArith2[i].arg1,ai,a,an))
5216            || (iiConvert(bt,dArith2[i].arg2,bi,b,bn))
5217            || (dArith2[i].p(res,an,bn)));
5218            // everything done, clean up temp. variables
5219            if (failed)
5220            {
5221              // leave loop, goto error handling
5222              break;
5223            }
5224            else
5225            {
5226              // everything ok, clean up and return
5227              an->CleanUp();
5228              bn->CleanUp();
5229              Free((ADDRESS)an,sizeof(sleftv));
5230              Free((ADDRESS)bn,sizeof(sleftv));
5231              a->CleanUp();
5232              b->CleanUp();
5233              return FALSE;
5234            }
5235          }
5236        }
5237        i++;
5238      }
5239      an->CleanUp();
5240      bn->CleanUp();
5241      Free((ADDRESS)an,sizeof(sleftv));
5242      Free((ADDRESS)bn,sizeof(sleftv));
5243    }
5244    // error handling ---------------------------------------------------
5245    const char *s=NULL;
5246    if (!errorreported)
5247    {
5248      if ((at==0) && (a->Fullname()!=sNoName))
5249      {
5250        s=a->Fullname();
5251      }
5252      else if ((bt==0) && (b->Fullname()!=sNoName))
5253      {
5254        s=b->Fullname();
5255      }
5256      if (s!=NULL)
5257        Werror("`%s` is not defined",s);
5258      else
5259      {
5260        i=iiTabIndex(dArithTab2,JJTAB2LEN,op);
5261        s = iiTwoOps(op);
5262        if (proccall)
5263        {
5264          Werror("%s(`%s`,`%s`) is not supported"
5265                ,s,Tok2Cmdname(at),Tok2Cmdname(bt));
5266        }
5267        else
5268        {
5269          Werror("`%s` %s `%s` is not supported"
5270                ,Tok2Cmdname(at),s,Tok2Cmdname(bt));
5271        }
5272        if (BVERBOSE(V_SHOW_USE))
5273        {
5274          while (dArith2[i].cmd==op)
5275          {
5276            if(((at==dArith2[i].arg1)||(bt==dArith2[i].arg2))
5277            && (dArith2[i].res!=0)
5278            && (dArith2[i].p!=jjWRONG2))
5279            {
5280              if (proccall)
5281                Werror("expected %s(`%s`,`%s`)"
5282                  ,s,Tok2Cmdname(dArith2[i].arg1),Tok2Cmdname(dArith2[i].arg2));
5283              else
5284                Werror("expected `%s` %s `%s`"
5285                  ,Tok2Cmdname(dArith2[i].arg1),s,Tok2Cmdname(dArith2[i].arg2));
5286            }
5287            i++;
5288          }
5289        }
5290      }
5291    }
5292    res->rtyp = UNKNOWN;
5293  }
5294  a->CleanUp();
5295  b->CleanUp();
5296#endif
5297  return TRUE;
5298}
5299
5300/*==================== operations with 1 arg. ===============================*/
5301
5302BOOLEAN iiExprArith1(leftv res, leftv a, int op)
5303{
5304#ifndef GENTABLE
5305  memset(res,0,sizeof(sleftv));
5306
5307  if (!errorreported)
5308  {
5309#ifdef SIQ
5310    if (siq>0)
5311    {
5312      //Print("siq:%d\n",siq);
5313      command d=(command)Alloc0(sizeof(ip_command));
5314      memcpy(&d->arg1,a,sizeof(sleftv));
5315      d->op=op;
5316      d->argc=1;
5317      res->data=(char *)d;
5318      res->rtyp=COMMAND;
5319      return FALSE;
5320    }
5321#endif
5322    int at=a->Typ();
5323    BOOLEAN failed=FALSE;
5324
5325    iiOp=op;
5326    int i=iiTabIndex(dArithTab1,JJTAB1LEN,op);
5327    while (dArith1[i].cmd==op)
5328    {
5329      if (at==dArith1[i].arg)
5330      {
5331        res->rtyp=dArith1[i].res;
5332        if (dArith1[i].res<0)
5333        {
5334          res->rtyp=-res->rtyp;
5335          #ifdef PROC_BUG
5336          dArith1[i].p(res,a);
5337          #else
5338          res->data=(char *)((Proc1)dArith1[i].p)((char *)a->Data());
5339          #endif
5340        }
5341        else if (dArith1[i].p(res,a))
5342        {
5343          break;// leave loop, goto error handling
5344        }
5345        if (a->Next()!=NULL)
5346        {
5347          res->next=(leftv)Alloc(sizeof(sleftv));
5348          failed=iiExprArith1(res->next,a->next,op);
5349        }
5350        a->CleanUp();
5351        return failed;
5352      }
5353      i++;
5354    }
5355    // implicite type conversion --------------------------------------------
5356    if (dArith1[i].cmd!=op)
5357    {
5358      int ai;
5359      leftv an = (leftv)Alloc0(sizeof(sleftv));
5360      i=iiTabIndex(dArithTab1,JJTAB1LEN,op);
5361      //Print("fuer %c , typ: %s\n",op,Tok2Cmdname(at));
5362      while (dArith1[i].cmd==op)
5363      {
5364        //Print("test %s\n",Tok2Cmdname(dArith1[i].arg));
5365        if ((ai=iiTestConvert(at,dArith1[i].arg))!=0)
5366        {
5367          res->rtyp=dArith1[i].res;
5368          if (dArith1[i].res<0)
5369          {
5370            res->rtyp=-res->rtyp;
5371            failed= iiConvert(at,dArith1[i].arg,ai,a,an);
5372            if (!failed)
5373            {
5374              #ifdef PROC_BUG
5375              dArith1[i].p(res,a);
5376              #else
5377              res->data=(char *)((Proc1)dArith1[i].p)((char *)an->Data());
5378              #endif
5379            }
5380          }
5381          else
5382          {
5383            failed= ((iiConvert(at,dArith1[i].arg,ai,a,an))
5384            || (dArith1[i].p(res,an)));
5385          }
5386          // everything done, clean up temp. variables
5387          if (failed)
5388          {
5389            // leave loop, goto error handling
5390            break;
5391          }
5392          else
5393          {
5394            // everything ok, clean up and return
5395            an->CleanUp();
5396            Free((ADDRESS)an,sizeof(sleftv));
5397            a->CleanUp();
5398            return FALSE;
5399          }
5400        }
5401        i++;
5402      }
5403      an->CleanUp();
5404      Free((ADDRESS)an,sizeof(sleftv));
5405    }
5406    // error handling
5407    if (!errorreported)
5408    {
5409      if ((at==0) && (a->Fullname()!=sNoName))
5410      {
5411        Werror("`%s` is not defined",a->Fullname());
5412      }
5413      else
5414      {
5415        i=iiTabIndex(dArithTab1,JJTAB1LEN,op);
5416        char *s = iiTwoOps(op);
5417        Werror("%s(`%s`) is not supported"
5418                ,s,Tok2Cmdname(at));
5419        if (BVERBOSE(V_SHOW_USE))
5420        {
5421          while (dArith1[i].cmd==op)
5422          {
5423            if ((dArith1[i].res!=0)
5424            && (dArith1[i].p!=jjWRONG))
5425              Werror("expected %s(`%s`)"
5426                ,s,Tok2Cmdname(dArith1[i].arg));
5427            i++;
5428          }
5429        }
5430      }
5431    }
5432    res->rtyp = UNKNOWN;
5433  }
5434  a->CleanUp();
5435#endif
5436  return TRUE;
5437}
5438
5439/*=================== operations with 3 args. ============================*/
5440
5441BOOLEAN iiExprArith3(leftv res, int op, leftv a, leftv b, leftv c)
5442{
5443#ifndef GENTABLE
5444  memset(res,0,sizeof(sleftv));
5445
5446  if (!errorreported)
5447  {
5448#ifdef SIQ
5449    if (siq>0)
5450    {
5451      //Print("siq:%d\n",siq);
5452      command d=(command)Alloc0(sizeof(ip_command));
5453      memcpy(&d->arg1,a,sizeof(sleftv));
5454      memcpy(&d->arg2,b,sizeof(sleftv));
5455      memcpy(&d->arg3,c,sizeof(sleftv));
5456      d->op=op;
5457      d->argc=3;
5458      res->data=(char *)d;
5459      res->rtyp=COMMAND;
5460      return FALSE;
5461    }
5462#endif
5463    int at=a->Typ();
5464    int bt=b->Typ();
5465    int ct=c->Typ();
5466
5467    iiOp=op;
5468    int i=0;
5469    while ((dArith3[i].cmd!=op)&&