# Changeset 8e682d in git

Ignore:
Timestamp:
Dec 3, 2014, 5:18:10 PM (10 years ago)
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b52fc4b2495505785981d640dcf7eb3e456778ef')
Children:
8af63a71de0015f0413204e1a7a8620a8d37cc7a
Parents:
41c9e458d3df4f560848651c9fc0710d72f703dc
Message:
`---`
Location:
Singular
Files:
5 edited

Unmodified
Removed
• ## Singular/LIB/grwalk.lib

 r41c9e4 } proc fwalk(ideal Go, int reduction, list #) proc fwalk(ideal Go, int reduction, int printout, list #) "SYNTAX: fwalk(ideal i,int reductioin); fwalk(ideal i, int reduction intvec v, intvec w); ideal G = fetch(xR, Go); G = system("Mfwalk", G, curr_weight, target_weight, reduction); G = system("Mfwalk", G, curr_weight, target_weight, reduction, printout); setring xR; ideal I = y3+xyz+y2z+xz3, 3+xy+x2y+y2z; int reduction = 1; fwalk(I,reduction); int printout = 1; fwalk(I,reduction,printout); }
• ## Singular/LIB/rwalk.lib

 r41c9e4 rwalk(ideal,int,int[,intvec,intvec]);   standard basis of ideal via Random Walk algorithm rwalk(ideal,int[,intvec,intvec]);       standard basis of ideal via Random Perturbation Walk algorithm frwalk(ideal,int[,intvec,intvec]);      standard basis of ideal via Random Fractal Walk algorithm frandwalk(ideal,int[,intvec,intvec]);      standard basis of ideal via Random Fractal Walk algorithm "; * Fractal Walk with random element * ************************************/ proc frandwalk(ideal Go, int radius, list #) "SYNTAX: frwalk(ideal i, int radius); frwalk(ideal i, int radius, intvec v, intvec w); proc frandwalk(ideal Go, int radius, int reduction, int printout, list #) "SYNTAX: frwalk(ideal i, int radius, int reduction, int printout); frwalk(ideal i, int radius, int reduction, int printout, intvec v, intvec w); TYPE:    ideal PURPOSE: compute the standard basis of the ideal, calculated via ideal G = fetch(xR, Go); int pert_deg = 2; G = system("Mfrwalk", G, curr_weight, target_weight, radius); G = system("Mfrwalk", G, curr_weight, target_weight, radius, reduction, printout); setring xR; ideal I = y3+xyz+y2z+xz3, 3+xy+x2y+y2z; int reduction = 0; frandwalk(I,2,0); } frandwalk(I,2,0,1); }
• ## Singular/extra.cc

 r41c9e4 if (strcmp(sys_cmd, "Mfwalk") == 0) { const short t[]={4,IDEAL_CMD,INTVEC_CMD,INTVEC_CMD,INT_CMD}; const short t[]={5,IDEAL_CMD,INTVEC_CMD,INTVEC_CMD,INT_CMD,INT_CMD}; if (!iiCheckTypes(h,t,1)) return TRUE; if (((intvec*) h->next->Data())->length() != currRing->N && Werror("system(\"Mfwalk\" ...) intvecs not of length %d\n", currRing->N); return TRUE; } ideal arg1 = (ideal) h->Data(); intvec* arg2 = (intvec*) h->next->Data(); intvec* arg3 = (intvec*) h->next->next->Data(); int arg4 = (int)(long) h->next->next->next->Data(); ideal result = (ideal) Mfwalk(arg1, arg2, arg3, arg4); res->rtyp = IDEAL_CMD; res->data =  result; return FALSE; } else #endif /*==================== Mfrwalk =================*/ #ifdef HAVE_WALK if (strcmp(sys_cmd, "Mfrwalk") == 0) { const short t[]={5,IDEAL_CMD,INTVEC_CMD,INTVEC_CMD,INT_CMD,INT_CMD}; if (!iiCheckTypes(h,t,1)) return TRUE; if (((intvec*) h->next->Data())->length() != currRing->N && ((intvec*) h->next->next->Data())->length() != currRing->N) { Werror("system(\"Mfrwalk\" ...) intvecs not of length %d\n",currRing->N); return TRUE; } int arg4 = (int)(long) h->next->next->next->Data(); int arg5 = (int)(long) h->next->next->next->next->Data(); ideal result = (ideal) Mfrwalk(arg1, arg2, arg3, arg4, arg5); ideal result = (ideal) Mfwalk(arg1, arg2, arg3, arg4, arg5); res->rtyp = IDEAL_CMD; res->data =  result; return FALSE; } else #endif /*==================== Mfrwalk =================*/ #ifdef HAVE_WALK if (strcmp(sys_cmd, "Mfrwalk") == 0) { const short t[]={6,IDEAL_CMD,INTVEC_CMD,INTVEC_CMD,INT_CMD,INT_CMD,INT_CMD}; if (!iiCheckTypes(h,t,1)) return TRUE; /* if (((intvec*) h->next->Data())->length() != currRing->N && ((intvec*) h->next->next->Data())->length() != currRing->N) { Werror("system(\"Mfrwalk\" ...) intvecs not of length %d\n",currRing->N); return TRUE; } */ if((((intvec*) h->next->Data())->length() != currRing->N && ((intvec*) h->next->next->Data())->length() != currRing->N ) && (((intvec*) h->next->Data())->length() != (currRing->N)*(currRing->N) && ((intvec*) h->next->next->Data())->length() != (currRing->N)*(currRing->N) )) { Werror("system(\"Mfrwalk\" ...) intvecs not of length %d or %d\n", currRing->N,(currRing->N)*(currRing->N)); return TRUE; } ideal arg1 = (ideal) h->Data(); intvec* arg2 = (intvec*) h->next->Data(); intvec* arg3 = (intvec*) h->next->next->Data(); int arg4 = (int)(long) h->next->next->next->Data(); int arg5 = (int)(long) h->next->next->next->next->Data(); int arg6 = (int)(long) h->next->next->next->next->next->Data(); ideal result = (ideal) Mfrwalk(arg1, arg2, arg3, arg4, arg5, arg6); res->rtyp = IDEAL_CMD; res->data =  result;