|
|
@ -101,7 +101,7 @@ int main(void) |
|
|
|
for (n0 = 0; n0 < 40; n0 = 2 * n0 + 1) { |
|
|
|
for (n0 = 0; n0 < 40; n0 = 2 * n0 + 1) { |
|
|
|
for (n1 = 0; n1 < 10; n1 = 2 * n1 + 1) { |
|
|
|
for (n1 = 0; n1 < 10; n1 = 2 * n1 + 1) { |
|
|
|
double best_error = 1000000000; |
|
|
|
double best_error = 1000000000; |
|
|
|
double bestpar0 = 1; |
|
|
|
double bestpar0 = n0 ? 1 : 100000; |
|
|
|
double bestpar1 = 1; |
|
|
|
double bestpar1 = 1; |
|
|
|
int better, i; |
|
|
|
int better, i; |
|
|
|
|
|
|
|
|
|
|
@ -118,8 +118,7 @@ int main(void) |
|
|
|
double par0, par1; |
|
|
|
double par0, par1; |
|
|
|
better = 0; |
|
|
|
better = 0; |
|
|
|
for (par0 = bestpar0 * 0.8; par0 <= bestpar0 * 1.21; par0 += bestpar0 * 0.05) { |
|
|
|
for (par0 = bestpar0 * 0.8; par0 <= bestpar0 * 1.21; par0 += bestpar0 * 0.05) { |
|
|
|
double range = n0 ? 0.8 : 1.0; |
|
|
|
for (par1 = bestpar1 * 0.8; par1 <= bestpar1 * 1.21; par1 += bestpar1 * 0.05) { |
|
|
|
for (par1 = bestpar1 * range; par1 <= bestpar1 * (2.01 - range); par1 += bestpar1 * 0.05) { |
|
|
|
|
|
|
|
double error = 0; |
|
|
|
double error = 0; |
|
|
|
TimeFilter *tf = ff_timefilter_new(1, par0, par1); |
|
|
|
TimeFilter *tf = ff_timefilter_new(1, par0, par1); |
|
|
|
for (i = 0; i < SAMPLES; i++) { |
|
|
|
for (i = 0; i < SAMPLES; i++) { |
|
|
|