This is very specialized for implementing point block Jacobi; I don't think 
it is something you would want to use directly. 

     If you do want to use it, it simply returns the inverses of the block 
diagonals in column major form. You can then call MatSetValues() with with each 
of those blocks into another PETSc matrix. values[i*bs*bs] is the starting 
point of each block in the array.


> On Sep 19, 2016, at 4:55 AM, Cyrill Vonplanta <> 
> wrote:
> Dear PETSc-Users,
> I would like to use the inverted block diagonals of a a matrix. I have seen 
> the function MatInvertBlockDiagonal() but I don’t know how to create a matrix 
> out of them or an array of block matrizes.
> Does anyone have an example on how to use  **values to create a PETSc matrix?
> Thanks
> Cyrill

Reply via email to