From e80b2ec787b1e4e1a6c894f2c274e6814d46fd81 Mon Sep 17 00:00:00 2001
From: wdaniau <wdaniau@b657c933-2333-4658-acf2-d3c7c2708721>
Date: Wed, 2 Dec 2009 16:54:21 +0000
Subject: [PATCH] Increase factor default value from 40 to 150 in dbesin.f90
 and dbesjn.f90 to increase accuracy of computed values compared to dbesri and
 dbesrj

git-svn-id: https://lxsd.femto-st.fr/svn/fvn@60 b657c933-2333-4658-acf2-d3c7c2708721
---
 fvn_fnlib/dbesin.f90 | 4 +++-
 fvn_fnlib/dbesjn.f90 | 4 +++-
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/fvn_fnlib/dbesin.f90 b/fvn_fnlib/dbesin.f90
index 9926b4e..ea59887 100644
--- a/fvn_fnlib/dbesin.f90
+++ b/fvn_fnlib/dbesin.f90
@@ -12,6 +12,8 @@ function dbesin(n,x,factor,big)
     ! big : a real that determine the threshold for taking anti overflow counter measure
     ! default value is 1e10
     !
+    ! 2009 : increasing factor defult value to 150
+    !
     real(dp_kind) :: dbesin
     integer :: n
     real(dp_kind) :: x
@@ -25,7 +27,7 @@ function dbesin(n,x,factor,big)
     real(dp_kind), external :: dbesi0,dbesi1
 
     ! Initialization of optional parameters
-    tfactor=40
+    tfactor=150
     if(present(factor)) tfactor=factor
     tbig=1e10
     if(present(big)) tbig=big
diff --git a/fvn_fnlib/dbesjn.f90 b/fvn_fnlib/dbesjn.f90
index e092866..afc184c 100644
--- a/fvn_fnlib/dbesjn.f90
+++ b/fvn_fnlib/dbesjn.f90
@@ -12,6 +12,8 @@ function dbesjn(n,x,factor,big)
     ! big : a real that determine the threshold for taking anti overflow counter measure
     ! default value is 1e10
     !
+    ! 2009 : increasing factor defult value to 150
+    !
     real(dp_kind) :: dbesjn
     integer :: n
     real(dp_kind) :: x
@@ -26,7 +28,7 @@ function dbesjn(n,x,factor,big)
     logical :: iseven
 
     ! Initialization of optional parameters
-    tfactor=40
+    tfactor=150
     if(present(factor)) tfactor=factor
     tbig=1d10
     if(present(big)) tbig=big
-- 
2.16.4