[petsc-dev] PCApplyTranspose in BiCG

[Barry Smith]

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
>