diff liboctave/oct-fftw.cc @ 5275:23b37da9fd5b

[project @ 2005-04-08 16:07:35 by jwe]
author jwe
date Fri, 08 Apr 2005 16:07:37 +0000
parents 50140aa23b32
children 4c8a2e4e0717
line wrap: on
line diff
--- a/liboctave/oct-fftw.cc	Thu Apr 07 21:51:37 2005 +0000
+++ b/liboctave/oct-fftw.cc	Fri Apr 08 16:07:37 2005 +0000
@@ -57,11 +57,11 @@
   octave_fftw_planner (void);
 
   fftw_plan create_plan (int dir, const int rank, const dim_vector dims, 
-			 int howmany, int stride, int dist, 
+			 octave_idx_type howmany, octave_idx_type stride, octave_idx_type dist, 
 			 const Complex *in, Complex *out);
 
   fftw_plan create_plan (const int rank, const dim_vector dims, 
-			 int howmany, int stride, int dist, 
+			 octave_idx_type howmany, octave_idx_type stride, octave_idx_type dist, 
 			 const double *in, Complex *out);
 
 private:
@@ -74,16 +74,16 @@
   fftw_plan plan[2];
 
   // dist
-  int d[2];
+  octave_idx_type d[2];
 
   // stride
-  int s[2];
+  octave_idx_type s[2];
 
   // rank
   int r[2];
 
   // howmany
-  int h[2];
+  octave_idx_type h[2];
 
   // dims
   dim_vector n[2];
@@ -95,16 +95,16 @@
   fftw_plan rplan;
 
   // dist
-  int rd;
+  octave_idx_type rd;
 
   // stride
-  int rs;
+  octave_idx_type rs;
 
   // rank
   int rr;
 
   // howmany
-  int rh;
+  octave_idx_type rh;
 
   // dims
   dim_vector rn;
@@ -136,8 +136,8 @@
 
 fftw_plan
 octave_fftw_planner::create_plan (int dir, const int rank,
-				  const dim_vector dims, int howmany,
-				  int stride, int dist, 
+				  const dim_vector dims, octave_idx_type howmany,
+				  octave_idx_type stride, octave_idx_type dist, 
 				  const Complex *in, Complex *out)
 {
   int which = (dir == FFTW_FORWARD) ? 0 : 1;
@@ -207,7 +207,7 @@
  
 fftw_plan
 octave_fftw_planner::create_plan (const int rank, const dim_vector dims, 
-				  int howmany, int stride, int dist, 
+				  octave_idx_type howmany, octave_idx_type stride, octave_idx_type dist, 
 				  const double *in, Complex *out)
 {
   fftw_plan *cur_plan_p = &rplan;
@@ -274,7 +274,7 @@
 
 static inline void
 convert_packcomplex_1d (Complex *out, size_t nr, size_t nc,
-			int stride, int dist)
+			octave_idx_type stride, octave_idx_type dist)
 {
   OCTAVE_QUIT;
 
@@ -351,7 +351,7 @@
 
 int
 octave_fftw::fft (const double *in, Complex *out, size_t npts, 
-		  size_t nsamples, int stride, int dist)
+		  size_t nsamples, octave_idx_type stride, octave_idx_type dist)
 {
   dist = (dist < 0 ? npts : dist);
 
@@ -371,7 +371,7 @@
 
 int
 octave_fftw::fft (const Complex *in, Complex *out, size_t npts, 
-		  size_t nsamples, int stride, int dist)
+		  size_t nsamples, octave_idx_type stride, octave_idx_type dist)
 {
   dist = (dist < 0 ? npts : dist);
 
@@ -388,7 +388,7 @@
 
 int
 octave_fftw::ifft (const Complex *in, Complex *out, size_t npts, 
-		   size_t nsamples, int stride, int dist)
+		   size_t nsamples, octave_idx_type stride, octave_idx_type dist)
 {
   dist = (dist < 0 ? npts : dist);
 
@@ -412,14 +412,14 @@
 octave_fftw::fftNd (const double *in, Complex *out, const int rank, 
 		    const dim_vector &dv)
 {
-  int dist = 1;
+  octave_idx_type dist = 1;
   for (int i = 0; i < rank; i++)
     dist *= dv(i);
 
   // Fool with the position of the start of the output matrix, so that
   // creating other half of the matrix won't cause cache problems.
 
-  int offset = (dv.numel () / dv(0)) * ((dv(0) - 1) / 2); 
+  octave_idx_type offset = (dv.numel () / dv(0)) * ((dv(0) - 1) / 2); 
   
   fftw_plan plan = fftw_planner.create_plan (rank, dv, 1, 1, dist,
 					     in, out + offset);
@@ -438,7 +438,7 @@
 octave_fftw::fftNd (const Complex *in, Complex *out, const int rank, 
 		    const dim_vector &dv)
 {
-  int dist = 1;
+  octave_idx_type dist = 1;
   for (int i = 0; i < rank; i++)
     dist *= dv(i);
 
@@ -456,7 +456,7 @@
 octave_fftw::ifftNd (const Complex *in, Complex *out, const int rank, 
 		     const dim_vector &dv)
 {
-  int dist = 1;
+  octave_idx_type dist = 1;
   for (int i = 0; i < rank; i++)
     dist *= dv(i);