solvers.F90

Go to the documentation of this file.

!! Copyright (C) 2002-2006 M. Marques, A. Castro, A. Rubio, G. Bertsch
!!
!! This program is free software; you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation; either version 2, or (at your option)
!! any later version.
!!
!! This program is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with this program; if not, write to the Free Software
!! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
!! 02110-1301, USA.
!!
 
#include "global.h"
 
module solvers_oct_m
  use blas_oct_m
  use debug_oct_m
  use global_oct_m
  use iso_c_binding
  use, intrinsic :: iso_fortran_env
  use lalg_basic_oct_m
  use loct_oct_m
  use math_oct_m
  use messages_oct_m
  use profiling_oct_m
 
  implicit none
 
  private
  public ::                     &
    dconjugate_gradients,       &
    zconjugate_gradients,       &
    zqmr_sym_gen_dotu,          &
    dqmr_sym_gen_dotu,          &
    zqmr_gen_dotu,              &
    dqmr_gen_dotu,              &
    didrs,                      &
    zidrs
 
  abstract interface
    subroutine doperator_i(x, y, userdata)
      use, intrinsic :: iso_fortran_env
      use, intrinsic :: iso_c_binding, only: c_ptr
      implicit none
      real(real64), contiguous, intent(in)  :: x(:)
      real(real64), contiguous, intent(out) :: y(:)
      type(c_ptr),              intent(in)  :: userdata(:)
    end subroutine doperator_i
    subroutine zoperator_i(x, y, userdata)
      use, intrinsic :: iso_fortran_env
      use, intrinsic :: iso_c_binding, only: c_ptr
      implicit none
      complex(real64), contiguous, intent(in)  :: x(:)
      complex(real64), contiguous, intent(out) :: y(:)
      type(c_ptr),                 intent(in)  :: userdata(:)
    end subroutine zoperator_i
    function ddotp_i(x, y)
      import real64
      implicit none
      real(real64) :: ddotp_i
      real(real64), intent(in) :: x(:)
      real(real64), intent(in) :: y(:)
    end function ddotp_i
    function dnrm_i(x)
      import real64
      implicit none
      real(real64), intent(in) :: x(:)
      real(real64) :: dnrm_i
    end function dnrm_i
    function zdotp_i(x, y)
      import real64
      implicit none
      complex(real64), intent(in) :: x(:)
      complex(real64), intent(in) :: y(:)
      complex(real64) :: zdotp_i
    end function zdotp_i
    function znrm_i(x)
      import real64
      implicit none
      complex(real64), intent(in) :: x(:)
      real(real64) :: znrm_i
    end function znrm_i
  end interface
 
 
 
 
  interface dconjugate_gradients
    module procedure dsym_conjugate_gradients, dbi_conjugate_gradients
  end interface dconjugate_gradients
 
  interface zconjugate_gradients
    module procedure zsym_conjugate_gradients, zbi_conjugate_gradients
  end interface zconjugate_gradients
 
contains
 
#include "undef.F90"
#include "complex.F90"
#include "solvers_inc.F90"
 
#include "undef.F90"
#include "real.F90"
#include "solvers_inc.F90"
 
end module solvers_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: