Great! Thank you! I will have a look at the code and I will also check how much 
function evaluations are needed in my case. 

Henrik


-- 
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
hbues...@eonerc.rwth-aachen.de

     -----Ursprüngliche Nachricht-----
     Von: Smith, Barry F. <bsm...@mcs.anl.gov>
     Gesendet: Donnerstag, 12. April 2018 16:42
     An: Buesing, Henrik <hbues...@eonerc.rwth-aachen.de>
     Cc: petsc-users <petsc-users@mcs.anl.gov>
     Betreff: Re: [petsc-users] Optimal coloring using -snes_fd_color with
     DMDACreate
     
     
        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 <hbues...@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
     > hbues...@eonerc.rwth-aachen.de
     > ------------------------------------------------------

Reply via email to