public void Theta3DoubleTest() { var data = new[] { 1, .8, 0.04246457514070379, .5, .8, 1.223823417425446, .5, .5, 1.484396862425167, }; for (int i = 0; i < data.Length; i += 3) { var z = data[i]; var q = data[i + 1]; var ex = data[i + 2]; var th = Theta.θ3(z, q); var err = Complex.Abs(th - ex) / Complex.Abs(th); Assert.IsTrue(Complex.Abs(err) < 1e-15); var t3 = Theta.θ3DoubleForNome(q); th = t3(z); err = Complex.Abs(th - ex) / Complex.Abs(th); Assert.IsTrue(Complex.Abs(err) < 1e-15); } }
public void ThetasAtZeroTest() { Complex th; double ex = 0.0; th = Theta.θ1(Complex.Zero, (Complex)(.5)); Assert.AreEqual(0, th.Imaginary); Assert.AreEqual(ex, th.Real); th = Theta.θ1(Complex.Zero, .5); Assert.AreEqual(0, th.Imaginary); Assert.AreEqual(ex, th.Real); th = Theta.θ1(0, .5); Assert.AreEqual(0, th.Imaginary); Assert.AreEqual(ex, th.Real); th = Theta.θ1ComplexForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.AreEqual(ex, th.Real); th = Theta.θ1DoubleForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.AreEqual(ex, th.Real); ex = 2.12893125051303; th = Theta.θ2(Complex.Zero, (Complex)(.5)); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ2(Complex.Zero, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ2(0, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ2ComplexForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ2DoubleForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); ex = 2.12893682721188; th = Theta.θ3(Complex.Zero, (Complex)(.5)); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ3(Complex.Zero, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ3(0, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ3ComplexForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ3DoubleForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); ex = 0.121124208002581; th = Theta.θ4(Complex.Zero, (Complex)(.5)); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ4(Complex.Zero, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ4(0, .5); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ4ComplexForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); th = Theta.θ4DoubleForNome(.5)(0); Assert.AreEqual(0, th.Imaginary); Assert.IsTrue(Math.Abs(ex - th.Real) < 1e-14); }