On Mon, Dec 26, 2016 at 13:53:47 +0200, Tommi Höynälänmaa wrote:
I'm using C++ and liboctave for programming. When I use
Matrix::lssolve for matrices given later in this message my C++
program gives false results even though command "mxa \ vb" works fine
in the Octave interpreter. The reason for the error is that BLAS
routine dnrm2 gives false results with standalone C++/liboctave
programs. When I wrote a straightforward new implementation of dnrm2
and used it in LAPACK routine dgelsd and related routines it uses
(dlarfg) the C++ program works fine.
Do you have any ideas what causes this odd behaviour of dnrm2?
I'm sure people would be happy to help, or at least try to reproduce, if
you could provide your C++ code sample to make it a bit easier.