source: git/kernel/GBEngine/syz.h @ e21795

spielwiese
Last change on this file since e21795 was 9b88e6, checked in by Hans Schoenemann <hannes@…>, 8 years ago
format: (long)0 -> 0L etc.
  • Property mode set to 100644
File size: 4.7 KB
Line 
1#ifndef SYZ_H
2#define SYZ_H
3/****************************************
4*  Computer Algebra System SINGULAR     *
5****************************************/
6/*
7* ABSTRACT: Resolutions
8*/
9#include <misc/mylimits.h>
10#include <kernel/structs.h>
11#include <polys/monomials/ring.h>
12#include <kernel/ideals.h>
13
14// Logarithm of estimate of maximal number of new components
15#define SYZ_SHIFT_MAX_NEW_COMP_ESTIMATE 8
16// Logarithm of "distance" between a new component and prev component
17#define SYZ_SHIFT_BASE_LOG (BIT_SIZEOF_LONG - 1 - SYZ_SHIFT_MAX_NEW_COMP_ESTIMATE)
18#define SYZ_SHIFT_BASE (1L << SYZ_SHIFT_BASE_LOG)
19struct sSObject{
20                 poly  p;
21                 poly  p1,p2; /*- the pair p comes from -*/
22                 poly  lcm;   /*- the lcm of p1,p2 -*/
23                 poly  syz;   /*- the syzygy associated to p1,p2 -*/
24                 int   ind1,ind2; /*- the indeces of p1,p2 -*/
25                 poly  isNotMinimal;
26                 int   syzind;
27                 int   order;
28                 int   length;
29                 int   reference;
30               };
31typedef struct sSObject SObject;
32typedef SObject * SSet;
33typedef SSet * SRes;
34
35class ssyStrategy;
36typedef ssyStrategy * syStrategy;
37class ssyStrategy{
38  public:
39  int ** truecomponents;
40  long** ShiftedComponents;
41  int ** backcomponents;
42  int ** Howmuch;
43  int ** Firstelem;
44  int ** elemLength;
45  intvec ** weights;
46  intvec ** hilb_coeffs;
47  resolvente res;              //polynomial data for internal use only
48  resolvente orderedRes;       //polynomial data for internal use only
49  SRes resPairs;               //polynomial data for internal use only
50  intvec * Tl;
51  intvec * resolution;
52  intvec * cw;
53  intvec * betti;
54  kBucket_pt bucket;
55  kBucket_pt syz_bucket;
56  ring syRing;
57  resolvente fullres;
58  resolvente minres;
59  unsigned long ** sev;
60  int length;
61  int regularity;
62  short list_length;
63  short references;
64};
65
66void sySchreyersSyzygiesM(polyset F,int Fmax,polyset* Shdl,int* Smax,
67   BOOLEAN noSort);
68
69void sySchreyersSyzygiesB(polyset F,int Fmax,polyset* Shdl,int* Smax,
70   BOOLEAN noSort);
71
72resolvente sySchreyerResolvente(ideal arg, int maxlength, int * length,
73   BOOLEAN isMonomial=FALSE, BOOLEAN notReplace=FALSE);
74
75syStrategy sySchreyer(ideal arg, int maxlength);
76
77resolvente syResolvente(ideal arg, int maxlength, int * length,
78                        intvec *** weights, BOOLEAN minim);
79
80syStrategy syResolution(ideal arg, int maxlength,intvec * w, BOOLEAN minim);
81
82void syMinimizeResolvente(resolvente res, int length, int first);
83
84intvec * syBetti(resolvente res,int length, int * regularity,
85         intvec* weights=NULL,BOOLEAN tomin=TRUE, int * row_shift=NULL);
86
87ideal syMinBase(ideal arg);
88
89BOOLEAN syTestOrder(ideal i);
90
91void syReOrderResolventFB(resolvente res,int length, int initial=1);
92
93resolvente syLaScala1(ideal arg,int * length);
94syStrategy syLaScala3(ideal arg,int * length);
95
96syStrategy syLaScala(ideal arg, int& maxlength, intvec* weights = NULL);
97
98syStrategy syHilb(ideal arg,int * length);
99syStrategy syKosz(ideal arg,int * length);
100
101void syKillComputation(syStrategy syzstr, ring r=currRing);
102intvec * syBettiOfComputation(syStrategy syzstr, BOOLEAN minim=TRUE,int * row_shift=NULL, intvec *weights=NULL);
103
104
105int sySize(syStrategy syzstr);
106int syDim(syStrategy syzstr);
107syStrategy syCopy(syStrategy syzstr);
108void syPrint(syStrategy syzstr,   const char *currRingName /* = currRingHdl->id */);
109
110syStrategy syMinimize(syStrategy syzstr);
111void syKillEmptyEntres(resolvente res,int length);
112
113extern int *  currcomponents;
114extern long *  currShiftedComponents;
115
116void syDeletePair(SObject * so);
117void syInitializePair(SObject * so);
118void syCopyPair(SObject * argso, SObject * imso);
119void syCompactifyPairSet(SSet sPairs, int sPlength, int first);
120void syCompactify1(SSet sPairs, int* sPlength, int first);
121SRes syInitRes(ideal arg,int * length, intvec * Tl, intvec * cw=NULL);
122void syResetShiftedComponents(syStrategy syzstr, int index,int hilb=0);
123void syEnlargeFields(syStrategy syzstr,int index);
124void syEnterPair(syStrategy syzstr, SObject * so, int * sPlength,int index);
125SSet syChosePairs(syStrategy syzstr, int *index, int *howmuch, int * actdeg);
126int syInitSyzMod(syStrategy syzstr, int index, int init=17);
127long syReorderShiftedComponents(long * sc, int n);
128void syGaussForOne(ideal arg,int gen,int ModComp,int from=-1,int till=0);
129void syEnterPair(SSet sPairs, SObject * so, int * sPlength,int index);
130void syEnterPair(syStrategy syzstr, SObject * so, int * sPlength,int index);
131syStrategy syKosz(ideal arg,int * length);
132
133resolvente syReorder(resolvente res,int length,
134        syStrategy syzstr,BOOLEAN toCopy=TRUE,resolvente totake=NULL);
135resolvente syReorder(resolvente res,int length,
136        syStrategy syzstr,BOOLEAN toCopy/*=TRUE*/,resolvente totake/*=NULL*/);
137
138#endif
Note: See TracBrowser for help on using the repository browser.