public void JaggedSetValues2()
{
Array jaggedGaussian = JaggedArray.ConvertToNew(
jagged, typeof(int), typeof(Gaussian),
delegate(object elt) { return(new Gaussian(1.0, 1.0)); });
Assert.Equal(
JaggedArray.GetLength(jaggedGaussian, typeof(Gaussian)),
JaggedArray.GetLength(jagged, typeof(int)));
JaggedArray.ConvertElements2(
jaggedGaussian, jagged, typeof(Gaussian),
delegate(object elt1, object elt2)
{
Gaussian g = (Gaussian)elt1;
double d = (int)elt2;
return(Gaussian.FromMeanAndPrecision(d + g.GetMean(), 2.0));
});
int i = 0;
foreach (Gaussian g in JaggedArray.ElementIterator(jaggedGaussian, typeof(Gaussian)))
{
Assert.Equal(g.GetMean(), 1.0 + arrExpected[i++]);
}
}
public void Array2DArrayEstimatorTest()
{
// Create the estimator
// element type is double[]
//Estimator<GaussianArray2DArray> est = null;
int length1 = ga2aArrayOfDist.GetLength(0);
int length2 = ga2aArrayOfDist.GetLength(1);
Array2DEstimator <GaussianArrayEstimator, GaussianArray2DArray, GaussianArray, double[]> est =
new Array2DEstimator <GaussianArrayEstimator, GaussianArray2DArray, GaussianArray, double[]>(
Utilities.Util.ArrayInit(length1, length2, (i, j) =>
new GaussianArrayEstimator(Utilities.Util.ArrayInit(ga2aArrayOfDist[i, j].Length, k => new GaussianEstimator()))
));
// Add some samples to the estimator
int nSamples = 5;
// Create a sample an mean variable of the right structure
double[, ][] mean = (double[, ][]) JaggedArray.ConvertToNew(
ga2aArrayOfDist, typeof(double), delegate(object elt) { return(0.0); });
double[, ][] sample = (double[, ][]) JaggedArray.ConvertToNew(
ga2aArrayOfDist, typeof(double), delegate(object elt) { return(0.0); });
// Create samples and add them to the estimator. Accumulate sum of samples at the same time
for (int nSamp = 0; nSamp < nSamples; nSamp++)
{
JaggedArray.ConvertElements2(
sample, ga2aArrayOfDist, delegate(object smp, object dst) { return(((Gaussian)dst).Sample()); });
est.Add(sample);
JaggedArray.ConvertElements2(
mean, sample, delegate(object mn, object smp) { return((double)mn + (double)smp); });
}
// Hand calculate the sample mean
JaggedArray.ConvertElements(
mean, delegate(object mn) { return(((double)mn) / ((double)nSamples)); });
// Let the estimator do the work
GaussianArray2DArray result = new GaussianArray2DArray(length1, length2, (i, j) => new GaussianArray(ga2aArrayOfDist[i, j].Length));
result = est.GetDistribution(result);
// The results should be identical to a very high precision
for (int i = 0; i < dim1; i++)
{
for (int j = 0; j < dim2; j++)
{
for (int k = 0; k < result[i, j].Count; k++)
{
Assert.True(System.Math.Abs(result[i, j][k].GetMean() - mean[i, j][k]) < 1e-9);
}
}
}
}
public void JaggedConvertToNewTypeGeneric()
{
Array jaggedGaussian = JaggedArray.ConvertToNew <int, Gaussian>(
jagged, delegate(int elt) { return(new Gaussian(0.5, 1.0)); });
Assert.Equal(
JaggedArray.GetLength(jaggedGaussian, typeof(Gaussian)),
JaggedArray.GetLength(jagged, typeof(int)));
foreach (Gaussian g in JaggedArray.ElementIterator <Gaussian>(jaggedGaussian))
{
Assert.Equal(0.5, g.GetMean());
}
}
public void JaggedVisitor()
{
Array jaggedDir = JaggedArray.ConvertToNew(
jagged, typeof(int), typeof(Dirichlet),
delegate(object elt) { return(new Dirichlet(2.0, 3.0, 4.0)); });
double[] newPseudoCount = { 3.0, 4.0, 5.0 };
JaggedArray.VisitElements(
jaggedDir, typeof(Dirichlet),
delegate(object elt) { ((Dirichlet)elt).PseudoCount.SetTo(newPseudoCount); });
foreach (Dirichlet d in JaggedArray.ElementIterator(jaggedDir, typeof(Dirichlet)))
{
Assert.Equal(3.0, d.PseudoCount[0]);
}
}
public void JaggedSetValues()
{
Array jaggedGaussian = JaggedArray.ConvertToNew(
jagged, typeof(int), typeof(Gaussian),
delegate(object elt) { return(new Gaussian(0.0, 1.0)); });
Assert.Equal(
JaggedArray.GetLength(jaggedGaussian, typeof(Gaussian)),
JaggedArray.GetLength(jagged, typeof(int)));
JaggedArray.ConvertElements(
jaggedGaussian, typeof(Gaussian),
delegate(object elt) { return(Gaussian.FromMeanAndPrecision(1.0, 2.0)); });
foreach (Gaussian g in JaggedArray.ElementIterator(jaggedGaussian, typeof(Gaussian)))
{
Assert.Equal(1.0, g.GetMean());
}
}
public void JaggedVisitor2()
{
Dirichlet d1 = new Dirichlet(2.0, 3.0, 4.0);
Dirichlet d2 = new Dirichlet(3.0, 4.0, 5.0);
Dirichlet d3 = d1 * d2;
Array jaggedDir1 = JaggedArray.ConvertToNew(
jagged, typeof(int), typeof(Dirichlet),
delegate(object elt) { return(new Dirichlet(d1)); });
Array jaggedDir2 = JaggedArray.ConvertToNew(
jagged, typeof(int), typeof(Dirichlet),
delegate(object elt) { return(new Dirichlet(d2)); });
JaggedArray.VisitElements2(
jaggedDir1, jaggedDir2, typeof(Dirichlet), typeof(Dirichlet),
delegate(object elt1, object elt2) { ((Dirichlet)elt1).PseudoCount.SetTo((((Dirichlet)elt1) * ((Dirichlet)elt2)).PseudoCount); });
foreach (Dirichlet d in JaggedArray.ElementIterator(jaggedDir1, typeof(Dirichlet)))
{
Assert.Equal(d.PseudoCount, d3.PseudoCount);
}
}
