math: add more tests for special cases of Bessel functions Y0, Y1, Yn

alexd765 · griesemer · commit ffb3b3698c69 · 2017-02-22T17:52:15.000Z
Test finite negative x with Y0(-1), Y1(-1), Yn(2,-1), Yn(-3,-1). Also test the special case Yn(0,0). Fixes golang#19130. Change-Id: I95f05a72e1c455ed8ddf202c56f4266f03f370fd Reviewed-on: https://go-review.googlesource.com/37310 Reviewed-by: Robert Griesemer <gri@golang.org>
diff --git a/src/math/all_test.go b/src/math/all_test.go
@@ -1742,30 +1742,35 @@ var vfy0SC = []float64{
 	0,
 	Inf(1),
 	NaN(),
+	-1,
 }
 var y0SC = []float64{
 	NaN(),
 	Inf(-1),
 	0,
 	NaN(),
+	NaN(),
 }
 var y1SC = []float64{
 	NaN(),
 	Inf(-1),
 	0,
 	NaN(),
+	NaN(),
 }
 var y2SC = []float64{
 	NaN(),
 	Inf(-1),
 	0,
 	NaN(),
+	NaN(),
 }
 var yM3SC = []float64{
 	NaN(),
 	Inf(1),
 	0,
 	NaN(),
+	NaN(),
 }
 
 // arguments and expected results for boundary cases
@@ -2716,6 +2721,9 @@ func TestYn(t *testing.T) {
 			t.Errorf("Yn(-3, %g) = %g, want %g", vfy0SC[i], f, yM3SC[i])
 		}
 	}
+	if f := Yn(0, 0); !alike(Inf(-1), f) {
+		t.Errorf("Yn(0, 0) = %g, want %g", f, Inf(-1))
+	}
 }
 
 // Check that math functions of high angle values

Original file line number	Diff line number	Diff line change
`@@ -1742,30 +1742,35 @@ var vfy0SC = []float64{`
`1742`	`1742`	`0,`
`1743`	`1743`	`Inf(1),`
`1744`	`1744`	`NaN(),`
	`1745`	`+ -1,`
`1745`	`1746`	`}`
`1746`	`1747`	`var y0SC = []float64{`
`1747`	`1748`	`NaN(),`
`1748`	`1749`	`Inf(-1),`
`1749`	`1750`	`0,`
`1750`	`1751`	`NaN(),`
	`1752`	`+ NaN(),`
`1751`	`1753`	`}`
`1752`	`1754`	`var y1SC = []float64{`
`1753`	`1755`	`NaN(),`
`1754`	`1756`	`Inf(-1),`
`1755`	`1757`	`0,`
`1756`	`1758`	`NaN(),`
	`1759`	`+ NaN(),`
`1757`	`1760`	`}`
`1758`	`1761`	`var y2SC = []float64{`
`1759`	`1762`	`NaN(),`
`1760`	`1763`	`Inf(-1),`
`1761`	`1764`	`0,`
`1762`	`1765`	`NaN(),`
	`1766`	`+ NaN(),`
`1763`	`1767`	`}`
`1764`	`1768`	`var yM3SC = []float64{`
`1765`	`1769`	`NaN(),`
`1766`	`1770`	`Inf(1),`
`1767`	`1771`	`0,`
`1768`	`1772`	`NaN(),`
	`1773`	`+ NaN(),`
`1769`	`1774`	`}`
`1770`	`1775`
`1771`	`1776`	`// arguments and expected results for boundary cases`
`@@ -2716,6 +2721,9 @@ func TestYn(t *testing.T) {`
`2716`	`2721`	`t.Errorf("Yn(-3, %g) = %g, want %g", vfy0SC[i], f, yM3SC[i])`
`2717`	`2722`	`}`
`2718`	`2723`	`}`
	`2724`	`+ if f := Yn(0, 0); !alike(Inf(-1), f) {`
	`2725`	`+ t.Errorf("Yn(0, 0) = %g, want %g", f, Inf(-1))`
	`2726`	`+ }`
`2719`	`2727`	`}`
`2720`	`2728`
`2721`	`2729`	`// Check that math functions of high angle values`