Homec4science
Diffusion COSOlver (improved_sugama) Browse

cosolver/src/test_mod.f90
cd1286e2d410improved_sugama

Last Change
Raw File

test_mod.f90
View Options

MODULE TEST
! TEST MODULE FOR COSOLVER
USE PREC_CONST
USE FMZM
USE COEFF
USE MODEL
USE SPEED_FUNCTIONS
USE COEFF
USE COULOMB
USE BASIS_TRANSFORMATION
USE AUXVAL
CONTAINS
SUBROUTINE maketest
! test routine
IMPLICIT NONE
TYPE(FM),SAVE :: X1,X2,X0
TYPE(FM) :: X3
INTEGER :: OUT = 18
CHARACTER(200) :: ST1,ST2
CHARACTER*100 :: T4STR_
CHARACTER*100 :: T5STR_
CALL FM_ENTER_USER_ROUTINE
CALL INIT_PREC_CONST
! set model to FM type for sampling tests
CALL set_model_fm
! write to file
! SET fm precision to give at least 60 significant digits
CALL FM_SET(50)
!
!! CALL GET_NNT4(15,15,15)
CALL GET_NNT5
CALL memory
!! CALL evaluate_auxval
!-----------------------------------------------------
!
! TEST COEFF SPEED FUNCTIONS
! !-----------------------------------------------------
! !
! ! TEST SPEED FUNCTIONS AND COEFF
! WRITE(OUT,*) 'Perform Sampling 1 ...'
! WRITE(OUT,*)' Computing Test-Particle speed function nusTabpjgd ...'
! WRITE(*,*) 'Perform Sampling 1 ...'
! WRITE(*,*)' Computing Test-Particle speed function nusTabpjgd ...'
! WRITE(*,*) 0,0,0,0
! X1 = nusTabpjgd(0,0,0,0)
! CALL print_fm(0,X1)
! WRITE(*,*) 1,0,1,0
! X1 = nusTabpjgd(1,0,1,0)
! CALL print_fm(0,X1)
! WRITE(*,*) 0,1,0,0
! X1 = nusTabpjgd(0,1,0,0)
! CALL print_fm(0,X1)
! WRITE(*,*) 2,2,2,0
! X1 = nusTabpjgd(2,2,2,0)
! CALL print_fm(0,X1)
! WRITE(*,*) 5,3,5,2
! X1 = nusTabpjgd(5,3,5,2)
! CALL print_fm(0,X1)
! WRITE(OUT,*)' Computing Back-Reaction speed function nusFabpjgd ...'
! WRITE(*,*)' Computing Back-Reaction speed function nusFabpjgd ...'
! WRITE(*,*) 0,0,0,0
! X2= nusFabpjgd(0,0,0,0)
! CALL print_fm(0,X2)
! WRITE(*,*) 1,0,1,0
! X2= nusFabpjgd(1,0,1,0)
! CALL print_fm(0,X2)
! WRITE(*,*) 0,1,0,0
! X2 = nusFabpjgd(0,1,0,0)
! CALL print_fm(0,X2)
! WRITE(*,*) 2,2,2,0
! X2 = nusFabpjgd(2,2,2,0)
! CALL print_fm(0,X2)
! WRITE(*,*) 5,3,5,2
! X2 = nusFabpjgd(5,3,5,2)
! CALL print_fm(0,X2)
! ! ! --------------------------------------
! CALL GET_NNT4(40,40,0)
! WRITE(OUT,*) 'Perform Sampling 2 ...'
! WRITE(*,*) 'Perform Sampling 2 ...'
! WRITE(OUT,*)'Compute Test-Particle coefficient A_ab^Tlkpjm ...'
! WRITE(*,*)'Compute Test-Particle coefficient A_ab^Tlkpjm ...'
! WRITE(OUT,*) 1,1,1,1,0
! X1= AabTLKPJM(1,1,1,1,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 2,2,1,1,0
! X1= AabTLKPJM(2,2,1,1,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,3,2,2,0
! X1 = AabTLKPJM(10,3,2,2,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,6,6,4,0
! X1= AabTLKPJM(10,6,6,4,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*)'Compute Back-Reaction coefficient B_ab^Flkpjm...'
! WRITE(*,*)'Compute Back-Reaction coefficient B_ab^Flkpjm...'
! WRITE(OUT,*) 1,1,1,1,0
! X1= BabFLKPJM(1,1,1,1,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 2,2,1,1,0
! X1= BabFLKPJM(2,2,1,1,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,3,2,2,0
! X1= BabFLKPJM(10,3,2,2,0)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,10,6,4,0
! X1= BabFLKPJM(10,10,6,4,0)
! CALL print_fm(OUT,X1)
!-----------------------------------------------------
!
! TEST BASIS TRANSFORMATION COEFF.
! !
! CALL GET_NNT4(10,10,0)
! T4STR_ = GET_T4(0,0,3,0)
! WRITE(*,*) T4STR_
! X2 = READ_T4_FROM_CSV(9,9,0,0,0)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! X2 = READ_T4_FROM_CSV(5,0,5,0,1)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! X2 = READ_T4_FROM_CSV(0,2,4,0,-1)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! X2 = READ_T4_FROM_CSV(0,2,4,0,-1)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! X2 = READ_T4_FROM_CSV(55,60,55,60,0)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
!
!
! T5STR_ = GET_T5(3,0,1,0,0)
! write(*,*) T5STR_
! T5STR_ = GET_T5(4,0,1,0,0)
! write(*,*) T5STR_
! T5STR_ = GET_T5(23,11,0,1,22)
! write(*,*) T5STR_
X0 = GET_T5inv(2,3,2,4,2)
WRITE(*,*) TO_DP(X0)
X0 = GET_T5inv(2,3,2,4,1)
WRITE(*,*) TO_DP(X0)
X0 = GET_T5inv(1,2,1,2,0)
WRITE(*,*) TO_DP(X0)
!----------------------------------------------------
!
! TEST LORENTZ MODEL SPEED VELOCITY INTEGRAL
! X1 = Injlkmp(10,3,3,2,0,1)
! CALL PRINT_FM(0,X1)
!----------------------------------------------------
!
!
! TEST IMPROVED SUGAMA
!
! Here we check that N_{ab}^{Sij} = M_ab^{S0i} N_{ab}^{S0j}/ M_{ab}^{00}
! DO i=0,9
! DO j=0,9
! X1 = NabSlk(i,j)
! X2 = MabSlk(0,i)
! X3 = NabSlk(0,j)
! CALL print_fm(11,X1)
! CALL print_fm(11,X2*X3/MabLlk(0,0))
! ! X1 = MabLlk(i,j)
! ! X2 = MabSlk(i,j)
! ! call print_fm(0,X1)
! ! call print_fm(0,X2)
! ! ENDDO
! ! X0 = NabSlk(i,j)
! ! CALL print_fm(0,X0)
! ! !
! ! X1 = NabLlk(i,0)
! ! X2 = NabLlk(0,j)
! ! X3 = X1*X2/NabLlk(0,0)
! ! CALL print_fm(0,X3)
! ENDDO
! ENDDO
! X1 = MabLlk(0,0)
! CALL print_fm(0,X1)
! !-----------------------------------------------------
! ! TEST SUGAMA SPEED FUNCTIONS
! X0 = curlyVab(0,0)
! CALL print_fm(0,X0)
! X0 = ALLX2L(0,0,0,1)
! CALL print_fm(0,X0)
! X0 = ALLX2L(0,0,1,1)
! CALL print_fm(0,X0)
! X0 = curlyIbals(0,0)
! CALL print_fm(0,X0)
! X0 = curlyIbals(0,1)
! CALL print_fm(0,X0)
! X0 = curlyIabls(0,0)
! CALL print_fm(0,X0)
! X0 = curlyIabls(0,1)
! CALL print_fm(0,X0)
! X0 = Erfkab(0)
! CALL print_fm(0,X0)
! X0 = Erfkab(1)
! CALL print_fm(0,X0)
! X0 = nuabDk(0)
! call print_fm(0,X0)
! X0 = nuabDk(1)
! call print_fm(0,X0)
! X0 = nuabDk(2)
! call print_fm(0,X0)
! X0 = nuabDk(3)
! call print_fm(0,X0)
! X0 = nuabDk(4)
! call print_fm(0,X0)
! X0 = nuabSDk(0)
! call print_fm(0,X0)
! X0 = nuabSDk(1)
! call print_fm(0,X0)
! X0 = nuabSDk(2)
! call print_fm(0,X0)
! X0 = nuabparak(0)
! call print_fm(0,X0)
! X0 = nuabparak(1)
! call print_fm(0,X0)
! X0 = nuabparak(2)
! call print_fm(0,X0)
! X0 = nuabparak(3)
! call print_fm(0,X0)
! ! X0 = nusSabpjd(0,0,0)
! ! CALL PRINT_FM(0,X0)
! ! X0 = nusSabpjd(0,0,1)
! ! CALL PRINT_FM(0,X0)
! ! X0 = nusSabpjd(1,2,0)
! ! CALL PRINT_FM(0,X0)
! ! X0 = nusSabpjd(2,5,1)
! ! CALL PRINT_FM(0,X0)
! !-----------------------------------------------------
! ! TEST OVERFLOWS
! X0 = ALLX2L(0,10,10,0)
! X0 = X2prodLegendre(2,2)
! CALL PRINT_FM(0,X0)
! X0 = X2prodLegendre(4,2)
! CALL PRINT_FM(0,X0)
! X0 = X2prodLegendre(4,6)
! CALL PRINT_FM(0,X0)
! TEST FM arithmetic
CALL FM_EXIT_USER_ROUTINE
END SUBROUTINE maketest
!-----------------------------------------------------
! SUBROUTINE SAMPLING1
! ! Perform Sampling 1
! !
! WRITE(OUT,*) 'Perform Sampling 1 ...'
! WRITE(OUT,*)' Computing Test-Particle speed function nusTabpjgd ...'
! WRITE(*,*) 'Perform Sampling 1 ...'
! WRITE(*,*)' Computing Test-Particle speed function nusTabpjgd ...'
! WRITE(OUT,*) 1,1,1,1
! CALL nusTabpjgd(1,1,1,1,X1)ssdsad
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 2,3,1,1
! CALL nusTabpjgd(2,3,1,1,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,8,5,3
! CALL nusTabpjgd(10,8,5,3,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*)' Compuing Back-Reaction speed function nusFabpjgd ...'
! WRITE(*,*)' Computing Back-Reaction speed function nusFabpjgd ...'
! WRITE(OUT,*) 1,1,1,1
! CALL nusFabpjgd(1,1,1,1,X2)
! CALL print_fm(OUT,X2)
! WRITE(OUT,*) 3,2,1,1
! CALL nusFabpjgd(3,2,1,1,X2)
! CALL print_fm(OUT,X2)
! WRITE(OUT,*) 10,8,4,3
! CALL nusFabpjgd(10,8,4,3,X2)
! CALL print_fm(OUT,X2)
! END SUBROUTINE SAMPLING1
! ! ---------------------------------------------
! SUBROUTINE SAMPLING2
! ! Perform Sampling 2
! !
! WRITE(OUT,*) 'Perform Sampling 2 ...'
! WRITE(*,*) 'Perform Sampling 2 ...'
! WRITE(OUT,*)'Compute Test-Particle coefficient A_ab^Tlkpjm ...'
! WRITE(*,*)'Compute Test-Particle coefficient A_ab^Tlkpjm ...'
! WRITE(OUT,*) 1,1,1,1,0
! CALL AabTLKPJM(1,1,1,1,0,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,3,2,2,0
! CALL AabTLKPJM(10,3,2,2,0,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,6,6,4,0
! CALL AabTLKPJM(10,6,6,4,0,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*)'Compute Back-Reaction coefficient B_ab^Flkpjm...'
! WRITE(*,*)'Compute Back-Reaction coefficient B_ab^Flkpjm...'
! WRITE(OUT,*) 1,1,1,1,0
! CALL BabFLKPJM(1,1,1,1,0,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,3,2,2,0
! CALL BabFLKPJM(10,3,2,2,0,X1)
! CALL print_fm(OUT,X1)
! WRITE(OUT,*) 10,10,6,4,0
! CALL BabFLKPJM(10,10,6,4,0,X1)
! CALL print_fm(OUT,X1)
! END SUBROUTINE SAMPLING2
! ! ---------------------------------------------
! SUBROUTINE SAMPLING_T4
! ! This routines the routine READ_T4_FROM_CSV, and inquire a seeded list of T4 and T4^{-1} coeff
! !
! USE FMZM
! IMPLICIT NONE
! !
! CALL FM_ENTER_USER_ROUTINE
! CALL READ_T4_FROM_CSV(5,0,5,0,1,X2)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! CALL READ_T4_FROM_CSV(0,2,4,0,-1,X2)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! CALL READ_T4_FROM_CSV(55,60,55,60,0,X2)
! CALL FM_FORM('F90.60',X2,ST1)
! WRITE(*,*) ADJUSTL(TRIM(ST1))
! CALL FM_EXIT_USER_ROUTINE
! END SUBROUTINE SAMPLING_T4
! ---------------------------------------------
END MODULE TEST
c4science · Help