procedures_changed_boxes_size_ratio_factory.f90 Source File

module procedures_changed_boxes_size_ratio_factory

use, intrinsic :: iso_fortran_env, only: DP => REAL64
use json_module, only: json_file
use procedures_errors, only: error_exit
use procedures_checks, only: check_data_found
use classes_periodic_box, only: Abstract_Periodic_Box
use procedures_environment_inquirers, only: periodicity_is_xyz, periodicity_is_xy
use classes_changed_box_size_ratio, only: Abstract_Changed_Box_Size_Ratio, &
    XYZ_Changed_Box_Size_Ratio, XY_Changed_Box_Size_Ratio, Null_Changed_Box_Size_Ratio

implicit none

private
public :: create, destroy

interface create
    module procedure :: create_line
    module procedure :: create_element
end interface create

interface destroy
    module procedure :: destroy_element
    module procedure :: destroy_line
end interface destroy

contains

    subroutine create_line(changed_boxes_size_ratio, periodic_boxes, boxes_size_can_change)
        class(Abstract_Changed_Box_Size_Ratio), allocatable, intent(out) :: &
            changed_boxes_size_ratio(:)
        class(Abstract_Periodic_Box), intent(in) :: periodic_boxes(:)
        logical, intent(in) :: boxes_size_can_change(:)

        if (all(boxes_size_can_change)) then
            if (all(periodicity_is_xyz(periodic_boxes))) then
                allocate(XYZ_Changed_Box_Size_Ratio :: &
                    changed_boxes_size_ratio(size(periodic_boxes)))
            else if (all(periodicity_is_xy(periodic_boxes))) then
                allocate(XY_Changed_Box_Size_Ratio :: &
                    changed_boxes_size_ratio(size(periodic_boxes)))
            else
                call error_exit("procedures_changed_boxes_size_ratio_factory: create_line: "//&
                        "box periodicity is unknown.")
            end if
        else
            allocate(Null_Changed_Box_Size_Ratio :: changed_boxes_size_ratio(size(periodic_boxes)))
        end if
    end subroutine create_line

    subroutine destroy_line(changed_boxes_size_ratio)
        class(Abstract_Changed_Box_Size_Ratio), allocatable, intent(inout) :: &
            changed_boxes_size_ratio(:)

        if (allocated(changed_boxes_size_ratio)) deallocate(changed_boxes_size_ratio)
    end subroutine destroy_line

    subroutine create_element(changed_box_size_ratio, periodic_box, box_size_can_change, &
        generating_data, prefix)
        class(Abstract_Changed_Box_Size_Ratio), allocatable, intent(out) :: changed_box_size_ratio
        class(Abstract_Periodic_Box), intent(in) :: periodic_box
        logical, intent(in) :: box_size_can_change
        type(json_file), intent(inout) :: generating_data
        character(len=*), intent(in) :: prefix

        call allocate(changed_box_size_ratio, periodic_box, box_size_can_change)
        call construct(changed_box_size_ratio, generating_data, prefix)
    end subroutine create_element

    subroutine allocate(changed_box_size_ratio, periodic_box, box_size_can_change)
        class(Abstract_Changed_Box_Size_Ratio), allocatable, intent(out) :: changed_box_size_ratio
        class(Abstract_Periodic_Box), intent(in) :: periodic_box
        logical, intent(in) :: box_size_can_change

        if (box_size_can_change) then
            if (periodicity_is_xyz(periodic_box)) then
                allocate(XYZ_Changed_Box_Size_Ratio :: changed_box_size_ratio)
            else if (periodicity_is_xy(periodic_box)) then
                allocate(XY_Changed_Box_Size_Ratio :: changed_box_size_ratio)
            else
                call error_exit("procedures_changed_boxes_size_ratio_factory: allocate: "//&
                        "box periodicity is unknown.")
            end if
        else
            allocate(Null_Changed_Box_Size_Ratio :: changed_box_size_ratio)
        end if
    end subroutine allocate

    subroutine construct(changed_box_size_ratio, generating_data, prefix)
        class(Abstract_Changed_Box_Size_Ratio), intent(inout) :: changed_box_size_ratio
        type(json_file), intent(inout) :: generating_data
        character(len=*), intent(in) :: prefix

        real(DP) :: initial_delta
        logical :: is_xyz_or_xy
        character(len=:), allocatable :: data_field
        logical :: data_found

        select type (changed_box_size_ratio)
            type is (XYZ_Changed_Box_Size_Ratio)
                is_xyz_or_xy = .true.
            type is (XY_Changed_Box_Size_Ratio)
                is_xyz_or_xy = .true.
            type is (Null_Changed_Box_Size_Ratio)
                is_xyz_or_xy = .false.
            class default
                call error_exit("procedures_changed_boxes_size_ratio_factory: construct: "//&
                    "changed_box_size: type unknown.")
        end select
        if (is_xyz_or_xy) then
            data_field = prefix//"initial delta"
            call generating_data%get(data_field, initial_delta, data_found)
            call check_data_found(data_field, data_found)
            call changed_box_size_ratio%set(initial_delta)
        end if
    end subroutine construct

    subroutine destroy_element(changed_box_size_ratio)
        class(Abstract_Changed_Box_Size_Ratio), allocatable, intent(inout) :: changed_box_size_ratio

        if (allocated(changed_box_size_ratio)) deallocate(changed_box_size_ratio)
    end subroutine destroy_element

end module procedures_changed_boxes_size_ratio_factory