next up previous contents
Next: 3.6 Matrix-Free Implementation Up: 3. JTpack90 Previous: 3.4 ELL ``Class'' Module   Contents

3.5 CG Module

module JT_CG_module
  implicit none
  interface JT_CG
     module procedure CG_Full
     module procedure CG_ELL
  end interface
  private
  public :: JT_CG

contains

  subroutine CG_Full (status,b,control,x,cpu,rnormt,errt,rnorm,err,a,ap)
    use JT_Full_module

#include ``CG-use.F90''

    real(kind=JT_real_type), intent(in), dimension(:,:) :: a
    real(kind=JT_real_type), intent(inout), dimension(:,:) :: ap

#include ``CG-guts.F90''

  end subroutine CG_Full

  subroutine CG_ELL (status,b,control,x,cpu,rnormt,errt,rnorm,err,a,ap)
    use JT_ELL_module

#include ``CG-use.F90''

    type(JT_ELL_matrix), intent(in) :: a
    type(JT_ELL_matrix), intent(inout) :: ap

#include ``CG-guts.F90''

  end subroutine CG_ELL

end module JT_CG_module



next up previous contents
Next: 3.6 Matrix-Free Implementation Up: 3. JTpack90 Previous: 3.4 ELL ``Class'' Module   Contents
John A. Turner