From f26a262db02dd7dcc101270a2a4c0ecb573d5cfa Mon Sep 17 00:00:00 2001
From: daniau <daniau@b657c933-2333-4658-acf2-d3c7c2708721>
Date: Mon, 18 Feb 2008 09:29:01 +0000
Subject: [PATCH] git-svn-id: https://lxsd.femto-st.fr/svn/fvn@37
 b657c933-2333-4658-acf2-d3c7c2708721

---
 fvn_test/Makefile          |  2 +-
 fvn_test/test_bsin.f90     | 33 +++++++++++++++++++++++++++++++++
 fvn_test/test_bsjn.f90     | 33 +++++++++++++++++++++++++++++++++
 fvn_test/test_bskn.f90     | 35 +++++++++++++++++++++++++++++++++++
 fvn_test/test_bsyn.f90     | 35 +++++++++++++++++++++++++++++++++++
 fvn_test/test_specfunc.f90 | 12 ++++++++++++
 6 files changed, 149 insertions(+), 1 deletion(-)
 create mode 100644 fvn_test/test_bsin.f90
 create mode 100644 fvn_test/test_bsjn.f90
 create mode 100644 fvn_test/test_bskn.f90
 create mode 100644 fvn_test/test_bsyn.f90

diff --git a/fvn_test/Makefile b/fvn_test/Makefile
index f3cb1f5..52612b8 100644
--- a/fvn_test/Makefile
+++ b/fvn_test/Makefile
@@ -4,7 +4,7 @@ include $(BTREE)/Make.inc
 programs = test_fac test_matinv test_specfunc  \
 test_det test_matcon test_matev test_sparse test_inter1d \
 test_inter2d test_inter3d test_akima test_lsp test_muller \
-test_integ
+test_integ test_bsyn test_bsjn test_bskn test_bsin
 
 prog:$(programs)
 
diff --git a/fvn_test/test_bsin.f90 b/fvn_test/test_bsin.f90
new file mode 100644
index 0000000..6a8ea70
--- /dev/null
+++ b/fvn_test/test_bsin.f90
@@ -0,0 +1,33 @@
+program test_bsin
+use fvn_fnlib
+implicit none
+integer :: n,npoints,i
+real(4) :: xmin,xmax,xstep,x,y
+real(4),dimension(200,6) :: bes
+!    interface
+!    real(4) function besjn(n,x,factor,big)
+!        integer :: n
+!        real(4) :: x
+!        integer,optional :: factor
+!        real(4),optional :: big
+!    end function besjn
+!end interface
+!bsyn
+!y=besjn(2,0.5)
+!stop
+open(2,file='bsin.dat')
+
+xmin=-5.
+xmax=5.
+npoints=200
+xstep=(xmax-xmin)/dble(npoints)
+do i=1,npoints
+      x=xmin+i*xstep
+      do n=0,5
+        bes(i,n)=bsin(n,x)
+      end do
+      write(2,'(7e22.14)') x,bes(i,0:5)
+end do
+close(2)
+
+end program
\ No newline at end of file
diff --git a/fvn_test/test_bsjn.f90 b/fvn_test/test_bsjn.f90
new file mode 100644
index 0000000..457fda9
--- /dev/null
+++ b/fvn_test/test_bsjn.f90
@@ -0,0 +1,33 @@
+program test_bsjn
+use fvn_fnlib
+implicit none
+integer :: n,npoints,i
+real(4) :: xmin,xmax,xstep,x,y
+real(4),dimension(200,6) :: bes
+!    interface
+!    real(4) function besjn(n,x,factor,big)
+!        integer :: n
+!        real(4) :: x
+!        integer,optional :: factor
+!        real(4),optional :: big
+!    end function besjn
+!end interface
+!bsyn
+!y=besjn(2,0.5)
+!stop
+open(2,file='bsjn.dat')
+
+xmin=-20.
+xmax=20.
+npoints=200
+xstep=(xmax-xmin)/dble(npoints)
+do i=1,npoints
+      x=xmin+i*xstep
+      do n=0,5
+        bes(i,n)=bsjn(n,x)
+      end do
+      write(2,'(7e22.14)') x,bes(i,0:5)
+end do
+close(2)
+
+end program
\ No newline at end of file
diff --git a/fvn_test/test_bskn.f90 b/fvn_test/test_bskn.f90
new file mode 100644
index 0000000..b412b36
--- /dev/null
+++ b/fvn_test/test_bskn.f90
@@ -0,0 +1,35 @@
+program test_bskn
+use fvn_fnlib
+implicit none
+integer :: n,npoints,i
+real(4) :: xmin,xmax,xstep,x,y
+real(4),dimension(200,6) :: bes
+!    interface
+!    real(4) function besjn(n,x,factor,big)
+!        integer :: n
+!        real(4) :: x
+!        integer,optional :: factor
+!        real(4),optional :: big
+!    end function besjn
+!end interface
+!bsyn
+!y=besjn(2,0.5)
+!stop
+open(2,file='bskn.dat')
+
+xmin=0.
+xmax=10.
+npoints=200
+xstep=(xmax-xmin)/dble(npoints)
+do i=1,npoints
+      x=xmin+i*xstep
+      do n=0,5
+        bes(i,n)=bskn(n,x)
+      end do
+      write(2,'(7e22.14)') x,bes(i,0:5)
+end do
+close(2)
+
+! with gnuplot -> set yrange [-1:1]
+
+end program
\ No newline at end of file
diff --git a/fvn_test/test_bsyn.f90 b/fvn_test/test_bsyn.f90
new file mode 100644
index 0000000..87f22fe
--- /dev/null
+++ b/fvn_test/test_bsyn.f90
@@ -0,0 +1,35 @@
+program test_bsyn
+use fvn_fnlib
+implicit none
+integer :: n,npoints,i
+real(4) :: xmin,xmax,xstep,x,y
+real(4),dimension(200,6) :: bes
+!    interface
+!    real(4) function besjn(n,x,factor,big)
+!        integer :: n
+!        real(4) :: x
+!        integer,optional :: factor
+!        real(4),optional :: big
+!    end function besjn
+!end interface
+!bsyn
+!y=besjn(2,0.5)
+!stop
+open(2,file='bsyn.dat')
+
+xmin=0.
+xmax=20.
+npoints=200
+xstep=(xmax-xmin)/dble(npoints)
+do i=1,npoints
+      x=xmin+i*xstep
+      do n=0,5
+        bes(i,n)=bsyn(n,x)
+      end do
+      write(2,'(7e22.14)') x,bes(i,0:5)
+end do
+close(2)
+
+! with gnuplot -> set yrange [-1:1]
+
+end program
\ No newline at end of file
diff --git a/fvn_test/test_specfunc.f90 b/fvn_test/test_specfunc.f90
index cb9fc10..c5e763a 100644
--- a/fvn_test/test_specfunc.f90
+++ b/fvn_test/test_specfunc.f90
@@ -164,6 +164,18 @@ write(2,'(2e22.14)') x,y
 end do
 close(2)
 
+! dawson
+xmin=-5.
+xmax=5.
+xstep=(xmax-xmin)/dble(npoints)
+open(2,file='daws.dat')
+do i=1,npoints
+x=xmin+i*xstep
+y=daws(x)
+write(2,'(2e22.14)') x,y
+end do
+close(2)
+
 
 
 end program test_specfunc
-- 
2.16.4