[PD-cvs] SF.net SVN: pure-data:[10316] trunk/externals/bsaylor/partconv~.c

[zmoelnig at users.sourceforge.net]

Revision: 10316
          http://pure-data.svn.sourceforge.net/pure-data/?rev=10316&view=rev
Author:   zmoelnig
Date:     2008-09-30 07:36:25 +0000 (Tue, 30 Sep 2008)

Log Message:
-----------
#ifdef's code that depends on fftw3: you have to define HAVE_FFTW3_H in order to compile this objecdt correctly

Modified Paths:
--------------
    trunk/externals/bsaylor/partconv~.c

Modified: trunk/externals/bsaylor/partconv~.c
===================================================================
--- trunk/externals/bsaylor/partconv~.c	2008-09-29 08:59:55 UTC (rev 10315)
+++ trunk/externals/bsaylor/partconv~.c	2008-09-30 07:36:25 UTC (rev 10316)
@@ -31,11 +31,13 @@
 
 #include <math.h>
 #include <string.h>
-#include <fftw3.h>
+#ifdef HAVE_FFTW3_H
+# include <fftw3.h>
+#endif /* FFTW3 */
 #include "m_pd.h"
 
 #ifdef __VEC__
-#include <altivec.h>
+# include <altivec.h>
 #endif
 
 #define MAXPARTS 256	// max number of partitions
@@ -48,9 +50,11 @@
 
 struct sumbuffer {
 	int index;
+#ifdef HAVE_FFTW3_H
 	fftwf_complex *fd;
 	float *td;
 	fftwf_plan plan;
+#endif /* FFTW3 */
 	int readpos;
 	struct sumbuffer *next, *prev;
 };
@@ -67,6 +71,7 @@
 	int ir_prepared;
 	int pd_blocksize;
 
+#ifdef HAVE_FFTW3_H
 	// partitions of impulse response
 	fftwf_plan irpart_plan;
 	float *irpart_td[MAXPARTS];
@@ -78,7 +83,9 @@
 	int inbufpos;
 	float *input_td;
 	fftwf_complex *input_fd;
+#endif /* FFTW3 */
 
+
 	// circular array/list of buffers for accumulating results of convolution
 	struct sumbuffer sumbufs[MAXPARTS+2];
 	int nsumbufs;  // number of sumbufs
@@ -127,6 +134,7 @@
 	t_float *in = (t_float *)(w[2]);
 	t_float *out = (t_float *)(w[3]);
 	int n = (int)(w[4]);
+#ifdef HAVE_FFTW3_H
 	int i;
 	int j;
 	int k;	// bin
@@ -140,7 +148,7 @@
 	v4sf *cursumbuf_fd;
 	v4sf *input_fd;
 	v4sf *irpart_fd;
-#else
+#elif defined HAVE_FFTW3_H
 	fftwf_complex *cursumbuf_fd;
 	fftwf_complex *input_fd;
 	fftwf_complex *irpart_fd;
@@ -165,9 +173,7 @@
 		x->curpart = 0;
 		memcpy(x->input_td, x->inbuf, x->partsize * sizeof(float));  // copy 'gathering' input buffer into 'transform' buffer
 		memset(&(x->input_td[x->partsize]), 0, (x->paddedsize - x->partsize) * sizeof(float));  // pad
-
 		fftwf_execute(x->input_plan);  // transform the input
-
 		// everything has been read out of prev sumbuf, so clear it
 		memset(x->sumbuf->prev->td, 0,  x->paddedsize * sizeof(float));
 
@@ -184,7 +190,7 @@
 	for (p = x->curpart; p < endpart; p++) {
 		// multiply the input block by the partition, accumulating the result in the appropriate sumbuf
 #ifdef USE_SSE
-#include "sse-conv.inc.c"
+# include "sse-conv.inc.c"
 #else
 		cursumbuf_fd = x->sumbufs[(x->sumbuf->index + p) % x->nsumbufs].fd;
 		input_fd = x->input_fd;
@@ -226,6 +232,10 @@
 	x->sumbuf->prev->readpos += n;
 
 	x->curcall++;
+#else /* !FFTW3 */
+  while(n-->0)
+    *out++=*in++;
+#endif
 
 	return (w+5);
 }
@@ -234,8 +244,8 @@
 
 static void partconv_free(t_partconv *x)
 {
+#ifdef HAVE_FFTW3_H
 	int i;
-
 	fftwf_free(x->inbuf);
 	for (i = 0; i < x->nparts; i++)
 		fftwf_free(x->irpart_td[i]);
@@ -245,6 +255,7 @@
 		fftwf_free(x->sumbufs[i].fd);
 		fftwf_destroy_plan(x->sumbufs[i].plan);
 	}
+#endif /* FFTW3 */
 }
 
 static void partconv_set(t_partconv *x, t_symbol *s)
@@ -267,7 +278,7 @@
 		pd_error(x, "%s: bad template", x->arrayname->s_name);
 		return;
 	}
-
+#ifdef HAVE_FFTW3_H
 	// if the IR has already been prepared, free everything first
 	if (x->ir_prepared == 1) {
 		partconv_free(x);
@@ -328,6 +339,7 @@
 
 	post("partconv~: using %s in %d partitions with FFT-size %d", x->arrayname->s_name, x->nparts, x->fftsize);
 	x->ir_prepared = 1;
+#endif /* FFTW3 */
 }
 
 static void partconv_dsp(t_partconv *x, t_signal **sp)
@@ -347,7 +359,7 @@
 	outlet_new(&x->x_obj, gensym("signal"));
 
 	if (argc != 2) {
-		post("argc = %d", argc);
+    //		post("argc = %d", argc);
 		error("partconv~: usage: [partconv~ <arrayname> <partsize>]\n\t- partition size must be a power of 2 >= blocksize");
 		return NULL;
 	}
@@ -373,6 +385,11 @@
 	x->ir_prepared = 0;
 	x->pd_blocksize = sys_getblksize();
 
+#ifndef HAVE_FFTW3_H
+  pd_error(x, "partconv~: compiled without FFTW3 support! this is a dummy!");
+#endif
+
+
 	return (x);
 }
 


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.