Dear all I get a segmentation fault when I call TSDestroy. Without TSDestroy the code runs fine. I have included portion of my code below.
subroutine runts(ctx) use userdata use comdata use mtsdata implicit none #include <petsc/finclude/petsc.h90> type(tsdata) :: ctx ! Local variables integer,parameter :: h = 100 ! File id for history file TS :: ts Vec :: u PetscErrorCode :: ierr external :: RHSFunction, Monitor call VecDuplicate(ctx%p%v_res, u, ierr); CHKERRQ(ierr) ! Copy initial condition into u call VecCopy(ctx%p%v_u, u, ierr); CHKERRQ(ierr) call TSCreate(PETSC_COMM_WORLD, ts, ierr); CHKERRQ(ierr) call TSSetProblemType(ts, TS_NONLINEAR, ierr); CHKERRQ(ierr) call TSSetRHSFunction(ts, PETSC_NULL_OBJECT, RHSFunction, ctx, ierr); CHKERRQ(ierr) call TSSetInitialTimeStep(ts, 0.0, dtg, ierr); CHKERRQ(ierr) call TSSetType(ts, TSRK, ierr); CHKERRQ(ierr); call TSSetDuration(ts, itmax, tfinal, ierr); CHKERRQ(ierr); call TSSetExactFinalTime(ts, TS_EXACTFINALTIME_MATCHSTEP, ierr); CHKERRQ(ierr); call TSMonitorSet(ts, Monitor, ctx, PETSC_NULL_OBJECT, ierr); CHKERRQ(ierr) call TSSetSolution(ts, u, ierr); CHKERRQ(ierr) call TSSetFromOptions(ts, ierr); CHKERRQ(ierr) call TSSetUp(ts, ierr); CHKERRQ(ierr) call TSSolve(ts, u, ierr); CHKERRQ(ierr) call VecCopy(u, ctx%p%v_u, ierr); CHKERRQ(ierr) call VecDestroy(u, ierr); CHKERRQ(ierr) call TSDestroy(ts, ierr); CHKERRQ(ierr) end subroutine runts Thanks praveen
