diff --git a/pub/ftest/mi.f2t b/pub/ftest/mi.f2t new file mode 100755 index 0000000000000000000000000000000000000000..d8b1633e102ae43c229e04882fe106f86d5434aa --- /dev/null +++ b/pub/ftest/mi.f2t @@ -0,0 +1,9 @@ +#!/usr/bin/env frida +fm 101 1 h +0 oy sin(10*x) +0 oy .8*x +2 mi 1 +exit_unless(abs(y[3,0,0]-(0))<1e-15,"mi_test11") +exit_unless(abs(y[3,0,50]-(sin(4)))<1e-7,"mi_test12") +exit_unless(abs(y[3,0,100]-(sin(8)))<1e-7,"mi_test13") +exit(0) diff --git a/pub/ftest/mpaf.f2t b/pub/ftest/mpaf.f2t index df2d09da28e649eb0745a4d3611965cb78f9e349..7703c7138f6714e7aa66c562f4bbe5c29684b80f 100755 --- a/pub/ftest/mpaf.f2t +++ b/pub/ftest/mpaf.f2t @@ -1,5 +1,5 @@ #!/usr/bin/env frida -fm 10 1 h +fm 11 1 h 0 oy i^2 1 mpaf 3 exit_unless(abs(x[2,0,2]-.7)<1e-15,"mpaf_x") diff --git a/pub/lib/import.cpp b/pub/lib/import.cpp index 21485c6ecef20ae23e36847c735df64fcf51c041..921a052ff347e88ef28ef7f3cf533eb1ace4ebb0 100644 --- a/pub/lib/import.cpp +++ b/pub/lib/import.cpp @@ -244,8 +244,12 @@ void NImport::make_grid() PSpec s(new CSpec); s->x.resize(ni); - for (int i = 0; i < ni; ++i) - s->x[i] = ((double)i) / ni; + if (ni==1) { + s->x[0] = 0; + } else { + for (int i = 0; i < ni; ++i) + s->x[i] = ((double)i) / (ni-1); + } s->y.clear(); s->y.resize(ni, 0.);