Blitz logo

Blitz Support :

From: Othmar Frey (ofrey_at_[hidden])
Date: 2004-09-02 05:14:24


Hi.

I am using (or rather trying to use) the fftw library in combination
with Blitz++.

I've learned that the arrays of complex data passed to fftw must be
16-byte aligned
to exploit the "Single Instruction Multiple Data" support of fftw.

Is there a way to make a Vector<complex<double> > be 16-byte aligned in
memory?

Until now, I used the fftw++ wrapper class by John C. Bowman in
combination with his
Array class (http://www.math.ualberta.ca/~bowman/Array.h) then created a
blitz::Vector
from preexisting memory by passing a pointer to an array of Bowman's
Array class.
But I am not that happy with this kind of "solution".

Any suggestions?

Best regards

Othmar

-- 
-----------------------------------------------------------------------
Othmar Frey                             Remote Sensing Laboratories RSL
ofrey_at_[hidden]                      University of Zurich
http://www.geo.unizh.ch/~ofrey          Winterthurerstrasse 190	
phone +41 44 635 51 64                  CH-8057 Zurich
fax   +41 44 635 68 46                  Switzerland
-----------------------------------------------------------------------