1a5eb4965SSatish Balay #ifdef PETSC_RCS_HEADER 2*005c665bSBarry Smith static char vcid[] = "$Id: snesj2.c,v 1.7 1997/07/09 20:59:37 balay Exp bsmith $"; 391627157SBarry Smith #endif 491627157SBarry Smith 5639f9d9dSBarry Smith #include "src/mat/matimpl.h" /*I "mat.h" I*/ 6639f9d9dSBarry Smith #include "src/snes/snesimpl.h" /*I "snes.h" I*/ 791627157SBarry Smith 891627157SBarry Smith 95615d1e5SSatish Balay #undef __FUNC__ 105615d1e5SSatish Balay #define __FUNC__ "SNESDefaultComputeJacobianWithColoring" 1191627157SBarry Smith /*@C 1291627157SBarry Smith SNESDefaultComputeJacobianWithColoring 1391627157SBarry Smith 1491627157SBarry Smith Input Parameters: 1591627157SBarry Smith . snes - nonlinear solver object 1691627157SBarry Smith . x1 - location at which to evaluate Jacobian 1791627157SBarry Smith . ctx - MatFDColoring contex 1891627157SBarry Smith 1991627157SBarry Smith Output Parameters: 2091627157SBarry Smith . J - Jacobian matrix 2191627157SBarry Smith . B - Jacobian preconditioner 2291627157SBarry Smith . flag - flag indicating if the matrix nonzero structure has changed 2391627157SBarry Smith 2491627157SBarry Smith .keywords: SNES, finite differences, Jacobian 2591627157SBarry Smith 2691627157SBarry Smith .seealso: SNESSetJacobian(), SNESTestJacobian() 2791627157SBarry Smith @*/ 2891627157SBarry Smith int SNESDefaultComputeJacobianWithColoring(SNES snes,Vec x1,Mat *JJ,Mat *B,MatStructure *flag,void *ctx) 2991627157SBarry Smith { 3091627157SBarry Smith MatFDColoring color = (MatFDColoring) ctx; 3143a90d84SBarry Smith int ierr; 3288389382SBarry Smith 33*005c665bSBarry Smith ierr = MatFDColoringApply(*B,color,x1,flag,snes); CHKERRQ(ierr); 3491627157SBarry Smith return 0; 3591627157SBarry Smith } 3691627157SBarry Smith 3791627157SBarry Smith 3891627157SBarry Smith 39