Date: 2013-07-02T14:51:01+00:00

Actual source code: tscreate.c
petsc-3.4.2 2013-07-02

      2: #include <petsc-private/tsimpl.h>      /*I "petscts.h"  I*/

      4: const char *const TSConvergedReasons_Shifted[] = {
      5:   "DIVERGED_STEP_REJECTED",
      6:   "DIVERGED_NONLINEAR_SOLVE",
      7:   "CONVERGED_ITERATING",
      8:   "CONVERGED_TIME",
      9:   "CONVERGED_ITS",
     10:   "TSConvergedReason","TS_",0};
     11: const char *const*TSConvergedReasons = TSConvergedReasons_Shifted + 2;

     13: #undef  __FUNCT__
     15: /*@C
     16:   TSCreate - This function creates an empty timestepper. The problem type can then be set with TSSetProblemType() and the
     17:        type of solver can then be set with TSSetType().

     19:   Collective on MPI_Comm

     21:   Input Parameter:
     22: . comm - The communicator

     24:   Output Parameter:
     25: . ts   - The TS

     27:   Level: beginner

     29: .keywords: TS, create
     30: .seealso: TSSetType(), TSSetUp(), TSDestroy(), MeshCreate(), TSSetProblemType()
     31: @*/
     32: PetscErrorCode  TSCreate(MPI_Comm comm, TS *ts)
     33: {
     34:   TS             t;

     39:   *ts = NULL;
     40: #if !defined(PETSC_USE_DYNAMIC_LIBRARIES)
     41:   TSInitializePackage();
     42: #endif

     44:   PetscHeaderCreate(t, _p_TS, struct _TSOps, TS_CLASSID, "TS", "Time stepping", "TS", comm, TSDestroy, TSView);
     45:   PetscMemzero(t->ops, sizeof(struct _TSOps));

     47:   /* General TS description */
     48:   t->problem_type      = TS_NONLINEAR;
     49:   t->vec_sol           = NULL;
     50:   t->numbermonitors    = 0;
     51:   t->snes              = NULL;
     52:   t->setupcalled       = 0;
     53:   t->data              = NULL;
     54:   t->user              = NULL;
     55:   t->ptime             = 0.0;
     56:   t->time_step         = 0.1;
     57:   t->time_step_orig    = 0.1;
     58:   t->max_time          = 5.0;
     59:   t->steps             = 0;
     60:   t->max_steps         = 5000;
     61:   t->ksp_its           = 0;
     62:   t->snes_its          = 0;
     63:   t->work              = NULL;
     64:   t->nwork             = 0;
     65:   t->max_snes_failures = 1;
     66:   t->max_reject        = 10;
     67:   t->errorifstepfailed = PETSC_TRUE;
     68:   t->rhsjacobian.time  = -1e20;
     69:   t->ijacobian.time    = -1e20;
     70:   t->equation_type     = TS_EQ_UNSPECIFIED;

     72:   t->atol             = 1e-4;
     73:   t->rtol             = 1e-4;
     74:   t->cfltime          = PETSC_MAX_REAL;
     75:   t->cfltime_local    = PETSC_MAX_REAL;
     76:   t->exact_final_time = TS_EXACTFINALTIME_STEPOVER;

     78:   *ts = t;
     79:   return(0);
     80: }