[petsc-users] [petsc-maint #107932] run a simple example ex19 in src/snes/examples/tutorials

[Jed Brown]

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>