[petsc-dev] CUDA guys, stuff still broken in master/next
    Karl Rupp 
    rupp at iue.tuwien.ac.at
       
    Sat Feb 21 13:50:48 CST 2015
    
    
  
Hi,
 >> I don't think that it is a good idea to unconditionally set the 
vector type to VECCUSP here (and at a few other locations in bjacobi.c). 
Particularly as this causes a test to leak memory which does not even 
use CUSP vectors. We may not have the 'correct' vector information 
available at this point, so we can only guess which type is correct.
>
>    Would a call to MatGetVecs() here tell you the type of the vectors you need? Yeah you create and then immediately destroy a vector but that is ok.
yeah, that should work. Particularly as VecSetType() currently destroys 
the just created vector anyway.
Thanks and best regards,
Karli
    
    
More information about the petsc-dev
mailing list