Re: opinions please for eigenvalues

Karl Meerbergen via ublas <[email protected]> Sat, 2 Jan 2021 21:12:28 +0000
Newsgroups gmane.comp.lib.boost.ublas
Message-ID <[email protected]>
--===============2234136174171983472==
Content-Language: en-US
Content-Type: multipart/alternative;
	boundary="_000_E3C5B1E406AF439E8BC9D62518D31BA1kuleuvenbe_"

--_000_E3C5B1E406AF439E8BC9D62518D31BA1kuleuvenbe_
Content-Type: text/plain; charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

Hi John,

The most reliable option is to link with Lapack and use bindings. I know th=
ere are people who only want c++ include files and no other libraries. If y=
ou are not scared of linking with fortran codes, this is the best choice. T=
here are uBLAS bindings for LAPACK in the boost.sandbox.

Best regards,

Karl

On 10 Dec 2020, at 13:10, Cem Bassoy via ublas <[email protected]<mailt=
o:[email protected]>> wrote:

Hi John,

thanks for using Boost!

Boost.uBlas is a C++ template class library that primarily focuses on BLAS =
level 1, 2, 3 functionality for dense, packed and sparse matrices, see http=
s://www.boost.org/doc/libs/1_74_0/libs/numeric/ublas/doc/overview.html#func=
tionality Recently, we also have introduced tensors.

Boost.uBlas is convient in terms familiarity and convenience if you have al=
ready experience with Boost and its setup. You will have most basic matrix =
operations such as matrix-matrix-multiplications. However, more complex lin=
ear algebra routines are not yet included.

Eigen is a framework itself and provides you with more functionality for no=
w.

Good luck with your software project.

Best,
CB



Am Mo., 30. Nov. 2020 um 18:38 Uhr schrieb John Mairs via ublas <ublas@list=
s.boost.org<mailto:[email protected]>>:
Hi,

Our software project has been using boost for years and we are very comfort=
able with it.  We have a new requirement to calculate eigenvalues from 6x6 =
matrix.  This is our first foray into linear algebra.  At first glance it a=
ppears boost ublas supports the 'vocabulary' types of matrix and vector and=
 basic operations between them.  I also read about boost bindings between u=
blas and other popular linear algebra libraries.

I took a python class on linear algebra and remember near the end of the co=
urse we had to write our own algorithm for computing eigen values so its pr=
obably in the realm of possible todo this ourselves using ublas only.

Or we could use something non boost like eigen:
http://eigen.tuxfamily.org/index.php?title=3DMain_Page

Any strong opinions folks?

John
_______________________________________________
ublas mailing list
[email protected]<mailto:[email protected]>
https://lists.boost.org/mailman/listinfo.cgi/ublas
Sent to: [email protected]<mailto:[email protected]>
_______________________________________________
ublas mailing list
[email protected]<mailto:[email protected]>
https://lists.boost.org/mailman/listinfo.cgi/ublas
Sent to: [email protected]


--_000_E3C5B1E406AF439E8BC9D62518D31BA1kuleuvenbe_
Content-Type: text/html; charset="us-ascii"
Content-ID: <[email protected]>
Content-Transfer-Encoding: quoted-printable

<html>
<head>
<meta http-equiv=3D"Content-Type" content=3D"text/html; charset=3Dus-ascii"=
>
</head>
<body style=3D"word-wrap: break-word; -webkit-nbsp-mode: space; line-break:=
 after-white-space;" class=3D"">
Hi John,
<div class=3D""><br class=3D"">
</div>
<div class=3D"">The most reliable option is to link with Lapack and use bin=
dings. I know there are people who only want c&#43;&#43; include files and =
no other libraries. If you are not scared of linking with fortran codes, th=
is is the best choice. There are uBLAS bindings
 for LAPACK in the boost.sandbox.</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Best regards,</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Karl</div>
