On 10.04.2012 14:19, Jed Brown wrote: > > There is no space allocated at (44,44); note that counting is 0-based. > This consistency check has been in PETSc for many years, so I don't > know why you didn't trip it in earlier runs. Do you have some options > in PETSC_OPTIONS or .petscrc?
Jed, This is not possible: >> A = PetscBinaryRead('A.dat','complex',true); >> find(diag(A)==0) ans = Empty matrix: 0-by-1 >> A(45,45) ans = (1,1) 7.999999999999997e+002 +6.316546778308559e-009i Sorry, I forgot to mention that this is double complex matrix and vector. Have you read it as real double? > > If it is of any help files are here: > http://dl.dropbox.com/u/60982984/A.dat > http://dl.dropbox.com/u/60982984/RHS.dat > > Is this possible to look up changes in one file only? > > > $ hg log src/ksp/ksp/impls/tfqmr/ > > But the problem was clearly not in that file. -- Regards, Alexander -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120411/00aef2b3/attachment.html>