Commit faf18836 authored by Jean-Eric Campagne's avatar Jean-Eric Campagne
Browse files

bug in the beam formula identified

parent 5a826df2
...@@ -395,12 +395,13 @@ r_8 Beam::Value(r_8 const xp[], r_8 const* Par) { ...@@ -395,12 +395,13 @@ r_8 Beam::Value(r_8 const xp[], r_8 const* Par) {
} else { } else {
#if TYPE_BEAM == 1 #if TYPE_BEAM == 1
lob = exp(-0.220544*x*x); //gaussian equivalent of J1 formula normalise a 1 //BUG JEC 20/9/19 the magic coeff was determined on spherical_bessel fonction and not the Bessel function lob = exp(-0.220544*x*x); //gaussian equivalent of J1 formula normalise a 1
lob = exp(-0.280092*x*x); //gaussian equivalent of J1 formula normalise a 1
#else #else
if(x==0.) { if(x==0.) {
lob = 1.; lob = 1.;
} else { } else {
lob = 3. * boost::math::cyl_bessel_j(1, x)/x; lob = 2. * boost::math::bessel_j(1, x)/x;
lob *= lob; lob *= lob;
} }
#endif #endif
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment