cmiss_petsc.f90

Go to the documentation of this file.

MODULE cmisspetsc
  
  USE base_routines
  USE cmisspetsctypes
  USE kinds
  USE iso_varying_string
  USE strings
  USE types
  
#include "macros.h"

  IMPLICIT NONE
 
  PRIVATE

#include "petscversion.h"
#include "petsc/finclude/petsc.h"
  
  !Module parameters

  !Insert mode types
  insertmode, PARAMETER :: petsc_insert_values = insert_values 
  insertmode, PARAMETER :: petsc_add_values = add_values 
  
  !Scatter mode types
  scattermode, PARAMETER :: petsc_scatter_forward = scatter_forward
  scattermode, PARAMETER :: petsc_scatter_reverse = scatter_reverse
  
  !KSP types
  ksptype, PARAMETER :: petsc_ksprichardson = ksprichardson 
  ksptype, PARAMETER :: petsc_kspchebyshev = kspchebyshev 
  ksptype, PARAMETER :: petsc_kspcg = kspcg 
  ksptype, PARAMETER :: petsc_kspcgne = kspcgne 
  ksptype, PARAMETER :: petsc_kspnash = kspnash
  ksptype, PARAMETER :: petsc_kspstcg = kspstcg
  ksptype, PARAMETER :: petsc_kspgltr = kspgltr
  ksptype, PARAMETER :: petsc_kspgmres = kspgmres 
  ksptype, PARAMETER :: petsc_kspfgmres = kspfgmres 
  ksptype, PARAMETER :: petsc_ksplgmres = ksplgmres
  ksptype, PARAMETER :: petsc_kspdgmres = kspdgmres
  ksptype, PARAMETER :: petsc_ksppgmres = ksppgmres
  ksptype, PARAMETER :: petsc_ksptcqmr = ksptcqmr 
  ksptype, PARAMETER :: petsc_kspbcgs = kspbcgs 
  ksptype, PARAMETER :: petsc_kspibcgs = kspibcgs
  ksptype, PARAMETER :: petsc_kspfbcgs = kspfbcgs
  ksptype, PARAMETER :: petsc_kspfbcgsr = kspfbcgsr
  ksptype, PARAMETER :: petsc_kspbcgsl = kspbcgsl
  ksptype, PARAMETER :: petsc_kspcgs = kspcgs 
  ksptype, PARAMETER :: petsc_ksptfqmr = ksptfqmr 
  ksptype, PARAMETER :: petsc_kspcr = kspcr 
  ksptype, PARAMETER :: petsc_ksplsqr = ksplsqr
  ksptype, PARAMETER :: petsc_ksppreonly = ksppreonly 
  ksptype, PARAMETER :: petsc_kspqcg = kspqcg
  ksptype, PARAMETER :: petsc_kspbicg = kspbicg 
  ksptype, PARAMETER :: petsc_kspminres = kspminres 
  ksptype, PARAMETER :: petsc_kspsymmlq = kspsymmlq
  ksptype, PARAMETER :: petsc_ksplcd = ksplcd 
  ksptype, PARAMETER :: petsc_ksppython = ksppython
  ksptype, PARAMETER :: petsc_kspgcr = kspgcr 
  
  !KSPConvergedReason types
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_rtol = ksp_converged_rtol
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_atol = ksp_converged_atol
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_its = ksp_converged_its
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_iterating = ksp_converged_iterating
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_cg_neg_curve = ksp_converged_cg_neg_curve
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_cg_constrained = ksp_converged_cg_constrained
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_step_length = ksp_converged_step_length
  kspconvergedreason, PARAMETER :: petsc_ksp_converged_happy_breakdown = ksp_converged_happy_breakdown
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_null = ksp_diverged_null
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_its = ksp_diverged_its
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_dtol = ksp_diverged_dtol
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_breakdown = ksp_diverged_breakdown
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_breakdown_bicg = ksp_diverged_breakdown_bicg
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_nonsymmetric = ksp_diverged_nonsymmetric
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_indefinite_pc = ksp_diverged_indefinite_pc
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_nanorinf = ksp_diverged_nanorinf
  kspconvergedreason, PARAMETER :: petsc_ksp_diverged_indefinite_mat = ksp_diverged_indefinite_mat

  !KSPNorm types
  kspnormtype, PARAMETER :: petsc_ksp_norm_none = ksp_norm_none
  kspnormtype, PARAMETER :: petsc_ksp_norm_preconditioned = ksp_norm_preconditioned
  kspnormtype, PARAMETER :: petsc_ksp_norm_unpreconditioned = ksp_norm_unpreconditioned
  kspnormtype, PARAMETER :: petsc_ksp_norm_natural = ksp_norm_natural

  !MatAssembly types
  matassemblytype, PARAMETER :: petsc_mat_flush_assembly = mat_flush_assembly
  matassemblytype, PARAMETER :: petsc_mat_final_assembly = mat_final_assembly

  !MatDuplicate types
  matduplicateoption, PARAMETER :: petsc_mat_do_not_copy_values = mat_do_not_copy_values
  matduplicateoption, PARAMETER :: petsc_mat_copy_values = mat_copy_values
  matduplicateoption, PARAMETER :: petsc_mat_share_nonzero_pattern = mat_share_nonzero_pattern

  !MatFactor types
  matfactortype, PARAMETER :: petsc_mat_factor_none = mat_factor_none
  matfactortype, PARAMETER :: petsc_mat_factor_lu = mat_factor_lu
  matfactortype, PARAMETER :: petsc_mat_factor_cholesky = mat_factor_cholesky
  matfactortype, PARAMETER :: petsc_mat_factor_ilu = mat_factor_ilu
  matfactortype, PARAMETER :: petsc_mat_factor_icc = mat_factor_icc

  !MatInfo types
  matinfo, PARAMETER :: petsc_mat_info_size = mat_info_size
  matinfo, PARAMETER :: petsc_mat_info_block_size = mat_info_block_size
  matinfo, PARAMETER :: petsc_mat_info_nz_allocated = mat_info_nz_allocated
  matinfo, PARAMETER :: petsc_mat_info_nz_used = mat_info_nz_used
  matinfo, PARAMETER :: petsc_mat_info_nz_unneeded = mat_info_nz_unneeded
  matinfo, PARAMETER :: petsc_mat_info_memory = mat_info_memory
  matinfo, PARAMETER :: petsc_mat_info_assemblies = mat_info_assemblies
  matinfo, PARAMETER :: petsc_mat_info_mallocs = mat_info_mallocs
  matinfo, PARAMETER :: petsc_mat_info_fill_ratio_given = mat_info_fill_ratio_given
  matinfo, PARAMETER :: petsc_mat_info_fill_ratio_needed = mat_info_fill_ratio_needed
  matinfo, PARAMETER :: petsc_mat_info_factor_mallocs = mat_info_factor_mallocs  

  !MatInfoType types
  matinfotype, PARAMETER :: petsc_mat_local = mat_local
  matinfotype, PARAMETER :: petsc_mat_global_max = mat_global_max
  matinfotype, PARAMETER :: petsc_mat_global_sum = mat_global_sum
  
  !MatOption types
  matoption, PARAMETER :: petsc_mat_row_oriented = mat_row_oriented 
  matoption, PARAMETER :: petsc_mat_new_nonzero_locations = mat_new_nonzero_locations 
  matoption, PARAMETER :: petsc_mat_symmetric = mat_symmetric 
  matoption, PARAMETER :: petsc_mat_structurally_symmetric = mat_structurally_symmetric 
  matoption, PARAMETER :: petsc_mat_new_diagonals = mat_new_diagonals
  matoption, PARAMETER :: petsc_mat_ignore_off_proc_entries = mat_ignore_off_proc_entries 
  matoption, PARAMETER :: petsc_mat_new_nonzero_location_err = mat_new_nonzero_location_err 
  matoption, PARAMETER :: petsc_mat_new_nonzero_allocation_err = mat_new_nonzero_allocation_err 
  matoption, PARAMETER :: petsc_mat_use_hash_table = mat_use_hash_table 
  matoption, PARAMETER :: petsc_mat_keep_nonzero_pattern = mat_keep_nonzero_pattern 
  matoption, PARAMETER :: petsc_mat_ignore_zero_entries = mat_ignore_zero_entries 
  matoption, PARAMETER :: petsc_mat_use_inodes = mat_use_inodes 
  matoption, PARAMETER :: petsc_mat_hermitian = mat_hermitian 
  matoption, PARAMETER :: petsc_mat_symmetry_eternal = mat_symmetry_eternal 
  matoption, PARAMETER :: petsc_mat_dummy = mat_dummy
  matoption, PARAMETER :: petsc_mat_ignore_lower_triangular = mat_ignore_lower_triangular 
  matoption, PARAMETER :: petsc_mat_error_lower_triangular = mat_error_lower_triangular
  matoption, PARAMETER :: petsc_mat_getrow_uppertriangular = mat_getrow_uppertriangular
  matoption, PARAMETER :: petsc_mat_unused_nonzero_location_err = mat_unused_nonzero_location_err
  matoption, PARAMETER :: petsc_num_mat_options = mat_option_max
  matoption, PARAMETER :: petsc_mat_spd = mat_spd 
  matoption, PARAMETER :: petsc_mat_no_off_proc_entries = mat_no_off_proc_entries
  matoption, PARAMETER :: petsc_mat_no_off_proc_zero_rows = mat_no_off_proc_zero_rows
  
  !Matrix Solver Package types
  matsolverpackage, PARAMETER :: petsc_mat_solver_superlu = matsolversuperlu
  matsolverpackage, PARAMETER :: petsc_mat_solver_superlu_dist = matsolversuperlu_dist
  matsolverpackage, PARAMETER :: petsc_mat_solver_umfpack = matsolverumfpack
  matsolverpackage, PARAMETER :: petsc_mat_solver_cholmod = matsolvercholmod
  matsolverpackage, PARAMETER :: petsc_mat_solver_essl = matsolveressl
  matsolverpackage, PARAMETER :: petsc_mat_solver_lusol = matsolverlusol
  matsolverpackage, PARAMETER :: petsc_mat_solver_mumps = matsolvermumps
  matsolverpackage, PARAMETER :: petsc_mat_solver_pastix = matsolverpastix
  matsolverpackage, PARAMETER :: petsc_mat_solver_matlab = matsolvermatlab
  matsolverpackage, PARAMETER :: petsc_mat_solver_petsc = matsolverpetsc
  matsolverpackage, PARAMETER :: petsc_mat_solver_bas = matsolverbas
  matsolverpackage, PARAMETER :: petsc_mat_solver_cusparse = matsolvercusparse
  matsolverpackage, PARAMETER :: petsc_mat_solver_bstrm = matsolverbstrm
  matsolverpackage, PARAMETER :: petsc_mat_solver_sbstrm = matsolversbstrm
  
  !MatStructure types
  matstructure, PARAMETER :: petsc_different_nonzero_pattern = different_nonzero_pattern
  matstructure, PARAMETER :: petsc_subset_nonzero_pattern = subset_nonzero_pattern
  matstructure, PARAMETER :: petsc_same_nonzero_pattern = same_nonzero_pattern

  !MatReuse types
  matreuse, PARAMETER :: petsc_mat_initial_matrix = mat_initial_matrix
  matreuse, PARAMETER :: petsc_mat_reuse_matrix = mat_reuse_matrix
  matreuse, PARAMETER :: petsc_mat_ignore_matrix = mat_ignore_matrix

  !MatColoring types
  matcoloringtype, PARAMETER :: petsc_matcoloring_natural = matcoloringnatural
  matcoloringtype, PARAMETER :: petsc_matcoloring_sl = matcoloringsl
  matcoloringtype, PARAMETER :: petsc_matcoloring_lf = matcoloringlf
  matcoloringtype, PARAMETER :: petsc_matcoloring_id = matcoloringid
  matcoloringtype, PARAMETER :: petsc_matcoloring_greedy = matcoloringgreedy
  matcoloringtype, PARAMETER :: petsc_matcoloring_jp = matcoloringjp

  !Norm types
  normtype, PARAMETER :: petsc_norm_1 = norm_1
  normtype, PARAMETER :: petsc_norm_2 = norm_2
  normtype, PARAMETER :: petsc_norm_infinity = norm_infinity
  
  !PC types
  pctype, PARAMETER ::  petsc_pcnone = pcnone
  pctype, PARAMETER ::  petsc_pcjacobi = pcjacobi
  pctype, PARAMETER ::  petsc_pcsor = pcsor
  pctype, PARAMETER ::  petsc_pclu = pclu
  pctype, PARAMETER ::  petsc_pcshell = pcshell
  pctype, PARAMETER ::  petsc_pcbjacobi = pcbjacobi
  pctype, PARAMETER ::  petsc_pcmg = pcmg
  pctype, PARAMETER ::  petsc_pceisenstat = pceisenstat
  pctype, PARAMETER ::  petsc_pcilu = pcilu
  pctype, PARAMETER ::  petsc_pcicc = pcicc
  pctype, PARAMETER ::  petsc_pcasm = pcasm
  pctype, PARAMETER ::  petsc_pcksp = pcksp
  pctype, PARAMETER ::  petsc_pccomposite = pccomposite
  pctype, PARAMETER ::  petsc_pcredundant = pcredundant
  pctype, PARAMETER ::  petsc_pcspai = pcspai
  pctype, PARAMETER ::  petsc_pcnn = pcnn
  pctype, PARAMETER ::  petsc_pccholesky = pccholesky
  pctype, PARAMETER ::  petsc_pcpbjacobi = pcpbjacobi
  pctype, PARAMETER ::  petsc_pcmat = pcmat
  pctype, PARAMETER ::  petsc_pchypre = pchypre
  pctype, PARAMETER ::  petsc_pcparms = pcparms
  pctype, PARAMETER ::  petsc_pcfieldsplit = pcfieldsplit
  pctype, PARAMETER ::  petsc_pctfs = pctfs
  pctype, PARAMETER ::  petsc_pcml = pcml
  pctype, PARAMETER ::  petsc_pcgalerkin = pcgalerkin
  pctype, PARAMETER ::  petsc_pcexotic = pcexotic
  pctype, PARAMETER ::  petsc_pcsupportgraph = pcsupportgraph
  pctype, PARAMETER ::  petsc_pccp = pccp
  pctype, PARAMETER ::  petsc_pcbfbt = pcbfbt
  pctype, PARAMETER ::  petsc_pclsc = pclsc
  pctype, PARAMETER ::  petsc_pcpython = pcpython
  pctype, PARAMETER ::  petsc_pcpfmg = pcpfmg
  pctype, PARAMETER ::  petsc_pcsyspfmg = pcsyspfmg
  pctype, PARAMETER ::  petsc_pcredistribute = pcredistribute
  pctype, PARAMETER ::  petsc_pcsvd = pcsvd
  pctype, PARAMETER ::  petsc_pcgamg = pcgamg
  pctype, PARAMETER ::  petsc_pcgasm = pcgasm
  pctype, PARAMETER ::  petsc_pcsacusp = pcsacusp
  pctype, PARAMETER ::  petsc_pcsacusppoly = pcsacusppoly
  pctype, PARAMETER ::  petsc_pcbicgstabcusp = pcbicgstabcusp
  pctype, PARAMETER ::  petsc_pcainvcusp = pcainvcusp
  pctype, PARAMETER ::  petsc_pcbddc = pcbddc

  !SNES types
  snestype, PARAMETER :: petsc_snesnewtonls = snesnewtonls
  snestype, PARAMETER :: petsc_snesnewtontr = snesnewtontr
  snestype, PARAMETER :: petsc_snespython = snespython
  snestype, PARAMETER :: petsc_snestest = snestest
  snestype, PARAMETER :: petsc_snesnrichardson = snesnrichardson
  snestype, PARAMETER :: petsc_snesksponly = snesksponly
  snestype, PARAMETER :: petsc_snesvinewtonrsls = snesvinewtonrsls
  snestype, PARAMETER :: petsc_snesvinewtonssls = snesvinewtonssls
  snestype, PARAMETER :: petsc_snesngmres = snesngmres
  snestype, PARAMETER :: petsc_snesqn = snesqn
  snestype, PARAMETER :: petsc_snesshell = snesshell
  snestype, PARAMETER :: petsc_snesncg = snesncg
  snestype, PARAMETER :: petsc_snesfas = snesfas
  snestype, PARAMETER :: petsc_snesms = snesms

  !SNES converged types
  snesconvergedreason, PARAMETER :: petsc_snes_converged_fnorm_abs = snes_converged_fnorm_abs
  snesconvergedreason, PARAMETER :: petsc_snes_converged_fnorm_relative = snes_converged_fnorm_relative
  snesconvergedreason, PARAMETER :: petsc_snes_converged_snorm_relative = snes_converged_snorm_relative
  snesconvergedreason, PARAMETER :: petsc_snes_converged_its = snes_converged_its
  snesconvergedreason, PARAMETER :: petsc_snes_converged_tr_delta = snes_converged_tr_delta
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_function_domain = snes_diverged_function_domain
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_function_count = snes_diverged_function_count
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_linear_solve = snes_diverged_linear_solve
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_fnorm_nan = snes_diverged_fnorm_nan
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_max_it = snes_diverged_max_it
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_line_search = snes_diverged_line_search
  snesconvergedreason, PARAMETER :: petsc_snes_diverged_local_min = snes_diverged_local_min
  snesconvergedreason, PARAMETER :: petsc_snes_converged_iterating = snes_converged_iterating
  
  !SNES line search type
  sneslinesearchtype, PARAMETER :: petsc_snes_linesearch_basic = sneslinesearchbasic
  sneslinesearchtype, PARAMETER :: petsc_snes_linesearch_bt = sneslinesearchbt
  sneslinesearchtype, PARAMETER :: petsc_snes_linesearch_l2 = sneslinesearchl2
  sneslinesearchtype, PARAMETER :: petsc_snes_linesearch_cp = sneslinesearchcp
  sneslinesearchtype, PARAMETER :: petsc_snes_linesearch_shell = sneslinesearchshell
  
  !SNES line search order  
  sneslinesearchorder, PARAMETER :: petsc_snes_linesearch_order_linear = snes_linesearch_order_linear
  sneslinesearchorder, PARAMETER :: petsc_snes_linesearch_order_quadratic = snes_linesearch_order_quadratic
  sneslinesearchorder, PARAMETER :: petsc_snes_linesearch_order_cubic = snes_linesearch_order_cubic

  !SNES norm schedules
  snesnormschedule, PARAMETER :: petsc_snes_norm_default = snes_norm_default
  snesnormschedule, PARAMETER :: petsc_snes_norm_none = snes_norm_none
  snesnormschedule, PARAMETER :: petsc_snes_norm_always = snes_norm_always
  snesnormschedule, PARAMETER :: petsc_snes_norm_initial_only = snes_norm_initial_only
  snesnormschedule, PARAMETER :: petsc_snes_norm_final_only = snes_norm_final_only
  snesnormschedule, PARAMETER :: petsc_snes_norm_initial_final_only = snes_norm_initial_final_only

  !SNES QN types
  snesqntype, PARAMETER :: petsc_snes_qn_lbfgs = snes_qn_lbfgs 
  snesqntype, PARAMETER :: petsc_snes_qn_broyden = snes_qn_broyden
  snesqntype, PARAMETER :: petsc_snes_qn_badbroyden = snes_qn_badbroyden
  
  !SNES QN restart types
  snesqnrestarttype, PARAMETER :: petsc_snes_qn_restart_default = snes_qn_restart_default
  snesqnrestarttype, PARAMETER :: petsc_snes_qn_restart_none = snes_qn_restart_none 
  snesqnrestarttype, PARAMETER :: petsc_snes_qn_restart_powell = snes_qn_restart_powell
  snesqnrestarttype, PARAMETER :: petsc_snes_qn_restart_periodic = snes_qn_restart_periodic 

  !SNES QN scaling types
  snesqnscaletype, PARAMETER :: petsc_snes_qn_scale_default = snes_qn_scale_default
  snesqnscaletype, PARAMETER :: petsc_snes_qn_scale_none = snes_qn_scale_none 
  snesqnscaletype, PARAMETER :: petsc_snes_qn_scale_shanno = snes_qn_scale_shanno 
  snesqnscaletype, PARAMETER :: petsc_snes_qn_scale_linesearch = snes_qn_scale_linesearch  
  snesqnscaletype, PARAMETER :: petsc_snes_qn_scale_jacobian = snes_qn_scale_jacobian   

  !TS types
  tstype, PARAMETER :: petsc_ts_euler = tseuler
  tstype, PARAMETER :: petsc_ts_beuler = tsbeuler
  tstype, PARAMETER :: petsc_ts_pseudo = tspseudo
  tstype, PARAMETER :: petsc_ts_cn = tscn
  tstype, PARAMETER :: petsc_ts_sundials = tssundials
  tstype, PARAMETER :: petsc_ts_rk = tsrk
  tstype, PARAMETER :: petsc_ts_python = tspython
  tstype, PARAMETER :: petsc_ts_theta = tstheta
  tstype, PARAMETER :: petsc_ts_alpha = tsalpha
  tstype, PARAMETER :: petsc_ts_gl = tsgl
  tstype, PARAMETER :: petsc_ts_ssp = tsssp
  tstype, PARAMETER :: petsc_ts_arkimex = tsarkimex
  tstype, PARAMETER :: petsc_ts_rosw = tsrosw
  tstype, PARAMETER :: petsc_ts_eimex = tseimex

  !TS convergence flags
  tsconvergedreason, PARAMETER :: petsc_ts_converged_iterating = ts_converged_iterating
  tsconvergedreason, PARAMETER :: petsc_ts_converged_time = ts_converged_time
  tsconvergedreason, PARAMETER :: petsc_ts_converged_its = ts_converged_its
  tsconvergedreason, PARAMETER :: petsc_ts_diverged_nonlinear_solve = ts_diverged_nonlinear_solve
  tsconvergedreason, PARAMETER :: petsc_ts_diverged_step_rejected = ts_diverged_step_rejected
  
  !TS problem types
  tsproblemtype, PARAMETER :: petsc_ts_linear = ts_linear
  tsproblemtype, PARAMETER :: petsc_ts_nonlinear = ts_nonlinear
  
  !TS Sundials types
  tssundialstype, PARAMETER :: petsc_sundials_adams = sundials_adams
  tssundialstype, PARAMETER :: petsc_sundials_bdf = sundials_bdf

  !TS Sundials Gram Schmidt Type
  tssundialsgramschmidttype, PARAMETER :: petsc_sundials_modified_gs = sundials_modified_gs
  tssundialsgramschmidttype, PARAMETER :: petsc_sundials_classical_gs = sundials_classical_gs
  
  !Module types

  !Module variables

  LOGICAL, SAVE :: petschandleerror

  !Interfaces

  INTERFACE

    !PETSc miscellanous routines

    SUBROUTINE petscfinalize(ierr)
      petscint ierr
    END SUBROUTINE petscfinalize

    SUBROUTINE petscinitialize(file,ierr)
      CHARACTER(LEN=*) file
      petscint ierr
    END SUBROUTINE petscinitialize

    SUBROUTINE petscpopsignalhandler(ierr)
      petscint ierr
    END SUBROUTINE petscpopsignalhandler
    
    SUBROUTINE petsclogview(viewer,ierr)
      petscviewer viewer
      petscint ierr
    END SUBROUTINE petsclogview

    !IS routines
    
    SUBROUTINE isdestroy(indexset,ierr)
      is indexset
      petscint ierr
    END SUBROUTINE isdestroy

    !IS coloring routines

    SUBROUTINE iscoloringdestroy(iscoloring,ierr)
      iscoloring iscoloring
      petscint ierr
    END SUBROUTINE iscoloringdestroy

    !IS local to global mapping routines
    
    SUBROUTINE islocaltoglobalmappingapply(islocaltoglobalmapping,n,idxin,idxout,ierr)
      islocaltoglobalmapping islocaltoglobalmapping
      petscint n
      petscint idxin(*)
      petscint idxout(*)
      petscint ierr
    END SUBROUTINE islocaltoglobalmappingapply
    
    SUBROUTINE islocaltoglobalmappingapplyis(islocaltoglobalmapping,isin,isout,ierr)
      islocaltoglobalmapping islocaltoglobalmapping
      is isin
      is isout
      petscint ierr
    END SUBROUTINE islocaltoglobalmappingapplyis
    
    SUBROUTINE islocaltoglobalmappingcreate(comm,blockSize,n,indices,mode,islocaltoglobalmapping,ierr)
      mpi_comm comm
      petscint blocksize
      petscint n
      petscint indices(*)
      petsccopymode mode
      islocaltoglobalmapping islocaltoglobalmapping
      petscint ierr
    END SUBROUTINE islocaltoglobalmappingcreate
    
    SUBROUTINE islocaltoglobalmappingdestroy(islocaltoglobalmapping,ierr)
      islocaltoglobalmapping islocaltoglobalmapping
      petscint ierr
    END SUBROUTINE islocaltoglobalmappingdestroy

    !KSP routines
    
    SUBROUTINE kspcreate(comm,ksp,ierr)
      mpi_comm comm
      ksp ksp
      petscint ierr
    END SUBROUTINE kspcreate
    
    SUBROUTINE kspdestroy(ksp,ierr)
      ksp ksp
      petscint ierr
    END SUBROUTINE kspdestroy
    
    SUBROUTINE kspgetconvergedreason(ksp,reason,ierr)
      ksp ksp
      kspconvergedreason reason
      petscint ierr
    END SUBROUTINE kspgetconvergedreason
    
    SUBROUTINE kspgetiterationnumber(ksp,its,ierr)
      ksp ksp
      petscint its
      petscint ierr
    END SUBROUTINE kspgetiterationnumber
    
    SUBROUTINE kspgetpc(ksp,pc,ierr)
      ksp ksp
      pc pc
      petscint ierr
    END SUBROUTINE kspgetpc
    
    SUBROUTINE kspgetresidualnorm(ksp,rnorm,ierr)
      ksp ksp
      petscreal rnorm
      petscint ierr
    END SUBROUTINE kspgetresidualnorm
    
    SUBROUTINE kspgmressetrestart(ksp,restart,ierr)
      ksp ksp
      petscint restart
      petscint ierr
    END SUBROUTINE kspgmressetrestart
    
    SUBROUTINE kspsetfromoptions(ksp,ierr)
      ksp ksp
      petscint ierr
    END SUBROUTINE kspsetfromoptions
    
    SUBROUTINE kspsetinitialguessnonzero(ksp,flg,ierr)
      ksp ksp
      petscbool flg
      petscint ierr
    END SUBROUTINE kspsetinitialguessnonzero
    
    SUBROUTINE kspsetoperators(ksp,amat,pmat,ierr)
      ksp ksp
      mat amat
      mat pmat
      petscint ierr
    END SUBROUTINE kspsetoperators
    
    SUBROUTINE kspsetreusepreconditioner(ksp,flag,ierr)
      ksp ksp
      petscbool flag
      petscint ierr
    END SUBROUTINE kspsetreusepreconditioner
    
    SUBROUTINE kspsettolerances(ksp,rtol,atol,dtol,maxits,ierr)
      ksp ksp
      petscreal rtol
      petscreal atol
      petscreal dtol
      petscint maxits
      petscint ierr
    END SUBROUTINE kspsettolerances
    
    SUBROUTINE kspsettype(ksp,method,ierr)
      ksp ksp
      ksptype method
      petscint ierr
    END SUBROUTINE kspsettype
    
    SUBROUTINE kspsetup(ksp,ierr)
      ksp ksp
      petscint ierr
    END SUBROUTINE kspsetup
    
    SUBROUTINE kspsolve(ksp,b,x,ierr)
      ksp ksp
      vec b
      vec x
      petscint ierr
    END SUBROUTINE kspsolve

    !Matrix routines
    
    SUBROUTINE matassemblybegin(A,assemblytype,ierr)
      mat a
      matassemblytype assemblytype
      petscint ierr
    END SUBROUTINE matassemblybegin
    
    SUBROUTINE matassemblyend(A,assemblytype,ierr)
      mat a
      matassemblytype assemblytype
      petscint ierr
    END SUBROUTINE matassemblyend

    SUBROUTINE matcreate(comm,A,ierr)
      mpi_comm comm
      mat a
      petscint ierr
    END SUBROUTINE matcreate

    SUBROUTINE matcreateaij(comm,localm,localn,globalm,globaln,diagnumbernzperrow,diagnumbernzeachrow,offdiagnumbernzperrow, &
      & offdiagnumbernzeachrow,a,ierr)
      mpi_comm comm
      petscint localm
      petscint localn
      petscint globalm
      petscint globaln
      petscint diagnumbernzperrow
      petscint diagnumbernzeachrow(*)
      petscint offdiagnumbernzperrow
      petscint offdiagnumbernzeachrow(*)
      mat a
      petscint ierr
    END SUBROUTINE matcreateaij
        
    SUBROUTINE matcreatedense(comm,localm,localn,globalm,globaln,matrixdata,A,ierr)
      mpi_comm comm
      petscint localm
      petscint localn
      petscint globalm
      petscint globaln
      petscscalar matrixdata(*)
      mat a
      petscint ierr
    END SUBROUTINE matcreatedense
        
    SUBROUTINE matcreateseqaij(comm,m,n,numbernzperrow,numbernzeachrow,A,ierr)
      mpi_comm comm
      petscint m
      petscint n
      petscint numbernzperrow
      petscint numbernzeachrow(*)
      mat a
      petscint ierr
    END SUBROUTINE matcreateseqaij
          
    SUBROUTINE matcreateseqdense(comm,m,n,matrixdata,A,ierr)
      mpi_comm comm
      petscint m
      petscint n
      petscscalar matrixdata(*)
      mat a
      petscint ierr
    END SUBROUTINE matcreateseqdense

    SUBROUTINE matdensegetarrayf90(A,array,ierr)
      mat a
      petscscalar, POINTER :: array(:,:)
      petscint ierr
    END SUBROUTINE matdensegetarrayf90

    SUBROUTINE matdenserestorearrayf90(A,array,ierr)
      mat a
      petscscalar, POINTER :: array(:,:)
      petscint ierr
    END SUBROUTINE matdenserestorearrayf90

    SUBROUTINE matdestroy(A,ierr)
      mat a
      petscint ierr
    END SUBROUTINE matdestroy
    
    SUBROUTINE matgetinfo(A,flag,info,ierr)
      mat a
      matinfotype flag
      matinfo info(*)
      petscint ierr
    END SUBROUTINE matgetinfo
    
    SUBROUTINE matgetownershiprange(A,firstrow,lastrow,ierr)
      mat a
      petscint firstrow
      petscint lastrow
      petscint ierr
    END SUBROUTINE matgetownershiprange

    SUBROUTINE matgetrow(A,row,ncols,cols,values,ierr)
      mat a
      petscint row
      petscint ncols
      petscint cols(*)
      petscscalar values(*)
      petscint ierr
    END SUBROUTINE matgetrow
    
    SUBROUTINE matgetvalues(A,m,idxm,n,idxn,values,ierr)
      mat a
      petscint m
      petscint idxm(*)
      petscint n
      petscint idxn(*)
      petscscalar values(*)
      petscint ierr
    END SUBROUTINE matgetvalues
    
    SUBROUTINE matmumpsseticntl(A,icntl,ival,ierr)
      mat a
      petscint icntl
      petscint ival
      petscint ierr
    END SUBROUTINE matmumpsseticntl

    SUBROUTINE matmumpssetcntl(A,icntl,val,ierr)
      mat a
      petscint icntl
      petscreal val
      petscint ierr
    END SUBROUTINE matmumpssetcntl

    SUBROUTINE matrestorerow(A,row,ncols,cols,values,ierr)
      mat a
      petscint row
      petscint ncols
      petscint cols(*)
      petscscalar values(*)
      petscint ierr
    END SUBROUTINE matrestorerow
    
    SUBROUTINE matseqaijgetarrayf90(A,array,ierr)
      mat a
      petscscalar, POINTER :: array(:,:)
      petscint ierr
    END SUBROUTINE matseqaijgetarrayf90

    SUBROUTINE matseqaijgetmaxrownonzeros(A,maxNumberNonZeros,ierr)
      mat a
      petscint maxnumbernonzeros
      petscint ierr
    END SUBROUTINE matseqaijgetmaxrownonzeros

    SUBROUTINE matseqaijrestorearrayf90(A,array,ierr)
      mat a
      petscscalar, POINTER :: array(:,:)
      petscint ierr
    END SUBROUTINE matseqaijrestorearrayf90

    SUBROUTINE matsetlocaltoglobalmapping(A,ctx,ierr)
      mat a
      islocaltoglobalmapping ctx
      petscint ierr
    END SUBROUTINE matsetlocaltoglobalmapping

    SUBROUTINE matsetoption(A,option,flag,ierr)
      mat a
      matoption option
      petscbool flag
      petscint ierr
    END SUBROUTINE matsetoption

    SUBROUTINE matsetsizes(A,localm,localn,globalM,globalN,ierr)
      mat a
      petscint localm
      petscint localn
      petscint globalm
      petscint globaln
      petscint ierr
    END SUBROUTINE matsetsizes
    
    SUBROUTINE matsettype(A,matrixtype,ierr)
      mat a
      mattype matrixtype
      petscint ierr
    END SUBROUTINE matsettype

    SUBROUTINE matsetvalue(A,row,col,value,insertmode,ierr)
      mat a
      petscint row
      petscint col
      petscscalar value
      insertmode insertmode
      petscint ierr
    END SUBROUTINE matsetvalue
    
    SUBROUTINE matsetvalues(A,m,mindices,n,nindices,values,insertmode,ierr)
      mat a
      petscint m
      petscint mindices(*)
      petscint n
      petscint nindices(*)
      petscscalar values(*)
      insertmode insertmode
      petscint ierr
    END SUBROUTINE matsetvalues
    
    SUBROUTINE matsetvaluelocal(A,row,col,value,insertmode,ierr)
      mat a
      petscint row
      petscint col
      petscscalar value
      insertmode insertmode
      petscint ierr
    END SUBROUTINE matsetvaluelocal
    
    SUBROUTINE matsetvalueslocal(A,m,mindices,n,nindices,values,insertmode,ierr)
      mat a
      petscint m
      petscint mindices(*)
      petscint n
      petscint nindices(*)
      petscscalar values(*)
      insertmode insertmode
      petscint ierr
    END SUBROUTINE matsetvalueslocal
    
    SUBROUTINE matview(A,v,ierr)
      mat a
      petscviewer v
      petscint ierr
    END SUBROUTINE matview

    SUBROUTINE matzeroentries(A,ierr)
      mat a
      petscint ierr
    END SUBROUTINE matzeroentries

    !Mat coloring routines

    SUBROUTINE matcoloringapply(coloring,isColoring,ierr)
      matcoloring coloring
      iscoloring iscoloring
      petscint ierr
    END SUBROUTINE matcoloringapply
    
    SUBROUTINE matcoloringcreate(A,coloring,ierr)
      mat a
      matcoloring coloring
      petscint ierr
    END SUBROUTINE matcoloringcreate
    
    SUBROUTINE matcoloringdestroy(coloring,ierr)
      matcoloring coloring
      petscint ierr
    END SUBROUTINE matcoloringdestroy
    
    SUBROUTINE matcoloringsetfromoptions(coloring,ierr)
      matcoloring coloring
      petscint ierr
    END SUBROUTINE matcoloringsetfromoptions
    
    SUBROUTINE matcoloringsettype(coloring,coloringType,ierr)
      matcoloring coloring
      matcoloringtype coloringtype
      petscint ierr
    END SUBROUTINE matcoloringsettype

    !Mat FD coloring routines

    SUBROUTINE matfdcoloringcreate(A,iscoloring,fdcoloring,ierr)
      mat a
      iscoloring iscoloring
      matfdcoloring fdcoloring
      petscint ierr
    END SUBROUTINE matfdcoloringcreate
    
    SUBROUTINE matfdcoloringdestroy(fdcoloring,ierr)
      matfdcoloring fdcoloring
      petscint ierr
    END SUBROUTINE matfdcoloringdestroy
    
    SUBROUTINE matfdcoloringsetfromoptions(fdcoloring,ierr)
      matfdcoloring fdcoloring
      petscint ierr
    END SUBROUTINE matfdcoloringsetfromoptions

    SUBROUTINE matfdcoloringsetfunction(fdcoloring,ffunction,ctx,ierr)
      USE types
      matfdcoloring fdcoloring
      EXTERNAL ffunction
      TYPE(solver_type), POINTER :: ctx
      petscint ierr
    END SUBROUTINE matfdcoloringsetfunction
    
    SUBROUTINE matfdcoloringsetparameters(fdcoloring,rerror,umin,ierr)
      matfdcoloring fdcoloring
      petscscalar rerror
      petscscalar umin
      petscint ierr
    END SUBROUTINE matfdcoloringsetparameters
    
    SUBROUTINE matfdcoloringsetup(A,iscoloring,fdcoloring,ierr)
      mat a
      iscoloring iscoloring
      matfdcoloring fdcoloring
      petscint ierr
    END SUBROUTINE matfdcoloringsetup
    
    !Pre-conditioner routines

    SUBROUTINE pcfactorgetmatrix(pc,A,ierr)
      pc pc
      mat a
      petscint ierr
    END SUBROUTINE pcfactorgetmatrix

    SUBROUTINE pcfactorsetmatsolverpackage(pc,solverpackage,ierr)
      pc pc
      matsolverpackage solverpackage
      petscint ierr
    END SUBROUTINE pcfactorsetmatsolverpackage

    SUBROUTINE pcfactorsetupmatsolverpackage(pc,ierr)
      pc pc
      petscint ierr
    END SUBROUTINE pcfactorsetupmatsolverpackage

    SUBROUTINE pcsetfromoptions(pc,ierr)
      pc pc
      petscint ierr
    END SUBROUTINE pcsetfromoptions
    
    SUBROUTINE pcsetreusepreconditioner(pc,flag,ierr)
      pc pc
      petscbool flag
      petscint ierr
    END SUBROUTINE pcsetreusepreconditioner
    
    SUBROUTINE pcsettype(pc,method,ierr)
      pc pc
      pctype method
      petscint ierr
    END SUBROUTINE pcsettype
    
    !SNES routines

    SUBROUTINE snescreate(comm,snes,ierr)
      mpi_comm comm
      snes snes
      petscint ierr
    END SUBROUTINE snescreate

    SUBROUTINE snesdestroy(snes,ierr)
      snes snes
      petscint ierr
    END SUBROUTINE snesdestroy

    SUBROUTINE snesgetapplicationcontext(snes,ctx,ierr)
      USE types
      snes snes
      TYPE(solver_type), POINTER :: ctx
      petscint ierr
    END SUBROUTINE snesgetapplicationcontext

    SUBROUTINE snesgetconvergedreason(snes,reason,ierr)
      snes snes
      snesconvergedreason reason
      petscint ierr
    END SUBROUTINE snesgetconvergedreason
    
    SUBROUTINE snesgetfunction(snes,f,ffunction,ctx,ierr)
      USE types
      snes snes
      vec f
      EXTERNAL ffunction
      petscint ctx
      petscint ierr
    END SUBROUTINE snesgetfunction

    SUBROUTINE snesgetiterationnumber(snes,iter,ierr)
      snes snes
      petscint iter
      petscint ierr
    END SUBROUTINE snesgetiterationnumber

    SUBROUTINE snesgetjacobian(snes,A,B,Jfunction,ctx,ierr)
      USE types
      snes snes
      mat a
      mat b      
      EXTERNAL jfunction
      petscint ctx
      petscint ierr
    END SUBROUTINE snesgetjacobian

    SUBROUTINE snesgetksp(snes,ksp,ierr)
      snes snes
      ksp ksp
      petscint ierr
    END SUBROUTINE snesgetksp

    SUBROUTINE snesgetlinesearch(snes,linesearch,ierr)
      snes snes
      sneslinesearch linesearch
      petscint ierr
    END SUBROUTINE snesgetlinesearch

    SUBROUTINE snesgetsolutionupdate(snes,solutionUpdate,ierr)
      snes snes
      vec solutionupdate
      petscint ierr
    END SUBROUTINE snesgetsolutionupdate

    SUBROUTINE snesmonitorset(snes,mfunction,mctx,monitordestroy,ierr)
      USE types
      snes snes
      EXTERNAL mfunction
      TYPE(solver_type), POINTER :: mctx
      EXTERNAL monitordestroy
      petscint ierr
    END SUBROUTINE snesmonitorset

    SUBROUTINE snesqnsetrestarttype(snes,rtype,ierr)
      snes snes
      snesqnrestarttype rtype
      petscint ierr
    END SUBROUTINE snesqnsetrestarttype

    SUBROUTINE snesqnsetscaletype(snes,stype,ierr)
      snes snes
      snesqnscaletype stype
      petscint ierr
    END SUBROUTINE snesqnsetscaletype

    SUBROUTINE snesqnsettype(snes,qtype,ierr)
      snes snes
      snesqntype qtype
      petscint ierr
    END SUBROUTINE snesqnsettype

    SUBROUTINE snessetapplicationcontext(snes,ctx,ierr)
      USE types
      snes snes
      TYPE(solver_type), POINTER :: ctx
      petscint ierr
    END SUBROUTINE snessetapplicationcontext

    SUBROUTINE snessetconvergencetest(snes,cfunction,ctx,destroyFunction,ierr)
      USE types
      snes snes
      EXTERNAL cfunction
      TYPE(solver_type), POINTER :: ctx
      EXTERNAL destroyfunction
      petscint ierr
    END SUBROUTINE snessetconvergencetest

    SUBROUTINE snessetfromoptions(snes,ierr)
      snes snes
      petscint ierr
    END SUBROUTINE snessetfromoptions

    SUBROUTINE snessetfunction(snes,f,ffunction,ctx,ierr)
      USE types
      snes snes
      vec f
      EXTERNAL ffunction
      TYPE(solver_type), POINTER :: ctx
      petscint ierr
    END SUBROUTINE snessetfunction

    SUBROUTINE snessetjacobian(snes,A,B,Jfunction,ctx,ierr)
      USE types
      snes snes
      mat a
      mat b      
      EXTERNAL jfunction
      TYPE(solver_type), POINTER :: ctx
      petscint ierr
    END SUBROUTINE snessetjacobian

    SUBROUTINE snessetksp(snes,ksp,ierr)
      snes snes
      ksp ksp
      petscint ierr
    END SUBROUTINE snessetksp

    SUBROUTINE snessetnormschedule(snes,normschedule,ierr)
      snes snes
      snesnormschedule normschedule
      petscint ierr
    END SUBROUTINE snessetnormschedule

    SUBROUTINE snessettolerances(snes,abstol,rtol,stol,maxit,maxf,ierr)
      snes snes
      petscreal abstol
      petscreal rtol
      petscreal stol
      petscint maxit
      petscint maxf
      petscint ierr
    END SUBROUTINE snessettolerances
    
    SUBROUTINE snessettrustregiontolerance(snes,trtol,ierr)
      snes snes
      petscreal trtol
      petscint ierr
    END SUBROUTINE snessettrustregiontolerance

    SUBROUTINE snessettype(snes,method,ierr)
      snes snes
      snestype method
      petscint ierr
    END SUBROUTINE snessettype

    SUBROUTINE snessolve(snes,b,x,ierr)
      snes snes
      vec b
      vec x
      petscint ierr
    END SUBROUTINE snessolve

    !SNES line search routines
    
    SUBROUTINE sneslinesearchbtsetalpha(linesearch,alpha,ierr)
      sneslinesearch linesearch
      petscreal alpha
      petscint ierr
    END SUBROUTINE sneslinesearchbtsetalpha

     SUBROUTINE sneslinesearchcomputenorms(linesearch,ierr)
      sneslinesearch linesearch
      petscint ierr
    END SUBROUTINE sneslinesearchcomputenorms

    SUBROUTINE sneslinesearchgetnorms(linesearch,xnorm,fnorm,ynorm,ierr)
      sneslinesearch linesearch
      petscreal xnorm
      petscreal fnorm
      petscreal ynorm
      petscint ierr
    END SUBROUTINE sneslinesearchgetnorms

    SUBROUTINE sneslinesearchgetvecs(linesearch,x,f,y,w,g,ierr)
      sneslinesearch linesearch
      vec x
      vec f
      vec y
      vec w
      vec g
      petscint ierr
    END SUBROUTINE sneslinesearchgetvecs

    SUBROUTINE sneslinesearchsetcomputenorms(linesearch,flag,ierr)
      sneslinesearch linesearch
      petscbool flag
      petscint ierr
    END SUBROUTINE sneslinesearchsetcomputenorms

    SUBROUTINE sneslinesearchsetmonitor(linesearch,flag,ierr)
      sneslinesearch linesearch
      petscbool flag
      petscint ierr
    END SUBROUTINE sneslinesearchsetmonitor

    SUBROUTINE sneslinesearchsetnorms(snes,xnorm,fnorm,ynorm,ierr)
      snes snes
      petscreal xnorm
      petscreal fnorm
      petscreal ynorm
      petscint ierr
    END SUBROUTINE sneslinesearchsetnorms

    SUBROUTINE sneslinesearchsetorder(linesearch,linesearchorder,ierr)
      sneslinesearch linesearch
      sneslinesearchorder linesearchorder
      petscint ierr
    END SUBROUTINE sneslinesearchsetorder

    SUBROUTINE sneslinesearchsettolerances(linesearch,steptol,maxstep,rtol,atol,ltol,maxIt,ierr)
      sneslinesearch linesearch
      petscreal steptol
      petscreal maxstep
      petscreal rtol
      petscreal atol
      petscreal ltol
      petscint maxit
      petscint ierr
    END SUBROUTINE sneslinesearchsettolerances
    
    SUBROUTINE sneslinesearchsettype(linesearch,linesearchtype,ierr)
      sneslinesearch linesearch
      sneslinesearchtype linesearchtype
      petscint ierr
    END SUBROUTINE sneslinesearchsettype
    
    !Time stepping routines
    
    SUBROUTINE tscreate(comm,ts,ierr)
      mpi_comm comm
      ts ts
      petscint ierr
    END SUBROUTINE tscreate

    SUBROUTINE tsdestroy(ts,ierr)
      ts ts
      petscint ierr
    END SUBROUTINE tsdestroy

    SUBROUTINE tsgetsolution(ts,currentsolution,ierr)
      ts ts
      vec currentsolution
      petscint ierr
    END SUBROUTINE tsgetsolution
    
    SUBROUTINE tsmonitorset(ts,mfunction,mctx,monitordestroy,ierr)
      USE types
      ts ts
      EXTERNAL mfunction
      TYPE(solver_type), POINTER :: mctx
      EXTERNAL monitordestroy
      petscint ierr
    END SUBROUTINE tsmonitorset

    SUBROUTINE tssetduration(ts,maxsteps,maxtime,ierr)
      ts ts
      petscint maxsteps
      petscreal maxtime
      petscint ierr
    END SUBROUTINE tssetduration

    SUBROUTINE tssetexactfinaltime(ts,eftopt,ierr)
      ts ts
      petscbool eftopt
      petscint ierr
    END SUBROUTINE tssetexactfinaltime

    SUBROUTINE tssetfromoptions(ts,ierr)
      ts ts
      petscint ierr
    END SUBROUTINE tssetfromoptions

    SUBROUTINE tssetinitialtimestep(ts,initial_time,time_step,ierr)
      ts ts
      petscreal initial_time
      petscreal time_step
      petscint ierr
    END SUBROUTINE tssetinitialtimestep
    
    SUBROUTINE tssetproblemtype(ts,probtype,ierr)
      ts ts
      tsproblemtype probtype
      petscint ierr
    END SUBROUTINE tssetproblemtype
    
    SUBROUTINE tssetrhsfunction(ts,r,rhsfunc,ctx,ierr)
      USE types
      ts ts
      vec r
      EXTERNAL rhsfunc
      TYPE(cellmlpetsccontexttype), POINTER :: ctx
      petscint ierr
    END SUBROUTINE tssetrhsfunction

    SUBROUTINE tssetsolution(ts,initialsolution,ierr)
      ts ts
      vec initialsolution
      petscint ierr
    END SUBROUTINE tssetsolution

    SUBROUTINE tssettimestep(ts,time_step,ierr)
      ts ts
      petscreal time_step
      petscint ierr
    END SUBROUTINE tssettimestep
    
    SUBROUTINE tssettype(ts,tstype,ierr)
      ts ts
      tstype tstype
      petscint ierr
    END SUBROUTINE tssettype
    
    SUBROUTINE tssolve(ts,x,ftime,ierr)
      ts ts
      vec x
      petscreal ftime
      petscint ierr
    END SUBROUTINE tssolve

    SUBROUTINE tsstep(ts,steps,ptime,ierr)
      ts ts
      petscint steps
      petscreal ptime
      petscint ierr
    END SUBROUTINE tsstep

    SUBROUTINE tssundialssettolerance(ts,abstol,reltol,ierr)
      ts ts
      petscreal abstol
      petscreal reltol
      petscint ierr
    END SUBROUTINE tssundialssettolerance

    SUBROUTINE tssundialssettype(ts,sundialstype,ierr)
      ts ts
      tssundialstype sundialstype
      petscint ierr
    END SUBROUTINE tssundialssettype

    !Vector routines

    SUBROUTINE vecassemblybegin(x,ierr)
      vec x
      petscint ierr
    END SUBROUTINE vecassemblybegin

    SUBROUTINE vecassemblyend(x,ierr)
      vec x
      petscint ierr
    END SUBROUTINE vecassemblyend

    SUBROUTINE veccopy(x,y,ierr)
      vec x
      vec y
      petscint ierr
    END SUBROUTINE veccopy

    SUBROUTINE veccreate(comm,x,ierr)
      mpi_comm comm
      vec x
      petscint ierr
    END SUBROUTINE veccreate

    SUBROUTINE veccreateghost(comm,localm,globalm,nghost,ghosts,x,ierr)
      mpi_comm comm
      petscint localm
      petscint globalm
      petscint nghost
      petscint ghosts(*)
      vec x
      petscint ierr
    END SUBROUTINE veccreateghost

    SUBROUTINE veccreateghostwitharray(comm,localm,globalm,nghost,ghosts,array,x,ierr)
      mpi_comm comm
      petscint localm
      petscint globalm
      petscint nghost
      petscint ghosts(*)
      petscscalar array(*)
      vec x
      petscint ierr
    END SUBROUTINE veccreateghostwitharray

    SUBROUTINE veccreatempi(comm,localm,globalm,x,ierr)
      mpi_comm comm
      petscint localm
      petscint globalm
      vec x
      petscint ierr
    END SUBROUTINE veccreatempi

    SUBROUTINE veccreatempiwitharray(comm,localn,globaln,array,x,ierr)
      mpi_comm comm
      petscint localn
      petscint globaln
      petscscalar array(*)
      vec x
      petscint ierr
    END SUBROUTINE veccreatempiwitharray

    SUBROUTINE veccreateseq(comm,m,x,ierr)
      mpi_comm comm
      petscint m
      vec x
      petscint ierr
    END SUBROUTINE veccreateseq

    SUBROUTINE veccreateseqwitharray(comm,n,array,x,ierr)
      mpi_comm comm
      petscint n
      petscscalar array(*)
      vec x
      petscint ierr
    END SUBROUTINE veccreateseqwitharray

    SUBROUTINE vecdestroy(x,ierr)
      vec x
      petscint ierr
    END SUBROUTINE vecdestroy

    SUBROUTINE vecduplicate(old,new,ierr)
      vec old,new
      petscint ierr
    END SUBROUTINE vecduplicate

    SUBROUTINE vecdot(x,y,val,ierr)
      vec x
      vec y
      petscscalar val
      petscint ierr
    END SUBROUTINE vecdot

    SUBROUTINE vecgetarrayf90(x,vec_data,ierr)
      vec x
      petscscalar, POINTER :: vec_data(:)
      petscint ierr
    END SUBROUTINE vecgetarrayf90

    SUBROUTINE vecgetarrayreadf90(x,vec_data,ierr)
      vec x
      petscscalar, POINTER :: vec_data(:)
      petscint ierr
    END SUBROUTINE vecgetarrayreadf90

    SUBROUTINE vecgetlocalsize(x,size,ierr)
      vec x
      petscint size
      petscint ierr
    END SUBROUTINE vecgetlocalsize

    SUBROUTINE vecgetownershiprange(x,low,high,ierr)
      vec x
      petscint low
      petscint high
      petscint ierr
    END SUBROUTINE vecgetownershiprange

    SUBROUTINE vecgetsize(x,size,ierr)
      vec x
      petscint size
      petscint ierr
    END SUBROUTINE vecgetsize

    SUBROUTINE vecgetvalues(x,n,indices,values,ierr)
      vec x
      petscint n
      petscint indices(*)
      petscscalar values(*)
      petscint ierr
    END SUBROUTINE vecgetvalues

    SUBROUTINE vecghostgetlocalform(g,l,ierr)
      vec g
      vec l
      petscint ierr
    END SUBROUTINE vecghostgetlocalform

    SUBROUTINE vecghostrestorelocalform(g,l,ierr)
      vec g
      vec l
      petscint ierr
    END SUBROUTINE vecghostrestorelocalform

   SUBROUTINE vecghostupdatebegin(x,insertmode,scattermode,ierr)
      vec x
      insertmode insertmode
      scattermode scattermode
      petscint ierr
    END SUBROUTINE vecghostupdatebegin

    SUBROUTINE vecghostupdateend(x,insertmode,scattermode,ierr)
      vec x
      insertmode insertmode
      scattermode scattermode
      petscint ierr
    END SUBROUTINE vecghostupdateend

    SUBROUTINE vecnorm(x,ntype,val,ierr)
      vec x
      normtype ntype
      petscreal val
      petscint ierr
    END SUBROUTINE vecnorm
    
    SUBROUTINE vecrestorearrayf90(x,vec_data,ierr)
      vec x
      petscscalar, POINTER :: vec_data(:)
      petscint ierr
    END SUBROUTINE vecrestorearrayf90
    
    SUBROUTINE vecrestorearrayreadf90(x,vec_data,ierr)
      vec x
      petscscalar, POINTER :: vec_data(:)
      petscint ierr
    END SUBROUTINE vecrestorearrayreadf90
    
    SUBROUTINE vecscale(x,alpha,ierr)
      vec x
      petscscalar alpha
      petscint ierr
    END SUBROUTINE vecscale

    SUBROUTINE vecset(x,value,ierr)
      vec x
      petscscalar value
      petscint ierr
    END SUBROUTINE vecset

    SUBROUTINE vecsetfromoptions(x,ierr)
      vec x
      petscint ierr
    END SUBROUTINE vecsetfromoptions

    SUBROUTINE vecsetlocaltoglobalmapping(v,ctx,ierr)
      vec v
      islocaltoglobalmapping ctx
      petscint ierr
    END SUBROUTINE vecsetlocaltoglobalmapping

    SUBROUTINE vecsetsizes(x,localm,globalm,ierr)
      vec x
      petscint localm,globalm
      petscint ierr
    END SUBROUTINE vecsetsizes

    SUBROUTINE vecsetvalues(x,n,indices,values,insertmode,ierr)
      vec x
      petscint n
      petscint indices(*)
      petscscalar values(*)
      insertmode insertmode
      petscint ierr
    END SUBROUTINE vecsetvalues

    SUBROUTINE vecsetvalueslocal(x,n,indices,values,insertmode,ierr)
      vec x
      petscint n
      petscint indices(*)
      petscscalar values(*)
      insertmode insertmode
      petscint ierr
    END SUBROUTINE vecsetvalueslocal

    SUBROUTINE vecview(x,v,ierr)
      vec x
      petscviewer v
      petscint ierr
    END SUBROUTINE vecview

  END INTERFACE

  INTERFACE petsc_snesgetjacobian
    MODULE PROCEDURE petsc_snesgetjacobiansolver
    MODULE PROCEDURE petsc_snesgetjacobianspecial
  END INTERFACE petsc_snesgetjacobian

  INTERFACE petsc_snessetjacobian
    MODULE PROCEDURE petsc_snessetjacobiansolver
  END INTERFACE petsc_snessetjacobian

  !Miscelaneous routines and constants

  PUBLIC petsc_true,petsc_false
  
  PUBLIC petsc_null_bool,petsc_null_character,petsc_null_function,petsc_null_integer,petsc_null_double,petsc_null_object, &
    & petsc_null_scalar,petsc_null_real

  PUBLIC petsc_default_integer,petsc_default_real

  PUBLIC petsc_decide

  PUBLIC petsc_comm_world,petsc_comm_self

  PUBLIC petsc_viewer_stdout_world,petsc_viewer_stdout_self,petsc_viewer_draw_world,petsc_viewer_draw_self

  PUBLIC petsc_initialise,petsc_finalise

  PUBLIC petsc_errorhandlingseton,petsc_errorhandlingsetoff

  PUBLIC petsc_logview
  
  !Value insert constants

  PUBLIC petsc_add_values,petsc_insert_values

  PUBLIC petsc_scatter_forward,petsc_scatter_reverse

  !Norm constants

  PUBLIC petsc_norm_1,petsc_norm_2,petsc_norm_infinity
  
  !IS routines

  PUBLIC petsc_isinitialise,petsc_isfinalise

  PUBLIC petsc_isdestroy
  
  !IS coloring routines

  PUBLIC petsc_iscoloringinitialise,petsc_iscoloringfinalise

  PUBLIC petsc_iscoloringdestroy

  !IS local to global mapping routines

  PUBLIC petsc_islocaltoglobalmappinginitialise,petsc_islocaltoglobalmappingfinalise

  PUBLIC petsc_islocaltoglobalmappingapply,petsc_islocaltoglobalmappingapplyis,petsc_islocaltoglobalmappingcreate, &
    & petsc_islocaltoglobalmappingdestroy
    
  !KSP routines and constants
 
  PUBLIC petsc_ksprichardson,petsc_kspchebyshev,petsc_kspcg,petsc_kspcgne,petsc_kspnash,petsc_kspstcg,petsc_kspgltr, &
    & petsc_kspgmres,petsc_kspfgmres,petsc_ksplgmres,petsc_kspdgmres,petsc_ksppgmres,petsc_ksptcqmr,petsc_kspbcgs, &
    & petsc_kspibcgs,petsc_kspfbcgs,petsc_kspfbcgsr,petsc_kspbcgsl,petsc_kspcgs,petsc_ksptfqmr,petsc_kspcr,petsc_ksplsqr, &
    & petsc_ksppreonly,petsc_kspqcg,petsc_kspbicg,petsc_kspminres,petsc_kspsymmlq,petsc_ksplcd,petsc_ksppython,petsc_kspgcr

  PUBLIC petsc_ksp_converged_rtol,petsc_ksp_converged_atol,petsc_ksp_converged_its,petsc_ksp_converged_iterating, &
    & petsc_ksp_converged_cg_neg_curve,petsc_ksp_converged_cg_constrained,petsc_ksp_converged_step_length, &
    & petsc_ksp_converged_happy_breakdown,petsc_ksp_diverged_null,petsc_ksp_diverged_its,petsc_ksp_diverged_dtol, &
    & petsc_ksp_diverged_breakdown,petsc_ksp_diverged_breakdown_bicg,petsc_ksp_diverged_nonsymmetric, &
    & petsc_ksp_diverged_indefinite_pc,petsc_ksp_diverged_nanorinf,petsc_ksp_diverged_indefinite_mat

  PUBLIC petsc_ksp_norm_none,petsc_ksp_norm_preconditioned,petsc_ksp_norm_unpreconditioned,petsc_ksp_norm_natural

  PUBLIC petsc_kspinitialise,petsc_kspfinalise

  PUBLIC petsc_kspcreate,petsc_kspdestroy,petsc_kspgetconvergedreason,petsc_kspgetiterationnumber,petsc_kspgetpc, &
    & petsc_kspgetresidualnorm,petsc_kspgmressetrestart,petsc_kspsetfromoptions,petsc_kspsetinitialguessnonzero, &
    & petsc_kspsetoperators,petsc_kspsetreusepreconditioner,petsc_kspsettolerances,petsc_kspsettype,petsc_kspsetup, &
    & petsc_kspsolve

  !Matrix routines and constants

  PUBLIC petsc_mat_flush_assembly,petsc_mat_final_assembly

  PUBLIC petsc_mat_do_not_copy_values,petsc_mat_copy_values,petsc_mat_share_nonzero_pattern
  
  PUBLIC petsc_mat_info_size,petsc_mat_info_block_size,petsc_mat_info_nz_allocated,petsc_mat_info_nz_used, &
    & petsc_mat_info_nz_unneeded,petsc_mat_info_memory,petsc_mat_info_assemblies,petsc_mat_info_mallocs, &
    & petsc_mat_info_fill_ratio_given,petsc_mat_info_fill_ratio_needed,petsc_mat_info_factor_mallocs

  PUBLIC petsc_mat_local,petsc_mat_global_max,petsc_mat_global_sum

  PUBLIC petsc_mat_row_oriented,petsc_mat_new_nonzero_locations,petsc_mat_symmetric,petsc_mat_structurally_symmetric, &
    & petsc_mat_new_diagonals,petsc_mat_ignore_off_proc_entries,petsc_mat_new_nonzero_location_err, &
    & petsc_mat_new_nonzero_allocation_err,petsc_mat_use_hash_table,petsc_mat_keep_nonzero_pattern, &
    & petsc_mat_ignore_zero_entries,petsc_mat_use_inodes,petsc_mat_hermitian,petsc_mat_symmetry_eternal, &
    & petsc_mat_dummy,petsc_mat_ignore_lower_triangular,petsc_mat_error_lower_triangular,petsc_mat_getrow_uppertriangular, &
    & petsc_mat_unused_nonzero_location_err,petsc_mat_spd,petsc_mat_no_off_proc_entries,petsc_mat_no_off_proc_zero_rows

  PUBLIC petsc_mat_solver_superlu,petsc_mat_solver_superlu_dist,petsc_mat_solver_umfpack,petsc_mat_solver_cholmod, &
    & petsc_mat_solver_essl,petsc_mat_solver_lusol,petsc_mat_solver_mumps,petsc_mat_solver_pastix,petsc_mat_solver_matlab, &
    & petsc_mat_solver_petsc,petsc_mat_solver_bas,petsc_mat_solver_cusparse,petsc_mat_solver_bstrm,petsc_mat_solver_sbstrm

  PUBLIC petsc_different_nonzero_pattern,petsc_subset_nonzero_pattern,petsc_same_nonzero_pattern

  PUBLIC petsc_mat_initial_matrix,petsc_mat_reuse_matrix,petsc_mat_ignore_matrix
 
  PUBLIC petsc_matinitialise,petsc_matfinalise

  PUBLIC petsc_matassemblybegin,petsc_matassemblyend,petsc_matcreate,petsc_matcreateaij,petsc_matcreatedense, &
    & petsc_matcreateseqaij,petsc_matcreateseqdense,petsc_matdensegetarrayf90,petsc_matdenserestorearrayf90, &
    & petsc_matdestroy,petsc_matgetinfo,petsc_matgetownershiprange,petsc_matgetrow,petsc_matgetvalues, &
    & petsc_matmumpsseticntl,petsc_matmumpssetcntl,petsc_matrestorerow,petsc_matseqaijgetarrayf90, &
    & petsc_matseqaijgetmaxrownonzeros,petsc_matseqaijrestorearrayf90,petsc_matsetlocaltoglobalmapping, &
    & petsc_matsetoption,petsc_matsetsizes,petsc_matsetvalue,petsc_matsetvalues,petsc_matsetvaluelocal, &
    & petsc_matsetvalueslocal,petsc_matview,petsc_matzeroentries

  !Matrix coloring routines and constants
  PUBLIC petsc_matcoloring_natural,petsc_matcoloring_sl,petsc_matcoloring_lf,petsc_matcoloring_id,petsc_matcoloring_greedy, &
    & petsc_matcoloring_jp

  PUBLIC petsc_matcoloringinitialise,petsc_matcoloringfinalise
  
  PUBLIC petsc_matcoloringapply,petsc_matcoloringcreate,petsc_matcoloringdestroy,petsc_matcoloringsetfromoptions, &
    & petsc_matcoloringsettype

  !Matrix FD coloring routines and constants

  PUBLIC petsc_matfdcoloringinitialise,petsc_matfdcoloringfinalise
  
  PUBLIC petsc_matfdcoloringcreate,petsc_matfdcoloringdestroy,petsc_matfdcoloringsetfromoptions,petsc_matfdcoloringsetfunction, &
    & petsc_matfdcoloringsetparameters,petsc_matfdcoloringsetup

  !Pre-conditioner routines and constants
  
  PUBLIC petsc_pcnone,petsc_pcjacobi,petsc_pcsor,petsc_pclu,petsc_pcshell,petsc_pcbjacobi,petsc_pcmg,petsc_pceisenstat, &
    & petsc_pcilu,petsc_pcicc,petsc_pcasm,petsc_pcksp,petsc_pccomposite,petsc_pcredundant,petsc_pcspai,petsc_pcnn, &
    & petsc_pccholesky,petsc_pcpbjacobi,petsc_pcmat,petsc_pchypre,petsc_pcparms,petsc_pcfieldsplit,petsc_pctfs,petsc_pcml, &
    & petsc_pcgalerkin,petsc_pcexotic,petsc_pcsupportgraph,petsc_pccp,petsc_pcbfbt,petsc_pclsc,petsc_pcpython,petsc_pcpfmg, &
    & petsc_pcsyspfmg,petsc_pcredistribute,petsc_pcsvd,petsc_pcgamg,petsc_pcgasm,petsc_pcsacusp,petsc_pcsacusppoly, &
    & petsc_pcbicgstabcusp,petsc_pcainvcusp,petsc_pcbddc
  
  PUBLIC petsc_pcinitialise,petsc_pcfinalise

  PUBLIC petsc_pcfactorgetmatrix,petsc_pcfactorsetmatsolverpackage,petsc_pcfactorsetupmatsolverpackage,petsc_pcsetfromoptions, &
    & petsc_pcsetreusepreconditioner,petsc_pcsettype

  !SNES routines and constants
  
  PUBLIC petsc_snesnewtonls,petsc_snesnewtontr,petsc_snespython,petsc_snestest,petsc_snesnrichardson,petsc_snesksponly, &
    & petsc_snesvinewtonrsls,petsc_snesvinewtonssls,petsc_snesngmres,petsc_snesqn,petsc_snesshell,petsc_snesncg,petsc_snesfas, &
    & petsc_snesms
  
  PUBLIC petsc_snes_converged_fnorm_abs,petsc_snes_converged_fnorm_relative,petsc_snes_converged_snorm_relative, &
    & petsc_snes_converged_its,petsc_snes_converged_tr_delta,petsc_snes_diverged_function_domain, &
    & petsc_snes_diverged_function_count,petsc_snes_diverged_linear_solve,petsc_snes_diverged_fnorm_nan, &
    & petsc_snes_diverged_max_it,petsc_snes_diverged_line_search,petsc_snes_diverged_local_min,petsc_snes_converged_iterating

  PUBLIC petsc_snes_norm_default,petsc_snes_norm_none,petsc_snes_norm_always,petsc_snes_norm_initial_only, &
    & petsc_snes_norm_final_only,petsc_snes_norm_initial_final_only

  PUBLIC petsc_snes_qn_lbfgs,petsc_snes_qn_broyden,petsc_snes_qn_badbroyden
  
  PUBLIC petsc_snes_qn_restart_none,petsc_snes_qn_restart_powell,petsc_snes_qn_restart_periodic
  
  PUBLIC petsc_snes_qn_scale_default,petsc_snes_qn_scale_none,petsc_snes_qn_scale_shanno,petsc_snes_qn_scale_linesearch, &
    & petsc_snes_qn_scale_jacobian

  PUBLIC petsc_snesinitialise,petsc_snesfinalise

  PUBLIC petsc_snescomputejacobiandefault,petsc_snescomputejacobiandefaultcolor

  PUBLIC petsc_snescreate,petsc_snesdestroy,petsc_snesgetapplicationcontext,petsc_snesgetconvergedreason,petsc_snesgetfunction, &
    & petsc_snesgetiterationnumber,petsc_snesgetjacobian,petsc_snesgetksp,petsc_snesgetlinesearch,petsc_snesgetsolutionupdate, &
    & petsc_snesmonitorset,petsc_snesqnsetrestarttype,petsc_snesqnsetscaletype,petsc_snesqnsettype, &
    & petsc_snessetapplicationcontext,petsc_snessetconvergencetest,petsc_snessetfromoptions,petsc_snessetfunction, &
    & petsc_snessetjacobian,petsc_snessetksp,petsc_snessetnormschedule,petsc_snessettolerances,petsc_snessettrustregiontolerance, &
    & petsc_snessettype,petsc_snessolve

  !SNES line search routines and constants

  PUBLIC petsc_snes_linesearch_basic,petsc_snes_linesearch_bt,petsc_snes_linesearch_l2,petsc_snes_linesearch_cp, &
    & petsc_snes_linesearch_shell
  
  PUBLIC petsc_snes_linesearch_order_linear,petsc_snes_linesearch_order_quadratic,petsc_snes_linesearch_order_cubic

  PUBLIC petsc_sneslinesearchinitialise,petsc_sneslinesearchfinalise

  PUBLIC petsc_sneslinesearchbtsetalpha,petsc_sneslinesearchcomputenorms,petsc_sneslinesearchgetnorms,petsc_sneslinesearchgetvecs, &
    & petsc_sneslinesearchsetcomputenorms,petsc_sneslinesearchsetmonitor,petsc_sneslinesearchsetnorms, &
    & petsc_sneslinesearchsetorder,petsc_sneslinesearchsettolerances,petsc_sneslinesearchsettype
  
  !TS routines and constants
  
  PUBLIC petsc_ts_euler,petsc_ts_beuler,petsc_ts_pseudo,petsc_ts_cn,petsc_ts_sundials,petsc_ts_rk,petsc_ts_python, &
    & petsc_ts_theta,petsc_ts_alpha,petsc_ts_gl,petsc_ts_ssp,petsc_ts_arkimex,petsc_ts_rosw,petsc_ts_eimex

  PUBLIC petsc_ts_converged_iterating,petsc_ts_converged_time,petsc_ts_converged_its,petsc_ts_diverged_nonlinear_solve, &
    & petsc_ts_diverged_step_rejected

  PUBLIC petsc_ts_linear,petsc_ts_nonlinear

  PUBLIC petsc_sundials_adams,petsc_sundials_bdf

  PUBLIC petsc_sundials_modified_gs,petsc_sundials_classical_gs

  PUBLIC petsc_tsinitialise,petsc_tsfinalise

  PUBLIC petsc_tscreate,petsc_tsdestroy,petsc_tsgetsolution,petsc_tsmonitorset,petsc_tssetduration,petsc_tssetexactfinaltime, &
    & petsc_tssetfromoptions,petsc_tssetinitialtimestep,petsc_tssetproblemtype,petsc_tssetrhsfunction,petsc_tssetsolution, &
    & petsc_tssettimestep,petsc_tssettype,petsc_tssolve,petsc_tsstep,petsc_tssundialssettolerance,petsc_tssundialssettype

  !Vector routines and constants

  PUBLIC petsc_vecinitialise,petsc_vecfinalise

  PUBLIC petsc_vecassemblybegin,petsc_vecassemblyend,petsc_veccopy,petsc_veccreate,petsc_veccreateghost, &
    & petsc_veccreateghostwitharray,petsc_veccreatempi,petsc_veccreatempiwitharray,petsc_veccreateseq, &
    & petsc_veccreateseqwitharray,petsc_vecdestroy,petsc_vecduplicate,petsc_vecdot,petsc_vecgetarrayf90, &
    & petsc_vecgetarrayreadf90,petsc_vecgetlocalsize,petsc_vecgetownershiprange,petsc_vecgetsize,petsc_vecgetvalues, &
    & petsc_vecnorm,petsc_vecrestorearrayf90,petsc_vecrestorearrayreadf90,petsc_vecscale,petsc_vecset,petsc_vecsetfromoptions, &
    & petsc_vecsetlocaltoglobalmapping,petsc_vecsetsizes,petsc_vecsetvalues,petsc_vecview
  
