Jed, Thank you. TFQMR works well now, however TCQMR doesn't. Although latter was working before I updated petsc-dev. Might be same issue?
On 18.04.2012 07:30, Jed Brown wrote: > Thanks for the details, we should have identified the issue earlier. > When VecDot changed from using custom code to BLAS xDOT, the > conjugation definition was accidentally changed. If you pull after > this patch, TFQMR will work for your problem again. > > http://petsc.cs.iit.edu/petsc/petsc-dev/rev/8f37c498b7b4 > > On Wed, Apr 11, 2012 at 00:38, Alexander Grayver > <agrayver at gfz-potsdam.de <mailto:agrayver at gfz-potsdam.de>> wrote: > > 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 > > -- Regards, Alexander -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120418/6b34b559/attachment.html>
