My Project
Loading...
Searching...
No Matches
walk_ip.cc
Go to the documentation of this file.
1/****************************************
2* Computer Algebra System SINGULAR *
3****************************************/
4/*
5* ABSTRACT: frwalk: interpreter link
6*/
7
8#include "kernel/mod2.h"
9#include "Singular/tok.h"
10#include "misc/options.h"
11#include "Singular/ipid.h"
12#include "misc/intvec.h"
13#include "omalloc/omalloc.h"
14#include "kernel/polys.h"
15#include "kernel/ideals.h"
16#include "polys/matpol.h"
19#include "Singular/subexpr.h"
21#include "kernel/GBEngine/syz.h"
22#include "coeffs/numbers.h"
23#include "Singular/lists.h"
24#include "Singular/attrib.h"
25#include "Singular/ipconv.h"
28#include "polys/weight.h"
32//#include "kernel/gnumpfl.h"
33//#include "kernel/mpr_base.h"
34//#include "kernel/ffields.h"
35#include "polys/clapsing.h"
37#include "Singular/ipshell.h"
41#include "polys/prCopy.h"
42
43///////////////////////////////////////////////////////////////////
44//walkProc
45///////////////////////////////////////////////////////////////////
46//Description: The main function for the Walk-Algorithm. Checks the
47//input-data, and calls walk64 (see walkMain.cc). Returns the new
48//groebner basis or something else if an error occoured.
49///////////////////////////////////////////////////////////////////
50//Uses: omAlloc0,walkConsistency,rGetGlobalOrderWeightVec,
51//omFreeSize,sizeof,IDIDEAL,walk64,rSetHdl,idrMoveR,Werror,idInit
52///////////////////////////////////////////////////////////////////
53
56{
57 WalkState state = WalkOk;
60 si_opt_1 &= (~Sy_bit(OPT_REDSB)); //make sure option noredSB is set
61
64 idhdl sourceRingHdl = (idhdl)first->data;
65 ring sourceRing = IDRING(sourceRingHdl);
66 rChangeCurrRing( sourceRing );
67
68 if(state==WalkOk)
69 {
70 int * vperm = (int *)omAlloc0( (currRing->N+1)*sizeof( int ) );
71 state= walkConsistency( sourceRing, destRing, vperm );
72 omFreeSize( (ADDRESS)vperm, (currRing->N+1)*sizeof(int) );
73 }
74
77
80 if ( state == WalkOk )
81 {
82 idhdl ih = currRing->idroot->get( second->Name(), myynest );
83 if ( (ih != NULL) && (IDTYP(ih)==IDEAL_CMD) )
84 {
88 }
89 }
90 else
91 {
92 state=WalkNoIdeal;
93 }
94 }
95
96 if ( state == WalkOk )
97 {
98 // Now the settings are compatible with Walk
101 }
102
103 SI_RESTORE_OPT(save1,save2);//making sure options are as before function call
104
107
108 switch (state) {
109 case WalkOk:
111 break;
112
114 Werror("ring %s and current ring are incompatible\n",
115 first->Name() );
117 break;
118
120 WerrorS( "Order of basering not allowed,\n must be a combination of a,A,lp,dp,Dp,wp,Wp,M and C.\n");
122 break;
123
125 Werror( "Order of %s not allowed,\n must be a combination of a,A,lp,dp,Dp,wp,Wp,M and C.\n",first->Name());
128 break;
129
130 case WalkNoIdeal:
131 Werror( "Can't find ideal %s in ring %s.\n",
132 second->Name(), first->Name() );
134 break;
135
137 WerrorS( "Overflow occurred.\n");
139 break;
140
141 default:
143 }
144
145 return destIdeal;
146}
147
148///////////////////////////////////////////////////////////////////
149//fractalWalkProc
150///////////////////////////////////////////////////////////////////
151//Description: The main function for the Fractalwalk-Algorithm.
152//Responsible for contact between user and walk64. Checks the
153//input-data, and calls fractalWalk64. Returns the new groebner
154//basis or something else if an error occurred.
155///////////////////////////////////////////////////////////////////
156//Uses: omAlloc0,fractalWalkConsistency,omFreeSize,sizeof,IDIDEAL,
157//fractalWalk64,rSetHdl,idrMoveR,Werror,idInit
158///////////////////////////////////////////////////////////////////
159
160ideal
162{
163
164 //unperturbedStartVectorStrategy SHOULD BE SET BY THE USER THROUGH
165 //A THIRD ARGUMENT. TRUE MEANS THAT THE UNPERTURBED START
166 //VECTOR STRATEGY IS USED AND FALSE THAT THE START VECTOR IS
167 //MAXIMALLY PERTURBED
168
170
171 WalkState state = WalkOk;
174 si_opt_1 &= (~Sy_bit(OPT_REDSB)); //make sure option noredSB is set
175
178 idhdl sourceRingHdl = (idhdl)first->data;
180 ring sourceRing = currRing;
181
182 int * vperm = (int *)omAlloc0( (currRing->N+1)*sizeof( int ) );
183 state= fractalWalkConsistency( sourceRing, destRing, vperm );
184 omFreeSize( (ADDRESS)vperm, (currRing->N+1)*sizeof(int) );
185
188 if ( state == WalkOk ) {
189 idhdl ih = currRing->idroot->get( second->Name(), myynest );
190 if ( (ih != NULL) && (IDTYP(ih)==IDEAL_CMD) ) {
192 if(hasFlag((leftv)ih,FLAG_STD)){
194 }
195 }
196 else {
197 state=WalkNoIdeal;
198 }
199 }
200
201 if ( state == WalkOk ) {
202 // Now the settings are compatible with Walk
206 }
207
208 SI_RESTORE_OPT(save1,save2);//making sure options are as before functiocall
209
210 if ( state == WalkOk )
211 {
215 }
216
217
218 switch (state) {
219
220 case WalkOk:
222 return(destIdeal);
223 break;
224
226 Werror( "ring %s and current ring are incompatible\n",
227 first->Name() );
230 return destIdeal;
231 break;
232
234 Werror( "Order of basering not allowed,\n must be a combination of lp,dp,Dp,wp,Wp and C or just M.\n");
237 return destIdeal;
238 break;
239
241 Werror( "Order of %s not allowed,\n must be a combination of lp,dp,Dp,wp,Wp and C or just M.\n",
242 first->Name());
245 return destIdeal;
246 break;
247
248 case WalkNoIdeal:
249 Werror( "Can't find ideal %s in ring %s.\n",
250 second->Name(), first->Name() );
253 return destIdeal;
254 break;
255
257 Werror( "Overflow occurred in ring %s.\n", first->Name() );
260 return destIdeal;
261 break;
262
263 default:
265 destIdeal= idInit(1,1);
266 return destIdeal;
267 }
268
269
270 return NULL;
271}
272
273
274///////////////////////////////////////////////////////////////////
275//getiv64
276///////////////////////////////////////////////////////////////////
277//Description: retrieves the int64vec from input list l
278///////////////////////////////////////////////////////////////////
279//Assumes: that the first entry of l is an int64vec
280///////////////////////////////////////////////////////////////////
281//Uses: none
282///////////////////////////////////////////////////////////////////
283
285{
286 return (int64vec*)(l->m[0].data);
287}
288
289///////////////////////////////////////////////////////////////////
290
291
292///////////////////////////////////////////////////////////////////
293//getint64
294///////////////////////////////////////////////////////////////////
295//Description: retrieves the int64 from input list l
296///////////////////////////////////////////////////////////////////
297//Assumes: that the second entry of l is an int64
298///////////////////////////////////////////////////////////////////
299//Uses: none
300///////////////////////////////////////////////////////////////////
301
302// not used, bad impl.
303//int64 getint64(lists l)
304//{
305// return (int64)(long)(l->m[1].data);
306//}
307
308///////////////////////////////////////////////////////////////////
int BOOLEAN
Definition auxiliary.h:87
#define TRUE
Definition auxiliary.h:100
#define FALSE
Definition auxiliary.h:96
int l
Definition cfEzgcd.cc:100
Definition idrec.h:35
Class used for (list of) interpreter objects.
Definition subexpr.h:83
const char * Name()
Definition subexpr.h:120
void * data
Definition subexpr.h:88
Definition lists.h:24
void WerrorS(const char *s)
Definition feFopen.cc:24
VAR int myynest
Definition febase.cc:41
@ IDEAL_CMD
Definition grammar.cc:284
ideal idCopy(ideal A)
Definition ideals.h:60
#define hasFlag(A, F)
Definition ipid.h:112
#define IDIDEAL(a)
Definition ipid.h:133
#define IDRING(a)
Definition ipid.h:127
#define IDTYP(a)
Definition ipid.h:119
#define FLAG_STD
Definition ipid.h:106
void rSetHdl(idhdl h)
Definition ipshell.cc:5128
#define omFreeSize(addr, size)
#define omAlloc0(size)
#define NULL
Definition omList.c:12
VAR unsigned si_opt_1
Definition options.c:5
#define SI_SAVE_OPT(A, B)
Definition options.h:20
#define OPT_REDSB
Definition options.h:76
#define SI_RESTORE_OPT(A, B)
Definition options.h:23
void rChangeCurrRing(ring r)
Definition polys.cc:15
VAR ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
Definition polys.cc:13
Compatibility layer for legacy polynomial operations (over currRing)
ideal idrMoveR(ideal &id, ring src_r, ring dest_r)
Definition prCopy.cc:248
void Werror(const char *fmt,...)
Definition reporter.cc:189
idrec * idhdl
Definition ring.h:21
ideal idInit(int idsize, int rank)
initialise an ideal / module
#define BITSET
Definition structs.h:16
WalkState fractalWalk64(ideal sourceIdeal, ring destRing, ideal &destIdeal, BOOLEAN sourceIsSB, BOOLEAN unperturbedStartVectorStrategy)
Definition walkMain.cc:614
WalkState walk64(ideal I, int64vec *currw64, ring destRing, int64vec *destVec64, ideal &destIdeal, BOOLEAN sourceIsSB)
Definition walkMain.cc:222
WalkState
Definition walkMain.h:7
@ WalkIncompatibleDestRing
Definition walkMain.h:28
@ WalkIncompatibleRings
Definition walkMain.h:9
@ WalkOverFlowError
Definition walkMain.h:11
@ WalkOk
Definition walkMain.h:30
@ WalkNoIdeal
Definition walkMain.h:8
@ WalkIncompatibleSourceRing
Definition walkMain.h:29
WalkState fractalWalkConsistency(ring sring, ring dring, int *vperm)
Definition walkProc.cc:220
WalkState walkConsistency(ring sring, ring dring, int *vperm)
Definition walkProc.cc:53
ideal sortRedSB(ideal G)
int64vec * rGetGlobalOrderWeightVec(ring r)
int64vec * getiv64(lists l)
Definition walk_ip.cc:284
ideal fractalWalkProc(leftv first, leftv second)
Definition walk_ip.cc:161
ideal walkProc(leftv first, leftv second)
Definition walk_ip.cc:55