further about PCComputeExplicitOperator()

[Dave May]

Hi Barry,
   Does PCComputeExplicitOperator() really build B*A or A*B?
The above code appears to just be applying the preconditioner to the each
column of the identity matrix, with no reference to the original operator or
the preconditioner side.

I only wanted to clarify this fact as at one stage I wrote a function to
compute B*A or A*B as I had convinced myself that
PCComputeExplicitOperator() just assembled the inverse of the
preconditioner.

Cheers,
    Dave.



On Sat, May 3, 2008 at 2:50 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>  Note also:
> PCComputeExplicitOperator - Computes the explicit preconditioned operator.
> this means it computes B*A or A*B (depending on left or right
> preconditioning). This beasty
> (unless you use Jacobi preconditioning) is always dense and it makes no
> sense to store in
> sparse format except for fun.
>
>   Barry
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080503/471d78d6/attachment.htm>