STRUMPACK
STRUMPACK copied to clipboard
[Comment] Consult on StructuredMatrix factor
I would appreciate some insight on the proper use of STRUMPACK. I have a symmetric matrix C for which I would like to compute
o^T C^{-1} o
and
log (det(C))
With DenseMatrix this can be calculated in an efficient way using Cholesky by determining L^{-1}o and summing the determinant of the diagonal of L.
Looking at the API and examples, it was not clear how to use StructoredMatrix::factor() to accomplish the same thing.
I am using small matrices in preliminary testing but our real problem requires large matrices in distributed memory.
I would appreciate insight on how to proceed. Thanks!