sort.F90

Go to the documentation of this file.

!! Copyright (C) 2002-2006 M. Marques, A. Castro, A. Rubio, G. Bertsch
!! Copyright (C) 2021 S. Ohlmann
!!
!! 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 sort_oct_m
  use debug_oct_m
  use global_oct_m
 
  implicit none
 
  private
  public ::                     &
    sort
 
  ! ------------------------------------------------------------------------------
  interface sort
    module procedure dsort, isort, lsort
    module procedure dshellsort1, zshellsort1, ishellsort1
    module procedure dshellsort2, zshellsort2, ishellsort2
  end interface sort
 
  interface
    subroutine dsort1(size, array)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: size
      real(real64),   intent(inout) :: array(*)
    end subroutine dsort1
 
    subroutine dsort2(size, array, indices)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: size
      real(real64),   intent(inout) :: array(*)
      integer, intent(out)   :: indices(*)
    end subroutine dsort2
 
    subroutine isort1(size, array)
      implicit none
      integer, intent(in)    :: size
      integer, intent(inout) :: array(*)
    end subroutine isort1
 
    subroutine isort2(size, array, indices)
      implicit none
      integer, intent(in)    :: size
      integer, intent(inout) :: array(*)
      integer, intent(out)   :: indices(*)
    end subroutine isort2
 
    subroutine lsort1(size, array)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,    intent(in)    :: size
      integer(int64), intent(inout) :: array(*)
    end subroutine lsort1
 
    subroutine lsort2(size, array, indices)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,    intent(in)    :: size
      integer(int64), intent(inout) :: array(*)
      integer,    intent(out)   :: indices(*)
    end subroutine lsort2
  end interface
 
contains
 
  ! ---------------------------------------------------------
  subroutine dsort(a, ind)
    real(real64),      intent(inout) :: a(:)
    integer, optional, intent(out)   :: ind(:)
 
    push_sub(dsort)
 
    if (size(a) > 0) then
 
      if (.not. present(ind)) then
        call dsort1(size(a), a)
      else
        call dsort2(size(a), a, ind)
      end if
 
    end if
 
    pop_sub(dsort)
  end subroutine dsort
 
 
  ! ---------------------------------------------------------
  subroutine isort(a, ind)
    integer,           intent(inout) :: a(:)
    integer, optional, intent(out)   :: ind(:)
 
    push_sub(isort)
 
    if (size(a) > 0) then
 
      if (.not. present(ind)) then
        call isort1(size(a), a)
      else
        call isort2(size(a), a, ind)
      end if
 
    end if
 
    pop_sub(isort)
  end subroutine isort
 
  ! ---------------------------------------------------------
  subroutine lsort(a, ind)
    integer(int64),       intent(inout) :: a(:)
    integer, optional, intent(out)   :: ind(:)
 
    push_sub(lsort)
 
    if (size(a) > 0) then
 
      if (.not. present(ind)) then
        call lsort1(size(a), a)
      else
        call lsort2(size(a), a, ind)
      end if
 
    end if
 
    pop_sub(lsort)
  end subroutine lsort
 
 
#include "undef.F90"
#include "complex.F90"
#include "sort_inc.F90"
 
#include "undef.F90"
#include "real.F90"
#include "sort_inc.F90"
 
#include "undef.F90"
#include "integer.F90"
#include "sort_inc.F90"
 
end module sort_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: