[petsc-users] Customizeing MatSetValuesBlocked(...)

Jed Brown jedbrown at mcs.anl.gov
Thu Aug 9 10:28:56 CDT 2012


On Thu, Aug 9, 2012 at 9:26 AM, Jinquan Zhong <jzhong at scsolutions.com>wrote:

> ·         A’x=b, where A’=[A^-1 U B]= ****
>
> What does the notation above mean? ****
>
> **ü  **That means A’ consists A^-1 and B.  A’ is a sparse matrix at the
> order of 2 million.  B = sparse matrix at the order of 2 million to be
> mixed and mingled with A^-1.
>
>
"mingled" is not a very useful term. How about using an actual equation so
we can have some clue what you mean.


> ****
>
>  What are the dimensions and sparsity of each piece. ****
>
> **ü  **The dimension of A^-1 is about 225,000.  B is about 2 millions. ***
> *
>
> ** **
>
> If you have an "inverse" floating around, what sort of operator is it the
> inverse of.****
>
> **ü  **A is a double complex matrix we have generated from a practical
> application using FEM.
>
Then there is no way A^{-1} should be stored as a dense matrix.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120809/cb91bf52/attachment-0001.html>


More information about the petsc-users mailing list