lxsd / fvn

Download as
Browse Code ยป

Commit e93f80fc630e5ff34c9bfd35f652cd47b1cf184e

Authored by daniau
1 parent 2919a9e2da
Exists in master and in 3 other branches geevx, geevx_2020, multi

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

Showing 4 changed files with 4 additions and 4 deletions Inline Diff

fvn_test/test_bsin.f90
Diff comments   View file @ e93f80f
program test_bsin 1 1 program test_bsin
use fvn_fnlib 2 2 use fvn_fnlib
implicit none 3 3 implicit none
integer :: n,npoints,i 4 4 integer :: n,npoints,i
real(4) :: xmin,xmax,xstep,x,y 5 5 real(4) :: xmin,xmax,xstep,x,y
real(4),dimension(200,6) :: bes 6 6 real(4),dimension(200,0:5) :: bes
! interface 7 7 ! interface
! real(4) function besjn(n,x,factor,big) 8 8 ! real(4) function besjn(n,x,factor,big)
! integer :: n 9 9 ! integer :: n
! real(4) :: x 10 10 ! real(4) :: x
! integer,optional :: factor 11 11 ! integer,optional :: factor
! real(4),optional :: big 12 12 ! real(4),optional :: big
! end function besjn 13 13 ! end function besjn
!end interface 14 14 !end interface
!bsyn 15 15 !bsyn
!y=besjn(2,0.5) 16 16 !y=besjn(2,0.5)
!stop 17 17 !stop
open(2,file='bsin.dat') 18 18 open(2,file='bsin.dat')
19 19
xmin=-5. 20 20 xmin=-5.
xmax=5. 21 21 xmax=5.
npoints=200 22 22 npoints=200
xstep=(xmax-xmin)/dble(npoints) 23 23 xstep=(xmax-xmin)/dble(npoints)
do i=1,npoints 24 24 do i=1,npoints
x=xmin+i*xstep 25 25 x=xmin+i*xstep
do n=0,5 26 26 do n=0,5
bes(i,n)=bsin(n,x) 27 27 bes(i,n)=bsin(n,x)
end do 28 28 end do
write(2,'(7e22.14)') x,bes(i,0:5) 29 29 write(2,'(7e22.14)') x,bes(i,0:5)
end do 30 30 end do
close(2) 31 31 close(2)
32 32
fvn_test/test_bsjn.f90
Diff comments   View file @ e93f80f
program test_bsjn 1 1 program test_bsjn
use fvn_fnlib 2 2 use fvn_fnlib
implicit none 3 3 implicit none
integer :: n,npoints,i 4 4 integer :: n,npoints,i
real(4) :: xmin,xmax,xstep,x,y 5 5 real(4) :: xmin,xmax,xstep,x,y
real(4),dimension(200,6) :: bes 6 6 real(4),dimension(200,0:5) :: bes
! interface 7 7 ! interface
! real(4) function besjn(n,x,factor,big) 8 8 ! real(4) function besjn(n,x,factor,big)
! integer :: n 9 9 ! integer :: n
! real(4) :: x 10 10 ! real(4) :: x
! integer,optional :: factor 11 11 ! integer,optional :: factor
! real(4),optional :: big 12 12 ! real(4),optional :: big
! end function besjn 13 13 ! end function besjn
!end interface 14 14 !end interface
!bsyn 15 15 !bsyn
!y=besjn(2,0.5) 16 16 !y=besjn(2,0.5)
!stop 17 17 !stop
open(2,file='bsjn.dat') 18 18 open(2,file='bsjn.dat')
19 19
xmin=-20. 20 20 xmin=-20.
xmax=20. 21 21 xmax=20.
npoints=200 22 22 npoints=200
xstep=(xmax-xmin)/dble(npoints) 23 23 xstep=(xmax-xmin)/dble(npoints)
do i=1,npoints 24 24 do i=1,npoints
x=xmin+i*xstep 25 25 x=xmin+i*xstep
do n=0,5 26 26 do n=0,5
bes(i,n)=bsjn(n,x) 27 27 bes(i,n)=bsjn(n,x)
end do 28 28 end do
write(2,'(7e22.14)') x,bes(i,0:5) 29 29 write(2,'(7e22.14)') x,bes(i,0:5)
end do 30 30 end do
close(2) 31 31 close(2)
32 32
fvn_test/test_bskn.f90
Diff comments   View file @ e93f80f
program test_bskn 1 1 program test_bskn
use fvn_fnlib 2 2 use fvn_fnlib
implicit none 3 3 implicit none
integer :: n,npoints,i 4 4 integer :: n,npoints,i
real(4) :: xmin,xmax,xstep,x,y 5 5 real(4) :: xmin,xmax,xstep,x,y
real(4),dimension(200,6) :: bes 6 6 real(4),dimension(200,0:5) :: bes
! interface 7 7 ! interface
! real(4) function besjn(n,x,factor,big) 8 8 ! real(4) function besjn(n,x,factor,big)
! integer :: n 9 9 ! integer :: n
! real(4) :: x 10 10 ! real(4) :: x
! integer,optional :: factor 11 11 ! integer,optional :: factor
! real(4),optional :: big 12 12 ! real(4),optional :: big
! end function besjn 13 13 ! end function besjn
!end interface 14 14 !end interface
!bsyn 15 15 !bsyn
!y=besjn(2,0.5) 16 16 !y=besjn(2,0.5)
!stop 17 17 !stop
open(2,file='bskn.dat') 18 18 open(2,file='bskn.dat')
19 19
xmin=0. 20 20 xmin=0.
xmax=10. 21 21 xmax=10.
npoints=200 22 22 npoints=200
xstep=(xmax-xmin)/dble(npoints) 23 23 xstep=(xmax-xmin)/dble(npoints)
do i=1,npoints 24 24 do i=1,npoints
x=xmin+i*xstep 25 25 x=xmin+i*xstep
do n=0,5 26 26 do n=0,5
bes(i,n)=bskn(n,x) 27 27 bes(i,n)=bskn(n,x)
end do 28 28 end do
write(2,'(7e22.14)') x,bes(i,0:5) 29 29 write(2,'(7e22.14)') x,bes(i,0:5)
end do 30 30 end do
close(2) 31 31 close(2)
32 32
! with gnuplot -> set yrange [-1:1] 33 33 ! with gnuplot -> set yrange [-1:1]
fvn_test/test_bsyn.f90
Diff comments   View file @ e93f80f
program test_bsyn 1 1 program test_bsyn
use fvn_fnlib 2 2 use fvn_fnlib
implicit none 3 3 implicit none
integer :: n,npoints,i 4 4 integer :: n,npoints,i
real(4) :: xmin,xmax,xstep,x,y 5 5 real(4) :: xmin,xmax,xstep,x,y
real(4),dimension(200,6) :: bes 6 6 real(4),dimension(200,0:5) :: bes
! interface 7 7 ! interface
! real(4) function besjn(n,x,factor,big) 8 8 ! real(4) function besjn(n,x,factor,big)
! integer :: n 9 9 ! integer :: n
! real(4) :: x 10 10 ! real(4) :: x
! integer,optional :: factor 11 11 ! integer,optional :: factor
! real(4),optional :: big 12 12 ! real(4),optional :: big
! end function besjn 13 13 ! end function besjn
!end interface 14 14 !end interface
!bsyn 15 15 !bsyn
!y=besjn(2,0.5) 16 16 !y=besjn(2,0.5)
!stop 17 17 !stop
open(2,file='bsyn.dat') 18 18 open(2,file='bsyn.dat')
19 19
xmin=0. 20 20 xmin=0.
xmax=20. 21 21 xmax=20.
npoints=200 22 22 npoints=200
xstep=(xmax-xmin)/dble(npoints) 23 23 xstep=(xmax-xmin)/dble(npoints)
do i=1,npoints 24 24 do i=1,npoints
x=xmin+i*xstep 25 25 x=xmin+i*xstep
do n=0,5 26 26 do n=0,5
bes(i,n)=bsyn(n,x) 27 27 bes(i,n)=bsyn(n,x)
end do 28 28 end do
write(2,'(7e22.14)') x,bes(i,0:5) 29 29 write(2,'(7e22.14)') x,bes(i,0:5)
end do 30 30 end do
close(2) 31 31 close(2)
32 32
! with gnuplot -> set yrange [-1:1] 33 33 ! with gnuplot -> set yrange [-1:1]