diff options
Diffstat (limited to 'charf.c')
| -rw-r--r-- | charf.c | 6 |
1 files changed, 4 insertions, 2 deletions
@@ -138,7 +138,9 @@ void compute(EXPTYPE p, int D, VALUETYPE* f, VALUETYPE* records){ } - for(int k=1; k<=K; k++){ + differentiate(df[0],df[1],D,1); + differentiate(dMf[0],dMf[1],D,1); + for(int k=2; k<=K; k++){ /*Compute kth derivative of f and Mf from (k-1)th derivative*/ differentiate(df[(k+1)%2],df[k%2],D,1); differentiate(dMf[(k+1)%2],dMf[k%2],D,1); @@ -167,7 +169,7 @@ void compute(EXPTYPE p, int D, VALUETYPE* f, VALUETYPE* records){ double t = to_double(r); if( - over_threshold_charf(t,k) && + //over_threshold_charf(t,k) && is_greater(r,records[k])){ /*extra check for printing only because in error mode for some reason floats randomly seem to increase by tiny amounts*/ if(t>to_double(records[k])){ |