[petsc-users] Diagnosing Poisson Solver Behavior

[K. N. Ramachandran]

Hello Matt,

On Thu, Oct 15, 2015 at 12:18 PM, Matthew Knepley <knepley at gmail.com> wrote:

>
> The bad design decision is using a storage format designed for high speed
> MatMult() to use for data interchange. You should be
> using the API for value insertion rather than some storage format.
>

Ok, noted. Instead of using CreateMatSeqAIJWithArrays, I have also tried
setting up the Mat passed into the ComputeMatrix function, using
MatSetValue where the cols and values are referred to values in the CSR
format, but I still get the same error. It looks like the Mat is not
preallocated properly in the DM, even when I use DMDACreate3d.

Could you suggest what I can do to resolve this?

Thanks,
K.N.Ramachandran
Ph: 814-441-4279
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151015/f7312e00/attachment.html>