Subject: Re: [ublas] numeric bindings nrm2 differences
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2011-11-29 09:26:57
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 );
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.