1bdad233fSMatthew Knepley 2af0996ceSBarry Smith #include <petsc/private/tsimpl.h> /*I "petscts.h" I*/ 3bdad233fSMatthew Knepley 4193ac0bcSJed Brown const char *const TSConvergedReasons_Shifted[] = { 5193ac0bcSJed Brown "DIVERGED_STEP_REJECTED", 6193ac0bcSJed Brown "DIVERGED_NONLINEAR_SOLVE", 7193ac0bcSJed Brown "CONVERGED_ITERATING", 8193ac0bcSJed Brown "CONVERGED_TIME", 9193ac0bcSJed Brown "CONVERGED_ITS", 10193ac0bcSJed Brown "TSConvergedReason","TS_",0}; 11193ac0bcSJed Brown const char *const*TSConvergedReasons = TSConvergedReasons_Shifted + 2; 12193ac0bcSJed Brown 13bdad233fSMatthew Knepley #undef __FUNCT__ 14bdad233fSMatthew Knepley #define __FUNCT__ "TSCreate" 1560893bc3SSatish Balay /*@C 16bd6a702fSBarry Smith TSCreate - This function creates an empty timestepper. The problem type can then be set with TSSetProblemType() and the 17bd6a702fSBarry Smith type of solver can then be set with TSSetType(). 18bdad233fSMatthew Knepley 19bdad233fSMatthew Knepley Collective on MPI_Comm 20bdad233fSMatthew Knepley 21bdad233fSMatthew Knepley Input Parameter: 22bdad233fSMatthew Knepley . comm - The communicator 23bdad233fSMatthew Knepley 24bdad233fSMatthew Knepley Output Parameter: 25bdad233fSMatthew Knepley . ts - The TS 26bdad233fSMatthew Knepley 27bdad233fSMatthew Knepley Level: beginner 28bdad233fSMatthew Knepley 29bdad233fSMatthew Knepley .keywords: TS, create 3004ceabe5SMatthew G. Knepley .seealso: TSSetType(), TSSetUp(), TSDestroy(), TSSetProblemType() 31bdad233fSMatthew Knepley @*/ 320adebc6cSBarry Smith PetscErrorCode TSCreate(MPI_Comm comm, TS *ts) 330adebc6cSBarry Smith { 34bdad233fSMatthew Knepley TS t; 35dfbe8321SBarry Smith PetscErrorCode ierr; 36bdad233fSMatthew Knepley 37bdad233fSMatthew Knepley PetscFunctionBegin; 384482741eSBarry Smith PetscValidPointer(ts,1); 390298fd71SBarry Smith *ts = NULL; 40607a6623SBarry Smith ierr = TSInitializePackage();CHKERRQ(ierr); 41bdad233fSMatthew Knepley 42*73107ff1SLisandro Dalcin ierr = PetscHeaderCreate(t, TS_CLASSID, "TS", "Time stepping", "TS", comm, TSDestroy, TSView);CHKERRQ(ierr); 43bdad233fSMatthew Knepley 44bdad233fSMatthew Knepley /* General TS description */ 45089b2837SJed Brown t->problem_type = TS_NONLINEAR; 460298fd71SBarry Smith t->vec_sol = NULL; 47bdad233fSMatthew Knepley t->numbermonitors = 0; 480298fd71SBarry Smith t->snes = NULL; 49bdad233fSMatthew Knepley t->setupcalled = 0; 500298fd71SBarry Smith t->data = NULL; 510298fd71SBarry Smith t->user = NULL; 52cdbf8f93SLisandro Dalcin t->ptime = 0.0; 53cdbf8f93SLisandro Dalcin t->time_step = 0.1; 5431748224SBarry Smith t->time_step_orig = 0.1; 55bdad233fSMatthew Knepley t->max_time = 5.0; 568392e04aSShri Abhyankar t->steprollback = PETSC_FALSE; 57bdad233fSMatthew Knepley t->steps = 0; 58cdbf8f93SLisandro Dalcin t->max_steps = 5000; 595ef26d82SJed Brown t->ksp_its = 0; 605ef26d82SJed Brown t->snes_its = 0; 610298fd71SBarry Smith t->work = NULL; 62bdad233fSMatthew Knepley t->nwork = 0; 63193ac0bcSJed Brown t->max_snes_failures = 1; 64193ac0bcSJed Brown t->max_reject = 10; 65193ac0bcSJed Brown t->errorifstepfailed = PETSC_TRUE; 660e4ef248SJed Brown t->rhsjacobian.time = -1e20; 67e1244c69SJed Brown t->rhsjacobian.scale = 1.; 68b41af12eSJed Brown t->ijacobian.shift = 1.; 69e817cc15SEmil Constantinescu t->equation_type = TS_EQ_UNSPECIFIED; 70bdad233fSMatthew Knepley 711c3436cfSJed Brown t->atol = 1e-4; 721c3436cfSJed Brown t->rtol = 1e-4; 738d59e960SJed Brown t->cfltime = PETSC_MAX_REAL; 748d59e960SJed Brown t->cfltime_local = PETSC_MAX_REAL; 7549354f04SShri Abhyankar t->exact_final_time = TS_EXACTFINALTIME_STEPOVER; 76e144a568SJed Brown 77bdad233fSMatthew Knepley *ts = t; 78bdad233fSMatthew Knepley PetscFunctionReturn(0); 79bdad233fSMatthew Knepley } 80