lxsd / fvn

Download as
Browse Code ยป

Commit fdbb9bb0a5c5ee1811ed27aace4452bd92e8529e

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

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

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

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