poisson_no.F90

Go to the documentation of this file.

!! Copyright (C) 2002-2011 M. Marques, A. Castro, A. Rubio, G. Bertsch, M. Oliveira
!!
!! 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 poisson_no_oct_m
  use debug_oct_m
  use global_oct_m
  use mesh_oct_m
 
  implicit none
  private
  public ::                  &
    poisson_no_t,           &
    poisson_no_init,        &
    poisson_no_end,         &
    poisson_no_solve
 
  type poisson_no_t
    private
    !integer         :: all_nodes_comm
  end type poisson_no_t
contains
 
  subroutine poisson_no_init(this)
    type(poisson_no_t),  intent(out)   :: this
! may need to add these later for housekeeping in no poisson case. Otherwise delete these 2 lines and
!  type member above
!    integer,             intent(in)    :: all_nodes_comm
!    logical, optional,   intent(in)    :: init_world
 
 
    push_sub(poisson_no_init)
 
    !this%all_nodes_comm = all_nodes_comm
 
    pop_sub(poisson_no_init)
  end subroutine poisson_no_init
 
  !-----------------------------------------------------------------
  subroutine poisson_no_end(this)
    type(poisson_no_t), intent(inout) :: this
 
    push_sub(poisson_no_end)
 
! nothing to do - only integer objects in poisson_no_t
 
    pop_sub(poisson_no_end)
  end subroutine poisson_no_end
 
  !-----------------------------------------------------------------
 
  subroutine poisson_no_solve(this, mesh, pot, rho)
    type(poisson_no_t),             intent(in)    :: this
    type(mesh_t),                   intent(in)    :: mesh
    real(real64),                   intent(out)   :: pot(:)
    real(real64),                   intent(in)    :: rho(:)
 
    push_sub(poisson_no_solve)
 
    pot(1:mesh%np) = m_zero
 
    pop_sub(poisson_no_solve)
  end subroutine poisson_no_solve
 
end module poisson_no_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: