[petsc-users] MatCreateMPIAIJWithSplitArrays for symmetric matrix?

[Jed Brown]

On Mon, Nov 8, 2010 at 09:32, Daniel Langr <daniel.langr at gmail.com> wrote:

> Is there an equivalent of MatCreateMPIAIJWithSplitArrays function for
> symmetric matrices? Or does this function work with symmetric matrices?
>

There is MatCreateMPISBAIJWithArrays() in petsc-dev.


> I work with symmetric matrices that are extremely expensive to compute and
> I don't know their structure in advance. So to do preallocation, it is not
> practical to compute a matrix twice - firstly only to get d_nnz and o_nnz
> arrays and secondly to set values. Now, I solve this problem by storing a
> matrix to a file system, which is faster but still far from ideal. The best
> solution I see is to compute a matrix only once, store it in memory as split
> CSR arrays and then just call something like
> MatCreateMPISBAIJWithSplitArrays.
>

Where do the matrices come from?  How do you know how to allocate to create
the arrays?

Jed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101108/756558d8/attachment.htm>