lxsd / fvn

Download as
Browse Code ยป

Commit 47b3dde4dcf34f679f4df299670b8a0baa508652

Authored by wdaniau
1 parent 6afbb16eea
Exists in master and in 3 other branches geevx, geevx_2020, multi

Ajout d'un programme test de la fonction ze1 

git-svn-id: https://lxsd.femto-st.fr/svn/fvn@57 b657c933-2333-4658-acf2-d3c7c2708721

Showing 2 changed files with 41 additions and 1 deletions Side-by-side Diff

... ... @@ -4,7 +4,7 @@
4 4 programs = test_fac$(exext) test_matinv$(exext) test_specfunc$(exext) \
5 5 test_det$(exext) test_matcon$(exext) test_matev$(exext) test_sparse$(exext) test_inter1d$(exext) \
6 6 test_inter2d$(exext) test_inter3d$(exext) test_akima$(exext) test_lsp$(exext) test_muller$(exext) \
7   -test_integ$(exext) test_bsyn$(exext) test_bsjn$(exext) test_bskn$(exext) test_bsin$(exext) test_operators$(exext)
  7 +test_integ$(exext) test_bsyn$(exext) test_bsjn$(exext) test_bskn$(exext) test_bsin$(exext) test_operators$(exext) test_ze1$(exext)
8 8  
9 9 prog:$(programs)
10 10  
fvn_test/test_ze1.f90
Diff comments   View file @ 47b3dde
  1 +program test_ze1
  2 +use fvn_fnlib
  3 +implicit none
  4 +
  5 +complex(8) :: z,resz
  6 +complex(8), dimension(6) :: abramowitz_x,abramowitz_y
  7 +integer(4) :: i
  8 +
  9 +abramowitz_x(1)=dcmplx(5.,5.)
  10 +abramowitz_y(1)=dcmplx(0.906058d0,0.070209d0)
  11 +abramowitz_x(2)=dcmplx(-5.,5.)
  12 +abramowitz_y(2)=dcmplx(1.079407d0,0.143879d0)
  13 +abramowitz_x(3)=dcmplx(11.,5.)
  14 +abramowitz_y(3)=dcmplx(0.932672d0,0.026361d0)
  15 +abramowitz_x(4)=dcmplx(-11.,5.)
  16 +abramowitz_y(4)=dcmplx(1.084526d0,0.049336d0)
  17 +abramowitz_x(5)=dcmplx(15.,18.)
  18 +abramowitz_y(5)=dcmplx(0.972359d0,0.029448d0)
  19 +abramowitz_x(6)=dcmplx(-15.,18.)
  20 +abramowitz_y(6)=dcmplx(1.026183d0,0.036552d0)
  21 +
  22 +! Value of z*exp(z)*e1(z) from Abramowitz & Stegun
  23 +! to cover the different cases of the algo
  24 +!
  25 +! Case (5,5) modulus <10 with positive real part
  26 +! Case (-5,5) modulus <10 with negative real part
  27 +! Case (11,5) modulus between 10 and 20 with positive real part
  28 +! Case (-11,5) modulus between 10 and 20 with negative real part
  29 +! Case (15,18) modulus > 20 with positive real part
  30 +! Case (-15,18) modulus > 20 with negative real part
  31 +
  32 +do i=1,6
  33 + z=abramowitz_x(i)
  34 + resz=z*exp(z)*e1(z)
  35 + write(*,*) "Tabulated : ",abramowitz_y(i)
  36 + write(*,*) "Calculated :",resz
  37 + write(*,*)
  38 +end do
  39 +end program