fix deallocation of PC context

MODFLOW-USGS · Sep 26, 2023 · 4e0a8dc · 4e0a8dc
1 parent 765c080
commit 4e0a8dc
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 7 deletions.
diff --git a/src/Solution/PETSc/PetscImsPreconditioner.F90 b/src/Solution/PETSc/PetscImsPreconditioner.F90
@@ -74,12 +74,6 @@ subroutine pcshell_destroy(pc, ierr)
     PC :: pc !< the shell preconditioner
     PetscErrorCode :: ierr !< PETSc error code
     ! local
-    type(PcShellCtxType), pointer :: pc_ctx => null()
-
-    call PCShellGetContext(pc, pc_ctx, ierr)
-    CHKERRQ(ierr)
-    call VecDestroy(pc_ctx%diag, ierr)
-    CHKERRQ(ierr)
 
   end subroutine pcshell_destroy
 

diff --git a/src/Solution/PETSc/PetscSolver.F90 b/src/Solution/PETSc/PetscSolver.F90
@@ -423,7 +423,9 @@ subroutine petsc_destroy(this)
     CHKERRQ(ierr)
     deallocate (this%petsc_ctx)
 
-    allocate (this%pc_context)
+    call VecDestroy(this%pc_context%diag, ierr)
+    CHKERRQ(ierr)
+    deallocate (this%pc_context)
 
   end subroutine petsc_destroy