public void generate_test()
{
double kappa = 35;
var vm = new VonMisesDistribution(0, kappa);
var target = new VonMisesFisherDistribution(new double[] { -1, 0, 0 }, kappa);
double[][] actual = target.Generate(100);
var sphere = new UniformSphereDistribution(3);
for (int i = 0; i < actual.Length; i++)
{
double p = sphere.ProbabilityDensityFunction(actual[i]);
Assert.AreEqual(0.025330295910584444, p, 1e-8);
}
}
public void ConstructorTest1()
{
// If p=2 the distribution reduces to the von Mises distribution on the circle.
double kappa = 4.2;
var vm = new VonMisesDistribution(0, kappa);
var target = new VonMisesFisherDistribution(new double[] { -1, 0 }, kappa);
double s = Math.Sqrt(2) / 2;
double[] mean = target.Mean;
double a000 = target.ProbabilityDensityFunction(new double[] { +1, +0 });
double a045 = target.ProbabilityDensityFunction(new double[] { +s, +s });
double a090 = target.ProbabilityDensityFunction(new double[] { +0, +1 });
double a135 = target.ProbabilityDensityFunction(new double[] { -s, +s });
double a180 = target.ProbabilityDensityFunction(new double[] { -1, +0 });
double a225 = target.ProbabilityDensityFunction(new double[] { -s, -s });
double a270 = target.ProbabilityDensityFunction(new double[] { +0, -1 });
double a315 = target.ProbabilityDensityFunction(new double[] { +s, -s });
double a360 = target.ProbabilityDensityFunction(new double[] { +1, +0 });
double offset = -Math.PI;
double e000 = vm.ProbabilityDensityFunction(offset + 0);
double e045 = vm.ProbabilityDensityFunction(offset + Math.PI / 4);
double e090 = vm.ProbabilityDensityFunction(offset + Math.PI / 2);
double e135 = vm.ProbabilityDensityFunction(offset + Math.PI * (3 / 4.0));
double e180 = vm.ProbabilityDensityFunction(offset + Math.PI);
double e225 = vm.ProbabilityDensityFunction(offset + Math.PI * (5 / 4.0));
double e270 = vm.ProbabilityDensityFunction(offset + Math.PI * (3 / 2.0));
double e315 = vm.ProbabilityDensityFunction(offset + Math.PI * (7 / 4.0));
double e360 = vm.ProbabilityDensityFunction(offset + Math.PI * 2);
Assert.AreEqual(e000, a000, 1e-8);
Assert.AreEqual(e045, a045, 1e-8);
Assert.AreEqual(e090, a090, 1e-8);
Assert.AreEqual(e135, a135, 1e-8);
Assert.AreEqual(e180, a180, 1e-8);
Assert.AreEqual(e225, a225, 1e-8);
Assert.AreEqual(e270, a270, 1e-8);
Assert.AreEqual(e315, a315, 1e-8);
Assert.AreEqual(e360, a360, 1e-8);
}