<div><br class=3D"">
<blockquote type=3D"cite" class=3D"">
<div class=3D"">On 10 Dec 2020, at 13:10, Cem Bassoy via ublas &lt;<a href=
=3D"mailto:[email protected]" class=3D"">[email protected]</a>&gt; =
wrote:</div>
<br class=3D"Apple-interchange-newline">
<div class=3D"">
<div dir=3D"ltr" class=3D"">Hi John,<br class=3D"">
<br class=3D"">
thanks for using Boost!<br class=3D"">
<br class=3D"">
Boost.uBlas is a C&#43;&#43; template class library that primarily focuses =
on BLAS level 1, 2, 3 functionality for dense, packed and sparse matrices, =
see
<a href=3D"https://www.boost.org/doc/libs/1_74_0/libs/numeric/ublas/doc/ove=
rview.html#functionality" class=3D"">
https://www.boost.org/doc/libs/1_74_0/libs/numeric/ublas/doc/overview.html#=
functionality</a> Recently, we also have introduced tensors.<br class=3D"">
<br class=3D"">
<div class=3D"">Boost.uBlas is convient in terms familiarity and convenienc=
e if you have already experience with Boost and its setup. You will have mo=
st basic matrix operations such as matrix-matrix-multiplications. However, =
more complex linear algebra routines
 are not yet included.<br class=3D"">
</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Eigen is a framework itself and provides you with more func=
tionality for now.
<br class=3D"">
</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Good luck with your software project.</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Best,</div>
<div class=3D"">CB<br class=3D"">
</div>
<div class=3D""><br class=3D"">
</div>
</div>
<br class=3D"">
<br class=3D"">
<div class=3D"gmail_quote">
<div dir=3D"ltr" class=3D"gmail_attr">Am Mo., 30. Nov. 2020 um 18:38&nbsp;U=
hr schrieb John Mairs via ublas &lt;<a href=3D"mailto:[email protected]=
" class=3D"">[email protected]</a>&gt;:<br class=3D"">
</div>
<blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px 0.8ex;border-=
left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir=3D"ltr" class=3D"">Hi,
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Our software project has been using boost for years and we =
are very comfortable with it.&nbsp; We have a new requirement to calculate =
eigenvalues from 6x6 matrix.&nbsp; This is our first foray into linear alge=
bra.&nbsp; At first glance it appears boost ublas
 supports the 'vocabulary' types of matrix and vector and basic operations =
between them.&nbsp; I also read about boost bindings between ublas and othe=
r popular linear algebra libraries.&nbsp;&nbsp;</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">I took a python class on linear algebra and remember near&n=
bsp;the end of the course we had to write our own algorithm for computing e=
igen values&nbsp;so its probably in the realm of possible todo this ourselv=
es using ublas only.&nbsp;&nbsp;</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Or we could&nbsp;use something non boost like eigen:</div>
<div class=3D""><a href=3D"http://eigen.tuxfamily.org/index.php?title=3DMai=
n_Page" target=3D"_blank" class=3D"">http://eigen.tuxfamily.org/index.php?t=
itle=3DMain_Page</a><br class=3D"">
</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">Any strong opinions folks?</div>
<div class=3D""><br class=3D"">
</div>
<div class=3D"">John</div>
</div>
_______________________________________________<br class=3D"">
ublas mailing list<br class=3D"">
<a href=3D"mailto:[email protected]" target=3D"_blank" class=3D"">ublas=
@lists.boost.org</a><br class=3D"">
<a href=3D"https://lists.boost.org/mailman/listinfo.cgi/ublas" rel=3D"noref=
errer" target=3D"_blank" class=3D"">https://lists.boost.org/mailman/listinf=
o.cgi/ublas</a><br class=3D"">
Sent to: <a href=3D"mailto:[email protected]" target=3D"_blank" class=3D=
"">[email protected]</a><br class=3D"">
</blockquote>
</div>
_______________________________________________<br class=3D"">
ublas mailing list<br class=3D"">
<a href=3D"mailto:[email protected]" class=3D"">[email protected]</=
a><br class=3D"">
https://lists.boost.org/mailman/listinfo.cgi/ublas<br class=3D"">
Sent to: [email protected]<br class=3D"">
</div>
</blockquote>
</div>
<br class=3D"">
</body>
</html>

--_000_E3C5B1E406AF439E8BC9D62518D31BA1kuleuvenbe_--

--===============2234136174171983472==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline