[petsc-users] preconditionning with a Matshell

Nicolas Pozin nicolas.pozin at inria.fr
Wed Jul 29 10:49:09 CDT 2015

Hello everyone, 

I'm trying to solve a system iteratively through conjugate gradient with jacobi preconditionning. 
Thing is the system matrix is a matshell. 
When I run the code matgetdiagonal returns an error : "MatGetDiagonal() line ......... /petsc/3.4.3/src/mat/interface/matrix.c Mat type shell" 
Is matgetiagonal function enable to act on a matshell? 

if someone has a clue on it.... 

thanks a lot 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150729/0a8ffc06/attachment.html>

More information about the petsc-users mailing list