[petsc-dev] PCApplyTranspose in BiCG

[Jed Brown]

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120417/7c4f0e57/attachment.html>