fisher_matrix_calculator.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.
!
!----------------------------------------------------------------------------------------



program fisher_matrix_calculator

    use precision
    use cosmicfish_types
    use init_from_file
    use modelparams
    use amlutils
    use inifile
    use cosmicfish_camb_interface
    use fisher_calculator_cls
    use fisher_calculator_sn
    use fisher_calculator_rd
    use fisher_calculator_derived
    use fisher_manipulation
    use cosmicfish_utilities

    implicit none

    Type(cambparams)       :: P
    Type(cosmicfish_params) :: FP
    character(LEN=Ini_max_string_len) paramfile

    integer  :: num_param, num_derived
    real(dl), allocatable, dimension(:,:) :: Fisher_matrix, Derived_Matrix
    real(dl), allocatable, dimension(:)   :: fiducial_params

    ! feedback:
    call cosmicfish_write_header( name=' CosmicFish Fisher matrix calculator. ' )

    ! read the parameter file name:
    paramfile = ''
    if (getparamcount() /= 0)  paramfile = getparam(1)
    if (paramfile == '') stop 'No parameter input file'
    ! initialize parameters:
    call init_cosmicfish_from_file( p, fp, paramfile )
    ! get the number of parameters:
    call num_param_fisher(p, fp, num_param)
    ! allocation:
    allocate( fisher_matrix(num_param,num_param) )

    fisher_matrix = 0._dl

    ! Cls Fisher matrix:
    if ( fp%cosmicfish_want_cls ) then
        ! compute the Fisher matrix:
        call fisher_cls( p, fp, num_param, fisher_matrix, outroot=fp%outroot )
        ! save it to file:
        call save_fisher_to_file( p, fp, num_param, fisher_matrix, filename=fp%outroot//'fisher_matrix_cls.dat' )
        ! save the parameter names to file:
        call save_paramnames_to_file( p, fp, num_param, filename=fp%outroot//'fisher_matrix_cls.paramnames' )
    end if

    if ( fp%cosmicfish_want_SN  ) then
        ! compute the Fisher matrix:
        call fisher_sn( p, fp, num_param, fisher_matrix, outroot=fp%outroot )
        ! save it to file:
        call save_fisher_to_file( p, fp, num_param, fisher_matrix, filename=fp%outroot//'fisher_matrix_SN.dat' )
        ! save the parameter names to file:
        call save_paramnames_to_file( p, fp, num_param, filename=fp%outroot//'fisher_matrix_SN.paramnames' )
    end if

    if ( fp%cosmicfish_want_RD  ) then
        ! compute the Fisher matrix:
        call fisher_rd( p, fp, num_param, fisher_matrix, outroot=fp%outroot)
        ! save it to file:
        call save_fisher_to_file( p, fp, num_param, fisher_matrix, filename=fp%outroot//'fisher_matrix_RD.dat' )
        ! save the parameter names to file:
        call save_paramnames_to_file( p, fp, num_param, filename=fp%outroot//'fisher_matrix_RD.paramnames' )
    end if

    if ( fp%cosmicfish_want_Mpk ) write(*,*) 'Not yet implemented'

    if (fp%cosmicfish_want_derived ) then
        ! get number of derived parameters:
        call dimension_derived_parameters( p, fp, num_derived )
        ! allocate the derived parameters matrix:
        allocate( derived_matrix(num_param, num_derived) )
        ! compute the Fisher matrix:
        call fisher_derived( p, fp, num_param, num_derived, derived_matrix, outroot=fp%outroot )
        ! save it to file:
        call save_derived_fisher_to_file( p, fp, num_param, num_derived, derived_matrix, filename=fp%outroot//'fisher_matrix_derived.dat' )
        ! save the parameter names to file:
        call save_derived_paramnames_to_file( p, fp, num_param, num_derived, filename=fp%outroot//'fisher_matrix_derived.paramnames' )
    end if

end program fisher_matrix_calculator