[petsc-users] example on manually setup matrix coloring for finite difference Jacobian?

[Matthew Knepley]

On Wed, Jul 17, 2013 at 6:02 PM, Zou (Non-US), Ling <ling.zou at inl.gov>wrote:

> Hi All,
>
> If I have a 2-d problem with unstructured triangle mesh, I suppose I
> need handle the coloring manually in case I need the finite
> differencing Jacobian. Also I guess it is doable as I know the mesh
> connectivity. I noticed there are several examples handling structured
> mesh case. Is there any example to setup the matrix coloring for this
> kind of case?
>

We have no examples of this, and have not pushed it in our development. I
would also consider forming a simplified operator for preconditioning while
using the finite difference approximation for the action.

   Matt


> Best,
>
> Ling
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130717/27d36234/attachment.html>