Aldo,
This is a bug. The allocated number was not multipled by the block size
squared. Because of the this the compute unneeded number was garbage. We fixed
this a while ago in petsc-dev (yes it should have been fixed as a patch for 3.1
but that didn't happen.)
Barry
On Jun 30, 2
Dear all,
I have a doubt concerning the interpretation of the output from -vec_view.
I am using ILU(k) on a BAIJ matrix; the run is uni-processor.
> KSP Object:
> type: gmres
> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt
> Orthogonalization with no iterative refinement
>
