jiri/LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

increased storage for precomputed simplicial()

Browse Source
This commit is contained in:
Jiri Pittner
2025-11-20 15:28:17 +01:00
parent 7ba44efaef
commit 54642a71cc
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
miscfunc.cc
View File

@@ -767,7 +767,7 @@ return(bern[n]);
// //
//enlarge storage size if needed //enlarge storage size if needed
// //
#define SIMPLICIAL_MAXD 8 #define SIMPLICIAL_MAXD 16
#define SIMPLICIAL_MAXN 1024 #define SIMPLICIAL_MAXN 1024
unsigned long long simplicial(int d, int n) unsigned long long simplicial(int d, int n)
@@ -783,7 +783,7 @@ if(d==1) return n;
//resize precomputed part as needed //resize precomputed part as needed
if(n>maxn) if(n>maxn)
{ {
if(n>=SIMPLICIAL_MAXN) laerror("storage overflow in simplicial"); if(n>=SIMPLICIAL_MAXN) laerror("storage overflow in simplicial, increase MAXN");
int newdim=2*n; int newdim=2*n;
if(newdim<2*maxn) newdim=2*maxn; if(newdim<2*maxn) newdim=2*maxn;
if(newdim>=SIMPLICIAL_MAXN) newdim=SIMPLICIAL_MAXN-1; if(newdim>=SIMPLICIAL_MAXN) newdim=SIMPLICIAL_MAXN-1;
@@ -797,7 +797,7 @@ if(n>maxn)
if(d>maxd) if(d>maxd)
{ {
if(d>=SIMPLICIAL_MAXD) laerror("storage overflow in simplicial"); if(d>=SIMPLICIAL_MAXD) laerror("storage overflow in simplicial, increase MAXD");
for(int dd=maxd+1; dd<=d; ++dd) for(int dd=maxd+1; dd<=d; ++dd)
{ {
stored[dd-2][0]=0; stored[dd-2][0]=0;
@@ -817,7 +817,7 @@ int inverse_simplicial(int d, unsigned long long s)
if(s==0) return 0; if(s==0) return 0;
if(d==0 || s==1) return 1; if(d==0 || s==1) return 1;
if(d==1) return (int)s; if(d==1) return (int)s;
if(d>=SIMPLICIAL_MAXD) laerror("storage overflow in inverse_simplicial"); if(d>=SIMPLICIAL_MAXD) laerror("storage overflow in inverse_simplicial - increase MAXD");
static int maxd=0; static int maxd=0;
static Polynomial<double> polynomials[SIMPLICIAL_MAXD]; static Polynomial<double> polynomials[SIMPLICIAL_MAXD];