diff src/fftw-3.3.5/simd-support/vsx.c @ 127:7867fa7e1b6b

Current fftw source

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/src/fftw-3.3.5/simd-support/vsx.c	Tue Oct 18 13:40:26 2016 +0100
@@ -0,0 +1,69 @@
+/*
+ * Copyright (c) 2003, 2007-14 Matteo Frigo
+ * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
+ *
+ * VSX SIMD implementation added 2015 Erik Lindahl.
+ * Erik Lindahl places his modifications in the public domain.
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+ *
+ */
+
+
+#include "ifftw.h"
+
+#if HAVE_VSX
+
+#if HAVE_SYS_SYSCTL_H
+#  include <sys/sysctl.h>
+#endif
+
+#include <signal.h>
+#include <setjmp.h>
+
+static jmp_buf jb;
+
+static void sighandler(int x)
+{
+     longjmp(jb, 1);
+}
+
+static int really_have_vsx(void)
+{
+     void (*oldsig)(int);
+     oldsig = signal(SIGILL, sighandler);
+     if (setjmp(jb)) {
+	  signal(SIGILL, oldsig);
+	  return 0;
+     } else {
+          float mem[2];
+          __asm__ __volatile__ ("stxsdx 0,0,%0" :: "r" (mem) : "memory" );
+	  signal(SIGILL, oldsig);
+	  return 1;
+     }
+     return 0;
+}
+
+int X(have_simd_vsx)(void)
+{
+     static int init = 0, res;
+     if (!init) {
+	  res = really_have_vsx();
+	  init = 1;
+     }
+     return res;
+}
+
+#endif