Re: [ublas] numeric bindings nrm2 differences

Subject: Re: [ublas] numeric bindings nrm2 differences
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2011-11-29 09:26:57

Hi Philipp,

if you used ATLAS and defined BOOST_NUMERIC_BINDINGS_BLAS_CBLAS, there shouldn't be any problems with return conventions. So as a first step, I recommend to replace the call

double norm2 = boost::numeric::bindings::blas::nrm2( v );
by 
double norm2 = norm_2( v );

and see whether this gives you the expected results.

Then you can try to call 
double norm2 = cblas_dnrm2( n, x, incx );
"directly", and see whether the problem still persists. If it does, it might indicate a problem of the used ATLAS library.

Regards,
Thomas