From 0ed7a8b68dd73e2b0473b8127aeca99f3bac9061 Mon Sep 17 00:00:00 2001
From: Thomas Grill <xovo@users.sourceforge.net>
Date: Wed, 1 Apr 2009 21:13:09 +0000
Subject: cleaned up grill externals - replaced with svn:externals to
 svn.grrrr.org/ext/trunk/

svn path=/trunk/; revision=10951
---
 externals/grill/vasp/source/rdx2fft.cpp | 82 ---------------------------------
 1 file changed, 82 deletions(-)
 delete mode 100644 externals/grill/vasp/source/rdx2fft.cpp

(limited to 'externals/grill/vasp/source/rdx2fft.cpp')

diff --git a/externals/grill/vasp/source/rdx2fft.cpp b/externals/grill/vasp/source/rdx2fft.cpp
deleted file mode 100644
index b4ecf6c9..00000000
--- a/externals/grill/vasp/source/rdx2fft.cpp
+++ /dev/null
@@ -1,82 +0,0 @@
-#include <math.h>
-
-#define PI 3.1415926535897932384f
-
-//////////////////////////////////////////////////////////////////////////
-
-/* calculate bidirectional fourier transform of complex data radix 2 */
-/* adapted from subroutine FOUREA listed in                          */
-/* Programs for Digital Signal Processing                            */
-/* edited by Digital Signal Processing Committee                     */
-/* IEEE Acoustics Speech and Signal Processing Committee             */
-/* Chapter 1 Section 1.1 Page 1.1-4,5                                */
-/* direct -1 forward +1 reverse                                      */
-
-bool fft_bidir_complex_radix2(int size,float *real,float *imag,int direct)
-{
-  int i,j,m,mmax,istep;
-  float c,s,treal,timag,theta; 
-
-  /* compute transform */
-
-  j=1;
-  for(i=1;i<=size;i++)
-  {
-    if(i<j)
-    {
-      treal=real[j-1];
-      timag=imag[j-1];
-      real[j-1]=real[i-1];
-      imag[j-1]=imag[i-1];
-      real[i-1]=treal;
-      imag[i-1]=timag;
-    }
-    m=size/2;
-    while(j>m)
-    {
-      j-=m;
-      m=(m+1)/2;
-    }
-    j+=m;
-  }
-  mmax=1;
-  while(size>mmax)
-  {
-    istep=2*mmax;
-    for(m=1;m<=mmax;m++)
-    {
-      theta=PI*(float)direct*(float)(m-1)/(float)mmax;
-      c=(float)cos(theta);
-      s=(float)sin(theta);
-      for(i=m;i<=size;i+=istep)
-      {
-	j=i+mmax;
-	treal=real[j-1]*c-imag[j-1]*s;
-	timag=imag[j-1]*c+real[j-1]*s;
-	real[j-1]=real[i-1]-treal;
-	imag[j-1]=imag[i-1]-timag;
-	real[i-1]+=treal;
-	imag[i-1]+=timag;
-      }
-    }
-    mmax=istep;
-  }
-
-  return true;
-}
-
-#if 0
-/* calculate forward fourier transform of complex data radix 2 */
-
-bool fft_fwd_complex_radix2(int size,float *real,float *imag)
-{
-  return fft_bidir_complex_radix2(size,real,imag,-1);
-}
-
-/* calculate inverse fourier transform of complex data radix 2 */
-
-bool fft_inv_complex_radix2(int size,float *real,float *imag)
-{
-  return fft_bidir_complex_radix2(size,real,imag,1);
-}
-#endif
-- 
cgit v1.2.1