[petsc-users] Write Non-Zero Values of MPI Matrix on an MPI Vector

[Jed Brown]

Ali Berk Kahraman <aliberkkahraman at yahoo.com> writes:

> OK, here is the thing. I have a 2D cartesian regular grid. I am working 
> on wavelet method collocation method, which creates an irregular 
> adaptive grid by turning grid points on an off on the previously 
> mentioned cartesian grid. I store the grid and the values as sparse Mat 
> objects, where each entry to the matrix denotes the x and y location of 
> the value (x:row, y:column). 

This will have terrible load balance.  Better to have a legit data
structure (p4est/DMForest may be useful) to represent your adaptive
grid.