sll_m_poisson_2d_base.F90

Go to the documentation of this file.

 
module sll_m_poisson_2d_base
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#include "sll_working_precision.h"
 
   implicit none
 
   public :: &
      sll_c_poisson_2d_base, &
      sll_i_function_of_position
 
   private
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
   type, abstract :: sll_c_poisson_2d_base
 
   contains
 
      procedure(signature_compute_phi_from_rho_2d), deferred, pass(poisson) :: &
         compute_phi_from_rho
 
      procedure(signature_compute_e_from_rho_2d), deferred, pass(poisson) :: &
         compute_e_from_rho
 
      procedure(signature_norm_squared), deferred :: &
         l2norm_squared
      procedure(signature_update_dofs_function), deferred :: &
         compute_rhs_from_function
 
      procedure(signature_empty), deferred :: &
         free
 
   end type sll_c_poisson_2d_base
 
   abstract interface
 
      function sll_i_function_of_position(x)
         use sll_m_working_precision
         sll_real64             :: sll_i_function_of_position
         sll_real64, intent(in) :: x(:)
      end function sll_i_function_of_position
   end interface
 
#ifndef DOXYGEN_SHOULD_SKIP_THIS
 
   abstract interface
 
      ! solves -\Delta phi = rho in 2d
 
      subroutine signature_compute_phi_from_rho_2d(poisson, phi, rho)
 
         use sll_m_working_precision
         import sll_c_poisson_2d_base
 
         class(sll_c_poisson_2d_base)           :: poisson
         sll_real64, dimension(:, :), intent(in)  :: rho
         sll_real64, dimension(:, :), intent(out) :: phi
 
      end subroutine signature_compute_phi_from_rho_2d
 
   end interface
 
   abstract interface
      ! solves E = -\nabla Phi with -\Delta phi = rho in 2d
      subroutine signature_compute_e_from_rho_2d(poisson, E1, E2, rho)
 
         use sll_m_working_precision
         import sll_c_poisson_2d_base
 
         class(sll_c_poisson_2d_base)            :: poisson
         sll_real64, dimension(:, :), intent(in)  :: rho
         sll_real64, dimension(:, :), intent(out) :: e1
         sll_real64, dimension(:, :), intent(out) :: e2
 
      end subroutine signature_compute_e_from_rho_2d
   end interface
 
   !---------------------------------------------------------------------------!
   abstract interface
      function signature_norm_squared(poisson, coefs_dofs) result(r)
         use sll_m_working_precision
         import sll_c_poisson_2d_base
         class(sll_c_poisson_2d_base), intent(in) :: poisson
         sll_real64, intent(in)                    :: coefs_dofs(:, :) 
         sll_real64                                 :: r
      end function signature_norm_squared
   end interface
 
   !---------------------------------------------------------------------------!
   abstract interface
      subroutine signature_update_dofs_function(poisson, func, coefs_dofs)
         use sll_m_working_precision
         import sll_c_poisson_2d_base
         import sll_i_function_of_position
         class(sll_c_poisson_2d_base)            :: poisson
         procedure(sll_i_function_of_position)    :: func
         sll_real64, intent(out)                  :: coefs_dofs(:) 
      end subroutine signature_update_dofs_function
   end interface
 
   abstract interface
 
      ! Destructor
 
      subroutine signature_empty(poisson)
         import sll_c_poisson_2d_base
         class(sll_c_poisson_2d_base)   :: poisson
 
      end subroutine signature_empty
 
   end interface
 
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
 
end module sll_m_poisson_2d_base