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>
