On Apr 18, 2012, at 2:20 AM, Alexander Grayver wrote: > Jed, > > Thank you. > TFQMR works well now, however TCQMR doesn't. Although latter was working > before I updated petsc-dev. Might be same issue?
Please clarify. When was TCQMR working and when did it stop working? Was TCQMR working before you just updated in the last 24 hours? Did TCQMR work a long time ago (like TFQMR worked a long time ago?)? Barry > > 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> 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 >