[petsc-users] Optimal coloring using -snes_fd_color with DMDACreate

[Smith, Barry F.]

   I'm not familiar with that paper so I am not sure if it produces exactly the same colorings. It is based on the same principle. The coloring is done in 
/src/dm/impls/da/fdda.c DMCreateColoring_DA_3d_MPIAIJ() so you could compare. Note if they have more efficient special cases in the paper they could be added to the code without too much difficulty.

   Barry


> On Apr 12, 2018, at 8:26 AM, Buesing, Henrik <hbuesing at eonerc.rwth-aachen.de> wrote:
> 
> Dear all, 
>  
> When I use –snes_fd_color with my da created by DMDACreate3d, will it use an optimal coloring algorithm (see [1]) to estimate the Jacobian?
>  
> Thank you!
> Henrik
>  
>  
>  
> [1]
> Goldfarb/Toint: Optimal estimation of Jacobian and Hessian matrices that arise in finite difference calculations
> https://doi.org/10.1090/S0025-5718-1984-0744925-5
>  
>  
> --
> Dipl.-Math. Henrik Büsing
> Institute for Applied Geophysics and Geothermal Energy
> E.ON Energy Research Center
> RWTH Aachen University
> ------------------------------------------------------
> Mathieustr. 10            |    Tel +49 (0)241 80 49907
> 52074 Aachen, Germany     |    Fax +49 (0)241 80 49889
> ------------------------------------------------------
> http://www.eonerc.rwth-aachen.de/GGE
> hbuesing at eonerc.rwth-aachen.de
> ------------------------------------------------------