matter.F90

Go to the documentation of this file.

!! Copyright (C) 2021 Micael 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 matter_oct_m
  use debug_oct_m
  use electrons_oct_m
  use global_oct_m
  use ions_oct_m
  use messages_oct_m
  use mpi_oct_m
  use multisystem_oct_m
  use namespace_oct_m
  use profiling_oct_m
  implicit none
 
  private
  public ::               &
    matter_t
 
  type, extends(multisystem_t) :: matter_t
    private
    class(electrons_t), pointer :: electrons => null()
    class(ions_t),      pointer :: ions => null()
  contains
    final :: matter_finalizer
  end type matter_t
 
  interface matter_t
    procedure matter_constructor
  end interface matter_t
 
contains
 
  ! ---------------------------------------------------------------------------------------
  function matter_constructor(namespace, grp) result(matter)
    type(namespace_t), intent(in) :: namespace
    type(mpi_grp_t),   intent(in) :: grp
    class(matter_t),   pointer    :: matter
 
    push_sub(matter_constructor)
 
    safe_allocate(matter)
 
    matter%namespace = namespace
 
    matter%ions => ions_t(namespace_t("ions", parent=matter%namespace), grp)
    matter%electrons => electrons_t(namespace_t("electrons", parent=matter%namespace), grp)
 
    call matter%list%add(matter%ions)
    call matter%list%add(matter%electrons)
 
    pop_sub(matter_constructor)
  end function matter_constructor
 
  ! ---------------------------------------------------------
  subroutine matter_finalizer(this)
    type(matter_t), intent(inout) :: this
 
    push_sub(matter_finalizer)
 
    call multisystem_end(this)
    nullify(this%electrons)
    nullify(this%ions)
 
    pop_sub(matter_finalizer)
  end subroutine matter_finalizer
 
end module matter_oct_m