On Tue, Mar 13, 2012 at 10:38, Paul Mullowney <paulm at txcorp.com> wrote:
> In bjacobi.c, the following comments are given for PCCreate_BJacobi. > These things should definitely be fixed, but I said "jacobi" (without the "b"). > > Developer Notes: This preconditioner does not currently work with > CUDA/CUSP for a couple of > reasons. > > (1) It creates seq vectors as work vectors that should be > cusp > > (2) The use of VecPlaceArray() is not handled properly by CUSP (that is it > will not know where the ownership of the vector is so may use wrong values) > even if it did know the ownership it may induce extra copy ups and downs. > Satish suggests a VecTransplantArray() to handle two vectors sharing the > same pointer and handling the CUSP side as well instead of > VecGetArray()/VecPlaceArray(). > > > I've had some success with comment (2) above in regards to Singleblock > BJacobi. I haven't yet tried it for Multiblock or Multiproc. > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120313/bd23aa8d/attachment.htm>