CONTAINS

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_errorhandlingsetoff(err,error,*)

    !Argument Variables
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ErrorHandlingSetOff",err,error,*999)

    petschandleerror=.false.
    
    exits("Petsc_ErrorHandlingSetOff")
    RETURN
999 errorsexits("Petsc_ErrorHandlingSetOff",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_errorhandlingsetoff
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_errorhandlingseton(err,error,*)

    !Argument Variables
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ErrorHandlingSetOn",err,error,*999)

    petschandleerror=.true.
    
    exits("Petsc_ErrorHandlingSetOn")
    RETURN
999 errorsexits("Petsc_ErrorHandlingSetOn",err,error)
    RETURN 1
  END SUBROUTINE petsc_errorhandlingseton
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_finalise(err,error,*)

    !Argument Variables
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_Finalise",err,error,*999)

    CALL petscfinalize(err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PetscFinalize.",err,error,*999)
    ENDIF
    
    exits("Petsc_Finalise")
    RETURN
999 errorsexits("Petsc_Finalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_finalise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_initialise(file,err,error,*)

    !Argument Variables
    CHARACTER(LEN=*), INTENT(IN) :: file
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_Initialise",err,error,*999)

    petschandleerror=.true.
    CALL petscinitialize(file,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PetscInitialize.",err,error,*999)
    ENDIF
    ! Disable PETSc's signal handler as we have our own OpenCMISS signal handlers in cmiss_c.c
    CALL petscpopsignalhandler(err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PetscPopSignalHandler.",err,error,*999)
    ENDIF
    
    exits("Petsc_Initialise")
    RETURN
999 errorsexits("Petsc_Initialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_initialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_logview(viewer,err,error,*)

    !Argument Variables
    petscviewer, INTENT(IN) :: viewer 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_LogView",err,error,*999)

    CALL petsclogview(viewer,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PetscLogView.",err,error,*999)
    ENDIF
    
    exits("Petsc_LogView")
    RETURN
999 errorsexits("Petsc_LogView",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_logview
  
  !
  !================================================================================================================================
  !

  !Finalise the PETSc IS structure and destroy the IS
  SUBROUTINE petsc_isfinalise(is,err,error,*)

    !Argument Variables
    TYPE(petscistype), INTENT(INOUT) :: is
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISFinalise",err,error,*999)

    IF(is%is/=petsc_null_object) THEN
      CALL petsc_isdestroy(is,err,error,*999)
    ENDIF
    
    exits("Petsc_ISFinalise")
    RETURN
999 errorsexits("Petsc_ISFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_isfinalise
    
  !
  !================================================================================================================================
  !
  
  !Initialise the PETSc IS structure
  SUBROUTINE petsc_isinitialise(is,err,error,*)

    !Argument Variables
    TYPE(petscistype), INTENT(INOUT) :: is
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISInitialise",err,error,*999)

    is%is=petsc_null_object
    
    exits("Petsc_ISInitialise")
    RETURN
999 errorsexits("Petsc_ISInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_isinitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_isdestroy(is,err,error,*)

    !Argument Variables
    TYPE(petscistype), INTENT(INOUT) :: is
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISDestroy",err,error,*999)

    CALL isdestroy(is%is,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISDestroy.",err,error,*999)
    ENDIF
    is%is=petsc_null_object
        
    exits("Petsc_ISDestroy")
    RETURN
999 errorsexits("Petsc_ISDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_isdestroy
    
  !
  !
  !================================================================================================================================
  !

  !Finalise the PETSc ISColoring structure and destroy the ISColoring
  SUBROUTINE petsc_iscoloringfinalise(iscoloring,err,error,*)

    !Argument Variables
    TYPE(petsciscoloringtype), INTENT(INOUT) :: iscoloring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISColoringFinalise",err,error,*999)

    IF(iscoloring%iscoloring/=petsc_null_object) THEN
      CALL petsc_iscoloringdestroy(iscoloring,err,error,*999)
    ENDIF
    
    exits("Petsc_ISColoringFinalise")
    RETURN
999 errorsexits("Petsc_ISColoringFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_iscoloringfinalise
    
  !
  !================================================================================================================================
  !
  
  !Initialise the PETSc ISColoring structure
  SUBROUTINE petsc_iscoloringinitialise(iscoloring,err,error,*)

    !Argument Variables
    TYPE(petsciscoloringtype), INTENT(INOUT) :: iscoloring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISColoringInitialise",err,error,*999)

    iscoloring%iscoloring=petsc_null_object
    
    exits("Petsc_ISColoringInitialise")
    RETURN
999 errorsexits("Petsc_ISColoringInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_iscoloringinitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_iscoloringdestroy(iscoloring,err,error,*)

    !Argument Variables
    TYPE(petsciscoloringtype), INTENT(INOUT) :: iscoloring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISColoringDestroy",err,error,*999)

    CALL iscoloringdestroy(iscoloring%iscoloring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISColoringDestroy.",err,error,*999)
    ENDIF
    iscoloring%iscoloring=petsc_null_object
    
    exits("Petsc_ISColoringDestroy")
    RETURN
999 errorsexits("Petsc_ISColoringDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_iscoloringdestroy
  
  !  
  !================================================================================================================================
  !

  !Finalise the PETSc ISLocalToGlobalMapping structure and destroy the ISLocalToGlobalMapping
  SUBROUTINE petsc_islocaltoglobalmappingfinalise(isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    TYPE(petscislocaltogloabalmappingtype), INTENT(INOUT) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingFinalise",err,error,*999)

    IF(islocaltoglobalmapping%isLocalToGlobalMapping/=petsc_null_object) THEN
      CALL petsc_islocaltoglobalmappingdestroy(islocaltoglobalmapping,err,error,*999)
    ENDIF
    
    exits("Petsc_ISLocalToGlobalMappingFinalise")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappingfinalise
    
  !
  !================================================================================================================================
  !
  
  !Initialise the PETSc ISLocalToGlobalMapping structure
  SUBROUTINE petsc_islocaltoglobalmappinginitialise(isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    TYPE(petscislocaltogloabalmappingtype), INTENT(INOUT) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingInitialise",err,error,*999)

    islocaltoglobalmapping%isLocalToGlobalMapping=petsc_null_object
    
    exits("Petsc_ISLocalToGlobalMappingInitialise")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappinginitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_islocaltoglobalmappingapply(isLocalToGlobalMapping,n,idxIn,idxOut,err,error,*)

    !Argument Variables
    TYPE(petscislocaltogloabalmappingtype), INTENT(INOUT) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: idxIn(:)
    INTEGER(INTG), INTENT(OUT) :: idxOut(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingApply",err,error,*999)

    CALL islocaltoglobalmappingapply(islocaltoglobalmapping%isLocalToGlobalMapping,n,idxin,idxout,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISLocalToGlobalMappingApply.",err,error,*999)
    ENDIF
    islocaltoglobalmapping%isLocalToGlobalMapping=petsc_null_object
    
    exits("Petsc_ISLocalToGlobalMappingApply")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingApply",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappingapply
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_islocaltoglobalmappingapplyis(isLocalToGlobalMapping,isIn,isOut,err,error,*)

    !Argument Variables
    TYPE(petscislocaltogloabalmappingtype), INTENT(IN) :: isLocalToGlobalMapping
    TYPE(petscistype), INTENT(IN) :: isIn
    TYPE(petscistype), INTENT(OUT) :: isOut
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingApplyIS",err,error,*999)

    CALL islocaltoglobalmappingapplyis(islocaltoglobalmapping%isLocalToGlobalMapping,isin%is,isout%is,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISLocalToGlobalMappingApplyIS.",err,error,*999)
    ENDIF
    
    exits("Petsc_ISLocalToGlobalMappingApplyIS")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingApplyIS",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappingapplyis
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_islocaltoglobalmappingcreate(communicator,blockSize,n,indices,mode,isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: blockSize
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: indices(:)
    petsccopymode, INTENT(IN) :: mode 
    TYPE(petscislocaltogloabalmappingtype), INTENT(INOUT) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingCreate",err,error,*999)

    CALL islocaltoglobalmappingcreate(communicator,blocksize,n,indices,mode,islocaltoglobalmapping%isLocalToGlobalMapping,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISLocalToGlobalMappingCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_ISLocalToGlobalMappingCreate")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappingcreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_islocaltoglobalmappingdestroy(isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    TYPE(petscislocaltogloabalmappingtype), INTENT(INOUT) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_ISLocalToGlobalMappingDestroy",err,error,*999)

    CALL islocaltoglobalmappingdestroy(islocaltoglobalmapping%isLocalToGlobalMapping,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in ISLocalToGlobalMappingDestroy.",err,error,*999)
    ENDIF
    islocaltoglobalmapping%isLocalToGlobalMapping=petsc_null_object
    
    exits("Petsc_ISLocalToGlobalMappingDestroy")
    RETURN
999 errorsexits("Petsc_ISLocalToGlobalMappingDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_islocaltoglobalmappingdestroy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspfinalise(ksp,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPFinalise",err,error,*999)

    IF(ksp%ksp/=petsc_null_object) THEN
      CALL petsc_kspdestroy(ksp,err,error,*999)
    ENDIF
    
    exits("Petsc_KSPFinalise")
    RETURN
999 errorsexits("Petsc_KSPFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspfinalise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspinitialise(ksp,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPInitialise",err,error,*999)

    ksp%ksp=petsc_null_object
    
    exits("Petsc_KSPInitialise")
    RETURN
999 errorsexits("Petsc_KSPInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspinitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspcreate(communicator,ksp,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPCreate",err,error,*999)

    CALL kspcreate(communicator,ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPCreate")
    RETURN
999 errorsexits("Petsc_KSPCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspcreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspdestroy(ksp,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPDestroy",err,error,*999)

    CALL kspdestroy(ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPDestroy.",err,error,*999)
    ENDIF
    ksp%ksp=petsc_null_object
    
    exits("Petsc_KSPDestroy")
    RETURN
999 errorsexits("Petsc_KSPDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspdestroy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspgetconvergedreason(ksp,reason,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: reason
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPGetConvergedReason",err,error,*999)

    CALL kspgetconvergedreason(ksp%ksp,reason,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPGetConvergedReason.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPGetConvergedReason")
    RETURN
999 errorsexits("Petsc_KSPGetConvergedReason",err,error)
    RETURN 1
  END SUBROUTINE petsc_kspgetconvergedreason
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspgetiterationnumber(ksp,iterationNumber,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: iterationNumber
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPGetIterationNumber",err,error,*999)

    CALL kspgetiterationnumber(ksp%ksp,iterationnumber,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPGetIterationNumber.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPGetIterationNumber")
    RETURN
999 errorsexits("Petsc_KSPGetIterationNumber",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspgetiterationnumber
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspgetpc(ksp,pc,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    TYPE(petscpctype), INTENT(INOUT) :: pc
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPGetPC",err,error,*999)

    CALL kspgetpc(ksp%ksp,pc%pc,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPGetPC.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPGetPC")
    RETURN
999 errorsexits("Petsc_KSPGetPC",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspgetpc
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspgetresidualnorm(ksp,residualNorm,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    REAL(DP), INTENT(OUT) :: residualNorm
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPGetResidualNorm",err,error,*999)

    CALL kspgetresidualnorm(ksp%ksp,residualnorm,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPGetResidualNorm.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPGetResidualNorm")
    RETURN
999 errorsexits("Petsc_KSPGetResidualNorm",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspgetresidualnorm
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspgmressetrestart(ksp,restart,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: restart
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPGMRESSetRestart",err,error,*999)

    CALL kspgmressetrestart(ksp%ksp,restart,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPGMRESSetRestart.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPGMRESSetRestart")
    RETURN
999 errorsexits("Petsc_KSPGMRESSetRestart",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspgmressetrestart
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsetfromoptions(ksp,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetFromOptions",err,error,*999)

    CALL kspsetfromoptions(ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetFromOptions")
    RETURN
999 errorsexits("Petsc_KSPSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsetfromoptions
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsetinitialguessnonzero(ksp,flag,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    LOGICAL, INTENT(IN) :: flag
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetInitialGuessNonZero",err,error,*999)

    IF(flag) THEN
      CALL kspsetinitialguessnonzero(ksp%ksp,petsc_true,err)
    ELSE
      CALL kspsetinitialguessnonzero(ksp%ksp,petsc_false,err)
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetInitialGuessNonzero.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetInitialGuessNonZero")
    RETURN
999 errorsexits("Petsc_KSPSetInitialGuessNonZero",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsetinitialguessnonzero
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsetoperators(ksp,amat,pmat,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    TYPE(petscmattype), INTENT(INOUT) :: amat
    TYPE(petscmattype), INTENT(INOUT) :: pmat
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetOperators",err,error,*999)

    CALL kspsetoperators(ksp%ksp,amat%mat,pmat%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetFromOperators.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetOperators")
    RETURN
999 errorsexits("Petsc_KSPSetOperators",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsetoperators
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsetreusepreconditioner(ksp,flag,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    LOGICAL, INTENT(IN) :: flag
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetReusePreconditioner",err,error,*999)

    IF(flag) THEN
      CALL kspsetreusepreconditioner(ksp%ksp,petsc_true,err)
    ELSE
      CALL kspsetreusepreconditioner(ksp%ksp,petsc_false,err)
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetReusePreconditioner.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetReusePreconditioner")
    RETURN
999 errorsexits("Petsc_KSPSetReusePreconditioner",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsetreusepreconditioner
  
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsettolerances(ksp,rTol,aTol,dTol,maxIterations,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    REAL(DP), INTENT(IN) :: rTol
    REAL(DP), INTENT(IN) :: aTol
    REAL(DP), INTENT(IN) :: dTol
    INTEGER(INTG), INTENT(IN) :: maxIterations
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetTolerances",err,error,*999)

    CALL kspsettolerances(ksp%ksp,rtol,atol,dtol,maxiterations,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetTolerances.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetTolerances")
    RETURN
999 errorsexits("Petsc_KSPSetTolerances",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsettolerances
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsettype(ksp,method,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    ksptype, INTENT(IN) :: method 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetType",err,error,*999)

    CALL kspsettype(ksp%ksp,method,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetType")
    RETURN
999 errorsexits("Petsc_KSPSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsettype
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsetup(ksp,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSetUp",err,error,*999)

    CALL kspsetup(ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSetUp.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSetUp")
    RETURN
999 errorsexits("Petsc_KSPSetUp",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsetup
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_kspsolve(ksp,b,x,err,error,*)

    !Argument Variables
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    TYPE(petscvectype), INTENT(INOUT) :: b
    TYPE(petscvectype), INTENT(INOUT)  :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_KSPSolve",err,error,*999)

    CALL kspsolve(ksp%ksp,b%vec,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in KSPSolve.",err,error,*999)
    ENDIF
    
    exits("Petsc_KSPSolve")
    RETURN
999 errorsexits("Petsc_KSPSolve",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_kspsolve
    
  !
  !================================================================================================================================
  !

  !Finalise the PETSc Mat structure
  SUBROUTINE petsc_matfinalise(a,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFinalise",err,error,*999)

    IF(a%mat/=petsc_null_object) THEN
      CALL petsc_matdestroy(a,err,error,*999)
    ENDIF
    
    exits("Petsc_MatFinalise")
    RETURN
999 errorsexits("Petsc_MatFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfinalise
    
  !
  !================================================================================================================================
  !

  !Initialise the PETSc Mat structure
  SUBROUTINE petsc_matinitialise(a,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatInitialise",err,error,*999)

    a%mat=petsc_null_object
     
    exits("Petsc_MatInitialise")
    RETURN
999 errorsexits("Petsc_MatInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matinitialise
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matassemblybegin(A,assemblyType,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: A !The matrix to assemble
    matassemblytype, INTENT(IN) :: assemblytype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatAssemblyBegin",err,error,*999)

    CALL matassemblybegin(a%mat,assemblytype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatAssemblyBegin.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatAssemblyBegin")
    RETURN
999 errorsexits("Petsc_MatAssemblyBegin",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matassemblybegin
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matassemblyend(A,assemblyType,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: A
    matassemblytype, INTENT(IN) :: assemblytype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatAssemblyEnd",err,error,*999)

    CALL matassemblyend(a%mat,assemblytype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatAssemblyEnd.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatAssemblyEnd")
    RETURN
999 errorsexits("Petsc_MatAssemblyEnd",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matassemblyend
       
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matcreate(communicator,A,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    TYPE(petscmattype), INTENT(INOUT) :: A
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatCreate",err,error,*999)

    CALL matcreate(communicator,a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatCreate")
    RETURN
999 errorsexits("Petsc_MatCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matcreateaij(communicator,localM,localN,globalM,globalN,diagNumberNonZerosPerRow,diagNumberNonZerosEachRow, &
    & offdiagnumbernonzerosperrow,offdiagnumbernonzeroseachrow,a,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localM
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalM
    INTEGER(INTG), INTENT(IN) :: globalN
    INTEGER(INTG), INTENT(IN) :: diagNumberNonZerosPerRow
    INTEGER(INTG), INTENT(IN) :: diagNumberNonZerosEachRow(:)
    INTEGER(INTG), INTENT(IN) :: offDiagNumberNonZerosPerRow
    INTEGER(INTG), INTENT(IN) :: offDiagNumberNonZerosEachRow(:)
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatCreateAIJ",err,error,*999)

    CALL matcreateaij(communicator,localm,localn,globalm,globaln,diagnumbernonzerosperrow,diagnumbernonzeroseachrow, &
      & offdiagnumbernonzerosperrow,offdiagnumbernonzeroseachrow,a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatCreateAIJ.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatCreateAIJ")
    RETURN
999 errorsexits("Petsc_MatCreateAIJ",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcreateaij
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matcreatedense(communicator,localM,localN,globalM,globalN,matrixData,a,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localM
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalM
    INTEGER(INTG), INTENT(IN) :: globalN
    REAL(DP), INTENT(IN) :: matrixData(:)
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatCreateDense",err,error,*999)

    CALL matcreatedense(communicator,localm,localn,globalm,globaln,matrixdata,a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatCreateDense.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatCreateDense")
    RETURN
999 errorsexits("Petsc_MatCreateDense",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcreatedense
  
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matcreateseqaij(communicator,m,n,numberNonZerosPerRow,numberNonZerosEachRow,a,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: m
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: numberNonZerosPerRow
    INTEGER(INTG), INTENT(IN) :: numberNonZerosEachRow(:)
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatCreateSeqAIJ",err,error,*999)

    CALL matcreateseqaij(communicator,m,n,numbernonzerosperrow,numbernonzeroseachrow,a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatCreateSeqAIJ.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatCreateSeqAIJ")
    RETURN
999 errorsexits("Petsc_MatCreateSeqAIJ",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcreateseqaij
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matcreateseqdense(communicator,m,n,matrixData,a,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: m
    INTEGER(INTG), INTENT(IN) :: n
    REAL(DP), INTENT(IN) :: matrixData(*)
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatCreateSeqDense",err,error,*999)

    CALL matcreateseqdense(communicator,m,n,matrixdata,a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatCreateSeqDense.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatCreateSeqDense")
    RETURN
999 errorsexits("Petsc_MatCreateSeqDense",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcreateseqdense

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matdensegetarrayf90(a,array,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT), TARGET :: a
    REAL(DP), POINTER :: array(:,:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatDenseGetArrayF90",err,error,*999)

    IF(ASSOCIATED(array)) THEN
      CALL flagerror("Array is already associated.",err,error,*999)
    ELSE
      CALL matdensegetarrayf90(a%mat,array,err)
      IF(err/=0) THEN
        IF(petschandleerror) THEN
          chkerrq(err)
        ENDIF
        CALL flagerror("PETSc error in MatDenseGetArrayF90.",err,error,*999)
      ENDIF
    ENDIF
    
    exits("Petsc_MatDenseGetArrayF90")
    RETURN
999 errorsexits("Petsc_MatDenseGetArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matdensegetarrayf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matdenserestorearrayf90(a,array,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    REAL(DP), POINTER :: array(:,:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatDenseRestoreArrayF90",err,error,*999)

    CALL matdenserestorearrayf90(a%mat,array,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatDenseRestoreArrayF90.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatDenseRestoreArrayF90")
    RETURN
999 errorsexits("Petsc_MatDenseRestoreArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matdenserestorearrayf90
    
  !
  !================================================================================================================================
  !    

  SUBROUTINE petsc_matdestroy(a,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatDestroy",err,error,*999)

    CALL matdestroy(a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatDestroy.",err,error,*999)
    ENDIF
    a%mat=petsc_null_object
     
    exits("Petsc_MatDestroy")
    RETURN
999 errorsexits("Petsc_MatDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matdestroy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matgetinfo(a,flag,info,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    matinfotype, INTENT(IN) :: flag 
    matinfo, INTENT(OUT) :: info(mat_info_size) 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatGetInfo",err,error,*999)

    CALL matgetinfo(a%mat,flag,info,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatGetInfo.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatGetInfo")
    RETURN
999 errorsexits("Petsc_MatGetInfo",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matgetinfo
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matgetownershiprange(a,firstRow,lastRow,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: firstRow
    INTEGER(INTG), INTENT(OUT) :: lastRow
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatGetOwnershipRange",err,error,*999)

    CALL matgetownershiprange(a%mat,firstrow,lastrow,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatGetOwnershipRange.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatGetOwnershipRange")
    RETURN
999 errorsexits("Petsc_MatGetOwnershipRange",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matgetownershiprange
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matgetrow(A,rowNumber,numberOfColumns,columns,values,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: A
    INTEGER(INTG), INTENT(IN) :: rowNumber
    INTEGER(INTG), INTENT(OUT) :: numberOfColumns
    INTEGER(INTG), INTENT(OUT) :: columns(:)
    REAL(DP), INTENT(OUT) :: values(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatGetRow",err,error,*999)

    CALL matgetrow(a%mat,rownumber,numberofcolumns,columns,values,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatGetRow.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatGetRow")
    RETURN
999 errorsexits("Petsc_MatGetRow",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matgetrow
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matgetvalues(a,m,mIndices,n,nIndices,values,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: m
    INTEGER(INTG), INTENT(IN) :: mIndices(*)
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: nIndices(*)
    REAL(DP), INTENT(OUT) :: values(*)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatGetValues",err,error,*999)

    CALL matgetvalues(a%mat,m,mindices,n,nindices,values,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatGetValues.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatGetValues")
    RETURN
999 errorsexits("Petsc_MatGetValues",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matgetvalues
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matmumpsseticntl(factoredMatrix,icntl,ival,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: factoredMatrix
    INTEGER(INTG), INTENT(IN) :: icntl
    INTEGER(INTG), INTENT(IN) :: ival
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatMumpsSetIcntl",err,error,*999)

    CALL matmumpsseticntl(factoredmatrix%mat,icntl,ival,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatMumpsSetIcntl.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatMumpsSetIcntl")
    RETURN
999 errorsexits("Petsc_MatMumpsSetIcntl",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matmumpsseticntl

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matmumpssetcntl(factoredMatrix,icntl,val,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: factoredMatrix
    INTEGER(INTG), INTENT(IN) :: icntl
    REAL(DP), INTENT(IN) :: val
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatMumpsSetCntl",err,error,*999)

    CALL matmumpssetcntl(factoredmatrix%mat,icntl,val,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatMumpsSetCntl.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatMumpsSetCntl")
    RETURN
999 errorsexits("Petsc_MatMumpsSetCntl",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matmumpssetcntl

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matrestorerow(A,rowNumber,numberOfColumns,columns,values,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: A
    INTEGER(INTG), INTENT(IN) :: rowNumber
    INTEGER(INTG), INTENT(OUT) :: numberOfColumns
    INTEGER(INTG), INTENT(OUT) :: columns(:)
    REAL(DP), INTENT(OUT) :: values(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatRestoreRow",err,error,*999)

    CALL matrestorerow(a%mat,rownumber,numberofcolumns,columns,values,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatRestoreRow.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatRestoreRow")
    RETURN
999 errorsexits("Petsc_MatRestoreRow",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matrestorerow
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matseqaijgetarrayf90(a,array,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT), TARGET :: a
    REAL(DP), POINTER :: array(:,:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSeqAIJGetArrayF90",err,error,*999)

    IF(ASSOCIATED(array)) THEN
      CALL flagerror("Array is already associated.",err,error,*999)
    ELSE
      CALL matseqaijgetarrayf90(a%mat,array,err)
      IF(err/=0) THEN
        IF(petschandleerror) THEN
          chkerrq(err)
        ENDIF
        CALL flagerror("PETSc error in MatSeqAIJGetArrayF90.",err,error,*999)
      ENDIF
    ENDIF
    
    exits("Petsc_MatSeqAIJGetArrayF90")
    RETURN
999 errorsexits("Petsc_MatSeqAIJGetArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matseqaijgetarrayf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matseqaijgetmaxrownonzeros(a,maxNumberNonZeros,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: maxNumberNonZeros
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSeqAIJGetMaxRowNonzeros",err,error,*999)

    !CALL MatSeqAIJGetMaxRowNonzeros(A%mat,maxNumberNonZeros,err)
    maxnumbernonzeros=0
    CALL flagerror("Not implemented.",err,error,*999)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSeqAIJGetMaxRowNonzeros.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSeqAIJGetMaxRowNonzeros")
    RETURN
999 errorsexits("Petsc_MatSeqAIJGetMaxRowNonzeros",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matseqaijgetmaxrownonzeros
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matseqaijrestorearrayf90(a,array,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    REAL(DP), POINTER :: array(:,:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("MatSeqAIJRestoreArrayF90",err,error,*999)

    CALL matseqaijrestorearrayf90(a%mat,array,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSeqAIJRestoreArrayF90.",err,error,*999)
    ENDIF
    
    exits("MatSeqAIJRestoreArrayF90")
    RETURN
999 errorsexits("MatSeqAIJRestoreArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matseqaijrestorearrayf90
  
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetlocaltoglobalmapping(a,isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petscislocaltogloabalmappingtype), INTENT(IN) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetLocalToGlobalMapping",err,error,*999)

    CALL matsetlocaltoglobalmapping(a%mat,islocaltoglobalmapping%isLocalToGlobalMapping,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetLocalToGlobalMapping.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetLocalToGlobalMapping")
    RETURN
999 errorsexits("Petsc_MatSetLocalToGlobalMapping",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetlocaltoglobalmapping
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetoption(a,option,flag,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    matoption, INTENT(IN) :: option 
    LOGICAL, INTENT(IN) :: flag
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetOption",err,error,*999)

    IF(flag) THEN
      CALL matsetoption(a%mat,option,petsc_true,err)
    ELSE
      CALL matsetoption(a%mat,option,petsc_false,err)
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetOption.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetOption")
    RETURN
999 errorsexits("Petsc_MatSetOption",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetoption
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetsizes(a,localM,localN,globalM,globalN,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: localM
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalM
    INTEGER(INTG), INTENT(IN) :: globalN
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetSizes",err,error,*999)

    CALL matsetsizes(a%mat,localm,localn,globalm,globaln,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetSizes.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetSizes")
    RETURN
999 errorsexits("Petsc_MatSetSizes",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetsizes
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matsettype(a,matrixType,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    mattype, INTENT(IN) :: matrixtype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetType",err,error,*999)

    CALL matsettype(a%mat,matrixtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetType")
    RETURN
999 errorsexits("Petsc_MatSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsettype

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetvalue(a,row,col,value,insertMode,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: row
    INTEGER(INTG), INTENT(IN) :: col
    REAL(DP), INTENT(IN) :: value
    insertmode, INTENT(IN) :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetValue",err,error,*999)

    CALL matsetvalue(a%mat,row,col,value,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetValue.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetValue")
    RETURN
999 errorsexits("Petsc_MatSetValue",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetvalue
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetvalues(a,m,mIndices,n,nIndices,values,insertMode,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: m
    INTEGER(INTG), INTENT(IN) :: mIndices(*)
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: nIndices(*)
    REAL(DP), INTENT(IN) :: values(*)
    insertmode, INTENT(IN) :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetValues",err,error,*999)

    CALL matsetvalues(a%mat,m,mindices,n,nindices,values,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetValues.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetValues")
    RETURN
999 errorsexits("Petsc_MatSetValues",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetvalues
    
  !
  !================================================================================================================================
  !
  SUBROUTINE petsc_matsetvaluelocal(a,row,col,VALUE,insertMode,err,error,*)
    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: row
    INTEGER(INTG), INTENT(IN) :: col
    REAL(DP), INTENT(IN) :: value
    insertmode, INTENT(IN) :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetValueLocal",err,error,*999)

    CALL matsetvaluelocal(a%mat,row,col,value,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetValueLocal.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetValueLocal")
    RETURN
999 errorsexits("Petsc_MatSetValueLocal",err,error)
    RETURN 1
  END SUBROUTINE petsc_matsetvaluelocal

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matsetvalueslocal(a,m,mIndices,n,nIndices,values,insertMode,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(IN) :: m
    INTEGER(INTG), INTENT(IN) :: mIndices(:)
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: nIndices(:)
    REAL(DP), INTENT(IN) :: values(:)
    insertmode, INTENT(IN) :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatSetValuesLocal",err,error,*999)

    CALL matsetvalueslocal(a%mat,m,mindices,n,nindices,values,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatSetValuesLocal.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatSetValuesLocal")
    RETURN
999 errorsexits("Petsc_MatSetValuesLocal",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matsetvalueslocal
  
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matview(a,viewer,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    petscviewer, INTENT(IN) :: viewer 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatView",err,error,*999)

    CALL matview(a%mat,viewer,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatView.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatView")
    RETURN
999 errorsexits("Petsc_MatView",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matview
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matzeroentries(a,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatZeroEntries",err,error,*999)

    CALL matzeroentries(a%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatZeroEntries.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatZeroEntries")
    RETURN
999 errorsexits("Petsc_MatZeroEntries",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matzeroentries
    
  !
  !================================================================================================================================
  !
  
  !Finalise the PETSc MatColoring structure and destroy the MatColoring
  SUBROUTINE petsc_matcoloringfinalise(matColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringFinalise",err,error,*999)

    IF(matcoloring%matColoring/=petsc_null_object) THEN
      CALL petsc_matcoloringdestroy(matcoloring,err,error,*999)
    ENDIF
    
    exits("Petsc_MatColoringFinalise")
    RETURN
999 errorsexits("Petsc_MatColoringFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringfinalise
    
  !
  !================================================================================================================================
  !
  
  !Initialise the PETSc MatColoring structure
  SUBROUTINE petsc_matcoloringinitialise(matColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringInitialise",err,error,*999)

    matcoloring%matColoring=petsc_null_object
    
    exits("Petsc_MatColoringInitialise")
    RETURN
999 errorsexits("Petsc_MatColoringInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringinitialise
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matcoloringapply(matColoring,isColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    TYPE(petsciscoloringtype), INTENT(INOUT) :: isColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringApply",err,error,*999)

    CALL matcoloringapply(matcoloring%matColoring,iscoloring%isColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatColoringApply.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatColoringApply")
    RETURN
999 errorsexits("Petsc_MatColoringApply",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringapply
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matcoloringcreate(a,matColoring,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringCreate",err,error,*999)

    CALL matcoloringcreate(a%mat,matcoloring%matColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatColoringCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatColoringCreate")
    RETURN
999 errorsexits("Petsc_MatColoringCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringcreate
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matcoloringdestroy(matColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringDestroy",err,error,*999)

    CALL matcoloringdestroy(matcoloring%matColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatColoringDestroy.",err,error,*999)
    ENDIF
    matcoloring%matColoring=petsc_null_object
    
    exits("Petsc_MatColoringDestroy")
    RETURN
999 errorsexits("Petsc_MatColoringDestroy",err,error)
    exits("Petsc_MatColoringDestroy")
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringdestroy
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matcoloringsetfromoptions(matColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringSetFromOptions",err,error,*999)

    CALL matcoloringsetfromoptions(matcoloring%matColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatColoringSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatColoringSetFromOptions")
    RETURN
999 errorsexits("Petsc_MatColoringSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringsetfromoptions
    
  !
  !================================================================================================================================
  !
  
  SUBROUTINE petsc_matcoloringsettype(matColoring,coloringType,err,error,*)

    !Argument Variables
    TYPE(petscmatcoloringtype), INTENT(INOUT) :: matColoring
    matcoloringtype, INTENT(IN) :: coloringtype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatColoringSetType",err,error,*999)

    CALL matcoloringsettype(matcoloring%matColoring,coloringtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatColoringSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatColoringSetType")
    RETURN
999 errorsexits("Petsc_MatColoringSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matcoloringsettype
    
  !
  !================================================================================================================================
  !

  !Finalise the PETSc MatFDColoring structure and destroy the MatFDColoring
  SUBROUTINE petsc_matfdcoloringfinalise(matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringFinalise",err,error,*999)

    IF(matfdcoloring%matFDColoring/=petsc_null_object) THEN
      CALL petsc_matfdcoloringdestroy(matfdcoloring,err,error,*999)
    ENDIF
    
    exits("Petsc_MatFDColoringFinalise")
    RETURN
999 errorsexits("Petsc_MatFDColoringFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringfinalise
    
  !
  !================================================================================================================================
  !
  
  !Initialise the PETSc MatFDColoring structure
  SUBROUTINE petsc_matfdcoloringinitialise(matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringInitialise",err,error,*999)

    matfdcoloring%matFDColoring=petsc_null_object
    
    exits("Petsc_MatFDColoringInitialise")
    RETURN
999 errorsexits("Petsc_MatFDColoringInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringinitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matfdcoloringcreate(a,isColoring,matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petsciscoloringtype), INTENT(IN) :: isColoring
    TYPE(petscmatfdcoloringtype), INTENT(OUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringCreate",err,error,*999)

    CALL matfdcoloringcreate(a%mat,iscoloring%isColoring,matfdcoloring%matFDColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatFDColoringCreate")
    RETURN
999 errorsexits("Petsc_MatFDColoringCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringcreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matfdcoloringdestroy(matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringDestroy",err,error,*999)

    CALL matfdcoloringdestroy(matfdcoloring%matFDColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringDestroy.",err,error,*999)
    ENDIF
    matfdcoloring%matFDColoring=petsc_null_object
     
    exits("Petsc_MatFDColoringDestroy")
    RETURN
999 errorsexits("Petsc_MatFDColoringDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringdestroy
    
  !
  !================================================================================================================================2
  !

  SUBROUTINE petsc_matfdcoloringsetfromoptions(matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringSetFromOptions",err,error,*999)

    CALL matfdcoloringsetfromoptions(matfdcoloring%matFDColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatFDColoringSetFromOptions")
    RETURN
999 errorsexits("Petsc_MatFDColoringSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringsetfromoptions
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matfdcoloringsetparameters(matFDColoring,rError,uMin,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    REAL(DP) :: rError
    REAL(DP) :: uMin
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringSetParameters",err,error,*999)

    CALL matfdcoloringsetparameters(matfdcoloring%matFDColoring,rerror,umin,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringSetParameters.",err,error,*999)
    ENDIF

    exits("Petsc_MatFDColoringSetParameters")
    RETURN
999 errorsexits("Petsc_MatFDColoringSetParameters",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringsetparameters

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matfdcoloringsetfunction(matFDColoring,fFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    EXTERNAL ffunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringSetFunction",err,error,*999)

    CALL matfdcoloringsetfunction(matfdcoloring%matFDColoring,ffunction,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringSetFunction.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatFDColoringSetFunction")
    RETURN
999 errorsexits("Petsc_MatFDColoringSetFunction",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringsetfunction
        
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_matfdcoloringsetup(a,isColoring,matFDColoring,err,error,*)

    !Argument Variables
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petsciscoloringtype), INTENT(INOUT) :: isColoring
    TYPE(petscmatfdcoloringtype), INTENT(INOUT) :: matFDColoring
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_MatFDColoringSetup",err,error,*999)

    CALL matfdcoloringsetup(a%mat,iscoloring%isColoring,matfdcoloring%matFDColoring,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in MatFDColoringSetup.",err,error,*999)
    ENDIF
    
    exits("Petsc_MatFDColoringSetup")
    RETURN
999 errorsexits("Petsc_MatFDColoringSetup",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_matfdcoloringsetup
    
  !
  !================================================================================================================================
  !

  !Finalise the PETSc PC structure 
  SUBROUTINE petsc_pcfinalise(pc,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCFinalise",err,error,*999)

    IF(pc%pc/=petsc_null_object) THEN
      !Do nothing - should be destroyed when the KSP is destroyed.
    ENDIF
    
    exits("Petsc_PCFinalise")
    RETURN
999 errorsexits("Petsc_PCFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcfinalise
    
  !
  !================================================================================================================================
  !

  !Initialise the PETSc PC structure
  SUBROUTINE petsc_pcinitialise(pc,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCInitialise",err,error,*999)

    pc%pc=petsc_null_object
    
    exits("Petsc_PCInitialise")
    RETURN
999 errorsexits("Petsc_PCInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcinitialise
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcfactorgetmatrix(pc,factoredMatrix,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    TYPE(petscmattype), INTENT(OUT) :: factoredMatrix
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCFactorGetMatrix",err,error,*999)

    CALL pcfactorgetmatrix(pc%pc,factoredmatrix%mat,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCFactorGetMatrix",err,error,*999)
    ENDIF
    
    exits("Petsc_PCFactorGetMatrix")
    RETURN
999 errorsexits("Petsc_PCFactorGetMatrix",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcfactorgetmatrix

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcfactorsetmatsolverpackage(pc,solverPackage,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    matsolverpackage, INTENT(IN) :: solverpackage 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCFactorSetMatSolverPackage",err,error,*999)

    CALL pcfactorsetmatsolverpackage(pc%pc,solverpackage,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCFactorSetMatSolverPackage.",err,error,*999)
    ENDIF
    
    exits("Petsc_PCFactorSetMatSolverPackage")
    RETURN
999 errorsexits("Petsc_PCFactorSetMatSolverPackage",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcfactorsetmatsolverpackage

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcfactorsetupmatsolverpackage(pc,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCFactorSetUpMatSolverPackage",err,error,*999)

    CALL pcfactorsetupmatsolverpackage(pc%pc,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCFactorSetUpMatSolverPackage.",err,error,*999)
    ENDIF
    
    exits("Petsc_PCFactorSetUpMatSolverPackage")
    RETURN
999 errorsexits("Petsc_PCFactorSetUpMatSolverPackage",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcfactorsetupmatsolverpackage

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcsetreusepreconditioner(pc,flag,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    LOGICAL, INTENT(IN) :: flag
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCSetReusePreconditioner",err,error,*999)

    IF(flag) THEN
      CALL pcsetreusepreconditioner(pc%pc,petsc_true,err)
    ELSE
      CALL pcsetreusepreconditioner(pc%pc,petsc_false,err)
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCSetReusePreconditioner.",err,error,*999)
    ENDIF
    
    exits("Petsc_PCSetReusePreconditioner")
    RETURN
999 errorsexits("Petsc_PCSetReusePreconditioner",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcsetreusepreconditioner
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcsetfromoptions(pc,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCSetFromOptions",err,error,*999)

    CALL pcsetfromoptions(pc%pc,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_PCSetFromOptions")
    RETURN
999 errorsexits("Petsc_PCSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcsetfromoptions
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_pcsettype(pc,method,err,error,*)

    !Argument Variables
    TYPE(petscpctype), INTENT(INOUT) :: pc
    pctype, INTENT(IN) :: method 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_PCSetType",err,error,*999)

    CALL pcsettype(pc%pc,method,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in PCSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_PCSetType")
    RETURN
999 errorsexits("Petsc_PCSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_pcsettype
    
  !
  !================================================================================================================================
  !

  !Finalise the PETSc SNES structure and destroy the SNES
  SUBROUTINE petsc_snesfinalise(snes,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesFinalise",err,error,*999)

    IF(snes%snes/=petsc_null_object) THEN
      CALL petsc_snesdestroy(snes,err,error,*999)
    ENDIF
    
    exits("Petsc_SnesFinalise")
    RETURN
999 errorsexits("Petsc_SnesFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesfinalise
    
  !
  !================================================================================================================================
  !

  !Initialise the PETSc SNES structure
  SUBROUTINE petsc_snesinitialise(snes,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesInitialise",err,error,*999)

    snes%snes=petsc_null_object
     
    exits("Petsc_SnesInitialise")
    RETURN
999 errorsexits("Petsc_SnesInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesinitialise
    
!
  !================================================================================================================================
  !

  SUBROUTINE petsc_snescomputejacobiandefault(snes,x,j,b,ctx,err,error,*)

    !Argument variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscmattype), INTENT(INOUT) :: j
    TYPE(petscmattype), INTENT(INOUT) :: b
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(INOUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesComputeJacobianDefault",err,error,*999)
    
    CALL snescomputejacobiandefault(snes%snes,x%vec,j%mat,b%mat,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESComputeJacobianDefault.",err,error,*999)
    ENDIF

    exits("Petsc_SnesComputeJacobianDefault")    
    RETURN
999 errorsexits("Petsc_SnesComputeJacobianDefault",err,error)
    RETURN
    
  END SUBROUTINE petsc_snescomputejacobiandefault

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snescomputejacobiandefaultcolor(snes,x,j,b,ctx,err,error,*)

    !Argument variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscmattype), INTENT(INOUT) :: j
    TYPE(petscmattype), INTENT(INOUT) :: b
    TYPE(petscmatfdcoloringtype), POINTER :: ctx
    INTEGER(INTG), INTENT(INOUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesComputeJacobianDefaultColor",err,error,*999)
    
    CALL snescomputejacobiandefaultcolor(snes%snes,x%vec,j%mat,b%mat,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESComputeJacobianDefaultColor.",err,error,*999)
    ENDIF

    exits("Petsc_SnesComputeJacobianDefaultColor")
    RETURN
999 errorsexits("Petsc_SnesComputeJacobianDefaultColor",err,error)
    RETURN
    
  END SUBROUTINE petsc_snescomputejacobiandefaultcolor

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snescreate(communicator,snes,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesCreate",err,error,*999)

    CALL snescreate(communicator,snes%snes,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesCreate")
    RETURN
999 errorsexits("Petsc_SnesCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snescreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesdestroy(snes,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesDestroy",err,error,*999)

    CALL snesdestroy(snes%snes,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESDestroy.",err,error,*999)
    ENDIF
    snes%snes=petsc_null_object
     
    exits("Petsc_SnesDestroy")
    RETURN
999 errorsexits("Petsc_SnesDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesdestroy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetapplicationcontext(snes,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetApplicationContext",err,error,*999)

    IF(ASSOCIATED(ctx)) THEN
      CALL flagerror("Context is already associated.",err,error,*999)
    ELSE
      CALL snesgetapplicationcontext(snes%snes,ctx,err)
      IF(err/=0) THEN
        IF(petschandleerror) THEN
          chkerrq(err)
        ENDIF
        CALL flagerror("PETSc error in SNESGetApplicationContext.",err,error,*999)
      ENDIF
    ENDIF
    
    exits("Petsc_SnesGetApplicationContext")
    RETURN
999 errorsexits("Petsc_SnesGetApplicationContext",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetapplicationcontext

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetconvergedreason(snes,reason,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: reason
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetConvergedReason",err,error,*999)

    CALL snesgetconvergedreason(snes%snes,reason,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetConvergedReason.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetConvergedReason")
    RETURN
999 errorsexits("Petsc_SnesGetConvergedReason",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetconvergedreason
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetfunction(snes,f,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(OUT) :: f
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables
    enters("Petsc_SnesGetFunction",err,error,*999)

    CALL snesgetfunction(snes%snes,f%vec,petsc_null_function,petsc_null_integer,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetFunction.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetFunction")
    RETURN
999 errorsexits("Petsc_SnesGetFunction",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetfunction

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetiterationnumber(snes,iterationNumber,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: iterationNumber
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetIterationNumber",err,error,*999)

    CALL snesgetiterationnumber(snes%snes,iterationnumber,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetIterationNumber.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetIterationNumber")
    RETURN
999 errorsexits("Petsc_SnesGetIterationNumber",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetiterationnumber
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetjacobiansolver(snes,a,b,jFunction,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petscmattype), INTENT(INOUT) :: b
    EXTERNAL jfunction 
!     TYPE(SOLVER_TYPE), POINTER :: CTX !<The solver data to pass to the function
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetJacobianSolver",err,error,*999)

    CALL snesgetjacobian(snes%snes,a%mat,b%mat,jfunction,petsc_null_integer,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetJacobian.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetJacobianSolver")
    RETURN
999 errorsexits("Petsc_SnesGetJacobianSolver",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetjacobiansolver

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetjacobianspecial(snes,a,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscmattype), INTENT(INOUT) :: a
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetJacobianSpecial",err,error,*999)

    CALL snesgetjacobian(snes%snes,a%mat,petsc_null_object,petsc_null_function,petsc_null_integer,err)

    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetJacobian.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetJacobianSpecial")
    RETURN
999 errorsexits("Petsc_SnesGetJacobianSpecial",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetjacobianspecial

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetksp(snes,ksp,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetKSP",err,error,*999)

    CALL snesgetksp(snes%snes,ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetKSP.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesGetKSP")
    RETURN
999 errorsexits("Petsc_SnesGetKSP",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetksp

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetlinesearch(snes,lineSearch,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscsneslinesearchtype), INTENT(OUT) :: lineSearch
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesGetLineSearch",err,error,*999)

    CALL snesgetlinesearch(snes%snes,linesearch%snesLineSearch,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetLineSearch.",err,error,*999)
    ENDIF

    exits("Petsc_SnesGetLineSearch")
    RETURN
999 errorsexits("Petsc_SnesGetLineSearch",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetlinesearch

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesgetsolutionupdate(snes,solutionUpdate,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(INOUT) :: solutionUpdate
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesGetSolutionUpdate",err,error,*999)

    CALL snesgetsolutionupdate(snes%snes,solutionupdate%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESGetSolutionUpdate.",err,error,*999)
    ENDIF

    exits("Petsc_SnesGetSolutionUpdate")
    RETURN
999 errorsexits("Petsc_SnesGetSolutionUpdate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesgetsolutionupdate
    
!
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesmonitorset(snes,mFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    EXTERNAL :: mfunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesMonitorSet",err,error,*999)

    CALL snesmonitorset(snes%snes,mfunction,ctx,petsc_null_function,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESMonitorSet.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesMonitorSet")
    RETURN
999 errorsexits("Petsc_SnesMonitorSet",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesmonitorset

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesqnsetrestarttype(snes,rType,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    snesqnrestarttype, INTENT(IN) :: rtype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesQNSetRestartType",err,error,*999)

    CALL snesqnsetrestarttype(snes%snes,rtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESQNSetRestartType.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesQNSetRestartType")
    RETURN
999 errorsexits("Petsc_SnesQNSetRestartType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesqnsetrestarttype

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesqnsetscaletype(snes,sType,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    snesqnscaletype, INTENT(IN) :: stype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesQNSetScaleType",err,error,*999)

    CALL snesqnsetscaletype(snes%snes,stype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESQNSetScaleType.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesQNSetScaleType")
    RETURN
999 errorsexits("Petsc_SnesQNSetScaleType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesqnsetscaletype

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snesqnsettype(snes,qType,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    snesqntype, INTENT(IN) :: qtype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesQNSetType",err,error,*999)

    CALL snesqnsettype(snes%snes,qtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESQNSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesQNSetType")
    RETURN
999 errorsexits("Petsc_SnesQNSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snesqnsettype

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetapplicationcontext(snes,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetApplicationContext",err,error,*999)

    CALL snessetapplicationcontext(snes%snes,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetApplicationContext.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetApplicationContext")
    RETURN
999 errorsexits("Petsc_SnesSetApplicationContext",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetapplicationcontext

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetconvergencetest(snes,cFunction,ctx,err,error,*)
    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    EXTERNAL cfunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetConvergenceTest",err,error,*999)

    CALL snessetconvergencetest(snes%snes,cfunction,ctx,petsc_null_function,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetConvergenceTest.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetConvergenceTest")
    RETURN
999 errorsexits("Petsc_SnesSetConvergenceTest",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetconvergencetest

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetfromoptions(snes,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetFromOptions",err,error,*999)

    CALL snessetfromoptions(snes%snes,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetFromOptions")
    RETURN
999 errorsexits("Petsc_SnesSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetfromoptions
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetfunction(snes,f,fFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(INOUT) :: f
    EXTERNAL ffunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetFunction",err,error,*999)

    CALL snessetfunction(snes%snes,f%vec,ffunction,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetFunction.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetFunction")
    RETURN
999 errorsexits("Petsc_SnesSetFunction",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetfunction

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetjacobianmatfdcoloring(snes,a,b,jFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petscmattype), INTENT(INOUT) :: b
    EXTERNAL jfunction 
    TYPE(petscmatfdcoloringtype) :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetJacobianMatFDColoring",err,error,*999)

    CALL snessetjacobianbuffer(snes,a,b,jfunction,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetJacobian.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetJacobianMatFDColoring")
    RETURN
999 errorsexits("Petsc_SnesSetJacobianMatFDColoring",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetjacobianmatfdcoloring
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetjacobiansolver(snes,a,b,jFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscmattype), INTENT(INOUT) :: a
    TYPE(petscmattype), INTENT(INOUT) :: b
    EXTERNAL jfunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetJacobianSolver",err,error,*999)

    CALL snessetjacobian(snes%snes,a%mat,b%mat,jfunction,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetJacobian.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetJacobianSolver")
    RETURN
999 errorsexits("Petsc_SnesSetJacobianSolver",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetjacobiansolver

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetksp(snes,ksp,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscksptype), INTENT(INOUT) :: ksp
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetKsp",err,error,*999)

    CALL snessetksp(snes%snes,ksp%ksp,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetKSP.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetKsp")
    RETURN
999 errorsexits("Petsc_SnesSetKsp",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetksp

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessetnormschedule(snes,normSchedule,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    snesnormschedule, INTENT(IN) :: normschedule 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetNormSchedule",err,error,*999)

    CALL snessetnormschedule(snes%snes,normschedule,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetNormSchedule.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetNormSchedule")
    RETURN
999 errorsexits("Petsc_SnesSetNormSchedule",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessetnormschedule

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessettolerances(snes,absTol,rTol,sTol,maxIterations,maxFunctionEvals,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    REAL(DP), INTENT(IN) :: absTol
    REAL(DP), INTENT(IN) :: rTol
    REAL(DP), INTENT(IN) :: sTol
    INTEGER(INTG), INTENT(IN) :: maxIterations
    INTEGER(INTG), INTENT(IN) :: maxFunctionEvals
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetTolerances",err,error,*999)

    CALL snessettolerances(snes%snes,abstol,rtol,stol,maxiterations,maxfunctionevals,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetTolerances.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetTolerances")
    RETURN
999 errorsexits("Petsc_SnesSetTolerances",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessettolerances
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessettrustregiontolerance(snes,trTol,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    REAL(DP), INTENT(IN) :: trTol
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetTrustRegionTolerance",err,error,*999)

    CALL snessettrustregiontolerance(snes%snes,trtol,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetTrustRegionTolerance.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetTrustRegionTolerance")
    RETURN
999 errorsexits("Petsc_SnesSetTrustRegionTolerance",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessettrustregiontolerance
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessettype(snes,method,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    snestype, INTENT(IN) :: method 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSetType",err,error,*999)

    CALL snessettype(snes%snes,method,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSetType")
    RETURN
999 errorsexits("Petsc_SnesSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessettype

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_snessolve(snes,b,x,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    TYPE(petscvectype), INTENT(INOUT) :: b
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesSolve",err,error,*999)

    CALL snessolve(snes%snes,b%vec,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESSolve.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesSolve")
    RETURN
999 errorsexits("Petsc_SnesSolve",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_snessolve

  !
  !================================================================================================================================
  !

  !Finalise the PETSc SNES line search structure
  SUBROUTINE petsc_sneslinesearchfinalise(lineSearch,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesLineSearchFinalise",err,error,*999)

    !We don't actually call PETSc's SNESLineSearchDestroy as PETSc accesses and destroys the LineSearch when calling
    !SNESDestroy, so we'll just let PETSc clean it up.
    linesearch%snesLineSearch=petsc_null_object

    exits("Petsc_SnesLineSearchFinalise")
    RETURN
999 errorsexits("Petsc_SnesLineSearchFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchfinalise

  !
  !================================================================================================================================
  !

  !Initialise the PETSc SNES line search structure
  SUBROUTINE petsc_sneslinesearchinitialise(lineSearch,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchInitialise",err,error,*999)
    linesearch%snesLineSearch=petsc_null_object
    exits("Petsc_SnesLineSearchInitialise")

    RETURN
999 errorsexits("Petsc_SnesLineSearchInitialise",err,error)
    RETURN 1
  END SUBROUTINE petsc_sneslinesearchinitialise

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchbtsetalpha(lineSearch,alpha,err,error,*)

    !Argument variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    REAL(DP), INTENT(IN) :: alpha
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesLineSearchBTSetAlpha",err,error,*999)

    CALL sneslinesearchbtsetalpha(linesearch%snesLineSearch,alpha,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      END IF
      CALL flagerror("PETSc error in SNESLineSearchBTSetAlpha.",err,error,*999)
    END IF

    exits("Petsc_SnesLineSearchBTSetAlpha")
    RETURN
999 errorsexits("Petsc_SnesLineSearchBTSetAlpha",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchbtsetalpha

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchcomputenorms(lineSearch,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchComputeNorms",err,error,*999)

    CALL sneslinesearchcomputenorms(linesearch%snesLineSearch,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SnesLineSearchComputeNorms.",err,error,*999)
    ENDIF

    exits("Petsc_SnesLineSearchComputeNorms")
    RETURN
999 errorsexits("Petsc_SnesLineSearchComputeNorms",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchcomputenorms

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchgetnorms(lineSearch,xNorm,fNorm,yNorm,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    REAL(DP), INTENT(INOUT) :: xNorm
    REAL(DP), INTENT(INOUT) :: fNorm
    REAL(DP), INTENT(INOUT) :: yNorm
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchGetNorms",err,error,*999)

    CALL sneslinesearchgetnorms(linesearch%snesLineSearch,xnorm,fnorm,ynorm,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SnesLineSearchGetNorms.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesLineSearchGetNorms")
    RETURN
999 errorsexits("Petsc_SnesLineSearchGetNorms",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchgetnorms

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchgetvecs(lineSearch,x,f,y,w,g,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscvectype), INTENT(INOUT) :: f
    TYPE(petscvectype), INTENT(INOUT) :: y
    TYPE(petscvectype), INTENT(INOUT) :: w
    TYPE(petscvectype), INTENT(INOUT) :: g
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchGetVecs",err,error,*999)

    CALL sneslinesearchgetvecs(linesearch%snesLineSearch,x%vec,f%vec,y%vec,w%vec,g%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESLineSearchGetVecs.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesLineSearchGetVecs")
    RETURN
999 errorsexits("Petsc_SnesLineSearchGetVecs",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchgetvecs

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchsetcomputenorms(lineSearch,computeNorms,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    LOGICAL, INTENT(IN) :: computeNorms
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchSetComputeNorms",err,error,*999)

    IF(computenorms) THEN
      CALL sneslinesearchsetcomputenorms(linesearch%snesLineSearch,petsc_true,err)
    ELSE
      CALL sneslinesearchsetcomputenorms(linesearch%snesLineSearch,petsc_false,err)      
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESLineSearchSetComputeNorms.",err,error,*999)
    ENDIF

    exits("Petsc_SnesLineSearchSetComputeNorms")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetComputeNorms",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchsetcomputenorms

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchsetmonitor(lineSearch,monitorLinesearch,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    LOGICAL, INTENT(IN) :: monitorLinesearch
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchSetMonitor",err,error,*999)

    IF(monitorlinesearch) THEN
      CALL sneslinesearchsetmonitor(linesearch%snesLineSearch,petsc_true,err)
    ELSE
      CALL sneslinesearchsetmonitor(linesearch%snesLineSearch,petsc_false,err)
    ENDIF
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESLineSearchSetMonitor.",err,error,*999)
    ENDIF

    exits("Petsc_SnesLineSearchSetMonitor")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetMonitor",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchsetmonitor

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchsetnorms(snes,xNorm,fNorm,yNorm,err,error,*)

    !Argument Variables
    TYPE(petscsnestype), INTENT(INOUT) :: snes
    REAL(DP), INTENT(INOUT) :: xNorm
    REAL(DP), INTENT(INOUT) :: fNorm
    REAL(DP), INTENT(INOUT) :: yNorm
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SnesLineSearchSetNorms",err,error,*999)

    CALL sneslinesearchsetnorms(snes%snes,xnorm,fnorm,ynorm,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("petsc error in SnesLineSearchSetNorms.",err,error,*999)
    ENDIF
    
    exits("Petsc_SnesLineSearchSetNorms")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetNorms",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchsetnorms

  !
  !================================================================================================================================
  !
  SUBROUTINE petsc_sneslinesearchsetorder(lineSearch,lineSearchOrder,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    sneslinesearchorder, INTENT(IN) :: linesearchorder 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables
    TYPE(varying_string) :: localError

    enters("Petsc_SnesLineSearchSetOrder",err,error,*999)

    SELECT CASE(linesearchorder)
    CASE(petsc_snes_linesearch_order_linear)
      CALL sneslinesearchsetorder(linesearch%snesLineSearch,snes_linesearch_order_linear,err)
    CASE(petsc_snes_linesearch_order_quadratic)
      CALL sneslinesearchsetorder(linesearch%snesLineSearch,snes_linesearch_order_quadratic,err)
    CASE(petsc_snes_linesearch_order_cubic)
      CALL sneslinesearchsetorder(linesearch%snesLineSearch,snes_linesearch_order_cubic,err)
    CASE DEFAULT
      localerror="The specified line search order of "//trim(numbertovstring(linesearchorder,"*",err,error))//" is invalid."
      CALL flagerror(localerror,err,error,*999)
    END SELECT
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESLineSearchSetOrder.",err,error,*999)
    ENDIF

    exits("Petsc_SnesLineSearchSetOrder")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetOrder",err,error)
    RETURN 1
  END SUBROUTINE petsc_sneslinesearchsetorder

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_sneslinesearchsettolerances(lineSearch,steptol,maxstep,rtol,atol,ltol,maxIt,err,error,*)

    !Argument variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    REAL(DP), INTENT(IN) :: steptol
    REAL(DP), INTENT(IN) :: maxstep
    REAL(DP), INTENT(IN) :: rtol
    REAL(DP), INTENT(IN) :: atol
    REAL(DP), INTENT(IN) :: ltol
    INTEGER(INTG), INTENT(IN) :: maxIt
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesLineSearchSetTolerances",err,error,*999)

    CALL sneslinesearchsettolerances(linesearch%snesLineSearch,steptol,maxstep,rtol,atol,ltol,maxit,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      END IF
      CALL flagerror("PETSc error in SNESLineSearchSetTolerances.",err,error,*999)
    END IF

    exits("Petsc_SnesLineSearchSetTolerances")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetTolerances",err,error)
    RETURN 1

  END SUBROUTINE petsc_sneslinesearchsettolerances

  !
  !================================================================================================================================
  !
 
  SUBROUTINE petsc_sneslinesearchsettype(lineSearch,lineSearchType,err,error,*)

    !Argument Variables
    TYPE(petscsneslinesearchtype), INTENT(INOUT) :: lineSearch
    sneslinesearchtype, INTENT(IN) :: linesearchtype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error

    enters("Petsc_SnesLineSearchSetType",err,error,*999)

    CALL sneslinesearchsettype(linesearch%snesLineSearch,linesearchtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in SNESLineSearchSetType.",err,error,*999)
    ENDIF

    exits("Petsc_SnesLineSearchSetType")
    RETURN
999 errorsexits("Petsc_SnesLineSearchSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_sneslinesearchsettype

  !
  !================================================================================================================================
  !

  !Finalise the PETSc TS structure and destroy the TS
  SUBROUTINE petsc_tsfinalise(ts,err,error,*)

    !Argument Variables
    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSFinalise",err,error,*999)

    IF(ts%ts/=petsc_null_object) THEN
      CALL petsc_tsdestroy(ts,err,error,*999)
    ENDIF
    
    exits("Petsc_TSFinalise")
    RETURN
999 errorsexits("Petsc_TSFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsfinalise
    
  !
  !================================================================================================================================
  !

  !Initialise the PETSc TS structure
  SUBROUTINE petsc_tsinitialise(ts,err,error,*)

    !Argument Variables
    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSInitialise",err,error,*999)

    ts%ts=petsc_null_object
     
    exits("Petsc_TSInitialise")
    RETURN
999 errorsexits("Petsc_TSInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsinitialise
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tscreate(communicator,ts,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(INOUT) :: communicator 
    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSCreate",err,error,*999)

    CALL tscreate(communicator,ts%ts,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSCreate")
    RETURN
999 errorsexits("Petsc_TSCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tscreate
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tsdestroy(ts,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSDestroy",err,error,*999)

    CALL tsdestroy(ts%ts,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSDestroy.",err,error,*999)
    ENDIF
    ts%ts=petsc_null_object
    
    exits("Petsc_TSDestroy")
    RETURN
999 errorsexits("Petsc_TSDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsdestroy
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tsgetsolution(ts,currentSolution,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    TYPE(petscvectype), INTENT(INOUT) :: currentSolution
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSGetSolution",err,error,*999)

    CALL tsgetsolution(ts%ts,currentsolution%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSGetSolution.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSGetSolution")
    RETURN
999 errorsexits("Petsc_TSGetSolution",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsgetsolution
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_tsmonitorset(ts,mFunction,ctx,err,error,*)

    !Argument Variables
    TYPE(petsctstype), INTENT(INOUT) :: ts
    EXTERNAL :: mfunction 
    TYPE(solver_type), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSMonitorSet",err,error,*999)

    CALL tsmonitorset(ts%ts,mfunction,ctx,petsc_null_function,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSMonitorSet.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSMonitorSet")
    RETURN
999 errorsexits("Petsc_TSMonitorSet",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsmonitorset
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetduration(ts,maxSteps,maxTime,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(IN) :: maxSteps
    REAL(DP), INTENT(IN) :: maxTime
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetDuration",err,error,*999)

    CALL tssetduration(ts%ts,maxsteps,maxtime,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetDuration.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetDuration")
    RETURN
999 errorsexits("Petsc_TSSetDuration",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetduration
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetexactfinaltime(ts,exactFinalTime,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    LOGICAL, INTENT(IN) :: exactFinalTime
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetExactFinalTime",err,error,*999)

    IF(exactfinaltime) THEN
      CALL tssetexactfinaltime(ts%ts,petsc_true,err)
    ELSE
      CALL tssetexactfinaltime(ts%ts,petsc_false,err)
    ENDIF
    
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetExactFinalTime.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetExactFinalTime")
    RETURN
999 errorsexits("Petsc_TSSetExactFinalTime",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetexactfinaltime

  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetfromoptions(ts,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetFromOptions",err,error,*999)

    CALL tssetfromoptions(ts%ts,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetFromOptions")
    RETURN
999 errorsexits("Petsc_TSSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetfromoptions

  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetinitialtimestep(ts,initialTime,timeStep,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    REAL(DP), INTENT(IN) :: initialTime
    REAL(DP), INTENT(IN) :: timeStep
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetInitialTimeStep",err,error,*999)

    CALL tssetinitialtimestep(ts%ts,initialtime,timestep,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetInitialTimeStep.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetInitialTimeStep")
    RETURN
999 errorsexits("Petsc_TSSetInitialTimeStep",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetinitialtimestep

  !
  !================================================================================================================================
  !
      
  SUBROUTINE petsc_tssetproblemtype(ts,probType,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(IN) :: probType
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetProblemType",err,error,*999)

    CALL tssetproblemtype(ts%ts,probtype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetProblemType.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetProblemType")
    RETURN
999 errorsexits("Petsc_TSSetProblemType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetproblemtype
    
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetrhsfunction(ts,rates,rhsFunction,ctx,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    TYPE(petscvectype), INTENT(INOUT) :: rates
    EXTERNAL rhsfunction 
    TYPE(cellmlpetsccontexttype), POINTER :: ctx
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetRHSFunction",err,error,*999)

    CALL tssetrhsfunction(ts%ts,rates%vec,rhsfunction,ctx,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetRHSFunction.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetRHSFunction")
    RETURN
999 errorsexits("Petsc_TSSetRHSFunction",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetrhsfunction

  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssetsolution(ts,initialSolution,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    TYPE(petscvectype), INTENT(IN) :: initialSolution
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetSolution",err,error,*999)

    CALL tssetsolution(ts%ts,initialsolution%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetSolution.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetSolution")
    RETURN
999 errorsexits("Petsc_TSSetSolution",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssetsolution

  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssettimestep(ts,timeStep,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    REAL(DP), INTENT(IN) :: timeStep
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetTimeStep",err,error,*999)

    CALL tssettimestep(ts%ts,timestep,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetTimeStep.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetTimeStep")
    RETURN
999 errorsexits("Petsc_TSSetTimeStep",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssettimestep
  
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssettype(ts,method,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    tstype, INTENT(IN) :: method 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSetType",err,error,*999)

    CALL tssettype(ts%ts,method,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSetType")
    RETURN
999 errorsexits("Petsc_TSSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssettype
  
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssolve(ts,x,finalTime,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    TYPE(petscvectype), INTENT(INOUT) :: x
    REAL(DP), INTENT(OUT) :: finalTime
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSolve",err,error,*999)

    CALL tssolve(ts%ts,x%vec,finaltime,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSolve.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSolve")
    RETURN
999 errorsexits("Petsc_TSSolve",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssolve
  
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tsstep(ts,steps,pTime,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    INTEGER(INTG), INTENT(IN) :: steps
    REAL(DP), INTENT(IN) :: pTime
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSStep",err,error,*999)

    CALL tsstep(ts%ts,steps,ptime,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSStep.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSStep")
    RETURN
999 errorsexits("Petsc_TSStep",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tsstep

  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssundialssettype(ts,sundialsType,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    tssundialstype, INTENT(IN) :: sundialstype 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSundialsSetType",err,error,*999)

    CALL tssundialssettype(ts%ts,sundialstype,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSundialsSetType.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSundialsSetType")
    RETURN
999 errorsexits("Petsc_TSSundialsSetType",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssundialssettype
  !
  !================================================================================================================================
  !
    
  SUBROUTINE petsc_tssundialssettolerance(ts,absTol,relTol,err,error,*)

    TYPE(petsctstype), INTENT(INOUT) :: ts
    REAL(DP), INTENT(IN) :: absTol
    REAL(DP), INTENT(IN) :: relTol
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_TSSundialsSetTolerance",err,error,*999)

    CALL tssundialssettolerance(ts%ts,abstol,reltol,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in TSSundialsSetTolerance.",err,error,*999)
    ENDIF
    
    exits("Petsc_TSSundialsSetTolerance")
    RETURN
999 errorsexits("Petsc_TSSundialsSetTolerance",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_tssundialssettolerance
  !
  !================================================================================================================================
  !

  !Finalise the PETSc Vec structure and destroy the KSP
  SUBROUTINE petsc_vecfinalise(x,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecFinalise",err,error,*999)

    IF(x%vec/=petsc_null_object) THEN
      CALL petsc_vecdestroy(x,err,error,*999)
    ENDIF
    
    exits("Petsc_VecFinalise")
    RETURN
999 errorsexits("Petsc_VecFinalise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecfinalise
    
  !
  !================================================================================================================================
  !

  !Initialise the PETSc Vec structure
  SUBROUTINE petsc_vecinitialise(x,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecInitialise",err,error,*999)

    x%vec=petsc_null_object
    
    exits("Petsc_VecInitialise")
    RETURN
999 errorsexits("Petsc_VecInitialise",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecinitialise
  
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecassemblybegin(x,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecAssemblyBegin",err,error,*999)

    CALL vecassemblybegin(x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecAssemblyBegin.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecAssemblyBegin")
    RETURN
999 errorsexits("Petsc_VecAssemblyBegin",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecassemblybegin
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecassemblyend(x,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecAssemblyEnd",err,error,*999)

    CALL vecassemblyend(x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecAssemblyEnd.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecAssemblyEnd")
    RETURN
999 errorsexits("Petsc_VecAssemblyEnd",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecassemblyend
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccopy(x,y,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscvectype), INTENT(INOUT) :: y
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCopy",err,error,*999)

    CALL veccopy(x%vec,y%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCopy.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCopy")
    RETURN
999 errorsexits("Petsc_VecCopy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccopy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreate(communicator,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreate",err,error,*999)

    CALL veccreate(communicator,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreate.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreate")
    RETURN
999 errorsexits("Petsc_VecCreate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreateghost(communicator,localN,globalN,numGhosts,ghosts,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalN
    INTEGER(INTG), INTENT(IN) :: numGhosts
    INTEGER(INTG), INTENT(IN) :: ghosts(:)
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateGhost",err,error,*999)

    CALL veccreateghost(communicator,localn,globaln,numghosts,ghosts,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateGhost.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateGhost")
    RETURN
999 errorsexits("Petsc_VecCreateGhost",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreateghost
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreateghostwitharray(communicator,localN,globalN,numGhosts,ghosts,array,x,err,error,*)

   !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalN
    INTEGER(INTG), INTENT(IN) :: numGhosts
    INTEGER(INTG), INTENT(IN) :: ghosts(:)
    REAL(DP), INTENT(OUT) :: array(:)
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateGhostWithArray",err,error,*999)

    CALL veccreateghostwitharray(communicator,localn,globaln,numghosts,ghosts,array,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateGhostWithArray.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateGhostWithArray")
    RETURN
999 errorsexits("Petsc_VecCreateGhostWithArray",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreateghostwitharray
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreatempi(communicator,localN,globalN,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalN
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateMPI",err,error,*999)

    CALL veccreatempi(communicator,localn,globaln,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateMPI.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateMPI")
    RETURN
999 errorsexits("Petsc_VecCreateMPI",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreatempi
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreatempiwitharray(communicator,localN,globalN,array,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalN
    REAL(DP), INTENT(OUT) :: array(:)
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateMPIWithArray",err,error,*999)

    CALL veccreatempiwitharray(communicator,localn,globaln,array,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateMPIWithArray.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateMPIWithArray")
    RETURN
999 errorsexits("Petsc_VecCreateMPIWithArray",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreatempiwitharray
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreateseq(communicator,n,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: n
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateSeq",err,error,*999)

    CALL veccreateseq(communicator,n,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateSeq.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateSeq")
    RETURN
999 errorsexits("Petsc_VecCreateSeq",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreateseq
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_veccreateseqwitharray(communicator,n,array,x,err,error,*)

    !Argument Variables
    mpi_comm, INTENT(IN) :: communicator 
    INTEGER(INTG), INTENT(IN) :: n
    REAL(DP), INTENT(OUT) :: array(:)
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecCreateSeqWithArray",err,error,*999)

    CALL veccreateseqwitharray(communicator,n,array,x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecCreateSeqWithArray.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecCreateSeqWithArray")
    RETURN
999 errorsexits("Petsc_VecCreateSeqWithArray",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_veccreateseqwitharray
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecdestroy(x,err,error,*)

   !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecDestroy",err,error,*999)

    CALL vecdestroy(x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecDestroy.",err,error,*999)
    ENDIF
    x%vec=petsc_null_object
    
    exits("Petsc_VecDestroy")
    RETURN
999 errorsexits("Petsc_VecDestroy",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecdestroy
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecduplicate(x,y,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscvectype), INTENT(OUT) :: y
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecDuplicate",err,error,*999)

    CALL vecduplicate(x%vec,y%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecDuplicate.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecDuplicate")
    RETURN
999 errorsexits("Petsc_VecDuplicate",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecduplicate
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecdot(x,y,dotProduct,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(IN) :: x
    TYPE(petscvectype), INTENT(IN) :: y
    REAL(DP), INTENT(OUT) :: dotProduct
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables
    
    enters("Petsc_VecDot",err,error,*999)
    
    CALL vecdot(x%vec,y%vec,dotproduct,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecDot.",err,error,*999)
    ENDIF

    exits("Petsc_VecDot")
    RETURN
999 errorsexits("Petsc_VecDot",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecdot
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetarrayf90(x,array,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT), TARGET :: x
    REAL(DP), POINTER :: array(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetArrayF90",err,error,*999)

    IF(ASSOCIATED(array)) THEN
      CALL flagerror("Array is already associated",err,error,*999)
    ELSE
      CALL vecgetarrayf90(x%vec,array,err)
      IF(err/=0) THEN
        IF(petschandleerror) THEN
          chkerrq(err)
        ENDIF
        CALL flagerror("PETSc error in VecGetArrayF90.",err,error,*999)
      ENDIF
    ENDIF
    
    exits("Petsc_VecGetArrayF90")
    RETURN
999 errorsexits("Petsc_VecGetArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetarrayf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetarrayreadf90(x,array,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT), TARGET :: x
    REAL(DP), POINTER :: array(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetArrayReadF90",err,error,*999)

    IF(ASSOCIATED(array)) THEN
      CALL flagerror("Array is already associated",err,error,*999)
    ELSE
      CALL vecgetarrayreadf90(x%vec,array,err)
      IF(err/=0) THEN
        IF(petschandleerror) THEN
          chkerrq(err)
        ENDIF
        CALL flagerror("PETSc error in VecGetArrayReadF90.",err,error,*999)
      ENDIF
    ENDIF
    
    exits("Petsc_VecGetArrayReadF90")
    RETURN
999 errorsexits("Petsc_VecGetArrayReadF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetarrayreadf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetlocalsize(x,n,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: n
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetLocalSize",err,error,*999)

    CALL vecgetlocalsize(x%vec,n,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGetLocalSize.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGetLocalSize")
    RETURN
999 errorsexits("Petsc_VecGetLocalSize",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetlocalsize
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetownershiprange(x,low,high,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: low
    INTEGER(INTG), INTENT(OUT) :: high
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetOwnershipRange",err,error,*999)

    CALL vecgetownershiprange(x%vec,low,high,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
        ENDIF
      CALL flagerror("PETSc error in VecGetOwnershipRange.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGetOwnershipRange")
    RETURN
999 errorsexits("Petsc_VecGetOwnershipRange",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetownershiprange

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetsize(x,n,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: n
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetSize",err,error,*999)

    CALL vecgetsize(x%vec,n,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGetSize.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGetSize")
    RETURN
999 errorsexits("Petsc_VecGetSize",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetsize
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecgetvalues(x,n,indices,values,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: indices(:)
    REAL(DP), INTENT(OUT) :: values(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGetValues",err,error,*999)

    CALL vecgetvalues(x%vec,n,indices,values,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGetValues.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGetValues")
    RETURN
999 errorsexits("Petsc_VecGetValues",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecgetvalues
    
!
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecghostgetlocalform(g,l,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: g
    TYPE(petscvectype), INTENT(INOUT) :: l
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGhostGetLocalForm",err,error,*999)

    CALL vecghostgetlocalform(g%vec,l%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGhostGetLocalForm.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGhostGetLocalForm")
    RETURN
999 errorsexits("Petsc_VecGhostGetLocalForm",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecghostgetlocalform
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecghostrestorelocalform(g,l,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: g
    TYPE(petscvectype), INTENT(INOUT) :: l
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGhostRestoreLocalForm",err,error,*999)

    CALL vecghostrestorelocalform(g%vec,l%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGhostRestoreLocalForm.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGhostRestoreLocalForm")
    RETURN
999 errorsexits("Petsc_VecGhostRestoreLocalForm",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecghostrestorelocalform
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecghostupdatebegin(x,insertMode,scatterMode,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    insertmode, INTENT(IN) :: insertmode 
    scattermode, INTENT(IN) :: scattermode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGhostUpdateBegin",err,error,*999)

    CALL vecghostupdatebegin(x%vec,insertmode,scattermode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGhostUpdateBegin.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGhostUpdateBegin")
    RETURN
999 errorsexits("Petsc_VecGhostUpdateBegin",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecghostupdatebegin
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecghostupdateend(x,insertMode,scatterMode,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    insertmode, INTENT(IN) :: insertmode 
    scattermode, INTENT(IN) :: scattermode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecGhostUpdateEnd",err,error,*999)

    CALL vecghostupdateend(x%vec,insertmode,scattermode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecGhostUpdateEnd.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecGhostUpdateEnd")
    RETURN
999 errorsexits("Petsc_VecGhostUpdateEnd",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecghostupdateend
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecnorm(x,normType,norm,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(IN) :: x
    normtype, INTENT(IN) :: normtype 
    REAL(DP), INTENT(OUT) :: norm
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables
    
    enters("Petsc_VecNorm",err,error,*999)
    
    CALL vecnorm(x%vec,normtype,norm,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecNorm.",err,error,*999)
    ENDIF

    exits("Petsc_VecNorm")
    RETURN
999 errorsexits("Petsc_VecNorm",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecnorm

  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecrestorearrayf90(x,array,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    REAL(DP), POINTER :: array(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecRestoreArrayF90",err,error,*999)

    CALL vecrestorearrayf90(x%vec,array,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecRestoreArrayF90.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecRestoreArrayF90")
    RETURN
999 errorsexits("Petsc_VecRestoreArrayF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecrestorearrayf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecrestorearrayreadf90(x,array,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT), TARGET :: x
    REAL(DP), POINTER :: array(:)
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecRestoreArrayReadF90",err,error,*999)

    CALL vecrestorearrayreadf90(x%vec,array,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecRestoreArrayReadF90.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecRestoreArrayReadF90")
    RETURN
999 errorsexits("Petsc_VecRestoreArrayReadF90",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecrestorearrayreadf90
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecscale(x,alpha,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    REAL(DP), INTENT(IN) :: alpha
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecScale",err,error,*999)

    CALL vecscale(x%vec,alpha,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecScale.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecScale")
    RETURN
999 errorsexits("Petsc_VecScale",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecscale
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecset(x,VALUE,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    REAL(DP), INTENT(IN) :: value
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecSet",err,error,*999)

    CALL vecset(x%vec,value,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSet.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecSet")
    RETURN
999 errorsexits("Petsc_VecSet",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecset
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecsetfromoptions(x,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecSetFromOptions",err,error,*999)

    CALL vecsetfromoptions(x%vec,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSetFromOptions.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecSetFromOptions")
    RETURN
999 errorsexits("Petsc_VecSetFromOptions",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecsetfromoptions
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecsetlocaltoglobalmapping(x,isLocalToGlobalMapping,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    TYPE(petscislocaltogloabalmappingtype), INTENT(IN) :: isLocalToGlobalMapping
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecSetLocalToGlobalMapping",err,error,*999)

    CALL vecsetlocaltoglobalmapping(x%vec,islocaltoglobalmapping%isLocalToGlobalMapping,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSetLocalToGlobalMapping.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecSetLocalToGlobalMapping")
    RETURN
999 errorsexits("Petsc_VecSetLocalToGlobalMapping",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecsetlocaltoglobalmapping
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecsetsizes(x,localN,globalN,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(IN) :: localN
    INTEGER(INTG), INTENT(IN) :: globalN
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecSetSizes",err,error,*999)

    CALL vecsetsizes(x%vec,localn,globaln,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSetSizes.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecSetSizes")
    RETURN
999 errorsexits("Petsc_VecSetSizes",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecsetsizes
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecsetvalues(x,n,indices,values,insertMode,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: indices(*)
    REAL(DP), INTENT(IN) :: values(*)
    insertmode, INTENT(IN) :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecSetValues",err,error,*999)

    CALL vecsetvalues(x%vec,n,indices,values,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSetValues.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecSetValues")
    RETURN
999 errorsexits("Petsc_VecSetValues",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecsetvalues
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_setvalueslocal(x,n,indices,values,insertMode,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    INTEGER(INTG), INTENT(IN) :: n
    INTEGER(INTG), INTENT(IN) :: indices(*)
    REAL(DP), INTENT(IN) :: values(*)
    insertmode :: insertmode 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_SetValuesLocal",err,error,*999)

    CALL vecsetvalueslocal(x%vec,n,indices,values,insertmode,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecSetValuesLocal.",err,error,*999)
    ENDIF
    
    exits("Petsc_SetValuesLocal")
    RETURN
999 errorsexits("Petsc_SetValuesLocal",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_setvalueslocal
    
  !
  !================================================================================================================================
  !

  SUBROUTINE petsc_vecview(x,viewer,err,error,*)

    !Argument Variables
    TYPE(petscvectype), INTENT(INOUT) :: x
    petscviewer, INTENT(IN) :: viewer 
    INTEGER(INTG), INTENT(OUT) :: err
    TYPE(varying_string), INTENT(OUT) :: error
    !Local Variables

    enters("Petsc_VecView",err,error,*999)

    CALL vecview(x%vec,viewer,err)
    IF(err/=0) THEN
      IF(petschandleerror) THEN
        chkerrq(err)
      ENDIF
      CALL flagerror("PETSc error in VecView.",err,error,*999)
    ENDIF
    
    exits("Petsc_VecView")
    RETURN
999 errorsexits("Petsc_VecView",err,error)
    RETURN 1
    
  END SUBROUTINE petsc_vecview


  !
  !================================================================================================================================
  !

END MODULE cmisspetsc
    
SUBROUTINE snessetjacobianbuffer(snes,A,B,jFunction,matFDColoring,err)

  USE cmisspetsctypes
  USE kinds

  IMPLICIT NONE
  
  !Argument Variables
  TYPE(petscsnestype), INTENT(INOUT) :: snes
  TYPE(petscmattype), INTENT(INOUT) :: a
  TYPE(petscmattype), INTENT(INOUT) :: b
  EXTERNAL jfunction 
  TYPE(petscmatfdcoloringtype) :: matFDColoring
  INTEGER(INTG), INTENT(OUT) :: err

  CALL snessetjacobian(snes%snes,a%mat,b%mat,jfunction,matfdcoloring%matFDColoring,err)

END SUBROUTINE snessetjacobianbuffer