TSEPAdapt.hpp File Reference

Step adaptation. More...

#include <petsc/private/tsimpl.h>

Go to the source code of this file.

Classes
struct	TSAdapt_EP

Macros
#define	TSADAPTEP   "ep"

Functions
static PetscErrorCode	TSAdaptChoose_EP (TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept, PetscReal *wlte, PetscReal *wltea, PetscReal *wlter)
static PetscErrorCode	TSAdaptReset_EP (TSAdapt adapt)
static PetscErrorCode	TSAdaptDestroy_EP (TSAdapt adapt)
PetscErrorCode	TSAdaptCreate_EP (TSAdapt adapt)

Detailed Description

Step adaptation.

Problem implementation for mix element for large-strain elasticity

Todo:

Implementation of plasticity

Definition in file TSEPAdapt.hpp.

Macro Definition Documentation

TSADAPTEP

#define TSADAPTEP   "ep"

Examples

ep.cpp.

Definition at line 15 of file TSEPAdapt.hpp.

Function Documentation

TSAdaptChoose_EP()

static PetscErrorCode TSAdaptChoose_EP ( TSAdapt  adapt, TS  ts, PetscReal  h, PetscInt *  next_sc, PetscReal *  next_h, PetscBool *  accept, PetscReal *  wlte, PetscReal *  wltea, PetscReal *  wlter  )

static

Definition at line 19 of file TSEPAdapt.hpp.

                                                                           {
  PetscFunctionBegin;
 
  TSAdapt_EP *basic = static_cast<TSAdapt_EP *>(adapt->data);
 
  *next_sc = 0; /* Reuse the same order scheme */
  *wlte = -1;   /* Weighted local truncation error was not evaluated */
  *wltea = -1;  /* Weighted absolute local truncation error is not used */
  *wlter = -1;  /* Weighted relative local truncation error is not used */
 
  *accept = PETSC_TRUE;
  *next_h = h; /* Reuse the old step */
 
  SNES snes;
  ierr = TSGetSNES(ts, &snes);
  CHKERRG(ierr);
 
  SNESConvergedReason reason;
  ierr = SNESGetConvergedReason(snes, &reason);
  CHKERRG(ierr);
 
  int it;
  ierr = SNESGetIterationNumber(snes, &it);
  CHKERRG(ierr);
 
  if (reason < 0) {
    h *= 0.75;
    *next_h = h;
    PetscPrintf(
        PETSC_COMM_WORLD,
        "\tDiverged set step length: it = %d, h = %3.4g set h = %3.4g \n", it,
        h, *next_h);
  } else if (reason > 0) {
 
    int its_d = 6;
    ierr = PetscOptionsGetInt(PETSC_NULL, "", "-desired_nb_of_its", &its_d,
                              PETSC_NULL);
    CHKERRG(ierr);
 
    h *= 1;//sqrt(static_cast<double>(its_d) / static_cast<double>(it + 1));
    *next_h = PetscClipInterval(h, adapt->dt_min, adapt->dt_max);
    PetscPrintf(
        PETSC_COMM_WORLD,
        "\tConverged set step length: it = %d, h = %3.4g set h = %3.4g \n", it,
        h, *next_h);
  }
 
  PetscFunctionReturn(0);
}

TSAdaptCreate_EP()

PetscErrorCode TSAdaptCreate_EP

(

TSAdapt

adapt

)

Examples

ep.cpp.

Definition at line 87 of file TSEPAdapt.hpp.

                                               {
  TSAdapt_EP *ep;
  PetscFunctionBegin;
  ierr = PetscNewLog(adapt, &ep);
  CHKERRG(ierr);
  adapt->data = (void *)ep;
  adapt->ops->choose = TSAdaptChoose_EP;
  adapt->ops->reset = TSAdaptReset_EP;
  adapt->ops->destroy = TSAdaptDestroy_EP;
  PetscFunctionReturn(0);
}

TSAdaptDestroy_EP()

static PetscErrorCode TSAdaptDestroy_EP ( TSAdapt  adapt )

static

Definition at line 78 of file TSEPAdapt.hpp.

                                                       {
  PetscFunctionBegin;
  ierr = TSAdaptReset_EP(adapt);
  CHKERRG(ierr);
  ierr = PetscFree(adapt->data);
  CHKERRG(ierr);
  PetscFunctionReturn(0);
}

TSAdaptReset_EP()

static PetscErrorCode TSAdaptReset_EP ( TSAdapt  adapt )

static

Definition at line 72 of file TSEPAdapt.hpp.

                                                     {
  TSAdapt_EP *basic = static_cast<TSAdapt_EP *>(adapt->data);
  PetscFunctionBegin;
  PetscFunctionReturn(0);
}