Actual source code: linesearchimpl.h

  1: #pragma once

  3: #include <petscsnes.h>
  4: #include <petsc/private/petscimpl.h>

  6: PETSC_EXTERN PetscBool      SNESLineSearchRegisterAllCalled;
  7: PETSC_EXTERN PetscErrorCode SNESLineSearchRegisterAll(void);
  8: PETSC_EXTERN PetscLogEvent  SNESLINESEARCH_Apply;

 10: typedef struct _LineSearchOps *LineSearchOps;

 12: struct _LineSearchOps {
 13:   PetscErrorCode (*view)(SNESLineSearch, PetscViewer);
 14:   SNESLineSearchApplyFn *apply;
 15:   PetscErrorCode (*precheck)(SNESLineSearch, Vec, Vec, PetscBool *, void *);
 16:   SNESLineSearchVIProjectFn *viproject;
 17:   SNESLineSearchVINormFn    *vinorm;
 18:   PetscErrorCode (*postcheck)(SNESLineSearch, Vec, Vec, Vec, PetscBool *, PetscBool *, void *);
 19:   PetscErrorCode (*setfromoptions)(SNESLineSearch, PetscOptionItems *);
 20:   PetscErrorCode (*reset)(SNESLineSearch);
 21:   PetscErrorCode (*destroy)(SNESLineSearch);
 22:   PetscErrorCode (*setup)(SNESLineSearch);
 23:   PetscErrorCode (*snesfunc)(SNES, Vec, Vec);
 24: };

 26: #define MAXSNESLSMONITORS 5

 28: struct _p_LineSearch {
 29:   PETSCHEADER(struct _LineSearchOps);

 31:   SNES snes;

 33:   void *data;

 35:   PetscBool setupcalled;

 37:   Vec vec_sol;
 38:   Vec vec_sol_new;
 39:   Vec vec_func;
 40:   Vec vec_func_new;
 41:   Vec vec_update;

 43:   PetscInt nwork;
 44:   Vec     *work;

 46:   PetscReal lambda;

 48:   PetscBool            norms;
 49:   PetscReal            fnorm;
 50:   PetscReal            ynorm;
 51:   PetscReal            xnorm;
 52:   SNESLineSearchReason result;
 53:   PetscBool            keeplambda;

 55:   PetscReal damping;
 56:   PetscReal maxstep;
 57:   PetscReal steptol;
 58:   PetscInt  max_its;
 59:   PetscReal rtol;
 60:   PetscReal atol;
 61:   PetscReal ltol;
 62:   PetscInt  order;

 64:   PetscReal precheck_picard_angle;

 66:   void *precheckctx;
 67:   void *postcheckctx;

 69:   PetscViewer monitor;
 70:   PetscErrorCode (*monitorftns[MAXSNESLSMONITORS])(SNESLineSearch, void *); /* monitor routine */
 71:   PetscErrorCode (*monitordestroy[MAXSNESLSMONITORS])(void **);             /* monitor context destroy routine */
 72:   void    *monitorcontext[MAXSNESLSMONITORS];                               /* monitor context */
 73:   PetscInt numbermonitors;                                                  /* number of monitors */
 74: };