Skip to content

Commit

Permalink
Update I0 and I1 to use better approximations from Pavel Holoborodko
Browse files Browse the repository at this point in the history
  • Loading branch information
jzmaddock committed Jan 18, 2017
1 parent 78320ad commit cc8e300
Show file tree
Hide file tree
Showing 26 changed files with 1,641 additions and 60 deletions.
75 changes: 75 additions & 0 deletions doc/equations/bessel17.mml
@@ -0,0 +1,75 @@
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/TR/MathML2/dtd/xhtml-math11-f.dtd" [<!ENTITY mathml 'http://www.w3.org/1998/Math/MathML'>]>
<html xmlns="http://www.w3.org/1999/xhtml">
<head><title>bessel17</title>
<!-- MathML created with MathCast Equation Editor version 0.92 -->
</head>
<body>
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block">
<mrow>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mfenced>
<mrow>
<mi>x</mi>
</mrow>
</mfenced>
<mo>&#x2248;</mo>
<mfrac>
<mi>x</mi>
<mn>2</mn>
</mfrac>
<mfenced>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>[</mo>
<mfrac>
<mi>x</mi>
<mn>2</mn>
</mfrac>
<msup>
<mo>]</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>[</mo>
<mfrac>
<mi>x</mi>
<mn>2</mn>
</mfrac>
<msup>
<mo>]</mo>
<mn>4</mn>
</msup>
<mi>P</mi>
<mfenced>
<mrow>
<mo>[</mo>
<mfrac>
<mi>x</mi>
<mn>2</mn>
</mfrac>
<msup>
<mo>]</mo>
<mn>2</mn>
</msup>
</mrow>
</mfenced>
</mrow>
</mfenced>
<mspace width="1.0em"/>
<mo>;</mo>
<mspace width="1.0em"/>
<mi>x</mi>
<mo>&lt;</mo>
<mn>7.75</mn>
</mrow>
</math></body>
</html>
2 changes: 2 additions & 0 deletions doc/equations/bessel17.svg
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
50 changes: 50 additions & 0 deletions doc/equations/bessel18.mml
@@ -0,0 +1,50 @@
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/TR/MathML2/dtd/xhtml-math11-f.dtd" [<!ENTITY mathml 'http://www.w3.org/1998/Math/MathML'>]>
<html xmlns="http://www.w3.org/1999/xhtml">
<head><title>bessel18</title>
<!-- MathML created with MathCast Equation Editor version 0.92 -->
</head>
<body>
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block">
<mrow>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mfenced>
<mrow>
<mi>x</mi>
</mrow>
</mfenced>
<mo>&#x2248;</mo>
<mfrac>
<mrow>
<msup>
<mi>e</mi>
<mi>x</mi>
</msup>
</mrow>
<mrow>
<msqrt>
<mi>x</mi>
</msqrt>
</mrow>
</mfrac>
<mi>P</mi>
<mfenced>
<mrow>
<mfrac>
<mn>1</mn>
<mi>x</mi>
</mfrac>
</mrow>
</mfenced>
<mspace width="1.0em"/>
<mo>;</mo>
<mspace width="1.0em"/>
<mi>x</mi>
<mo>&gt;</mo>
<mn>7.75</mn>
</mrow>
</math></body>
</html>
2 changes: 2 additions & 0 deletions doc/equations/bessel18.svg
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
50 changes: 50 additions & 0 deletions doc/equations/bessel20.mml
@@ -0,0 +1,50 @@
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/TR/MathML2/dtd/xhtml-math11-f.dtd" [<!ENTITY mathml 'http://www.w3.org/1998/Math/MathML'>]>
<html xmlns="http://www.w3.org/1999/xhtml">
<head><title>bessel20</title>
<!-- MathML created with MathCast Equation Editor version 0.92 -->
</head>
<body>
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block">
<mrow>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mfenced>
<mrow>
<mi>x</mi>
</mrow>
</mfenced>
<mo>&#x2248;</mo>
<mfrac>
<mrow>
<msup>
<mi>e</mi>
<mi>x</mi>
</msup>
</mrow>
<mrow>
<msqrt>
<mi>x</mi>
</msqrt>
</mrow>
</mfrac>
<mi>P</mi>
<mfenced>
<mrow>
<mfrac>
<mn>1</mn>
<mi>x</mi>
</mfrac>
</mrow>
</mfenced>
<mspace width="1.0em"/>
<mo>;</mo>
<mspace width="1.0em"/>
<mi>x</mi>
<mo>&gt;</mo>
<mn>7.75</mn>
</mrow>
</math></body>
</html>

1 comment on commit cc8e300

@smr99
Copy link

@smr99 smr99 commented on cc8e300 Sep 26, 2017

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This commit includes a change to tools/config.hpp that -- on the face of it -- appears wrong. What is the motivation for commenting out this line:

#define BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS

??

Please sign in to comment.