[petsc-dev] on implementing a new PetscViewer

Marco Zocca zocca.marco at gmail.com
Sat Jul 18 06:48:10 CDT 2015

 you are right, I shouldn't think of Mat objects as data.

On the other hand, there should be a way to look e.g. at the sparsity pattern.

Thank you for the pointer to VecGetArrayRead (and again, the
VecGet/VecRestore pair can be described as a bracket, thus protecting
the user from accidental leaks), moreover it's great to know that you
are familiar with Hs.

I'm getting very polarized opinions about this project, which makes me
very eager to get it done and release it.

On 18 July 2015 at 01:58, Jed Brown <jed at jedbrown.org> wrote:
> Marco Zocca <zocca.marco at gmail.com> writes:
>> (so: no need for copying data back to memory, during normal operation
>> of the library)
> Great.
>> I was asking about the Viewer for e.g. Vecs and Mats to hook up e.g.
>> with plotting functionality on the Hs side, but since this would be
>> needed "offline", I might just dump to file and load from there.
> How do you intend to plot matrices (linear operators, not arrays of
> numbers)?
> For vectors, you might as well access the array via VecGetArrayRead
> (managed via a withVec* construct in Haskell).

More information about the petsc-dev mailing list