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