Error in LCD

[Barry Smith]

  Keita,

    Thanks for reporting the problem. It is our error, we changed the management of
work vectors so KSPSetUp() could be used before KSPSolve() but it was not updated
correctly in this file. 

   I have pushed a fix to the release repository and it will be in our next
patch. For now you simply change the two lines from

  ierr = VecDuplicateVecs(ksp->vec_rhs,restart+1,&lcd->P); CHKERRQ(ierr);
  ierr = VecDuplicateVecs(ksp->vec_rhs, restart + 1, &lcd->Q); CHKERRQ(ierr); 

to
   ierr = VecDuplicateVecs(ksp->work[0],restart+1,&lcd->P); CHKERRQ(ierr);
  ierr = VecDuplicateVecs(ksp->work[0], restart + 1, &lcd->Q); CHKERRQ(ierr); 

  Please report these bugs to petsc-maint at mcs.anl.gov instead, thanks

  Barry

  

On Mon, 29 Oct 2007, Keita Teranishi wrote:

> Hi,
> 
>  
> 
> I am having a trouble with LCD in KSP module.  I got following error messages as I executed src/ksp/ksp/example/tutorials/ex10.  There is no error message when running much simpler examples such as ex2.  
> 
>  
> 
> Regards,  
> 
>  
> 
> keita $ ex10 -f0 ./xerox2000c3.spd -f1 ./xerox2000c3.spd -ksp_type lcd
> 
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> 
> [0]PETSC ERROR: Null argument, when expecting valid pointer!
> 
> [0]PETSC ERROR: Null Object: Parameter # 1!
> 
> [0]PETSC ERROR: ------------------------------------------------------------------------
> 
> :           
> 
> :
> 
> :
> 
> :
> 
> [0]PETSC ERROR: ------------------------------------------------------------------------
> 
> [0]PETSC ERROR: VecDuplicateVecs() line 533 in src/vec/vec/interface/vector.c
> 
> [0]PETSC ERROR: KSPSetUp_LCD() line 48 in src/ksp/ksp/impls/lcd/lcd.c
> 
> [0]PETSC ERROR: KSPSetUp() line 198 in src/ksp/ksp/interface/itfunc.c
> 
> [0]PETSC ERROR: main() line 293 in src/ksp/ksp/examples/tutorials/ex10.c
> 
> [unset]: aborting job:
> 
> application called MPI_Abort(MPI_COMM_WORLD, 85) - process 0
> 
>  
> 
> ================================
>  Keita Teranishi
>  Math Software Group
>  Cray Inc.
>  keita at cray.com
> ================================
> 
>  
> 
>