[petsc-users] preconditionning with a Matshell

[Nicolas Pozin]

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 
regards 
Nicolas 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150729/0a8ffc06/attachment.html>