05p2_bilinear_parametrizations_2D.f90

Go to the documentation of this file.

!----------------------------------------------------------------------------------------
!
! This file is part of EFTCAMB.
!
! Copyright (C) 2013-2016 by the EFTCAMB authors
!
! The EFTCAMB 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 eftcamb/LICENSE at
! the top level of the EFTCAMB distribution.
!
!----------------------------------------------------------------------------------------



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


module eftcamb_bilinear_parametrizations_2d

    use precision
    use amlutils
    use eft_def
    use eftcamb_cache
    use eftcamb_abstract_parametrizations_2d

    implicit none

    private

    public bilinear_parametrization_2d

    ! ---------------------------------------------------------------------------------------------
    type, extends ( parametrized_function_2d ) :: bilinear_parametrization_2d

        real(dl) :: linear_value_1
        real(dl) :: linear_value_2

    contains

        ! utility functions:
        procedure :: set_param_number      => bilinearparametrized2dsetparamnumber            
        procedure :: init_parameters       => bilinearparametrized2dinitparams                
        procedure :: parameter_value       => bilinearparametrized2dparametervalues           
        procedure :: feedback              => bilinearparametrized2dfeedback                  

        ! evaluation procedures:
        procedure :: value                 => bilinearparametrized2dvalue                   
        procedure :: first_derivative_x    => bilinearparametrized2dfirstderivativex        
        procedure :: first_derivative_y    => bilinearparametrized2dfirstderivativey        
        procedure :: second_derivative_x   => bilinearparametrized2dsecondderivativex       
        procedure :: second_derivative_y   => bilinearparametrized2dsecondderivativey       
        procedure :: second_derivative_xy  => bilinearparametrized2dsecondderivativexy      

    end type bilinear_parametrization_2d

contains

    ! ---------------------------------------------------------------------------------------------
    ! Implementation of the bilinear parametrization for 2D functions.
    ! ---------------------------------------------------------------------------------------------

    ! ---------------------------------------------------------------------------------------------
    subroutine bilinearparametrized2dsetparamnumber( self )

        implicit none

        class(bilinear_parametrization_2d) :: self

        ! initialize the number of parameters:
        self%parameter_number = 2

    end subroutine bilinearparametrized2dsetparamnumber

    ! ---------------------------------------------------------------------------------------------
    subroutine bilinearparametrized2dinitparams( self, array )

        implicit none

        class(bilinear_parametrization_2d)                        :: self
        real(dl), dimension(self%parameter_number), intent(in)  :: array

        self%linear_value_1 = array(1)
        self%linear_value_2 = array(2)

    end subroutine bilinearparametrized2dinitparams

    ! ---------------------------------------------------------------------------------------------
    subroutine bilinearparametrized2dparametervalues( self, i, value )

        implicit none

        class(bilinear_parametrization_2d) :: self
        integer     , intent(in)         :: i
        real(dl)    , intent(out)        :: value

        select case (i)
            case(1)
                value = self%linear_value_1
            case(2)
                value = self%linear_value_2
            case default
                write(*,*) 'Illegal index for parameter_names.'
                write(*,*) 'Maximum value is:', self%parameter_number
                call mpistop('EFTCAMB error')
        end select

    end subroutine bilinearparametrized2dparametervalues

    ! ---------------------------------------------------------------------------------------------
    subroutine bilinearparametrized2dfeedback( self, print_params )

        implicit none

        class(bilinear_parametrization_2d)  :: self
        logical, optional                   :: print_params

        integer                             :: i
        real(dl)                            :: param_value
        character(len=EFT_names_max_length) :: param_name
        logical                             :: print_params_temp

        if ( present(print_params) ) then
            print_params_temp = print_params
        else
            print_params_temp = .true.
        end if

        write(*,*)     'Bilinear function: ', self%name
        if ( print_params_temp ) then
            do i=1, self%parameter_number
                call self%parameter_names( i, param_name  )
                call self%parameter_value( i, param_value )
                write(*,'(a23,a,F12.6)') param_name, '=', param_value
            end do
        end if

    end subroutine bilinearparametrized2dfeedback

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

    function bilinearparametrized2dvalue( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dvalue

        bilinearparametrized2dvalue = self%linear_value_1*x + self%linear_value_2*y

    end function bilinearparametrized2dvalue

    ! ---------------------------------------------------------------------------------------------
    ! first variable, of the bilinear function.

    function bilinearparametrized2dfirstderivativex( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dfirstderivativex

        bilinearparametrized2dfirstderivativex = self%linear_value_1

    end function bilinearparametrized2dfirstderivativex

    ! ---------------------------------------------------------------------------------------------
    ! second variable, of the bilinear function.

    function bilinearparametrized2dfirstderivativey( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dfirstderivativey

        bilinearparametrized2dfirstderivativey = self%linear_value_2

    end function bilinearparametrized2dfirstderivativey

    ! ---------------------------------------------------------------------------------------------
    function bilinearparametrized2dsecondderivativex( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dsecondderivativex

        bilinearparametrized2dsecondderivativex = 0._dl

    end function bilinearparametrized2dsecondderivativex

    ! ---------------------------------------------------------------------------------------------
    function bilinearparametrized2dsecondderivativey( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dsecondderivativey

        bilinearparametrized2dsecondderivativey = 0._dl

    end function bilinearparametrized2dsecondderivativey

    ! ---------------------------------------------------------------------------------------------
    function bilinearparametrized2dsecondderivativexy( self, x, y, eft_cache )

        implicit none

        class(bilinear_parametrization_2d)                 :: self
        real(dl), intent(in)                               :: x
        real(dl), intent(in)                               :: y
        type(eftcamb_timestep_cache), intent(in), optional :: eft_cache
        real(dl) :: bilinearparametrized2dsecondderivativexy

        bilinearparametrized2dsecondderivativexy = 0._dl

    end function bilinearparametrized2dsecondderivativexy

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

end module eftcamb_bilinear_parametrizations_2d

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