Re: [ublas] lapack bindings

Subject: Re: [ublas] lapack bindings
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2011-12-14 03:57:49

On 2011-12-13 21:09, David Mebane wrote:
>
>
> It compiles, but the lapack routines fail.  Could it have something to
> do with my version of lapack?
>
>
> Thanks,
> -David
>

Hello David,

on my machine this code produces the output shown below. I guess it has 
something to do with your toolchain/library environment. Could you 
provide more information on that?

Cheers,

Rutger

$ g++ test.cpp -llapack
$ ./a.out

real symmetric

sal = 
[5,5]((1,4,7,10,13),(4,5,8,11,14),(7,8,9,12,15),(10,11,12,13,16),(13,14,15,16,17))

al = 
[5,5]((1,0,0,0,0),(4,5,0,0,0),(7,8,9,0,0),(10,11,12,13,0),(13,14,15,16,17))

sau = 
[5,5]((1,4,7,10,13),(4,5,8,11,14),(7,8,9,12,15),(10,11,12,13,16),(13,14,15,16,17))

au = 
[5,5]((1,4,7,10,13),(0,5,8,11,14),(0,0,9,12,15),(0,0,0,13,16),(0,0,0,0,17))

bl = [5,2]((35,70),(42,84),(51,102),(62,124),(75,150))

bu = [5,2]((35,70),(42,84),(51,102),(62,124),(75,150))