Hello Justin,
did you try the following: Create a BlockSparsityPattern (2,2). Then
copy the sparsity patterns sparsity_M, sparsity_S, sparsity_T into the
corresponding blocks. Do the same with the BlockSparseMatrix.
Best Regards,
Markus
On 03.05.2012 03:42, Justin Droba wrote:
Hi,
I have three SparseMatrix objects already created, M, S, and T. I'd
like to combine them into one BlockSparseMatrix of the form A = [M, S;
T, M]. All of these matrices are the same size (square) and have the
same sparsity patterns. I've read through the documentation on
BlockSparseMatrix and I can't quite figure out how to do this. I'm
sure I'm missing something obvious (unless what I want to do is
impossible!), but can someone point me to the methods I need to use?
I'd like to do this so I can use SparseDirectUMFPACK to solve Ax = b
with the above A. The documentation says that SparseDirectUMFPACK
should work with BlockSparseMatrix objects too.
Thank you!
Best,
-J
--
-----------------------------------------------------------------------------
Justin Droba
Office: A111 Wells Hall
Ph.D. Candidate in Applied Mathematics
Michigan State University
Department of Mathematics
East Lansing, MI 48824
-----------------------------------------------------------------------------
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
