[ublas] Example of use of boost::numeric::bindings::lapack with symmetric or triangular matrices?

Subject: [ublas] Example of use of boost::numeric::bindings::lapack with symmetric or triangular matrices?
From: Douglas Bates (bates_at_[hidden])
Date: 2011-06-09 18:33:58

I am struggling to determine how to use the lapack bindings.  As a
simple example, how would use potrf declared in
boost/numeric/bindings/lapack/computational/potrf.hpp to evaluate the
Cholesky decomposition of a positive definite symmetric matrix.  I
can't decide what type of matrix to pass to it.  On input it is a
symmetric matrix but on output it is a triangular matrix.

On a more general level, are there code examples using the lapack
bindings that I could study?