[petsc-dev] PCApplyTranspose in BiCG

[Alexander Grayver]

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>