/dports/math/palp/palp-2.20/ |
H A D | Global.h | 22 #ifndef POLY_Dmax /* You can set POLY_Dmax at compilation: 24 #define POLY_Dmax 6 /* max dim of polytope */ macro 32 #if (POLY_Dmax <= 3) 38 #elif (POLY_Dmax == 4) 139 int nw, N, z[POLY_Dmax][AMBI_Dmax], m[POLY_Dmax], nz, index;} CWS; 172 GL_Long G[VERT_Nmax][POLY_Dmax][POLY_Dmax]; 273 Long NF[POLY_Dmax][VERT_Nmax], int t, int S, int N); 285 void IP_Simplex_Decomp(Long CM[][POLY_Dmax], int p, int d, 527 INCI v[POLY_Dmax+1][FACE_Nmax]; /* vertex info */ 528 INCI f[POLY_Dmax+1][FACE_Nmax]; /* V-on-dual info */ [all …]
|
H A D | LG.h | 3 #define W_Nmax (POLY_Dmax+1) 25 typedef struct {int d, N, z[POLY_Dmax][W_Nmax], m[POLY_Dmax], M, r, R;/* Ref */ 26 Long w[W_Nmax], B[W_Nmax][POLY_Dmax], A[W_Nmax], rI[POLY_Dmax]; 30 typedef struct {int D,E,sts; int h[POLY_Dmax][POLY_Dmax];} VaHo; 40 Long V_to_G_GI(Long *V,int d, Long G[][POLY_Dmax],Long GI[][POLY_Dmax]);
|
H A D | Polynf.c | 547 { SL_Long S[POLY_Dmax][POLY_Dmax], X[POLY_Dmax]; 728 { GL_Long AV[POLY_Dmax],AG[POLY_Dmax][POLY_Dmax],AGI[POLY_Dmax][POLY_Dmax], 802 { GL_Long G[POLY_Dmax][POLY_Dmax],B[POLY_Dmax][POLY_Dmax]; 849 GL_Long G[POLY_Dmax][POLY_Dmax],GI[POLY_Dmax][POLY_Dmax]; int f=0; 1002 Y[POLY_Dmax][POLY_Dmax], G[POLY_Dmax][POLY_Dmax], *g[POLY_Dmax]; 1017 G[POLY_Dmax][POLY_Dmax],*g[POLY_Dmax]; 1434 *GI[POLY_Dmax][POLY_Dmax], GX[(POLY_Dmax*(POLY_Dmax+3))/2][POLY_Dmax]; 1598 *GI[POLY_Dmax][POLY_Dmax], GX[(POLY_Dmax*(POLY_Dmax+3))/2][POLY_Dmax]; 1623 *GI[POLY_Dmax][POLY_Dmax], GX[(POLY_Dmax*(POLY_Dmax+3))/2][POLY_Dmax]; 1682 GL_Long G[POLY_Dmax][POLY_Dmax],Ginv[POLY_Dmax][POLY_Dmax]; [all …]
|
H A D | Subpoly.h | 42 #if (POLY_Dmax < 5) 61 #define NB_MAX POLY_Dmax*VERT_Nmax /* an uneducated guess */ 63 #if(POLY_Dmax<5) 136 EqList *E, Long ANF[POLY_Dmax][VERT_Nmax]); 156 void Make_All_Sublat(NF_List *_L, int n, int v, subl_int diag[POLY_Dmax], 170 #if (POLY_Dmax <6) 194 Long NF[POLY_Dmax][VERT_Nmax], int *MS); /* OUT */ 203 Long pNF[POLY_Dmax][VERT_Nmax]); /* 1 if reflexive */ 207 Long NF[POLY_Dmax][VERT_Nmax], int *MS); /* OUT */
|
H A D | Nef.h | 3 #define W_Nmax (POLY_Dmax+1) 9 #define Pos_Max (POLY_Dmax + 2) 14 Long VM[FIB_POINT_Nmax][POLY_Dmax]; 96 Long x[POLY_Dmax];
|
H A D | Vertex.c | 337 { int p[POLY_Dmax], i, j, J=0; Long g=0, W[POLY_Dmax], *G[POLY_Dmax]; in VZ_to_Base() 350 { int i, j, k; Long W[POLY_Dmax], G[POLY_Dmax][POLY_Dmax]; in OrthBase_red_by_V() 397 B[(POLY_Dmax*(POLY_Dmax+1))/2][POLY_Dmax], W[POLY_Dmax]; if(_P->np<2) in GLZ_Start_Simplex() 749 Long InvMat[POLY_Dmax][POLY_Dmax]){ in Compute_InvMat() argument 752 LRat ind[POLY_Dmax][POLY_Dmax], x[POLY_Dmax], y[POLY_Dmax], f, in Compute_InvMat() 753 PInvMat[POLY_Dmax][POLY_Dmax]; in Compute_InvMat() 800 Long MaxDist[EQUA_Nmax], InvMat[POLY_Dmax][POLY_Dmax], Den=1; in Complete_Poly() 802 int OrdFac[VERT_Nmax], BasFac[POLY_Dmax], position[POLY_Dmax]; in Complete_Poly() 959 Long InvMat[POLY_Dmax][POLY_Dmax], FaceInfo *_I){ in lastline() argument 1006 Long MaxDist[EQUA_Nmax], InvMat[POLY_Dmax][POLY_Dmax], Den=1; in QComplete_Poly() [all …]
|
H A D | Coord.c | 93 i+=2; assert(*nz<POLY_Dmax); in Read_CWS_Zinfo() 109 if (POLY_Dmax < (polyDim + codim - 1)){ in checkDimension() 116 else if (POLY_Dmax < polyDim +1){ in checkDimension() 322 if(_CW->N-_CW->nw > POLY_Dmax){ in Read_CWS() 475 int *m, Long *M, Long G[POLY_Dmax][POLY_Dmax]) /* out */ in CWS_2_SublatZ() argument 483 { Long X[POLY_Dmax]; for(i=0;i<P->n;i++) in Reduce_PPL_2_Sublat() 496 { int i, j, Amin[POLY_Dmax+1]; Long *x=_P->x[_P->np=0], xmin[POLY_Dmax], in Old_Make_CWS_Points() 497 xmax[POLY_Dmax], Xmax[AMBI_Dmax], xaux[POLY_Dmax], L, R; CWS *_C=Cin; in Old_Make_CWS_Points() 498 CWLatticeBasis B; Long G[POLY_Dmax][POLY_Dmax],M[POLY_Dmax];int m=Cin->nz; in Old_Make_CWS_Points() 598 xmax[POLY_Dmax], Xmax[AMBI_Dmax], X0[AMBI_Dmax], xaux[POLY_Dmax], L, R; in Make_CWS_Points() [all …]
|
H A D | E_Poly.c | 14 Long G[VERT_Nmax][POLY_Dmax][POLY_Dmax],int fd[VERT_Nmax], 43 Long yold[POLY_Dmax]; in DYNadd_for_completion() 62 Long MaxDist[EQUA_Nmax], InvMat[POLY_Dmax][POLY_Dmax], Den=1; in DYNComplete_Poly() 63 Long yDen[POLY_Dmax]; in DYNComplete_Poly() 65 BasFac[POLY_Dmax], one[POLY_Dmax], position[POLY_Dmax]; in DYNComplete_Poly() 66 LRat ind[POLY_Dmax][POLY_Dmax], x[POLY_Dmax], y[POLY_Dmax], f, in DYNComplete_Poly() 67 PInvMat[POLY_Dmax][POLY_Dmax]; in DYNComplete_Poly() 285 Long G[VERT_Nmax][POLY_Dmax][POLY_Dmax]; in PRINT_Fibrations() 583 int d[POLY_Dmax], i, j; in PrintDegrees() 673 int h[POLY_Dmax][POLY_Dmax] = {{0}, {0}}, S[VERT_Nmax]; in Output() [all …]
|
H A D | Mori.h | 136 void IP_Simplex_Fiber(Long PM[][POLY_Dmax], int p, int d, /* need PM[i]!=0 */ FibW *F, int Wmax, in… 146 Long SimplexVolume(Long *V[POLY_Dmax+1],int d);
|
H A D | poly.c | 225 int SymNum, VPMSymNum; Long NF[POLY_Dmax][VERT_Nmax]; in main() 235 Long vB[POLY_Dmax],Z; int j; in main() 241 Long VM[POLY_Dmax][VERT_Nmax]; in main() 246 Long ANF[POLY_Dmax][VERT_Nmax]; in main()
|
H A D | Subpoly.c | 6 #define UnAided_IP_CHECK (POLY_Dmax>4) 373 Long drop_point[POLY_Dmax]; in Drop_and_Keep() 517 Long drop_point[POLY_Dmax]; in Start_Make_All_Subpolys() 576 Long perm_vec_in[POLY_Dmax], perm_vec_out[POLY_Dmax], abs_sum, max_abs_sum=0; in IREgcd() 603 Long RV[POLY_Dmax], FE[POLY_Dmax], dist=0; in Make_RedVec() 632 Long RedVec[POLY_Dmax]; in Reduce_Poly() 910 subl_int s[POLY_Dmax][POLY_Dmax], t[POLY_Dmax][POLY_Dmax], ts_lincol; in Make_All_Sublat() 1041 subl_int diag[POLY_Dmax]; in MakePolyOnSublat() 1207 int i,j,drop_point[POLY_Dmax]; in irred() 1341 Long drop_point[POLY_Dmax]; in Find_RSP_Drop_and_Keep() [all …]
|
H A D | LG.c | 35 Long BM[W_Nmax][W_Nmax], *B[W_Nmax], Wa[POLY_Dmax], Za[POLY_Dmax], in Read_WZ_PP() 36 F[W_Nmax], G[POLY_Dmax][POLY_Dmax], GI[POLY_Dmax][POLY_Dmax],X; in Read_WZ_PP() 60 i+=2; assert(*nz<POLY_Dmax); WZ->m[*nz]=0; in Read_WZ_PP() 136 { Long *PB[W_Nmax], Z[POLY_Dmax][VERT_Nmax], M[POLY_Dmax], D[POLY_Dmax]; 451 W.N=N; assert((N<=W_Nmax)&&(N<POLY_Dmax+2)); assert(P!=NULL); in MakeRefWeights() 605 { int b=0,c=0,l,M[POLY_Dmax],I[POLY_Dmax],eq[W_Nmax],j,J=W->M,N=W->N,vac=0; in Count_b01() 682 Long zsum1=0, zsum2=0, mo=W->d, kgV=mo, omega[POLY_Dmax]; in WIndex_HTrace() 788 int j,k,s, J=W->M, M[POLY_Dmax],I[POLY_Dmax],X=W->m[0],x[POLY_Dmax]; in LGO_VaHo() 799 int n[W_Nmax], N=0,h[POLY_Dmax+1],hn=0; in LGO_VaHo() 839 { int u,b,l,L[POLY_Dmax],w[W_Nmax],o[W_Nmax],io[W_Nmax],mm=1; in LGO_VaHo() [all …]
|
H A D | class.c | 21 #if ( POLY_Dmax * POINT_Nmax > 83400000 ) 22 #error decrease POLY_Dmax or/and POINT_Nmax for compiling class 259 #if (POLY_Dmax != 4) in main()
|
H A D | Subdb.c | 456 { int I,J,MS; Long NF[POLY_Dmax][VERT_Nmax]; in Print_Missing_Mirror() 907 Long NF[POLY_Dmax][VERT_Nmax]; Init_FInfoList(&L); in Bin2a() 960 Long NF[POLY_Dmax][VERT_Nmax]; in Bin2aDBsl() 1019 #if (POLY_Dmax < 6) 1473 Long NF[POLY_Dmax][VERT_Nmax]; in Extract_from_Hodge_db() 1668 Long G[][POLY_Dmax],Long *D,int *x) in Aux_Print_SLpoly() argument 1713 { Long D[POLY_Dmax],G[POLY_Dmax][POLY_Dmax],PM[VERT_Nmax][VERT_Nmax]; in PH_Sublat_Polys() 1735 { Long Z[POLY_Dmax][VERT_Nmax], M[POLY_Dmax]; int r; in PH_Sublat_Polys() 1767 { Long D[POLY_Dmax],G[POLY_Dmax][POLY_Dmax]; in V_Sublat_Polys() 1776 { int i,j; subl_int diag[POLY_Dmax], U[POLY_Dmax][VERT_Nmax]; in V_Sublat_Polys() [all …]
|
H A D | Nefpart.c | 11 GL_Long G[][POLY_Dmax]); /* current=best */ 48 int M[POLY_Dmax][POLY_Dmax]; 56 Long X[POLY_Dmax][VERT_Nmax]; 61 GL_Long G[POLY_Dmax][POLY_Dmax]; 70 int A[POLY_Dmax]; 259 int s[POLY_Dmax], d=0, n_flag, i, j; in NForm_S() 427 int gp[POLY_Dmax] = {0}, gp_flag = 1, i, j=0; in Codim_Check()
|
H A D | Subadd.c | 53 void Eval_Poly_NF(int *d,int *v,int *f, Long VM[POLY_Dmax][VERT_Nmax], 183 assert(*d <= POLY_Dmax); in Read_Bin_Info() 763 #if POLY_Dmax>3 in Print_Weight_Info() 805 void VPrint(int *d,int *v,Long V[POLY_Dmax][VERT_Nmax]) in VPrint() 935 void AuxVnf2ucNF(Long nf[POLY_Dmax][VERT_Nmax], int *d, int *v, in AuxVnf2ucNF() 948 Long NF[POLY_Dmax][VERT_Nmax], int *MS) /* OUT */ in UCnf2vNF() 980 { Long V_NF[POLY_Dmax][VERT_Nmax], F_NF[POLY_Dmax][VERT_Nmax]; in VF_2_ucNF() 981 Long VM[POLY_Dmax][VERT_Nmax], VPM[VERT_Nmax][VERT_Nmax]; in VF_2_ucNF() 982 unsigned char auxUC[POLY_Dmax*VERT_Nmax]; in VF_2_ucNF() 1079 { Long tNF[POLY_Dmax][VERT_Nmax]; int i, j, tMS, NV, NUC; in Test_ucNF() [all …]
|
H A D | cws.c | 114 Long x[POLY_Dmax+1][POLY_Dmax]; 115 EqList q[POLY_Dmax]; 116 INCI qI[POLY_Dmax][EQUA_Nmax]; 117 int f0[POLY_Dmax]; 167 assert(n<POLY_Dmax); in PrintPoint() 176 assert(n<POLY_Dmax); in PrintQ() 318 Long yq[POLY_Dmax]; in RecConstructRgcWeights() 352 Long y[POLY_Dmax]; in WsIpCheck() 353 Long yq[POLY_Dmax]; in WsIpCheck() 834 assert((N<=W_Nmax)&&(N<POLY_Dmax+2)); in MakeMoonWeights() [all …]
|
H A D | MoriCone.c | 145 Long cdiv=0, sv, nok, *X[POLY_Dmax]; int d=P->n, C[VERT_Nmax]; in DivClassBasis() 304 if(r>POLY_Dmax){fprintf(outFILE,"need POLY_Dmax>=%d\n",UT->n);exit(0);} in Check_Mori() 350 if(r>POLY_Dmax){fprintf(outFILE,"need POLY_Dmax>=%d\n",UT->n);exit(0);} in Print_Mori() 474 if(P->n<(POLY_Dmax+1)){ in InterSectionRing() 560 void AuxPrintRays(int R[VERT_Nmax][POLY_Dmax],int nrp[VERT_Nmax],int nr){ in AuxPrintRays() argument 601 { int d=P->n,i,j,k=0,tnt=0, F[POLY_Dmax+1]; Inci64 CI=0; Matrix A,B; in Triang1dSFan() 636 { int d=P->n,v=d+2,i,j,k=0,r,F[POLY_Dmax+2],Z[POLY_Dmax+2],z=0, // GKZ polygon: in Triang2dSFan() 637 R[VERT_Nmax][POLY_Dmax],nrp[VERT_Nmax],nr,neg=0,tnt=0;// num.max.tri in Triang2dSFan() 1048 Inci64 MT[VERT_Nmax*POLY_Dmax], *_CT[ANfan]; in GKZsubdivide() 1600 if(r>POLY_Dmax){fprintf(outFILE,"need POLY_Dmax>=%d\n",UT->n);exit(0);} in Print_Mori_Old() [all …]
|
H A D | nef.c | 15 typedef struct { Long x[POLY_Dmax][W_Nmax]; int N, n; } AmbiLatticeBasis; 110 int Make_Mirror(EPoly *_EP, int h[][POLY_Dmax], int D, int dim); 286 if (POLY_Dmax < (_P->n + codim - 1)){ in main()
|
H A D | SingularInput.c | 202 Long *X[POLY_Dmax], Vijkl; in HyperSurfSingular() 203 int x[POLY_Dmax], n=0; i=j=0; in HyperSurfSingular()
|
H A D | mori.c | 181 if (POLY_Dmax < (_P->np - _P->n)){ in main()
|