Re: [eigen] vectorization of complex

[Benoit Jacob <jacob.benoit.1@xxxxxxxxx>]

(It would then be nice to replace the division+sqrt by a single inverse_sqrt call, but we don't have that yet).

2010/7/23 Benoit Jacob <jacob.benoit.1@xxxxxxxxx>

There isn't, but this will compile to optimized code (only 1 array traversal):

array /= array.abs();

Benoit

2010/7/22 Carlos Becker <carlosbecker@xxxxxxxxx>

Ok, works now, thanks!

Btw, is there any Eigen function available to do what we were talking about in the other thread? I mean to 'normalize' every complex number in a given vector/matrix w.r.t. its abs value..

Thanks.

Carlos

On Thu, Jul 22, 2010 at 4:20 PM, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:

hm I fixed that during the merge, maybe I forgot to include this change, this is a matter of a few seconds...


gael

On Thu, Jul 22, 2010 at 5:11 PM, Carlos Becker <carlosbecker@xxxxxxxxx> wrote:

I am getting an error here, I was compiling with your branch before and no problem, but when compiling with the current dev branch (after your merge) I get:

In file included from eigen/Eigen/Geometry:51,

from eigen/Eigen/Dense:6,

from main.cpp:16:

eigen/Eigen/src/Geometry/arch/Geometry_SSE.h: In static member function ‘static Eigen::Quaternion<double> Eigen::ei_quat_product<1, Derived, OtherDerived, double, 1>::run(const Eigen::QuaternionBase<Derived>&, const Eigen::QuaternionBase<Derived1>&)’:

eigen/Eigen/src/Geometry/arch/Geometry_SSE.h:82: error: no matching function for call to ‘ei_pset1(const double&)’

eigen/Eigen/src/Geometry/arch/Geometry_SSE.h:83: error: no matching function for call to ‘ei_pset1(const double&)’

eigen/Eigen/src/Geometry/arch/Geometry_SSE.h:84: error: no matching function for call to ‘ei_pset1(const double&)’

eigen/Eigen/src/Geometry/arch/Geometry_SSE.h:85: error: no matching function for call to ‘ei_pset1(const double&)’

I guess there is something missing, but not sure. I am not using Geometry explicitly anywhere.

Thanks

On Thu, Jul 22, 2010 at 4:03 PM, Gael Guennebaud <gael.guennebaud@xxxxxxxxx> wrote:

merge done.

as usual I expect some troubles on non linux/gcc platforms, we'll see.....

gael