public void Theta2DoubleTest()
{
var data = new[]
{
-0.712045573932278, 0.0115779438837883, 0.496601444431198,
-1.02921651633888, 0.0251305156397967, 0.409986239688389,
-0.270601468194827, 0.212623597863526, 1.35096326526392,
2.06182432952114, 0.018128943757431, -0.3458068199639,
-0.0161641166929001, 0.0974889810203627, 1.128018761735100,
};
for (int i = 0; i < data.Length; i += 3)
{
var z = data[i];
var q = data[i + 1];
var th = Theta.θ2(z, q);
var ex = data[i + 2];
var err = th - ex;
Assert.IsTrue(Complex.Abs(err) < 1e-14);
var tf = Theta.θ2DoubleForNome(q);
th = tf(z);
err = th - ex;
Assert.IsTrue(Complex.Abs(err) < 1e-14);
}
}
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);
}
