[0e1846] | 1 | /**************************************** |
---|
| 2 | * Computer Algebra System SINGULAR * |
---|
| 3 | ****************************************/ |
---|
| 4 | /* |
---|
[6ae4f5] | 5 | * ABSTRACT: handling of leftv |
---|
[0e1846] | 6 | */ |
---|
| 7 | |
---|
| 8 | #include <stdlib.h> |
---|
| 9 | #include <stdio.h> |
---|
| 10 | #include <string.h> |
---|
| 11 | #include <ctype.h> |
---|
| 12 | #include <unistd.h> |
---|
| 13 | |
---|
[762407] | 14 | #include "config.h" |
---|
[b1dfaf] | 15 | #include <kernel/mod2.h> |
---|
[5ff1d3] | 16 | |
---|
| 17 | #include <omalloc/omalloc.h> |
---|
| 18 | |
---|
[0fb34ba] | 19 | #include <misc/intvec.h> |
---|
| 20 | #include <misc/options.h> |
---|
[5ff1d3] | 21 | |
---|
| 22 | |
---|
| 23 | #include <coeffs/ffields.h> |
---|
| 24 | #include <coeffs/numbers.h> |
---|
[75f10d] | 25 | #include <coeffs/bigintmat.h> |
---|
[5ff1d3] | 26 | |
---|
[0fb34ba] | 27 | #include <polys/monomials/maps.h> |
---|
| 28 | #include <polys/matpol.h> |
---|
[5ff1d3] | 29 | #include <polys/monomials/ring.h> |
---|
[737a68] | 30 | #include <kernel/polys.h> |
---|
[5ff1d3] | 31 | |
---|
| 32 | #include <kernel/longrat.h> |
---|
| 33 | // #include <coeffs/longrat.h> |
---|
| 34 | |
---|
| 35 | #include <kernel/febase.h> |
---|
| 36 | #include <kernel/ideals.h> |
---|
[599326] | 37 | #include <kernel/kstd1.h> |
---|
| 38 | #include <kernel/timer.h> |
---|
[5ff1d3] | 39 | #include <kernel/syz.h> |
---|
| 40 | |
---|
| 41 | #include <Singular/tok.h> |
---|
| 42 | #include <Singular/ipid.h> |
---|
[599326] | 43 | #include <Singular/ipshell.h> |
---|
| 44 | #include <Singular/lists.h> |
---|
| 45 | #include <Singular/attrib.h> |
---|
[44ca2f] | 46 | #include <Singular/links/silink.h> |
---|
[599326] | 47 | #include <Singular/attrib.h> |
---|
| 48 | #include <Singular/subexpr.h> |
---|
[b117e8] | 49 | #include <Singular/blackbox.h> |
---|
[855a8f] | 50 | |
---|
[5ff1d3] | 51 | |
---|
| 52 | |
---|
[0e051f] | 53 | omBin sSubexpr_bin = omGetSpecBin(sizeof(_ssubexpr)); |
---|
[c232af] | 54 | omBin sleftv_bin = omGetSpecBin(sizeof(sleftv)); |
---|
| 55 | omBin procinfo_bin = omGetSpecBin(sizeof(procinfo)); |
---|
| 56 | omBin libstack_bin = omGetSpecBin(sizeof(libstack)); |
---|
[ec7aac] | 57 | static omBin size_two_bin = omGetSpecBin(2); |
---|
[c232af] | 58 | |
---|
[0e1846] | 59 | sleftv sLastPrinted; |
---|
| 60 | const char sNoName[]="_"; |
---|
| 61 | #ifdef SIQ |
---|
| 62 | BOOLEAN siq=FALSE; |
---|
| 63 | #endif |
---|
| 64 | |
---|
| 65 | int sleftv::listLength() |
---|
| 66 | { |
---|
| 67 | int n = 1; |
---|
| 68 | leftv sl = next; |
---|
| 69 | while (sl!=NULL) |
---|
| 70 | { |
---|
| 71 | n++; |
---|
| 72 | sl=sl->next; |
---|
| 73 | } |
---|
| 74 | return n; |
---|
| 75 | } |
---|
| 76 | |
---|
| 77 | void sleftv::Print(leftv store, int spaces) |
---|
| 78 | { |
---|
| 79 | int t=Typ(); |
---|
| 80 | if (errorreported) return; |
---|
| 81 | #ifdef SIQ |
---|
| 82 | if (rtyp==COMMAND) |
---|
| 83 | { |
---|
| 84 | command c=(command)data; |
---|
| 85 | char ch[2]; |
---|
| 86 | ch[0]=c->op;ch[1]='\0'; |
---|
[85e68dd] | 87 | const char *s=ch; |
---|
[2c89d2] | 88 | if (c->op>127) s=iiTwoOps(c->op); |
---|
[0e1846] | 89 | ::Print("##command %d(%s), %d args\n", |
---|
| 90 | c->op, s, c->argc); |
---|
| 91 | if (c->argc>0) |
---|
| 92 | c->arg1.Print(NULL,spaces+2); |
---|
| 93 | if(c->argc<4) |
---|
| 94 | { |
---|
| 95 | if (c->argc>1) |
---|
| 96 | c->arg2.Print(NULL,spaces+2); |
---|
| 97 | if (c->argc>2) |
---|
| 98 | c->arg3.Print(NULL,spaces+2); |
---|
| 99 | } |
---|
| 100 | PrintS("##end"); |
---|
| 101 | } |
---|
| 102 | else |
---|
| 103 | #endif |
---|
| 104 | { |
---|
[2ba9a6] | 105 | const char *n=Name(); |
---|
| 106 | char *s; |
---|
| 107 | void *d=Data(); |
---|
[49fbc3f] | 108 | if (errorreported) return; |
---|
[2ba9a6] | 109 | if ((store!=NULL)&&(store!=this)) |
---|
| 110 | store->CleanUp(); |
---|
[a3432c] | 111 | |
---|
[2ba9a6] | 112 | switch (t /*=Typ()*/) |
---|
[0e1846] | 113 | { |
---|
| 114 | case UNKNOWN: |
---|
| 115 | case DEF_CMD: |
---|
[490a23b] | 116 | PrintNSpaces(spaces); |
---|
| 117 | PrintS("`");PrintS(n);PrintS("`"); |
---|
[0e1846] | 118 | break; |
---|
[95197e] | 119 | case PACKAGE_CMD: |
---|
[4c9cf1] | 120 | PrintNSpaces(spaces); |
---|
| 121 | paPrint(n,(package)d); |
---|
| 122 | break; |
---|
[0e1846] | 123 | case NONE: |
---|
| 124 | return; |
---|
| 125 | case INTVEC_CMD: |
---|
| 126 | case INTMAT_CMD: |
---|
[6d3ffec] | 127 | ((intvec *)d)->show(t,spaces); |
---|
[0e1846] | 128 | break; |
---|
[75f10d] | 129 | case BIGINTMAT_CMD: |
---|
| 130 | ((bigintmat *)d)->pprint(80); |
---|
| 131 | break; |
---|
[0e1846] | 132 | case RING_CMD: |
---|
| 133 | case QRING_CMD: |
---|
[db34980] | 134 | { |
---|
[490a23b] | 135 | PrintNSpaces(spaces); |
---|
[db34980] | 136 | const ring r = (const ring)d; |
---|
[fb85f97] | 137 | rWrite(r, currRing == r); |
---|
[0e1846] | 138 | break; |
---|
[db34980] | 139 | } |
---|
[0e1846] | 140 | case MATRIX_CMD: |
---|
[8fa9ea] | 141 | iiWriteMatrix((matrix)d,n,2, currRing, spaces); |
---|
[0e1846] | 142 | break; |
---|
| 143 | case MODUL_CMD: |
---|
| 144 | case IDEAL_CMD: |
---|
[3e7095] | 145 | if ((TEST_V_QRING) &&(currQuotient!=NULL) |
---|
| 146 | &&(!hasFlag(this,FLAG_QRING))) |
---|
[6c4db17] | 147 | { |
---|
| 148 | jjNormalizeQRingId(this); |
---|
| 149 | d=Data(); |
---|
| 150 | } |
---|
| 151 | // no break: |
---|
| 152 | case MAP_CMD: |
---|
[8fa9ea] | 153 | iiWriteMatrix((matrix)d,n,1, currRing, spaces); |
---|
[0e1846] | 154 | break; |
---|
| 155 | case POLY_CMD: |
---|
| 156 | case VECTOR_CMD: |
---|
[3e7095] | 157 | if ((TEST_V_QRING) &&(currQuotient!=NULL) |
---|
| 158 | &&(!hasFlag(this,FLAG_QRING))) |
---|
[6c4db17] | 159 | { |
---|
| 160 | jjNormalizeQRingP(this); |
---|
| 161 | d=Data(); |
---|
| 162 | } |
---|
[490a23b] | 163 | PrintNSpaces(spaces); |
---|
[0e1846] | 164 | pWrite0((poly)d); |
---|
| 165 | break; |
---|
[dfc6b54] | 166 | case RESOLUTION_CMD: |
---|
[22579cf] | 167 | { |
---|
[e7463f] | 168 | syStrategy tmp=(syStrategy)d; |
---|
[9b01b1] | 169 | syPrint(tmp,IDID(currRingHdl)); |
---|
[a3432c] | 170 | break; |
---|
[22579cf] | 171 | } |
---|
[0e1846] | 172 | case STRING_CMD: |
---|
[490a23b] | 173 | PrintNSpaces(spaces); |
---|
| 174 | PrintS((char *)d); |
---|
[0e1846] | 175 | break; |
---|
| 176 | case INT_CMD: |
---|
[490a23b] | 177 | PrintNSpaces(spaces); |
---|
| 178 | ::Print("%d",(int)(long)d); |
---|
[0e1846] | 179 | break; |
---|
[2ba9a6] | 180 | case PROC_CMD: |
---|
| 181 | { |
---|
| 182 | procinfov pi=(procinfov)d; |
---|
[490a23b] | 183 | |
---|
| 184 | PrintNSpaces(spaces); |
---|
| 185 | PrintS("// libname : "); |
---|
| 186 | PrintS(piProcinfo(pi, "libname")); |
---|
[17681b] | 187 | PrintLn(); |
---|
[490a23b] | 188 | |
---|
| 189 | PrintNSpaces(spaces); |
---|
| 190 | PrintS("// procname : "); |
---|
| 191 | PrintS(piProcinfo(pi, "procname")); |
---|
[17681b] | 192 | PrintLn(); |
---|
[490a23b] | 193 | |
---|
| 194 | PrintNSpaces(spaces); |
---|
| 195 | PrintS("// type : "); |
---|
| 196 | PrintS(piProcinfo(pi, "type")); |
---|
[a3432c] | 197 | // ::Print("%-*.*s// ref : %s",spaces,spaces," ", |
---|
| 198 | // piProcinfo(pi, "ref")); |
---|
| 199 | break; |
---|
[2ba9a6] | 200 | } |
---|
[46d09b] | 201 | case POINTER_CMD: |
---|
| 202 | { package pack = (package)d; |
---|
[17681b] | 203 | PrintNSpaces(spaces); |
---|
[490a23b] | 204 | PrintS("// PointerTest\n"); |
---|
[17681b] | 205 | PrintNSpaces(spaces); |
---|
[490a23b] | 206 | ::Print("// %s\n",IDID(pack->idroot)); |
---|
[46d09b] | 207 | //::Print(((char *)(pack->idroot)->data), spaces); |
---|
| 208 | break; |
---|
| 209 | } |
---|
[0e1846] | 210 | case LINK_CMD: |
---|
| 211 | { |
---|
| 212 | si_link l=(si_link)d; |
---|
[17681b] | 213 | PrintNSpaces(spaces); |
---|
[490a23b] | 214 | ::Print("// type : %s\n", slStatus(l, "type")); |
---|
[17681b] | 215 | PrintNSpaces(spaces); |
---|
[490a23b] | 216 | ::Print("// mode : %s\n", slStatus(l, "mode")); |
---|
[17681b] | 217 | PrintNSpaces(spaces); |
---|
[490a23b] | 218 | ::Print("// name : %s\n", slStatus(l, "name")); |
---|
[17681b] | 219 | PrintNSpaces(spaces); |
---|
[490a23b] | 220 | ::Print("// open : %s\n", slStatus(l, "open")); |
---|
[17681b] | 221 | PrintNSpaces(spaces); |
---|
[490a23b] | 222 | ::Print("// read : %s\n", slStatus(l, "read")); |
---|
[17681b] | 223 | PrintNSpaces(spaces); |
---|
[490a23b] | 224 | ::Print("// write: %s", slStatus(l, "write")); |
---|
[d754b7] | 225 | break; |
---|
[0e1846] | 226 | } |
---|
| 227 | case NUMBER_CMD: |
---|
[84a6a7] | 228 | case BIGINT_CMD: |
---|
[0e1846] | 229 | s=String(d); |
---|
| 230 | if (s==NULL) return; |
---|
[17681b] | 231 | PrintNSpaces(spaces); |
---|
[0e1846] | 232 | PrintS(s); |
---|
[c232af] | 233 | omFree((ADDRESS)s); |
---|
[0e1846] | 234 | break; |
---|
| 235 | case LIST_CMD: |
---|
| 236 | { |
---|
| 237 | lists l=(lists)d; |
---|
[ecf019] | 238 | if (lSize(l)<0) |
---|
[17681b] | 239 | { |
---|
| 240 | PrintNSpaces(spaces); |
---|
[490a23b] | 241 | PrintS("empty list\n"); |
---|
[17681b] | 242 | } |
---|
[0e1846] | 243 | else |
---|
| 244 | { |
---|
| 245 | int i=0; |
---|
| 246 | for (;i<=l->nr;i++) |
---|
| 247 | { |
---|
| 248 | if (l->m[i].rtyp!=DEF_CMD) |
---|
| 249 | { |
---|
[17681b] | 250 | PrintNSpaces(spaces); |
---|
[490a23b] | 251 | ::Print("[%d]:\n",i+1); |
---|
[0e1846] | 252 | l->m[i].Print(NULL,spaces+3); |
---|
| 253 | } |
---|
| 254 | } |
---|
| 255 | } |
---|
| 256 | break; |
---|
| 257 | } |
---|
[b117e8] | 258 | |
---|
[0e1846] | 259 | default: |
---|
[c4dab4] | 260 | if (t>MAX_TOK) |
---|
| 261 | { |
---|
| 262 | blackbox * bb=getBlackboxStuff(t); |
---|
[b117e8] | 263 | PrintNSpaces(spaces); |
---|
[c4dab4] | 264 | if (bb!=NULL) { bb->blackbox_Print(bb,d); } |
---|
| 265 | else { ::Print("Print: blackbox %d(bb=NULL)",t); } |
---|
| 266 | } |
---|
| 267 | else |
---|
[0e1846] | 268 | ::Print("Print:unknown type %s(%d)", Tok2Cmdname(t),t); |
---|
| 269 | } /* end switch: (Typ()) */ |
---|
| 270 | } |
---|
| 271 | if (next!=NULL) |
---|
| 272 | { |
---|
| 273 | if (t==COMMAND) PrintLn(); |
---|
| 274 | else if (t!=LIST_CMD) PrintS(" "); |
---|
| 275 | next->Print(NULL,spaces); |
---|
| 276 | } |
---|
| 277 | else if (t!=LIST_CMD) |
---|
| 278 | { |
---|
| 279 | PrintLn(); |
---|
[a3432c] | 280 | } |
---|
[0e1846] | 281 | #ifdef SIQ |
---|
| 282 | if (rtyp!=COMMAND) |
---|
| 283 | #endif |
---|
| 284 | { |
---|
| 285 | if ((store!=NULL) |
---|
[e7463f] | 286 | && (store!=this)) |
---|
[0e1846] | 287 | { |
---|
[e7463f] | 288 | if((t/*Typ()*/!=LINK_CMD) |
---|
| 289 | && (t/*Typ()*/!=POINTER_CMD) |
---|
| 290 | && (t/*Typ()*/!=PACKAGE_CMD) |
---|
| 291 | && (t/*Typ()*/!=DEF_CMD) |
---|
| 292 | ) |
---|
[3074334] | 293 | { |
---|
[e7463f] | 294 | store->rtyp=t/*Typ()*/; |
---|
| 295 | store->data=CopyD(); |
---|
| 296 | if(attribute!=NULL) |
---|
| 297 | { |
---|
| 298 | store->attribute=CopyA(); |
---|
| 299 | } |
---|
| 300 | store->flag=flag; |
---|
[3074334] | 301 | } |
---|
[0e1846] | 302 | } |
---|
| 303 | } |
---|
| 304 | } |
---|
| 305 | |
---|
[16acb0] | 306 | void sleftv::CleanUp(ring r) |
---|
[0e1846] | 307 | { |
---|
[ea947e] | 308 | if ((name!=NULL) && (name!=sNoName) && (rtyp!=IDHDL) && (rtyp!=ALIAS_CMD)) |
---|
[0e1846] | 309 | { |
---|
| 310 | //::Print("free %x (%s)\n",name,name); |
---|
[c232af] | 311 | omFree((ADDRESS)name); |
---|
[0e1846] | 312 | } |
---|
[f2dff02] | 313 | //name=NULL; |
---|
| 314 | //flag=0; |
---|
[0e1846] | 315 | if (data!=NULL) |
---|
| 316 | { |
---|
[c4dab4] | 317 | if (rtyp==IDHDL) attribute=NULL; // is only a pointer to attribute of id |
---|
| 318 | else s_internalDelete(rtyp,data,r); |
---|
[ff7e529] | 319 | //data=NULL; // will be done by Init() at the end |
---|
[0e1846] | 320 | } |
---|
| 321 | if (attribute!=NULL) |
---|
| 322 | { |
---|
| 323 | switch (rtyp) |
---|
| 324 | { |
---|
[46d09b] | 325 | case POINTER_CMD: |
---|
[0e1846] | 326 | case PACKAGE_CMD: |
---|
| 327 | case IDHDL: |
---|
| 328 | case ANY_TYPE: |
---|
| 329 | case VECHO: |
---|
| 330 | case VPRINTLEVEL: |
---|
| 331 | case VCOLMAX: |
---|
| 332 | case VTIMER: |
---|
[7a7652] | 333 | case VRTIMER: |
---|
[0e1846] | 334 | case VOICE: |
---|
| 335 | case VMAXDEG: |
---|
| 336 | case VMAXMULT: |
---|
| 337 | case TRACE: |
---|
| 338 | case VSHORTOUT: |
---|
| 339 | case VNOETHER: |
---|
| 340 | case VMINPOLY: |
---|
[7a7652] | 341 | case LIB_CMD: |
---|
[0e1846] | 342 | case 0: |
---|
[ff7e529] | 343 | //attribute=NULL; // will be done by Init() at the end |
---|
[0e1846] | 344 | break; |
---|
| 345 | default: |
---|
| 346 | { |
---|
[c227a2f] | 347 | attribute->killAll(r); |
---|
[0e1846] | 348 | } |
---|
| 349 | } |
---|
| 350 | } |
---|
| 351 | Subexpr h; |
---|
| 352 | while (e!=NULL) |
---|
| 353 | { |
---|
| 354 | h=e->next; |
---|
[c232af] | 355 | omFreeBin((ADDRESS)e, sSubexpr_bin); |
---|
[0e1846] | 356 | e=h; |
---|
| 357 | } |
---|
[ff7e529] | 358 | //rtyp=NONE; // will be done by Init() at the end |
---|
[0e1846] | 359 | if (next!=NULL) |
---|
| 360 | { |
---|
[df4422] | 361 | leftv tmp_n; |
---|
| 362 | do |
---|
| 363 | { |
---|
| 364 | tmp_n=next->next; |
---|
| 365 | //next->name=NULL; |
---|
| 366 | next->next=NULL; |
---|
[16acb0] | 367 | next->CleanUp(r); |
---|
[c232af] | 368 | omFreeBin((ADDRESS)next, sleftv_bin); |
---|
[df4422] | 369 | next=tmp_n; |
---|
| 370 | } while (next!=NULL); |
---|
[0e1846] | 371 | } |
---|
[f2dff02] | 372 | Init(); |
---|
[0e1846] | 373 | } |
---|
| 374 | |
---|
[4cbe5d] | 375 | BOOLEAN sleftv::RingDependend() |
---|
| 376 | { |
---|
| 377 | int rt=Typ(); |
---|
[df5fc1] | 378 | if(::RingDependend(rt) && (rt!=QRING_CMD)) |
---|
[4cbe5d] | 379 | return TRUE; |
---|
| 380 | if (rt==LIST_CMD) |
---|
| 381 | return lRingDependend((lists)Data()); |
---|
[a3bc95e] | 382 | return FALSE; |
---|
[4cbe5d] | 383 | } |
---|
| 384 | |
---|
[a1a595] | 385 | static inline void * s_internalCopy(const int t, void *d) |
---|
[0e1846] | 386 | { |
---|
[dfc6b54] | 387 | switch (t) |
---|
[0e1846] | 388 | { |
---|
| 389 | case INTVEC_CMD: |
---|
| 390 | case INTMAT_CMD: |
---|
| 391 | return (void *)ivCopy((intvec *)d); |
---|
[75f10d] | 392 | case BIGINTMAT_CMD: |
---|
| 393 | return (void*)bimCopy((bigintmat *)d); |
---|
[0e1846] | 394 | case MATRIX_CMD: |
---|
[8fa9ea] | 395 | return (void *)mp_Copy((matrix)d, currRing); |
---|
[0e1846] | 396 | case IDEAL_CMD: |
---|
| 397 | case MODUL_CMD: |
---|
| 398 | return (void *)idCopy((ideal)d); |
---|
| 399 | case STRING_CMD: |
---|
[c232af] | 400 | return (void *)omStrDup((char *)d); |
---|
[46d09b] | 401 | case POINTER_CMD: |
---|
| 402 | return d; |
---|
[daeb6d] | 403 | case PACKAGE_CMD: |
---|
| 404 | return (void *)paCopy((package) d); |
---|
[0e1846] | 405 | case PROC_CMD: |
---|
[2ba9a6] | 406 | return (void *)piCopy((procinfov) d); |
---|
[0e1846] | 407 | case POLY_CMD: |
---|
| 408 | case VECTOR_CMD: |
---|
| 409 | return (void *)pCopy((poly)d); |
---|
| 410 | case INT_CMD: |
---|
| 411 | return d; |
---|
| 412 | case NUMBER_CMD: |
---|
| 413 | return (void *)nCopy((number)d); |
---|
[84a6a7] | 414 | case BIGINT_CMD: |
---|
[4ff07f] | 415 | return (void *)n_Copy((number)d, coeffs_BIGINT); |
---|
[0e1846] | 416 | case MAP_CMD: |
---|
[8fa9ea] | 417 | return (void *)maCopy((map)d, currRing); |
---|
[0e1846] | 418 | case LIST_CMD: |
---|
| 419 | return (void *)lCopy((lists)d); |
---|
[e6969d] | 420 | case LINK_CMD: |
---|
| 421 | return (void *)slCopy((si_link) d); |
---|
[2f12a6f] | 422 | case RING_CMD: |
---|
| 423 | case QRING_CMD: |
---|
| 424 | { |
---|
| 425 | ring r=(ring)d; |
---|
| 426 | r->ref++; |
---|
[e7463f] | 427 | //Print("+ ring %d, ref %d\n",r,r->ref); |
---|
[2f12a6f] | 428 | return d; |
---|
| 429 | } |
---|
[dfc6b54] | 430 | case RESOLUTION_CMD: |
---|
| 431 | return (void*)syCopy((syStrategy)d); |
---|
[0e1846] | 432 | case DEF_CMD: |
---|
| 433 | case NONE: |
---|
[8739c9] | 434 | case 0: /* type in error case */ |
---|
[0e1846] | 435 | break; /* error recovery: do nothing */ |
---|
| 436 | //case COMMAND: |
---|
| 437 | default: |
---|
[b117e8] | 438 | { |
---|
| 439 | if (t>MAX_TOK) |
---|
| 440 | { |
---|
| 441 | blackbox *b=getBlackboxStuff(t); |
---|
[c4dab4] | 442 | if (b!=NULL) return b->blackbox_Copy(b,d); |
---|
[b117e8] | 443 | return NULL; |
---|
| 444 | } |
---|
| 445 | else |
---|
[b70e54] | 446 | Warn("s_internalCopy: cannot copy type %s(%d)", |
---|
[a1a595] | 447 | Tok2Cmdname(t),t); |
---|
[b117e8] | 448 | } |
---|
[0e1846] | 449 | } |
---|
| 450 | return NULL; |
---|
| 451 | } |
---|
| 452 | |
---|
[c4dab4] | 453 | void s_internalDelete(const int t, void *d, const ring r) |
---|
| 454 | { |
---|
| 455 | switch (t) |
---|
| 456 | { |
---|
| 457 | case INTVEC_CMD: |
---|
| 458 | case INTMAT_CMD: |
---|
| 459 | { |
---|
| 460 | intvec *v=(intvec*)d; |
---|
| 461 | delete v; |
---|
| 462 | break; |
---|
| 463 | } |
---|
[75f10d] | 464 | case BIGINTMAT_CMD: |
---|
| 465 | { |
---|
| 466 | bigintmat *v=(bigintmat*)d; |
---|
| 467 | delete v; |
---|
| 468 | break; |
---|
| 469 | } |
---|
[c4dab4] | 470 | case MAP_CMD: |
---|
| 471 | { |
---|
| 472 | map m=(map)d; |
---|
[95197e] | 473 | omFreeBinAddr((ADDRESS)m->preimage); |
---|
[c4dab4] | 474 | m->preimage=NULL; |
---|
| 475 | /* no break: continue as IDEAL*/ |
---|
| 476 | } |
---|
| 477 | case MATRIX_CMD: |
---|
| 478 | case IDEAL_CMD: |
---|
| 479 | case MODUL_CMD: |
---|
| 480 | { |
---|
| 481 | ideal i=(ideal)d; |
---|
| 482 | id_Delete(&i,r); |
---|
| 483 | break; |
---|
| 484 | } |
---|
| 485 | case STRING_CMD: |
---|
| 486 | omFree(d); |
---|
| 487 | break; |
---|
| 488 | //case POINTER_CMD: |
---|
| 489 | // return d; |
---|
| 490 | //case PACKAGE_CMD: |
---|
| 491 | // return (void *)paCopy((package) d); |
---|
| 492 | case PROC_CMD: |
---|
| 493 | piKill((procinfo*)d); |
---|
| 494 | break; |
---|
| 495 | case POLY_CMD: |
---|
| 496 | case VECTOR_CMD: |
---|
| 497 | { |
---|
| 498 | poly p=(poly)d; |
---|
| 499 | p_Delete(&p,r); |
---|
| 500 | break; |
---|
| 501 | } |
---|
| 502 | case NUMBER_CMD: |
---|
| 503 | { |
---|
| 504 | number n=(number)d; |
---|
| 505 | n_Delete(&n,r); |
---|
| 506 | break; |
---|
| 507 | } |
---|
| 508 | case BIGINT_CMD: |
---|
| 509 | { |
---|
| 510 | number n=(number)d; |
---|
[4ff07f] | 511 | n_Delete(&n,coeffs_BIGINT); |
---|
[c4dab4] | 512 | break; |
---|
| 513 | } |
---|
| 514 | case LIST_CMD: |
---|
| 515 | { |
---|
| 516 | lists l=(lists)d; |
---|
| 517 | l->Clean(r); |
---|
| 518 | break; |
---|
| 519 | } |
---|
| 520 | case LINK_CMD: |
---|
| 521 | { |
---|
| 522 | si_link l=(si_link)d; |
---|
| 523 | slKill(l); |
---|
| 524 | break; |
---|
| 525 | } |
---|
| 526 | case RING_CMD: |
---|
| 527 | case QRING_CMD: |
---|
| 528 | { |
---|
| 529 | ring R=(ring)d; |
---|
| 530 | #ifdef TEST |
---|
[e7463f] | 531 | if ((R==currRing)&&(R->ref<=0)) |
---|
| 532 | Print("currRing? ref=%d\n",R->ref); |
---|
[c4dab4] | 533 | else |
---|
| 534 | #endif |
---|
| 535 | rKill(R); |
---|
| 536 | break; |
---|
| 537 | } |
---|
| 538 | case RESOLUTION_CMD: |
---|
| 539 | { |
---|
| 540 | syStrategy s=(syStrategy)d; |
---|
| 541 | if (s!=NULL) syKillComputation(s,r); |
---|
| 542 | break; |
---|
| 543 | } |
---|
| 544 | case COMMAND: |
---|
| 545 | { |
---|
| 546 | command cmd=(command)d; |
---|
| 547 | if (cmd->arg1.rtyp!=0) cmd->arg1.CleanUp(r); |
---|
| 548 | if (cmd->arg2.rtyp!=0) cmd->arg2.CleanUp(r); |
---|
| 549 | if (cmd->arg3.rtyp!=0) cmd->arg3.CleanUp(r); |
---|
| 550 | omFreeBin((ADDRESS)d, sip_command_bin); |
---|
| 551 | break; |
---|
| 552 | } |
---|
| 553 | case INT_CMD: |
---|
| 554 | case DEF_CMD: |
---|
| 555 | case ALIAS_CMD: |
---|
| 556 | case POINTER_CMD: |
---|
| 557 | case PACKAGE_CMD: |
---|
| 558 | case IDHDL: |
---|
| 559 | case NONE: |
---|
| 560 | case ANY_TYPE: |
---|
| 561 | case VECHO: |
---|
| 562 | case VPRINTLEVEL: |
---|
| 563 | case VCOLMAX: |
---|
| 564 | case VTIMER: |
---|
| 565 | case VRTIMER: |
---|
| 566 | case VOICE: |
---|
| 567 | case VMAXDEG: |
---|
| 568 | case VMAXMULT: |
---|
| 569 | case TRACE: |
---|
| 570 | case VSHORTOUT: |
---|
| 571 | case VNOETHER: |
---|
| 572 | case VMINPOLY: |
---|
| 573 | case LIB_CMD: |
---|
| 574 | case 0: /* type in error case */ |
---|
| 575 | break; /* error recovery: do nothing */ |
---|
| 576 | //case COMMAND: |
---|
| 577 | //case COMMAND: |
---|
| 578 | default: |
---|
| 579 | { |
---|
| 580 | if (t>MAX_TOK) |
---|
| 581 | { |
---|
| 582 | blackbox *b=getBlackboxStuff(t); |
---|
| 583 | if (b!=NULL) b->blackbox_destroy(b,d); |
---|
| 584 | break; |
---|
| 585 | } |
---|
| 586 | else |
---|
| 587 | Warn("s_internalDelete: cannot delete type %s(%d)", |
---|
| 588 | Tok2Cmdname(t),t); |
---|
| 589 | } |
---|
| 590 | } |
---|
| 591 | } |
---|
[f71681] | 592 | |
---|
[a1a595] | 593 | void * slInternalCopy(leftv source, const int t, void *d, Subexpr e) |
---|
| 594 | { |
---|
| 595 | if (t==STRING_CMD) |
---|
| 596 | { |
---|
| 597 | if ((e==NULL) |
---|
| 598 | || (source->rtyp==LIST_CMD) |
---|
| 599 | || ((source->rtyp==IDHDL)&&(IDTYP((idhdl)source->data)==LIST_CMD))) |
---|
| 600 | return (void *)omStrDup((char *)d); |
---|
| 601 | else if (e->next==NULL) |
---|
| 602 | { |
---|
| 603 | char *s=(char*)omAllocBin(size_two_bin); |
---|
| 604 | s[0]=*(char *)d; |
---|
| 605 | s[1]='\0'; |
---|
| 606 | return s; |
---|
| 607 | } |
---|
| 608 | #ifdef TEST |
---|
| 609 | else |
---|
| 610 | { |
---|
| 611 | Werror("not impl. string-op in `%s`",my_yylinebuf); |
---|
| 612 | return NULL; |
---|
| 613 | } |
---|
| 614 | #endif |
---|
| 615 | } |
---|
| 616 | return s_internalCopy(t,d); |
---|
| 617 | } |
---|
| 618 | |
---|
[0e1846] | 619 | void sleftv::Copy(leftv source) |
---|
| 620 | { |
---|
[f2dff02] | 621 | Init(); |
---|
[0e1846] | 622 | rtyp=source->Typ(); |
---|
| 623 | void *d=source->Data(); |
---|
| 624 | if(!errorreported) |
---|
| 625 | { |
---|
[a1a595] | 626 | data=s_internalCopy(rtyp,d); |
---|
[8c94ba] | 627 | if ((source->attribute!=NULL)||(source->e!=NULL)) |
---|
| 628 | attribute=source->CopyA(); |
---|
| 629 | flag=source->flag; |
---|
[24c186a] | 630 | if (source->next!=NULL) |
---|
| 631 | { |
---|
[c232af] | 632 | next=(leftv)omAllocBin(sleftv_bin); |
---|
[24c186a] | 633 | next->Copy(source->next); |
---|
| 634 | } |
---|
[0e1846] | 635 | } |
---|
| 636 | } |
---|
| 637 | |
---|
| 638 | void * sleftv::CopyD(int t) |
---|
| 639 | { |
---|
[98adcd] | 640 | if ((rtyp!=IDHDL)&&(rtyp!=ALIAS_CMD)&&(e==NULL)) |
---|
[0e1846] | 641 | { |
---|
[a3432c] | 642 | if (iiCheckRing(t)) return NULL; |
---|
[8b87135] | 643 | void *x = data; |
---|
| 644 | if (rtyp==VNOETHER) x = (void *)pCopy((currRing->ppNoether)); |
---|
| 645 | else if ((rtyp==VMINPOLY) && nCoeff_is_algExt(currRing->cf) && (!nCoeff_is_GF(currRing->cf))) |
---|
| 646 | { |
---|
| 647 | const ring A = currRing->cf->extRing; |
---|
[db34980] | 648 | |
---|
[8b87135] | 649 | assume( A != NULL ); |
---|
[dd668f] | 650 | assume( A->qideal != NULL ); |
---|
[db34980] | 651 | |
---|
[dd668f] | 652 | x=(void *)p_Copy(A->qideal->m[0], A); |
---|
[8b87135] | 653 | } |
---|
[0e1846] | 654 | data=NULL; |
---|
| 655 | return x; |
---|
| 656 | } |
---|
[a3432c] | 657 | void *d=Data(); // will also do a iiCheckRing |
---|
| 658 | if ((!errorreported) && (d!=NULL)) return slInternalCopy(this,t,d,e); |
---|
[0e1846] | 659 | return NULL; |
---|
| 660 | } |
---|
| 661 | |
---|
[2166ad3] | 662 | //void * sleftv::CopyD() |
---|
| 663 | //{ |
---|
[24c186a] | 664 | //if ((rtyp!=IDHDL)&&(e==NULL) |
---|
| 665 | //&&(rtyp!=VNOETHER)&&(rtyp!=LIB_CMD)&&(rtyp!=VMINPOLY)) |
---|
| 666 | //{ |
---|
| 667 | // void *x=data; |
---|
| 668 | // data=NULL; |
---|
| 669 | // return x; |
---|
| 670 | //} |
---|
[2166ad3] | 671 | // return CopyD(Typ()); |
---|
| 672 | //} |
---|
[0e1846] | 673 | |
---|
[8c94ba] | 674 | attr sleftv::CopyA() |
---|
| 675 | { |
---|
| 676 | attr *a=Attribute(); |
---|
| 677 | if ((a!=NULL) && (*a!=NULL)) |
---|
| 678 | return (*a)->Copy(); |
---|
| 679 | return NULL; |
---|
| 680 | } |
---|
[0e1846] | 681 | |
---|
[a79a128] | 682 | char * sleftv::String(void *d, BOOLEAN typed, int dim) |
---|
[0e1846] | 683 | { |
---|
| 684 | #ifdef SIQ |
---|
| 685 | if (rtyp==COMMAND) |
---|
| 686 | { |
---|
| 687 | ::Print("##command %d\n",((command)data)->op); |
---|
| 688 | if (((command)data)->arg1.rtyp!=0) |
---|
| 689 | ((command)data)->arg1.Print(NULL,2); |
---|
| 690 | if (((command)data)->arg2.rtyp!=0) |
---|
| 691 | ((command)data)->arg2.Print(NULL,2); |
---|
| 692 | if (((command)data)->arg3.rtyp==0) |
---|
| 693 | ((command)data)->arg3.Print(NULL,2); |
---|
| 694 | PrintS("##end\n"); |
---|
[c232af] | 695 | return omStrDup(""); |
---|
[0e1846] | 696 | } |
---|
| 697 | #endif |
---|
| 698 | if (d==NULL) d=Data(); |
---|
| 699 | if (!errorreported) |
---|
| 700 | { |
---|
| 701 | char *s; |
---|
[84a6a7] | 702 | int t=Typ(); |
---|
| 703 | switch (t /*Typ()*/) |
---|
[0e1846] | 704 | { |
---|
[4b2155] | 705 | case INT_CMD: |
---|
[7604db] | 706 | if (typed) |
---|
| 707 | { |
---|
[c232af] | 708 | s=(char *)omAlloc(MAX_INT_LEN+7); |
---|
[7447d8] | 709 | sprintf(s,"int(%d)",(int)(long)d); |
---|
[7604db] | 710 | } |
---|
| 711 | else |
---|
| 712 | { |
---|
[c232af] | 713 | s=(char *)omAlloc(MAX_INT_LEN+2); |
---|
[7447d8] | 714 | sprintf(s,"%d",(int)(long)d); |
---|
[7604db] | 715 | } |
---|
| 716 | return s; |
---|
| 717 | |
---|
[4b2155] | 718 | case STRING_CMD: |
---|
[7604db] | 719 | if (d == NULL) |
---|
| 720 | { |
---|
[c232af] | 721 | if (typed) return omStrDup("\"\""); |
---|
| 722 | return omStrDup(""); |
---|
[7604db] | 723 | } |
---|
| 724 | if (typed) |
---|
| 725 | { |
---|
[c232af] | 726 | s = (char*) omAlloc(strlen((char*) d) + 3); |
---|
[7604db] | 727 | sprintf(s,"\"%s\"", (char*) d); |
---|
| 728 | return s; |
---|
| 729 | } |
---|
| 730 | else |
---|
| 731 | { |
---|
[c232af] | 732 | return omStrDup((char*)d); |
---|
[7604db] | 733 | } |
---|
[c931a2] | 734 | |
---|
[4b2155] | 735 | case POLY_CMD: |
---|
| 736 | case VECTOR_CMD: |
---|
[7604db] | 737 | if (typed) |
---|
| 738 | { |
---|
| 739 | char* ps = pString((poly) d); |
---|
[c232af] | 740 | s = (char*) omAlloc(strlen(ps) + 10); |
---|
[84a6a7] | 741 | sprintf(s,"%s(%s)", (t /*Typ()*/ == POLY_CMD ? "poly" : "vector"), ps); |
---|
[4c9cf1] | 742 | omFree(ps); |
---|
[7604db] | 743 | return s; |
---|
| 744 | } |
---|
| 745 | else |
---|
[ea9fd60] | 746 | return pString((poly)d); |
---|
[7604db] | 747 | |
---|
[4b2155] | 748 | case NUMBER_CMD: |
---|
[a70441f] | 749 | StringSetS((char*) (typed ? "number(" : "")); |
---|
[4b2155] | 750 | if ((rtyp==IDHDL)&&(IDTYP((idhdl)data)==NUMBER_CMD)) |
---|
| 751 | { |
---|
| 752 | nWrite(IDNUMBER((idhdl)data)); |
---|
| 753 | } |
---|
| 754 | else if (rtyp==NUMBER_CMD) |
---|
| 755 | { |
---|
| 756 | number n=(number)data; |
---|
| 757 | nWrite(n); |
---|
| 758 | data=(char *)n; |
---|
| 759 | } |
---|
[8fa9ea] | 760 | else if((rtyp==VMINPOLY)&&(rField_is_GF(currRing))) |
---|
[4b2155] | 761 | { |
---|
[8fa9ea] | 762 | nfShowMipo(currRing->cf); |
---|
[4b2155] | 763 | } |
---|
| 764 | else |
---|
| 765 | { |
---|
| 766 | number n=nCopy((number)d); |
---|
| 767 | nWrite(n); |
---|
| 768 | nDelete(&n); |
---|
| 769 | } |
---|
[538512] | 770 | StringAppendS((char*) (typed ? ")" : "")); |
---|
| 771 | return StringEndS(); |
---|
[c931a2] | 772 | |
---|
[84a6a7] | 773 | case BIGINT_CMD: |
---|
| 774 | { |
---|
| 775 | StringSetS((char*) (typed ? "bigint(" : "")); |
---|
| 776 | number nl=(number)d; |
---|
[4ff07f] | 777 | n_Write(nl,coeffs_BIGINT); |
---|
[538512] | 778 | StringAppendS((char*) (typed ? ")" : "")); |
---|
| 779 | return StringEndS(); |
---|
[84a6a7] | 780 | } |
---|
| 781 | |
---|
[4b2155] | 782 | case MATRIX_CMD: |
---|
[8fa9ea] | 783 | s= iiStringMatrix((matrix)d,dim, currRing); |
---|
[7604db] | 784 | if (typed) |
---|
| 785 | { |
---|
[c232af] | 786 | char* ns = (char*) omAlloc(strlen(s) + 40); |
---|
[c931a2] | 787 | sprintf(ns, "matrix(ideal(%s),%d,%d)", s, |
---|
[7604db] | 788 | ((ideal) d)->nrows, ((ideal) d)->ncols); |
---|
[c232af] | 789 | omCheckAddr(ns); |
---|
[7604db] | 790 | return ns; |
---|
| 791 | } |
---|
| 792 | else |
---|
| 793 | { |
---|
[c232af] | 794 | return omStrDup(s); |
---|
[7604db] | 795 | } |
---|
| 796 | |
---|
[4b2155] | 797 | case MODUL_CMD: |
---|
| 798 | case IDEAL_CMD: |
---|
| 799 | case MAP_CMD: |
---|
[8fa9ea] | 800 | s= iiStringMatrix((matrix)d,dim, currRing); |
---|
[7604db] | 801 | if (typed) |
---|
| 802 | { |
---|
[c232af] | 803 | char* ns = (char*) omAlloc(strlen(s) + 10); |
---|
[84a6a7] | 804 | sprintf(ns, "%s(%s)", (t/*Typ()*/==MODUL_CMD ? "module" : "ideal"), s); |
---|
[c232af] | 805 | omCheckAddr(ns); |
---|
[7604db] | 806 | return ns; |
---|
| 807 | } |
---|
[c232af] | 808 | return omStrDup(s); |
---|
[7604db] | 809 | |
---|
[4b2155] | 810 | case INTVEC_CMD: |
---|
| 811 | case INTMAT_CMD: |
---|
| 812 | { |
---|
| 813 | intvec *v=(intvec *)d; |
---|
[a79a128] | 814 | s = v->String(dim); |
---|
[7604db] | 815 | if (typed) |
---|
| 816 | { |
---|
| 817 | char* ns; |
---|
[84a6a7] | 818 | if (t/*Typ()*/ == INTMAT_CMD) |
---|
[7604db] | 819 | { |
---|
[c232af] | 820 | ns = (char*) omAlloc(strlen(s) + 40); |
---|
[7604db] | 821 | sprintf(ns, "intmat(intvec(%s),%d,%d)", s, v->rows(), v->cols()); |
---|
| 822 | } |
---|
| 823 | else |
---|
| 824 | { |
---|
[c232af] | 825 | ns = (char*) omAlloc(strlen(s) + 10); |
---|
[7604db] | 826 | sprintf(ns, "intvec(%s)", s); |
---|
| 827 | } |
---|
[c232af] | 828 | omCheckAddr(ns); |
---|
| 829 | omFree(s); |
---|
[7604db] | 830 | return ns; |
---|
| 831 | } |
---|
| 832 | else |
---|
| 833 | return s; |
---|
[4b2155] | 834 | } |
---|
[75f10d] | 835 | case BIGINTMAT_CMD: |
---|
| 836 | { |
---|
| 837 | bigintmat *bim=(bigintmat*)d; |
---|
| 838 | s = bim->String(); |
---|
| 839 | if (typed) |
---|
| 840 | { |
---|
| 841 | char* ns = (char*) omAlloc0(strlen(s) + 40); |
---|
| 842 | sprintf(ns, "bigintmat(bigintvec(%s),%d,%d)", s, bim->rows(), bim->cols()); |
---|
| 843 | omCheckAddr(ns); |
---|
| 844 | return ns; |
---|
| 845 | } |
---|
| 846 | else |
---|
| 847 | return omStrDup(s); |
---|
[4c9cf1] | 848 | } |
---|
[7604db] | 849 | |
---|
[4b2155] | 850 | case RING_CMD: |
---|
| 851 | case QRING_CMD: |
---|
[7604db] | 852 | s = rString((ring)d); |
---|
[c931a2] | 853 | |
---|
[7604db] | 854 | if (typed) |
---|
| 855 | { |
---|
| 856 | char* ns; |
---|
[84a6a7] | 857 | if (t/*Typ()*/ == QRING_CMD) |
---|
[7604db] | 858 | { |
---|
[8fa9ea] | 859 | char* id = iiStringMatrix((matrix) ((ring) d)->qideal, dim, |
---|
[db34980] | 860 | currRing); |
---|
[c232af] | 861 | ns = (char*) omAlloc(strlen(s) + strlen(id) + 20); |
---|
[7604db] | 862 | sprintf(ns, "\"%s\";%sideal(%s)", s,(dim == 2 ? "\n" : " "), id); |
---|
| 863 | } |
---|
| 864 | else |
---|
| 865 | { |
---|
[c232af] | 866 | ns = (char*) omAlloc(strlen(s) + 4); |
---|
[7604db] | 867 | sprintf(ns, "\"%s\"", s); |
---|
| 868 | } |
---|
[c232af] | 869 | omFree(s); |
---|
| 870 | omCheckAddr(ns); |
---|
[7604db] | 871 | return ns; |
---|
| 872 | } |
---|
| 873 | return s; |
---|
[4b2155] | 874 | case RESOLUTION_CMD: |
---|
[077e9c] | 875 | { |
---|
[4b2155] | 876 | lists l = syConvRes((syStrategy)d); |
---|
[a79a128] | 877 | s = lString(l, typed, dim); |
---|
[4b2155] | 878 | l->Clean(); |
---|
[7604db] | 879 | return s; |
---|
[077e9c] | 880 | } |
---|
[7604db] | 881 | |
---|
[4b2155] | 882 | case PROC_CMD: |
---|
[0e1846] | 883 | { |
---|
[4b2155] | 884 | procinfo* pi = (procinfo*) d; |
---|
| 885 | if((pi->language == LANG_SINGULAR) && (pi->data.s.body!=NULL)) |
---|
[7604db] | 886 | s = (pi->data.s.body); |
---|
[4b2155] | 887 | else |
---|
[85e68dd] | 888 | s = (char *)""; |
---|
[7604db] | 889 | if (typed) |
---|
| 890 | { |
---|
[c232af] | 891 | char* ns = (char*) omAlloc(strlen(s) + 4); |
---|
[7604db] | 892 | sprintf(ns, "\"%s\"", s); |
---|
[c232af] | 893 | omCheckAddr(ns); |
---|
[7604db] | 894 | return ns; |
---|
| 895 | } |
---|
[c232af] | 896 | return omStrDup(s); |
---|
[0e1846] | 897 | } |
---|
[c931a2] | 898 | |
---|
[f82470] | 899 | case LINK_CMD: |
---|
[a79a128] | 900 | s = slString((si_link) d); |
---|
[7604db] | 901 | if (typed) |
---|
| 902 | { |
---|
[c232af] | 903 | char* ns = (char*) omAlloc(strlen(s) + 10); |
---|
[7604db] | 904 | sprintf(ns, "link(\"%s\")", s); |
---|
[95197e] | 905 | omFreeBinAddr(s); |
---|
[c232af] | 906 | omCheckAddr(ns); |
---|
[7604db] | 907 | return ns; |
---|
| 908 | } |
---|
| 909 | return s; |
---|
[c931a2] | 910 | |
---|
[4b2155] | 911 | case LIST_CMD: |
---|
[7604db] | 912 | return lString((lists) d, typed, dim); |
---|
[b117e8] | 913 | |
---|
| 914 | default: |
---|
[c4dab4] | 915 | if(t> MAX_TOK) |
---|
| 916 | { |
---|
| 917 | blackbox *bb=getBlackboxStuff(t); |
---|
| 918 | if (bb!=NULL) return bb->blackbox_String(bb,d); |
---|
| 919 | } |
---|
[0e1846] | 920 | } /* end switch: (Typ()) */ |
---|
| 921 | } |
---|
[c232af] | 922 | return omStrDup(""); |
---|
[0e1846] | 923 | } |
---|
| 924 | |
---|
[c931a2] | 925 | |
---|
[0e1846] | 926 | int sleftv::Typ() |
---|
| 927 | { |
---|
| 928 | if (e==NULL) |
---|
| 929 | { |
---|
[38cfbb] | 930 | switch (rtyp) |
---|
[0e1846] | 931 | { |
---|
[11ca48] | 932 | case IDHDL: |
---|
| 933 | return IDTYP((idhdl)data); |
---|
[ea947e] | 934 | case ALIAS_CMD: |
---|
| 935 | { |
---|
| 936 | idhdl h=(idhdl)data; |
---|
| 937 | return ((idhdl)h->data.ustring)->typ; |
---|
| 938 | } |
---|
[0e1846] | 939 | case VECHO: |
---|
| 940 | case VPRINTLEVEL: |
---|
| 941 | case VCOLMAX: |
---|
| 942 | case VTIMER: |
---|
[7a7652] | 943 | case VRTIMER: |
---|
[0e1846] | 944 | case VOICE: |
---|
| 945 | case VMAXDEG: |
---|
| 946 | case VMAXMULT: |
---|
| 947 | case TRACE: |
---|
| 948 | case VSHORTOUT: |
---|
| 949 | return INT_CMD; |
---|
| 950 | case VMINPOLY: |
---|
| 951 | return NUMBER_CMD; |
---|
| 952 | case VNOETHER: |
---|
| 953 | return POLY_CMD; |
---|
| 954 | //case COMMAND: |
---|
| 955 | // return COMMAND; |
---|
| 956 | default: |
---|
| 957 | return rtyp; |
---|
| 958 | } |
---|
| 959 | } |
---|
| 960 | int r=0; |
---|
| 961 | int t=rtyp; |
---|
[11ca48] | 962 | if (t==IDHDL) t=IDTYP((idhdl)data); |
---|
| 963 | else if (t==ALIAS_CMD) { idhdl h=(idhdl)IDDATA((idhdl)data); t=IDTYP(h); } |
---|
[0e1846] | 964 | switch (t) |
---|
| 965 | { |
---|
| 966 | case INTVEC_CMD: |
---|
| 967 | case INTMAT_CMD: |
---|
| 968 | r=INT_CMD; |
---|
| 969 | break; |
---|
[75f10d] | 970 | case BIGINTMAT_CMD: |
---|
| 971 | r=BIGINT_CMD; |
---|
| 972 | break; |
---|
[0e1846] | 973 | case IDEAL_CMD: |
---|
| 974 | case MATRIX_CMD: |
---|
| 975 | case MAP_CMD: |
---|
| 976 | r=POLY_CMD; |
---|
| 977 | break; |
---|
| 978 | case MODUL_CMD: |
---|
| 979 | r=VECTOR_CMD; |
---|
| 980 | break; |
---|
| 981 | case STRING_CMD: |
---|
| 982 | r=STRING_CMD; |
---|
| 983 | break; |
---|
[8db6c3] | 984 | default: |
---|
[0e1846] | 985 | { |
---|
[8db6c3] | 986 | blackbox *b=NULL; |
---|
| 987 | if (t>MAX_TOK) |
---|
[0e1846] | 988 | { |
---|
[8db6c3] | 989 | b=getBlackboxStuff(t); |
---|
[0e1846] | 990 | } |
---|
[8db6c3] | 991 | if ((t==LIST_CMD)||((b!=NULL)&&BB_LIKE_LIST(b))) |
---|
[0e1846] | 992 | { |
---|
[8db6c3] | 993 | lists l; |
---|
| 994 | if (rtyp==IDHDL) l=IDLIST((idhdl)data); |
---|
[db34980] | 995 | else if (rtyp==ALIAS_CMD) |
---|
| 996 | { |
---|
| 997 | idhdl h=(idhdl)data; |
---|
| 998 | l=(lists)(((idhdl)h->data.ustring)->data.ustring); |
---|
| 999 | } |
---|
[8db6c3] | 1000 | else l=(lists)data; |
---|
| 1001 | if ((0<e->start)&&(e->start<=l->nr+1)) |
---|
| 1002 | { |
---|
| 1003 | Subexpr tmp=l->m[e->start-1].e; |
---|
| 1004 | l->m[e->start-1].e=e->next; |
---|
| 1005 | r=l->m[e->start-1].Typ(); |
---|
| 1006 | e->next=l->m[e->start-1].e; |
---|
| 1007 | l->m[e->start-1].e=tmp; |
---|
| 1008 | } |
---|
| 1009 | else |
---|
| 1010 | { |
---|
| 1011 | //Warn("out of range: %d not in 1..%d",e->start,l->nr+1); |
---|
| 1012 | r=NONE; |
---|
| 1013 | } |
---|
[0e1846] | 1014 | } |
---|
[8db6c3] | 1015 | else |
---|
| 1016 | Werror("cannot index type %s(%d)",Tok2Cmdname(t),t); |
---|
[0e1846] | 1017 | break; |
---|
| 1018 | } |
---|
| 1019 | } |
---|
| 1020 | return r; |
---|
| 1021 | } |
---|
| 1022 | |
---|
| 1023 | int sleftv::LTyp() |
---|
| 1024 | { |
---|
| 1025 | lists l=NULL; |
---|
| 1026 | int r; |
---|
| 1027 | if (rtyp==LIST_CMD) |
---|
| 1028 | l=(lists)data; |
---|
| 1029 | else if ((rtyp==IDHDL)&& (IDTYP((idhdl)data)==LIST_CMD)) |
---|
| 1030 | l=IDLIST((idhdl)data); |
---|
| 1031 | else |
---|
| 1032 | return Typ(); |
---|
| 1033 | //if (l!=NULL) |
---|
| 1034 | { |
---|
| 1035 | if ((e!=NULL) && (e->next!=NULL)) |
---|
| 1036 | { |
---|
| 1037 | if ((0<e->start)&&(e->start<=l->nr+1)) |
---|
| 1038 | { |
---|
| 1039 | l->m[e->start-1].e=e->next; |
---|
| 1040 | r=l->m[e->start-1].LTyp(); |
---|
| 1041 | l->m[e->start-1].e=NULL; |
---|
| 1042 | } |
---|
| 1043 | else |
---|
| 1044 | { |
---|
| 1045 | //Warn("out of range: %d not in 1..%d",e->start,l->nr+1); |
---|
| 1046 | r=NONE; |
---|
| 1047 | } |
---|
| 1048 | return r; |
---|
| 1049 | } |
---|
| 1050 | return LIST_CMD; |
---|
| 1051 | } |
---|
| 1052 | return Typ(); |
---|
| 1053 | } |
---|
| 1054 | |
---|
| 1055 | void * sleftv::Data() |
---|
| 1056 | { |
---|
[17681b] | 1057 | if ((rtyp!=IDHDL) && iiCheckRing(rtyp)) |
---|
[0e1846] | 1058 | return NULL; |
---|
| 1059 | if (e==NULL) |
---|
| 1060 | { |
---|
| 1061 | switch (rtyp) |
---|
| 1062 | { |
---|
[11ca48] | 1063 | case ALIAS_CMD: |
---|
| 1064 | { |
---|
| 1065 | idhdl h=(idhdl)data; |
---|
| 1066 | return ((idhdl)h->data.ustring)->data.ustring; |
---|
| 1067 | } |
---|
[0e1846] | 1068 | case VECHO: return (void *)si_echo; |
---|
| 1069 | case VPRINTLEVEL:return (void *)printlevel; |
---|
| 1070 | case VCOLMAX: return (void *)colmax; |
---|
| 1071 | case VTIMER: return (void *)getTimer(); |
---|
[ea986d] | 1072 | case VRTIMER: return (void *)getRTimer(); |
---|
[0e1846] | 1073 | case VOICE: return (void *)(myynest+1); |
---|
| 1074 | case VMAXDEG: return (void *)Kstd1_deg; |
---|
| 1075 | case VMAXMULT: return (void *)Kstd1_mu; |
---|
| 1076 | case TRACE: return (void *)traceit; |
---|
[c5f4b9] | 1077 | case VSHORTOUT: return (void *)(currRing != NULL ? currRing->ShortOut : 0); |
---|
[8b87135] | 1078 | case VMINPOLY: |
---|
| 1079 | if ( (currRing != NULL) && nCoeff_is_algExt(currRing->cf) && !nCoeff_is_GF(currRing->cf)) |
---|
| 1080 | { |
---|
| 1081 | /* Q(a), Fp(a), but not GF(q) */ |
---|
| 1082 | const ring A = currRing->cf->extRing; |
---|
[db34980] | 1083 | |
---|
[8b87135] | 1084 | assume( A != NULL ); |
---|
[dd668f] | 1085 | assume( A->qideal != NULL ); |
---|
[db34980] | 1086 | |
---|
[dd668f] | 1087 | return (void *)A->qideal->m[0]; |
---|
[8b87135] | 1088 | } |
---|
| 1089 | else |
---|
| 1090 | return (void *)currRing->cf->nNULL; |
---|
[db34980] | 1091 | |
---|
[993ae2] | 1092 | case VNOETHER: return (void *) (currRing->ppNoether); |
---|
[0e1846] | 1093 | case IDHDL: |
---|
| 1094 | return IDDATA((idhdl)data); |
---|
[46d09b] | 1095 | case POINTER_CMD: |
---|
| 1096 | return IDDATA((idhdl)data); |
---|
[0e1846] | 1097 | case COMMAND: |
---|
| 1098 | //return NULL; |
---|
| 1099 | default: |
---|
| 1100 | return data; |
---|
| 1101 | } |
---|
| 1102 | } |
---|
| 1103 | /* e != NULL : */ |
---|
| 1104 | int t=rtyp; |
---|
| 1105 | void *d=data; |
---|
| 1106 | if (t==IDHDL) |
---|
| 1107 | { |
---|
| 1108 | t=((idhdl)data)->typ; |
---|
| 1109 | d=IDDATA((idhdl)data); |
---|
| 1110 | } |
---|
[11ca48] | 1111 | else if (t==ALIAS_CMD) |
---|
| 1112 | { |
---|
| 1113 | idhdl h=(idhdl)IDDATA((idhdl)data); |
---|
| 1114 | t=IDTYP(h); |
---|
| 1115 | d=IDDATA(h); |
---|
| 1116 | } |
---|
[0e1846] | 1117 | if (iiCheckRing(t)) |
---|
| 1118 | return NULL; |
---|
| 1119 | char *r=NULL; |
---|
[c931a2] | 1120 | int index=e->start; |
---|
[0e1846] | 1121 | switch (t) |
---|
| 1122 | { |
---|
| 1123 | case INTVEC_CMD: |
---|
| 1124 | { |
---|
| 1125 | intvec *iv=(intvec *)d; |
---|
[c931a2] | 1126 | if ((index<1)||(index>iv->length())) |
---|
[0e1846] | 1127 | { |
---|
| 1128 | if (!errorreported) |
---|
[c931a2] | 1129 | Werror("wrong range[%d] in intvec(%d)",index,iv->length()); |
---|
[0e1846] | 1130 | } |
---|
| 1131 | else |
---|
[c931a2] | 1132 | r=(char *)((*iv)[index-1]); |
---|
[0e1846] | 1133 | break; |
---|
| 1134 | } |
---|
| 1135 | case INTMAT_CMD: |
---|
| 1136 | { |
---|
| 1137 | intvec *iv=(intvec *)d; |
---|
[c931a2] | 1138 | if ((index<1) |
---|
| 1139 | ||(index>iv->rows()) |
---|
[0e1846] | 1140 | ||(e->next->start<1) |
---|
| 1141 | ||(e->next->start>iv->cols())) |
---|
| 1142 | { |
---|
| 1143 | if (!errorreported) |
---|
[c931a2] | 1144 | Werror("wrong range[%d,%d] in intmat(%dx%d)",index,e->next->start, |
---|
[0e1846] | 1145 | iv->rows(),iv->cols()); |
---|
| 1146 | } |
---|
| 1147 | else |
---|
[c931a2] | 1148 | r=(char *)(IMATELEM((*iv),index,e->next->start)); |
---|
[0e1846] | 1149 | break; |
---|
| 1150 | } |
---|
[75f10d] | 1151 | case BIGINTMAT_CMD: |
---|
| 1152 | { |
---|
| 1153 | bigintmat *m=(bigintmat *)d; |
---|
| 1154 | if ((index<1) |
---|
| 1155 | ||(index>m->rows()) |
---|
| 1156 | ||(e->next->start<1) |
---|
| 1157 | ||(e->next->start>m->cols())) |
---|
| 1158 | { |
---|
| 1159 | if (!errorreported) |
---|
| 1160 | Werror("wrong range[%d,%d] in bigintmat(%dx%d)",index,e->next->start, |
---|
| 1161 | m->rows(),m->cols()); |
---|
| 1162 | } |
---|
| 1163 | else |
---|
| 1164 | r=(char *)(BIMATELEM((*m),index,e->next->start)); |
---|
| 1165 | break; |
---|
| 1166 | } |
---|
[0e1846] | 1167 | case IDEAL_CMD: |
---|
| 1168 | case MODUL_CMD: |
---|
| 1169 | case MAP_CMD: |
---|
| 1170 | { |
---|
| 1171 | ideal I=(ideal)d; |
---|
[c931a2] | 1172 | if ((index<1)||(index>IDELEMS(I))) |
---|
[0e1846] | 1173 | { |
---|
| 1174 | if (!errorreported) |
---|
[c931a2] | 1175 | Werror("wrong range[%d] in ideal/module(%d)",index,IDELEMS(I)); |
---|
[0e1846] | 1176 | } |
---|
| 1177 | else |
---|
[c931a2] | 1178 | r=(char *)I->m[index-1]; |
---|
[0e1846] | 1179 | break; |
---|
| 1180 | } |
---|
| 1181 | case STRING_CMD: |
---|
| 1182 | { |
---|
[c931a2] | 1183 | // this was a memory leak |
---|
| 1184 | // we evalute it, cleanup and replace this leftv by it's evalutated form |
---|
| 1185 | // the evalutated form will be build in tmp |
---|
| 1186 | sleftv tmp; |
---|
[f2dff02] | 1187 | tmp.Init(); |
---|
[c931a2] | 1188 | tmp.rtyp=STRING_CMD; |
---|
[ec7aac] | 1189 | r=(char *)omAllocBin(size_two_bin); |
---|
[c931a2] | 1190 | if ((index>0)&& (index<=(int)strlen((char *)d))) |
---|
[0e1846] | 1191 | { |
---|
[c931a2] | 1192 | r[0]=*(((char *)d)+index-1); |
---|
[0e1846] | 1193 | r[1]='\0'; |
---|
| 1194 | } |
---|
| 1195 | else |
---|
| 1196 | { |
---|
| 1197 | r[0]='\0'; |
---|
| 1198 | } |
---|
[c931a2] | 1199 | tmp.data=r; |
---|
[0bf7dff] | 1200 | if ((rtyp==IDHDL)||(rtyp==STRING_CMD)) |
---|
[c931a2] | 1201 | { |
---|
| 1202 | tmp.next=next; next=NULL; |
---|
[17681b] | 1203 | //if (rtyp==STRING_CMD) { omFree((ADDRESS)data); } |
---|
| 1204 | //data=NULL; |
---|
| 1205 | d=NULL; |
---|
[c931a2] | 1206 | CleanUp(); |
---|
| 1207 | memcpy(this,&tmp,sizeof(tmp)); |
---|
| 1208 | } |
---|
| 1209 | // and, remember, r is also the result... |
---|
| 1210 | else |
---|
| 1211 | { |
---|
| 1212 | // ??? |
---|
| 1213 | // here we still have a memory leak... |
---|
[82063a] | 1214 | // example: list L="123","456"; |
---|
| 1215 | // L[1][2]; |
---|
| 1216 | // therefore, it should never happen: |
---|
| 1217 | assume(0); |
---|
| 1218 | // but if it happens: here is the temporary fix: |
---|
| 1219 | // omMarkAsStaticAddr(r); |
---|
[c931a2] | 1220 | } |
---|
[0e1846] | 1221 | break; |
---|
| 1222 | } |
---|
| 1223 | case MATRIX_CMD: |
---|
| 1224 | { |
---|
[c931a2] | 1225 | if ((index<1) |
---|
| 1226 | ||(index>MATROWS((matrix)d)) |
---|
[0e1846] | 1227 | ||(e->next->start<1) |
---|
| 1228 | ||(e->next->start>MATCOLS((matrix)d))) |
---|
| 1229 | { |
---|
| 1230 | if (!errorreported) |
---|
[c931a2] | 1231 | Werror("wrong range[%d,%d] in intmat(%dx%d)", |
---|
| 1232 | index,e->next->start, |
---|
| 1233 | MATROWS((matrix)d),MATCOLS((matrix)d)); |
---|
[0e1846] | 1234 | } |
---|
| 1235 | else |
---|
[c931a2] | 1236 | r=(char *)MATELEM((matrix)d,index,e->next->start); |
---|
[0e1846] | 1237 | break; |
---|
| 1238 | } |
---|
[8db6c3] | 1239 | default: |
---|
[0e1846] | 1240 | { |
---|
[8db6c3] | 1241 | blackbox *b=NULL; |
---|
| 1242 | if (t>MAX_TOK) |
---|
| 1243 | { |
---|
| 1244 | b=getBlackboxStuff(t); |
---|
| 1245 | } |
---|
| 1246 | if ((t==LIST_CMD)||((b!=NULL)&&(BB_LIKE_LIST(b)))) |
---|
[0e1846] | 1247 | { |
---|
[8db6c3] | 1248 | lists l=(lists)d; |
---|
| 1249 | if ((0<index)&&(index<=l->nr+1)) |
---|
[82063a] | 1250 | { |
---|
[8db6c3] | 1251 | if ((e->next!=NULL) |
---|
| 1252 | && (l->m[index-1].rtyp==STRING_CMD)) |
---|
| 1253 | // string[..].Data() modifies sleftv, so let's do it ourself |
---|
[82063a] | 1254 | { |
---|
[8db6c3] | 1255 | char *dd=(char *)l->m[index-1].data; |
---|
| 1256 | int j=e->next->start-1; |
---|
| 1257 | r=(char *)omAllocBin(size_two_bin); |
---|
| 1258 | if ((j>=0) && (j<(int)strlen(dd))) |
---|
| 1259 | { |
---|
| 1260 | r[0]=*(dd+j); |
---|
| 1261 | r[1]='\0'; |
---|
| 1262 | } |
---|
| 1263 | else |
---|
| 1264 | { |
---|
| 1265 | r[0]='\0'; |
---|
| 1266 | } |
---|
[82063a] | 1267 | } |
---|
| 1268 | else |
---|
| 1269 | { |
---|
[8db6c3] | 1270 | Subexpr tmp=l->m[index-1].e; |
---|
| 1271 | l->m[index-1].e=e->next; |
---|
| 1272 | r=(char *)l->m[index-1].Data(); |
---|
| 1273 | e->next=l->m[index-1].e; |
---|
| 1274 | l->m[index-1].e=tmp; |
---|
[82063a] | 1275 | } |
---|
| 1276 | } |
---|
[8db6c3] | 1277 | else //if (!errorreported) |
---|
| 1278 | Werror("wrong range[%d] in list(%d)",index,l->nr+1); |
---|
[0e1846] | 1279 | } |
---|
[8db6c3] | 1280 | else |
---|
| 1281 | Werror("cannot index type %s(%d)",Tok2Cmdname(t),t); |
---|
[0e1846] | 1282 | break; |
---|
| 1283 | } |
---|
| 1284 | } |
---|
| 1285 | return r; |
---|
| 1286 | } |
---|
| 1287 | |
---|
| 1288 | attr * sleftv::Attribute() |
---|
| 1289 | { |
---|
[8c94ba] | 1290 | if (e==NULL) return &attribute; |
---|
[13e2da] | 1291 | if ((rtyp==LIST_CMD) |
---|
| 1292 | ||((rtyp==IDHDL)&&(IDTYP((idhdl)data)==LIST_CMD))) |
---|
[0e1846] | 1293 | { |
---|
| 1294 | leftv v=LData(); |
---|
| 1295 | return &(v->attribute); |
---|
| 1296 | } |
---|
| 1297 | return NULL; |
---|
| 1298 | } |
---|
| 1299 | |
---|
| 1300 | leftv sleftv::LData() |
---|
| 1301 | { |
---|
| 1302 | if (e!=NULL) |
---|
| 1303 | { |
---|
| 1304 | lists l=NULL; |
---|
| 1305 | |
---|
| 1306 | if (rtyp==LIST_CMD) |
---|
| 1307 | l=(lists)data; |
---|
[db34980] | 1308 | else if ((rtyp==IDHDL)&& (IDTYP((idhdl)data)==LIST_CMD)) |
---|
[0e1846] | 1309 | l=IDLIST((idhdl)data); |
---|
[db34980] | 1310 | else if (rtyp==ALIAS_CMD) |
---|
| 1311 | { |
---|
| 1312 | idhdl h=(idhdl)data; |
---|
| 1313 | l= (lists)(((idhdl)h->data.ustring)->data.ustring); |
---|
| 1314 | } |
---|
[0e1846] | 1315 | if (l!=NULL) |
---|
| 1316 | { |
---|
| 1317 | if ((0>=e->start)||(e->start>l->nr+1)) |
---|
| 1318 | return NULL; |
---|
| 1319 | if (e->next!=NULL) |
---|
| 1320 | { |
---|
| 1321 | l->m[e->start-1].e=e->next; |
---|
| 1322 | leftv r=l->m[e->start-1].LData(); |
---|
| 1323 | l->m[e->start-1].e=NULL; |
---|
| 1324 | return r; |
---|
| 1325 | } |
---|
| 1326 | return &(l->m[e->start-1]); |
---|
| 1327 | } |
---|
| 1328 | } |
---|
| 1329 | return this; |
---|
| 1330 | } |
---|
| 1331 | |
---|
| 1332 | leftv sleftv::LHdl() |
---|
| 1333 | { |
---|
| 1334 | if (e!=NULL) |
---|
| 1335 | { |
---|
[577d19f] | 1336 | lists l=NULL; |
---|
| 1337 | |
---|
| 1338 | if (rtyp==LIST_CMD) |
---|
| 1339 | l=(lists)data; |
---|
| 1340 | if ((rtyp==IDHDL)&& (IDTYP((idhdl)data)==LIST_CMD)) |
---|
| 1341 | l=IDLIST((idhdl)data); |
---|
| 1342 | if (l!=NULL) |
---|
| 1343 | { |
---|
| 1344 | if ((0>=e->start)||(e->start>l->nr+1)) |
---|
| 1345 | return NULL; |
---|
| 1346 | if (e->next!=NULL) |
---|
| 1347 | { |
---|
| 1348 | l->m[e->start-1].e=e->next; |
---|
| 1349 | leftv r=l->m[e->start-1].LHdl(); |
---|
| 1350 | l->m[e->start-1].e=NULL; |
---|
| 1351 | return r; |
---|
| 1352 | } |
---|
| 1353 | return &(l->m[e->start-1]); |
---|
| 1354 | } |
---|
[0e1846] | 1355 | } |
---|
| 1356 | return this; |
---|
| 1357 | } |
---|
| 1358 | |
---|
| 1359 | BOOLEAN assumeStdFlag(leftv h) |
---|
| 1360 | { |
---|
| 1361 | if ((h->e!=NULL)&&(h->LTyp()==LIST_CMD)) |
---|
| 1362 | { |
---|
| 1363 | return assumeStdFlag(h->LData()); |
---|
| 1364 | } |
---|
| 1365 | if (!hasFlag(h,FLAG_STD)) |
---|
| 1366 | { |
---|
| 1367 | if (!TEST_VERB_NSB) |
---|
[95197e] | 1368 | { |
---|
| 1369 | if (TEST_V_ALLWARN) |
---|
| 1370 | Warn("%s is no standard basis in >>%s<<",h->Name(),my_yylinebuf); |
---|
| 1371 | else |
---|
| 1372 | Warn("%s is no standard basis",h->Name()); |
---|
| 1373 | } |
---|
[0e1846] | 1374 | return FALSE; |
---|
| 1375 | } |
---|
| 1376 | return TRUE; |
---|
| 1377 | } |
---|
| 1378 | |
---|
| 1379 | /*2 |
---|
[c232af] | 1380 | * transforms a name (as an string created by omAlloc or omStrDup) |
---|
[0e1846] | 1381 | * into an expression (sleftv), deletes the string |
---|
| 1382 | * utility for grammar and iparith |
---|
| 1383 | */ |
---|
[85e68dd] | 1384 | void syMake(leftv v,const char * id, idhdl packhdl) |
---|
[0e1846] | 1385 | { |
---|
| 1386 | /* resolv an identifier: (to DEF_CMD, if siq>0) |
---|
| 1387 | * 1) reserved id: done by scanner |
---|
[8f1a54] | 1388 | * 2) `basering` / 'Current` |
---|
[0e1846] | 1389 | * 3) existing identifier, local |
---|
[4c9cf1] | 1390 | * 4) ringvar, ringpar, local ring |
---|
[0e1846] | 1391 | * 5) existing identifier, global |
---|
| 1392 | * 6) monom (resp. number), local ring: consisting of: |
---|
[4c9cf1] | 1393 | * 6') ringvar, ringpar,global ring |
---|
[0e1846] | 1394 | * 6'') monom (resp. number), local ring |
---|
| 1395 | * 7) monom (resp. number), non-local ring |
---|
| 1396 | * 8) basering |
---|
| 1397 | * 9) `_` |
---|
| 1398 | * 10) everything else is of type 0 |
---|
| 1399 | */ |
---|
| 1400 | #ifdef TEST |
---|
| 1401 | if ((*id<' ')||(*id>(char)126)) |
---|
| 1402 | { |
---|
| 1403 | Print("wrong id :%s:\n",id); |
---|
| 1404 | } |
---|
| 1405 | #endif |
---|
[0e051f] | 1406 | idhdl save_ring=currRingHdl; |
---|
[f2dff02] | 1407 | v->Init(); |
---|
[bd4cb92] | 1408 | if(packhdl != NULL) |
---|
[17ccde] | 1409 | { |
---|
| 1410 | // Print("setting req_packhdl to %s\n",IDID(packhdl)); |
---|
[bd4cb92] | 1411 | v->req_packhdl = IDPACKAGE(packhdl); |
---|
[ff7e529] | 1412 | } |
---|
[a3bc95e] | 1413 | else v->req_packhdl = currPack; |
---|
| 1414 | // if (v->req_packhdl!=basePack) |
---|
| 1415 | // Print("search %s in %s\n",id,v->req_packhdl->libname); |
---|
[610881] | 1416 | idhdl h=NULL; |
---|
[0e1846] | 1417 | #ifdef SIQ |
---|
| 1418 | if (siq<=0) |
---|
| 1419 | #endif |
---|
| 1420 | { |
---|
| 1421 | if (!isdigit(id[0])) |
---|
| 1422 | { |
---|
| 1423 | if (strcmp(id,"basering")==0) |
---|
| 1424 | { |
---|
| 1425 | if (currRingHdl!=NULL) |
---|
| 1426 | { |
---|
[95197e] | 1427 | if (id!=IDID(currRingHdl)) omFreeBinAddr((ADDRESS)id); |
---|
[610881] | 1428 | h=currRingHdl; |
---|
| 1429 | goto id_found; |
---|
[0e1846] | 1430 | } |
---|
| 1431 | else |
---|
| 1432 | { |
---|
| 1433 | v->name = id; |
---|
| 1434 | return; /* undefined */ |
---|
| 1435 | } |
---|
| 1436 | } |
---|
[610881] | 1437 | else if (strcmp(id,"Current")==0) |
---|
[8f1a54] | 1438 | { |
---|
| 1439 | if (currPackHdl!=NULL) |
---|
| 1440 | { |
---|
[95197e] | 1441 | omFreeBinAddr((ADDRESS)id); |
---|
[610881] | 1442 | h=currPackHdl; |
---|
| 1443 | goto id_found; |
---|
[8f1a54] | 1444 | } |
---|
| 1445 | else |
---|
| 1446 | { |
---|
| 1447 | v->name = id; |
---|
| 1448 | return; /* undefined */ |
---|
| 1449 | } |
---|
| 1450 | } |
---|
[bd4cb92] | 1451 | if(v->req_packhdl!=currPack) |
---|
| 1452 | { |
---|
| 1453 | h=v->req_packhdl->idroot->get(id,myynest); |
---|
| 1454 | } |
---|
| 1455 | else |
---|
[9c9fb9] | 1456 | h=ggetid(id); |
---|
[0e1846] | 1457 | /* 3) existing identifier, local */ |
---|
| 1458 | if ((h!=NULL) && (IDLEV(h)==myynest)) |
---|
| 1459 | { |
---|
[95197e] | 1460 | if (id!=IDID(h)) omFreeBinAddr((ADDRESS)id); /*assume strlen(id) <1000 */ |
---|
[610881] | 1461 | goto id_found; |
---|
[0e1846] | 1462 | } |
---|
| 1463 | } |
---|
[0e051f] | 1464 | if (yyInRingConstruction) |
---|
| 1465 | { |
---|
| 1466 | currRingHdl=NULL; |
---|
| 1467 | } |
---|
[0e1846] | 1468 | /* 4. local ring: ringvar */ |
---|
[b9e18a3] | 1469 | if ((currRingHdl!=NULL) && (IDLEV(currRingHdl)==myynest) |
---|
[0e051f] | 1470 | /*&& (!yyInRingConstruction)*/) |
---|
[0e1846] | 1471 | { |
---|
| 1472 | int vnr; |
---|
[8fa9ea] | 1473 | if ((vnr=r_IsRingVar(id, currRing))>=0) |
---|
[0e1846] | 1474 | { |
---|
| 1475 | poly p=pOne(); |
---|
| 1476 | pSetExp(p,vnr+1,1); |
---|
| 1477 | pSetm(p); |
---|
| 1478 | v->data = (void *)p; |
---|
| 1479 | v->name = id; |
---|
| 1480 | v->rtyp = POLY_CMD; |
---|
| 1481 | return; |
---|
| 1482 | } |
---|
[4c9cf1] | 1483 | if(rField_is_Extension(currRing)&&((vnr=r_IsRingVar(id, currRing->cf->extRing))>=0)) |
---|
| 1484 | { |
---|
| 1485 | BOOLEAN ok=FALSE; |
---|
| 1486 | poly p = pmInit(id,ok); |
---|
| 1487 | if (ok) |
---|
| 1488 | { |
---|
| 1489 | v->data = pGetCoeff(p); |
---|
| 1490 | pGetCoeff(p)=NULL; |
---|
| 1491 | pLmFree(p); |
---|
| 1492 | v->rtyp = NUMBER_CMD; |
---|
| 1493 | v->name = id; |
---|
| 1494 | return; |
---|
| 1495 | } |
---|
| 1496 | } |
---|
[0e1846] | 1497 | } |
---|
| 1498 | /* 5. existing identifier, global */ |
---|
| 1499 | if (h!=NULL) |
---|
| 1500 | { |
---|
[95197e] | 1501 | if (id!=IDID(h)) omFreeBinAddr((ADDRESS)id); /*assume strlen(id) <1000 */ |
---|
[610881] | 1502 | goto id_found; |
---|
[0e1846] | 1503 | } |
---|
| 1504 | /* 6. local ring: number/poly */ |
---|
| 1505 | if ((currRingHdl!=NULL) && (IDLEV(currRingHdl)==myynest)) |
---|
| 1506 | { |
---|
| 1507 | BOOLEAN ok=FALSE; |
---|
[0e051f] | 1508 | /*poly p = (!yyInRingConstruction) ? pmInit(id,ok) : (poly)NULL;*/ |
---|
| 1509 | poly p = pmInit(id,ok); |
---|
[0e1846] | 1510 | if (ok) |
---|
| 1511 | { |
---|
| 1512 | if (p==NULL) |
---|
| 1513 | { |
---|
| 1514 | v->data = (void *)nInit(0); |
---|
| 1515 | v->rtyp = NUMBER_CMD; |
---|
[875ec3] | 1516 | #ifdef HAVE_PLURAL |
---|
| 1517 | // in this case we may have monomials equal to 0 in p_Read |
---|
| 1518 | v->name = id; |
---|
| 1519 | #else |
---|
[95197e] | 1520 | omFreeBinAddr((ADDRESS)id); |
---|
[875ec3] | 1521 | #endif |
---|
[0e1846] | 1522 | } |
---|
[0e051f] | 1523 | else if (pIsConstant(p)) |
---|
[0e1846] | 1524 | { |
---|
| 1525 | v->data = pGetCoeff(p); |
---|
| 1526 | pGetCoeff(p)=NULL; |
---|
[512a2b] | 1527 | pLmFree(p); |
---|
[0e1846] | 1528 | v->rtyp = NUMBER_CMD; |
---|
| 1529 | v->name = id; |
---|
| 1530 | } |
---|
| 1531 | else |
---|
| 1532 | { |
---|
| 1533 | v->data = p; |
---|
| 1534 | v->rtyp = POLY_CMD; |
---|
| 1535 | v->name = id; |
---|
| 1536 | } |
---|
| 1537 | return; |
---|
| 1538 | } |
---|
| 1539 | } |
---|
| 1540 | /* 7. non-local ring: number/poly */ |
---|
| 1541 | { |
---|
| 1542 | BOOLEAN ok=FALSE; |
---|
[1aa559b] | 1543 | poly p = ((currRing!=NULL) /* ring required */ |
---|
| 1544 | && (currRingHdl!=NULL) |
---|
[0e051f] | 1545 | /*&& (!yyInRingConstruction) - not in decl */ |
---|
[1aa559b] | 1546 | && (IDLEV(currRingHdl)!=myynest)) /* already in case 4/6 */ |
---|
[0e1846] | 1547 | ? pmInit(id,ok) : (poly)NULL; |
---|
| 1548 | if (ok) |
---|
| 1549 | { |
---|
| 1550 | if (p==NULL) |
---|
| 1551 | { |
---|
| 1552 | v->data = (void *)nInit(0); |
---|
| 1553 | v->rtyp = NUMBER_CMD; |
---|
[875ec3] | 1554 | #ifdef HAVE_PLURAL |
---|
| 1555 | // in this case we may have monomials equal to 0 in p_Read |
---|
| 1556 | v->name = id; |
---|
| 1557 | #else |
---|
[95197e] | 1558 | omFreeBinAddr((ADDRESS)id); |
---|
[875ec3] | 1559 | #endif |
---|
[0e1846] | 1560 | } |
---|
| 1561 | else |
---|
| 1562 | if (pIsConstant(p)) |
---|
| 1563 | { |
---|
| 1564 | v->data = pGetCoeff(p); |
---|
| 1565 | pGetCoeff(p)=NULL; |
---|
[512a2b] | 1566 | pLmFree(p); |
---|
[0e1846] | 1567 | v->rtyp = NUMBER_CMD; |
---|
| 1568 | v->name = id; |
---|
| 1569 | } |
---|
| 1570 | else |
---|
| 1571 | { |
---|
| 1572 | v->data = p; |
---|
| 1573 | v->rtyp = POLY_CMD; |
---|
| 1574 | v->name = id; |
---|
| 1575 | } |
---|
| 1576 | return; |
---|
| 1577 | } |
---|
| 1578 | } |
---|
| 1579 | /* 8. basering ? */ |
---|
| 1580 | if ((myynest>1)&&(currRingHdl!=NULL)) |
---|
| 1581 | { |
---|
| 1582 | if (strcmp(id,IDID(currRingHdl))==0) |
---|
| 1583 | { |
---|
[95197e] | 1584 | if (IDID(currRingHdl)!=id) omFreeBinAddr((ADDRESS)id); /*assume strlen |
---|
| 1585 | (id) <1000 */ |
---|
[610881] | 1586 | h=currRingHdl; |
---|
| 1587 | goto id_found; |
---|
[0e1846] | 1588 | } |
---|
| 1589 | } |
---|
[8f1a54] | 1590 | if((v->req_packhdl!=basePack) && (v->req_packhdl==currPack)) |
---|
[a3bc95e] | 1591 | { |
---|
[5bce7c] | 1592 | h=basePack->idroot->get(id,myynest); |
---|
| 1593 | if (h!=NULL) |
---|
[ff7e529] | 1594 | { |
---|
[95197e] | 1595 | if (id!=IDID(h)) omFreeBinAddr((ADDRESS)id); /*assume strlen(id) <1000 */ |
---|
[5bce7c] | 1596 | v->req_packhdl=basePack; |
---|
[610881] | 1597 | goto id_found; |
---|
[5bce7c] | 1598 | } |
---|
[a3bc95e] | 1599 | } |
---|
[0e1846] | 1600 | } |
---|
| 1601 | #ifdef SIQ |
---|
| 1602 | else |
---|
| 1603 | v->rtyp=DEF_CMD; |
---|
| 1604 | #endif |
---|
| 1605 | /* 9: _ */ |
---|
| 1606 | if (strcmp(id,"_")==0) |
---|
| 1607 | { |
---|
[95197e] | 1608 | omFreeBinAddr((ADDRESS)id); |
---|
[0e1846] | 1609 | v->Copy(&sLastPrinted); |
---|
| 1610 | } |
---|
| 1611 | else |
---|
| 1612 | { |
---|
| 1613 | /* 10: everything else */ |
---|
| 1614 | /* v->rtyp = UNKNOWN;*/ |
---|
| 1615 | v->name = id; |
---|
| 1616 | } |
---|
[0e051f] | 1617 | currRingHdl=save_ring; |
---|
[610881] | 1618 | return; |
---|
| 1619 | id_found: // we have an id (in h) found, to set the data in from h |
---|
[ea947e] | 1620 | if (IDTYP(h)!=ALIAS_CMD) |
---|
| 1621 | { |
---|
| 1622 | v->rtyp = IDHDL; |
---|
| 1623 | v->flag = IDFLAG(h); |
---|
| 1624 | v->attribute=IDATTR(h); |
---|
| 1625 | } |
---|
| 1626 | else |
---|
| 1627 | { |
---|
| 1628 | v->rtyp = ALIAS_CMD; |
---|
| 1629 | } |
---|
[610881] | 1630 | v->name = IDID(h); |
---|
[ea947e] | 1631 | v->data = (char *)h; |
---|
[0e051f] | 1632 | currRingHdl=save_ring; |
---|
[0e1846] | 1633 | } |
---|
| 1634 | |
---|
| 1635 | int sleftv::Eval() |
---|
| 1636 | { |
---|
| 1637 | BOOLEAN nok=FALSE; |
---|
| 1638 | leftv nn=next; |
---|
| 1639 | next=NULL; |
---|
| 1640 | if(rtyp==IDHDL) |
---|
| 1641 | { |
---|
| 1642 | int t=Typ(); |
---|
[2a2c07] | 1643 | if (t!=PROC_CMD) |
---|
| 1644 | { |
---|
| 1645 | void *d=CopyD(t); |
---|
| 1646 | data=d; |
---|
| 1647 | rtyp=t; |
---|
| 1648 | name=NULL; |
---|
| 1649 | e=NULL; |
---|
| 1650 | } |
---|
[0e1846] | 1651 | } |
---|
| 1652 | else if (rtyp==COMMAND) |
---|
| 1653 | { |
---|
| 1654 | command d=(command)data; |
---|
| 1655 | if(d->op==PROC_CMD) //assume d->argc==2 |
---|
| 1656 | { |
---|
| 1657 | char *what=(char *)(d->arg1.Data()); |
---|
| 1658 | idhdl h=ggetid(what); |
---|
| 1659 | if((h!=NULL)&&(IDTYP(h)==PROC_CMD)) |
---|
| 1660 | { |
---|
| 1661 | nok=d->arg2.Eval(); |
---|
| 1662 | if(!nok) |
---|
| 1663 | { |
---|
[5edb77] | 1664 | nok=iiMake_proc(h,req_packhdl,&d->arg2); |
---|
| 1665 | if (!nok) |
---|
| 1666 | memcpy(this,&iiRETURNEXPR,sizeof(sleftv)); |
---|
[0e1846] | 1667 | } |
---|
| 1668 | } |
---|
| 1669 | else nok=TRUE; |
---|
| 1670 | } |
---|
| 1671 | else if (d->op=='=') //assume d->argc==2 |
---|
| 1672 | { |
---|
[faeb82] | 1673 | if ((d->arg1.rtyp!=IDHDL)&&(d->arg1.rtyp!=DEF_CMD)) |
---|
[0e1846] | 1674 | { |
---|
[faeb82] | 1675 | nok=d->arg1.Eval(); |
---|
| 1676 | } |
---|
| 1677 | if (!nok) |
---|
| 1678 | { |
---|
[85e68dd] | 1679 | const char *n=d->arg1.name; |
---|
[fff984] | 1680 | nok=(n == NULL) || d->arg2.Eval(); |
---|
[0e1846] | 1681 | if (!nok) |
---|
| 1682 | { |
---|
[d754b7] | 1683 | int save_typ=d->arg1.rtyp; |
---|
[85e68dd] | 1684 | omCheckAddr((ADDRESS)n); |
---|
[faeb82] | 1685 | if (d->arg1.rtyp!=IDHDL) |
---|
[fff984] | 1686 | syMake(&d->arg1,n); |
---|
[85e68dd] | 1687 | omCheckAddr((ADDRESS)d->arg1.name); |
---|
[0e1846] | 1688 | if (d->arg1.rtyp==IDHDL) |
---|
| 1689 | { |
---|
[c232af] | 1690 | n=omStrDup(IDID((idhdl)d->arg1.data)); |
---|
[0e1846] | 1691 | killhdl((idhdl)d->arg1.data); |
---|
[a3bc95e] | 1692 | d->arg1.Init(); |
---|
[f077d2] | 1693 | //d->arg1.data=NULL; |
---|
[0e1846] | 1694 | d->arg1.name=n; |
---|
| 1695 | } |
---|
[f077d2] | 1696 | d->arg1.rtyp=DEF_CMD; |
---|
[0e1846] | 1697 | sleftv t; |
---|
[d754b7] | 1698 | if(save_typ!=PROC_CMD) save_typ=d->arg2.rtyp; |
---|
[df5fc1] | 1699 | if (::RingDependend(d->arg2.rtyp)) |
---|
[d754b7] | 1700 | nok=iiDeclCommand(&t,&d->arg1,0,save_typ,&currRing->idroot); |
---|
[0e1846] | 1701 | else |
---|
[46d09b] | 1702 | nok=iiDeclCommand(&t,&d->arg1,0,save_typ,&IDROOT); |
---|
[0e1846] | 1703 | memcpy(&d->arg1,&t,sizeof(sleftv)); |
---|
[85e68dd] | 1704 | omCheckAddr((ADDRESS)d->arg1.name); |
---|
[0e1846] | 1705 | nok=nok||iiAssign(&d->arg1,&d->arg2); |
---|
[9d72fe] | 1706 | omCheckIf(d->arg1.name != NULL, // OB: ???? |
---|
[85e68dd] | 1707 | omCheckAddr((ADDRESS)d->arg1.name)); |
---|
[0e1846] | 1708 | if (!nok) |
---|
| 1709 | { |
---|
| 1710 | memset(&d->arg1,0,sizeof(sleftv)); |
---|
[4f568c] | 1711 | this->CleanUp(); |
---|
| 1712 | rtyp=NONE; |
---|
[0e1846] | 1713 | } |
---|
| 1714 | } |
---|
| 1715 | } |
---|
| 1716 | else nok=TRUE; |
---|
| 1717 | } |
---|
[91ef34] | 1718 | else |
---|
[0e1846] | 1719 | { |
---|
[91ef34] | 1720 | int toktype=iiTokType(d->op); |
---|
| 1721 | if ((toktype==CMD_M) |
---|
| 1722 | ||( toktype==ROOT_DECL_LIST) |
---|
[1f03aba] | 1723 | ||( toktype==RING_DECL_LIST)) |
---|
[91ef34] | 1724 | { |
---|
| 1725 | if (d->argc <=3) |
---|
| 1726 | { |
---|
| 1727 | if (d->argc>=1) nok=d->arg1.Eval(); |
---|
| 1728 | if ((!nok) && (d->argc>=2)) |
---|
| 1729 | { nok=d->arg2.Eval(); d->arg1.next=&d->arg2; } |
---|
| 1730 | if ((!nok) && (d->argc==3)) |
---|
| 1731 | { nok=d->arg3.Eval(); d->arg2.next=&d->arg3; } |
---|
| 1732 | if (d->argc==0) |
---|
| 1733 | nok=nok||iiExprArithM(this,NULL,d->op); |
---|
| 1734 | else |
---|
| 1735 | nok=nok||iiExprArithM(this,&d->arg1,d->op); |
---|
| 1736 | d->arg1.next=NULL; |
---|
| 1737 | d->arg2.next=NULL; |
---|
| 1738 | d->arg3.next=NULL; |
---|
| 1739 | } |
---|
| 1740 | else |
---|
| 1741 | { |
---|
| 1742 | nok=d->arg1.Eval(); |
---|
| 1743 | nok=nok||iiExprArithM(this,&d->arg1,d->op); |
---|
| 1744 | } |
---|
| 1745 | } |
---|
| 1746 | else if (d->argc==1) |
---|
| 1747 | { |
---|
| 1748 | nok=d->arg1.Eval(); |
---|
| 1749 | nok=nok||iiExprArith1(this,&d->arg1,d->op); |
---|
| 1750 | } |
---|
| 1751 | else if(d->argc==2) |
---|
| 1752 | { |
---|
| 1753 | nok=d->arg1.Eval(); |
---|
| 1754 | nok=nok||d->arg2.Eval(); |
---|
| 1755 | nok=nok||iiExprArith2(this,&d->arg1,d->op,&d->arg2); |
---|
| 1756 | } |
---|
| 1757 | else if(d->argc==3) |
---|
| 1758 | { |
---|
| 1759 | nok=d->arg1.Eval(); |
---|
| 1760 | nok=nok||d->arg2.Eval(); |
---|
| 1761 | nok=nok||d->arg3.Eval(); |
---|
| 1762 | nok=nok||iiExprArith3(this,d->op,&d->arg1,&d->arg2,&d->arg3); |
---|
| 1763 | } |
---|
| 1764 | else if(d->argc!=0) |
---|
| 1765 | { |
---|
| 1766 | nok=d->arg1.Eval(); |
---|
| 1767 | nok=nok||iiExprArithM(this,&d->arg1,d->op); |
---|
| 1768 | } |
---|
| 1769 | else // d->argc == 0 |
---|
| 1770 | { |
---|
| 1771 | nok = iiExprArithM(this, NULL, d->op); |
---|
| 1772 | } |
---|
[0e1846] | 1773 | } |
---|
| 1774 | } |
---|
| 1775 | else if (((rtyp==0)||(rtyp==DEF_CMD)) |
---|
| 1776 | &&(name!=NULL)) |
---|
| 1777 | { |
---|
| 1778 | syMake(this,name); |
---|
| 1779 | } |
---|
| 1780 | #ifdef MDEBUG |
---|
| 1781 | switch(Typ()) |
---|
| 1782 | { |
---|
| 1783 | case NUMBER_CMD: |
---|
| 1784 | #ifdef LDEBUG |
---|
| 1785 | nTest((number)Data()); |
---|
[84a6a7] | 1786 | #endif |
---|
| 1787 | break; |
---|
| 1788 | case BIGINT_CMD: |
---|
| 1789 | #ifdef LDEBUG |
---|
[4ff07f] | 1790 | n_Test((number)Data(),coeffs_BIGINT); |
---|
[0e1846] | 1791 | #endif |
---|
| 1792 | break; |
---|
| 1793 | case POLY_CMD: |
---|
| 1794 | pTest((poly)Data()); |
---|
| 1795 | break; |
---|
| 1796 | case IDEAL_CMD: |
---|
| 1797 | case MODUL_CMD: |
---|
| 1798 | case MATRIX_CMD: |
---|
| 1799 | { |
---|
| 1800 | ideal id=(ideal)Data(); |
---|
[c232af] | 1801 | omCheckAddrSize(id,sizeof(*id)); |
---|
[0e1846] | 1802 | int i=id->ncols*id->nrows-1; |
---|
| 1803 | for(;i>=0;i--) pTest(id->m[i]); |
---|
| 1804 | } |
---|
| 1805 | break; |
---|
| 1806 | } |
---|
| 1807 | #endif |
---|
| 1808 | if (nn!=NULL) nok=nok||nn->Eval(); |
---|
| 1809 | next=nn; |
---|
| 1810 | return nok; |
---|
| 1811 | } |
---|
| 1812 | |
---|
[85e68dd] | 1813 | const char *iiSleftv2name(leftv v) |
---|
[c12cf96] | 1814 | { |
---|
| 1815 | return(v->name); |
---|
| 1816 | } |
---|
[f71681] | 1817 | |
---|
| 1818 | void * sattr::CopyA() |
---|
| 1819 | { |
---|
| 1820 | omCheckAddrSize(this,sizeof(sattr)); |
---|
| 1821 | return s_internalCopy(atyp,data); |
---|
| 1822 | } |
---|
| 1823 | |
---|