pm_clusTest.F90

Go to the documentation of this file.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!                                                                                                                            !!!!
!!!!    ParaMonte: Parallel Monte Carlo and Machine Learning Library.                                                           !!!!
!!!!                                                                                                                            !!!!
!!!!    Copyright (C) 2012-present, The Computational Data Science Lab                                                          !!!!
!!!!                                                                                                                            !!!!
!!!!    This file is part of the ParaMonte library.                                                                             !!!!
!!!!                                                                                                                            !!!!
!!!!    LICENSE                                                                                                                 !!!!
!!!!                                                                                                                            !!!!
!!!!       https://github.com/cdslaborg/paramonte/blob/main/LICENSE.md                                                          !!!!
!!!!                                                                                                                            !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
module pm_clusTest
 
    use pm_err, only: err_type
    use pm_kind, only: SK, IK, LK, RKG => RK
    use pm_distUnif, only: rngf, rngf_type, xoshiro256ssw_type, rngx_type => xoshiro256ssw_type
 
    implicit none
 
    character(*, SK), parameter :: MODULE_NAME = "@pm_clusTest"
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
    type :: range_type
        integer(IK) :: ndim(2) = [2, 2]
        integer(IK) :: nell(2) = [1, 20]
        integer(IK) :: nsam(2) = [5000, 5000]
        real(RKG)   :: mean(2) = [-1._RKG, 1._RKG]
        real(RKG)   :: std(2) = [0.001_RKG, 1._RKG]
    end type
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
    type :: mmvue_type
        type(range_type)                :: range
        integer(IK)                     :: ndim
        integer(IK)                     :: nell
        integer(IK)                     :: nsam
        integer(IK)                     :: nsim = 10000
        real(RKG)                       :: logVolUnitBall
        real(RKG)                       :: logDensity
        real(RKG)                       :: logSumVolNormedEff
        real(RKG)       , allocatable   :: std(:,:)
        real(RKG)       , allocatable   :: sample(:,:)
        real(RKG)       , allocatable   :: mean(:,:)
        real(RKG)       , allocatable   :: invGram(:,:,:)
        real(RKG)       , allocatable   :: choLowGramUpp(:,:,:)
        real(RKG)       , allocatable   :: cumPropVolNormed(:)
        real(RKG)       , allocatable   :: logVolNormed(:)
        real(RKG)       , allocatable   :: mahalSq(:,:)
        real(RKG)       , allocatable   :: invmul(:)
        integer(IK)     , allocatable   :: size(:)
        integer(IK)     , allocatable   :: membership(:)
        type(err_type)                  :: err
    contains
        procedure, pass :: write => mmvue_type_write
    end type
 
    interface mmvue_type
        module procedure :: mmvue_typer_rngf, mmvue_typer_rngx
    end interface
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
    interface mmvue_typer
 
    impure module function mmvue_typer_rngf(rng, range, nsim) result(self)
#if INTEL_COMPILER_ENABLED && DLL_ENABLED && (WINDOWS_ENABLED || DARWIN_ENABLED)
        !DEC$ ATTRIBUTES DLLEXPORT :: mmvue_typer_rngf
#endif
        type(range_type), intent(in), optional :: range
        integer(IK), intent(in), optional :: nsim
        type(rngf_type), intent(in) :: rng
        type(mmvue_type) :: self
    end function
 
    impure module function mmvue_typer_rngx(rng, range, nsim) result(self)
#if INTEL_COMPILER_ENABLED && DLL_ENABLED && (WINDOWS_ENABLED || DARWIN_ENABLED)
        !DEC$ ATTRIBUTES DLLEXPORT :: mmvue_typer_rngx
#endif
        type(range_type), intent(in), optional :: range
        integer(IK), intent(in), optional :: nsim
        type(rngx_type), intent(inout) :: rng
        type(mmvue_type) :: self
    end function
 
    end interface
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
    interface mmvue_type_write
 
    impure module subroutine mmvue_type_write(self, funit)
#if INTEL_COMPILER_ENABLED && DLL_ENABLED && (WINDOWS_ENABLED || DARWIN_ENABLED)
        !DEC$ ATTRIBUTES DLLEXPORT :: mmvue_type_write
#endif
        class(mmvue_type), intent(in) :: self
        integer(IK), intent(in) :: funit
    end subroutine
 
    end interface
 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
end module pm_clusTest