public void SolutionUtilityTest()
{
var p = new Problem(nameof(SolutionUtilityTest));
var a = (Proposition)"a";
a.Utility = 1;
var b = (Proposition)"b";
b.Utility = -1;
var c = (Proposition)"c";
c.Utility = 1.5f;
var d = (Proposition)"d";
d.Utility = 3.3f;
var zero = (Proposition)"zeroUtility";
p.Quantify(1, 4, a, b, c, d, zero);
for (int i = 0; i < 100; i++)
{
var s = p.Solve();
float U(Proposition prop)
{
if (s[prop])
{
return(prop.Utility);
}
return(0);
}
Assert.IsTrue(Math.Abs(s.Utility - (U(a) + U(b) + U(c) + U(d) + U(zero))) < 0.0001f);
}
}
public void UtilityMaximizationTest()
{
var p = new Problem(nameof(UtilityMaximizationTest));
var a = (Proposition)"a";
a.Utility = 1;
var b = (Proposition)"b";
b.Utility = -1;
var c = (Proposition)"c";
c.Utility = 1.5f;
var d = (Proposition)"d";
d.Utility = 3.3f;
var zero = (Proposition)"zeroUtility";
p.Quantify(1, 4, a, b, c, d, zero);
for (int i = 0; i < 100; i++)
{
var s = p.HighUtilitySolution(1000);
Console.WriteLine(s.Model);
Assert.IsTrue(Math.Abs((1 + 1.5 + 3.3) - s.Utility) < 0.00001f);
}
}
public void FamilyGeneratorTest()
{
var p = new Problem("family generator");
var FamilySize = 25;
var generationCount = 5;
var matriarch = 1;
var cast = Enumerable.Range(matriarch, FamilySize).ToArray();
var kids = Enumerable.Range(matriarch + 1, FamilySize - 1).ToArray();
var childGenerations = Enumerable.Range(1, generationCount - 1).ToArray();
var female = Predicate <int>("female");
var generation = Predicate <int, int>("generation");
var parent = Predicate <int, int>("parent");
// Interestingly, this doesn't seem to speed things up.
//Func<int, int, Proposition> parent = (child, par) =>
// child > par ? p.GetProposition(Call.FromArgs(Problem.Current, "parent", child, par)) : false;
// Make of person # who is person number -who
int Mate(int who) => - who;
// Family must be 40-60% female
p.Quantify((int)(FamilySize * .4), (int)(FamilySize * .6), kids, female);
// Matriarch is the generation 0 female
p.Assert(female(matriarch), Not(female(Mate(matriarch))));
p.Assert(generation(matriarch, 0));
foreach (var child in kids)
{
// Everyone has exactly one parent from within the family (and one from outside)
p.Unique(cast, par => parent(child, par));
foreach (var par in cast)
{
parent(child, par).InitialProbability = 0;
}
// Everyone has a generation number
p.Unique(childGenerations, g => generation(child, g));
p.Assert(
// Only matriarch and patriarch are generation 0
Not(generation(child, 0)),
// Heteronormativity
female(child) == Not(female(Mate(child))));
foreach (var par in cast)
{
foreach (var g in childGenerations)
{
// Child's generation is one more than parent's generation
p.Assert(generation(child, g) <= (parent(child, par) & generation(par, g - 1)));
}
}
}
// Every generation has at least one kid
foreach (var g in childGenerations)
{
p.Exists(kids, k => generation(k, g));
}
p.Optimize();
Console.WriteLine(p.Stats);
Console.WriteLine(p.PerformanceStatistics);
Console.WriteLine();
for (int i = 0; i < 100; i++)
{
var s = p.Solve();
Console.WriteLine(p.PerformanceStatistics);
void PrintTree(int who, int depth)
{
for (int j = 0; j < depth; j++)
{
Console.Write(" ");
}
Console.WriteLine($"{Gender(who)} {who}: + {Mate(who)}");
foreach (var child in kids.Where(k => s[parent(k, who)]))
{
PrintTree(child, depth + 1);
}
}
string Gender(int who)
{
return(s[female(who)] ? "female" : "male");
}
PrintTree(matriarch, 0);
}
}
