003_Fisher_manipulation.f90

Go to the documentation of this file.

!----------------------------------------------------------------------------------------
!
! This file is part of CosmicFish.
!
! Copyright (C) 2015-2016 by the CosmicFish authors
!
! The CosmicFish code is free software;
! You can use it, redistribute it, and/or modify it under the terms
! of the GNU General Public License as published by the Free Software Foundation;
! either version 3 of the License, or (at your option) any later version.
! The full text of the license can be found in the file LICENSE at
! the top level of the CosmicFish distribution.
!
!----------------------------------------------------------------------------------------


!----------------------------------------------------------------------------------------


module fisher_manipulation

    use precision
    use matrix_utilities
    use cosmicfish_utilities

    implicit none

    private

    public fisher_delete_column_rows, confidence_coefficient, confidence_interval_params, &
        & marginalise_fisher, fisher_pca, fisher_fom, fisher_protection_unconstrained,   &
        & fisher_protection_degenerate

contains

    ! ---------------------------------------------------------------------------------------------
    subroutine fisher_delete_column_rows( Fisher_in, Fisher_out, indexes )

        implicit none

        real(dl), dimension(:,:), intent(in)  :: Fisher_in
        real(dl), dimension(:,:), intent(out) :: Fisher_out
        integer , dimension(:)  , intent(in)  :: indexes

        fisher_out = 0._dl

        stop 'Fisher_delete_column_rows not yet implemented'

    end subroutine fisher_delete_column_rows

    ! ---------------------------------------------------------------------------------------------
    function confidence_coefficient( confidence_level )

        implicit none

        real(dl), intent(in) :: confidence_level
        real(dl) :: confidence_coefficient

        confidence_coefficient = sqrt(2._dl)*dierf(confidence_level)

    end function confidence_coefficient

    ! ---------------------------------------------------------------------------------------------
    subroutine confidence_interval_params( Fisher_in, confidence_level, error_out )

        implicit none

        real(dl), intent(in)   :: Fisher_in(:,:)
        real(dl), intent(in)   :: confidence_level
        real(dl), intent(out)  :: error_out(:)

        error_out = 0._dl

        stop 'confidence_interval_params not yet implemented'

    end subroutine confidence_interval_params

    ! ---------------------------------------------------------------------------------------------
    subroutine marginalise_fisher( Fisher_in, indexes, Fisher_out )

        implicit none

        real(dl), dimension(:,:), intent(in)  :: Fisher_in
        real(dl), dimension(:,:), intent(out) :: Fisher_out
        integer , dimension(:)  , intent(in)  :: indexes

        fisher_out = 0._dl

        stop 'marginalise_Fisher not yet implemented'

    end subroutine marginalise_fisher

    ! ---------------------------------------------------------------------------------------------
    subroutine fisher_pca( Fisher_in, eigenvalues, eigenvectors )

        implicit none

        real(dl), dimension(:,:), intent(in)  :: Fisher_in
        real(dl), dimension(:)  , intent(out) :: eigenvalues
        real(dl), dimension(:,:), intent(out) :: eigenvectors

        eigenvalues  = 0._dl
        eigenvectors = 0._dl

        stop 'Fisher_PCA not yet implemented'

    end subroutine fisher_pca

    ! ---------------------------------------------------------------------------------------------
    subroutine fisher_fom( Fisher_in, FOM )

        implicit none

        real(dl), dimension(:,:), intent(in)  :: Fisher_in
        real(dl), intent(out) :: FOM

        fom = 0._dl

        stop 'Fisher_FOM not yet implemented'

    end subroutine fisher_fom

    ! ---------------------------------------------------------------------------------------------
    subroutine fisher_protection_unconstrained( Fisher_in )

        implicit none

        real(dl), dimension(:,:), intent(inout)  :: Fisher_in

        integer :: fisher_dim, i, j
        logical :: is_zero

        fisher_dim = size( fisher_in, 1 )

        ! cycle over the Fisher matrix to see wether there are columns or row that are
        ! exactly zero.
        do i=1, fisher_dim

            is_zero = .true.

            do j=1, fisher_dim
                if ( fisher_in(i,j) /= 0._dl ) is_zero = .false.
            end do

            if ( is_zero ) then
                fisher_in(i,i) = 1.d-16
            end if

        end do

    end subroutine fisher_protection_unconstrained

    ! ---------------------------------------------------------------------------------------------
    subroutine fisher_protection_degenerate( Fisher_in )

        implicit none

        real(dl), dimension(:,:), intent(inout)  :: Fisher_in

        integer :: fisher_dim, i

        fisher_dim = size( fisher_in, 1 )

        ! cycle over the Fisher matrix
        do i=1, fisher_dim
            fisher_in(i,i) = fisher_in(i,i) + 1.d-16
        end do

    end subroutine fisher_protection_degenerate

    ! ---------------------------------------------------------------------------------------------

end module fisher_manipulation