1 files changed, 75 insertions, 0 deletions
diff --git a/graphics/asymptote/fftw++.cc b/graphics/asymptote/fftw++.cc
new file mode 100644
index 0000000000..7c1095f930
--- /dev/null
+++ b/graphics/asymptote/fftw++.cc
@@ -0,0 +1,75 @@
+#include <cstring>
+#include <sstream>
+#include "fftw++.h"
+
+using namespace std;
+
+namespace fftwpp {
+
+const double fftw::twopi=2.0*acos(-1.0);
+
+// User settings:
+unsigned int fftw::effort=FFTW_MEASURE;
+const char *fftw::WisdomName=".wisdom";
+unsigned int fftw::maxthreads=1;
+double fftw::testseconds=0.2; // Time limit for threading efficiency tests
+
+fftw_plan (*fftw::planner)(fftw *f, Complex *in, Complex *out)=Planner;
+
+const char *fftw::oddshift="Shift is not implemented for odd nx";
+const char *inout=
+  "constructor and call must be both in place or both out of place";
+
+fft1d::Table fft1d::threadtable;
+mfft1d::Table mfft1d::threadtable;
+rcfft1d::Table rcfft1d::threadtable;
+crfft1d::Table crfft1d::threadtable;
+mrcfft1d::Table mrcfft1d::threadtable;
+mcrfft1d::Table mcrfft1d::threadtable;
+fft2d::Table fft2d::threadtable;
+
+void LoadWisdom()
+{
+  static bool Wise=false;
+  if(!Wise) {
+    ifstream ifWisdom;
+    ifWisdom.open(fftw::WisdomName);
+    ostringstream wisdom;
+    wisdom << ifWisdom.rdbuf();
+    ifWisdom.close();
+    const string& s=wisdom.str();
+    fftw_import_wisdom_from_string(s.c_str());
+    Wise=true;
+  }
+}
+
+void SaveWisdom()
+{
+  ofstream ofWisdom;
+  ofWisdom.open(fftw::WisdomName);
+  char *wisdom=fftw_export_wisdom_to_string();
+  ofWisdom << wisdom;
+  fftw_free(wisdom);
+  ofWisdom.close();
+}
+
+fftw_plan Planner(fftw *F, Complex *in, Complex *out)
+{
+  LoadWisdom();
+  fftw::effort |= FFTW_WISDOM_ONLY;
+  fftw_plan plan=F->Plan(in,out);
+  fftw::effort &= !FFTW_WISDOM_ONLY;
+  if(!plan) {
+    plan=F->Plan(in,out);
+    SaveWisdom();
+  }
+  return plan;
+}
+
+ThreadBase::ThreadBase() {threads=fftw::maxthreads;}
+
+}
+
+namespace utils {
+unsigned int defaultmpithreads=1;
+}