public PositionSetEditSet getRandomPositionSet(int pointNum)
{
unchecked
{
int seed = (int)DateTime.Now.Ticks;
this.pointNum = pointNum;
LaplaceDistribution distributionX = new LaplaceDistribution(new StandardGenerator(seed++));
distributionX.Alpha = X_Alpha;
LaplaceDistribution distributionY = new LaplaceDistribution(new StandardGenerator(seed++));
distributionY.Alpha = Y_Alpha;
for (int i = 0; i < clusterPointNum; i++)
{
distributionX.Mu = RandomMaker.RapidBetween(minMu, maxMu);
distributionY.Mu = RandomMaker.RapidBetween(minMu, maxMu);
RandomPositionSet randomPositionSet =
new RandomPositionSet((int)(pointNum / clusterPointNum), 1000, distributionX, distributionY);
positionSetEditSet.AddPositionSet(randomPositionSet);
}
}
return(positionSetEditSet);
}
public static void lanch(IConvexHullEngine chEngine, int pointCount, float min, float max)
{
//求解凸包并作图
//IPositionSet ps = testData();
IPositionSet ps = new RandomPositionSet(pointCount, min, max);
//System.Console.Out.WriteLine("position set:");
//printPositionSet(ps);
IPositionSet cps = chEngine.ConvexHull(ps);
//System.Console.Out.WriteLine("Jarvis Match:");
//printPositionSet(cps);
PainterDialog.Clear();
PainterDialog.DrawPositionSet(ps);
PainterDialog.DrawConvexHull(cps);
PainterDialog.Show();
AnalyzeReport report = new AnalyzeReport();
//测试点集引用
cps.InitToTraverseSet();
ps.InitToTraverseSet();
bool correct = true;
while (correct && cps.NextPosition())
{
IPosition cp = cps.GetPosition();
bool find = false;
ps.InitToTraverseSet();
while (!find && ps.NextPosition())
{
IPosition p = ps.GetPosition();
if (cp == p)
{
find = true;
}
}
if (!find)
{
correct = false;
}
}
report.content += "引用测试：" + (correct ? "正确" : "错误") + "\n";
//测试凸包
if (correct)
{
report.content += "凸包正确性测试：\n";
IPositionSet cps_ref = (new QuickHull()).ConvexHull(ps);
// System.Console.Out.WriteLine("Quick Hull:");
// printPositionSet(cps_ref);
IPosition[] cpa = (IPosition[])(cps.ToArray());
IPosition[] cpa_ref = (IPosition[])(cps_ref.ToArray());
if (cpa.Length != cpa_ref.Length)
{
report.content += "数目不等";
correct = false;
}
if (correct)
{
int n = cpa.Length;
int m = 0;
int p = 0;
for (; p < n; p++)
{
if (cpa[0] == cpa_ref[p])
{
break;
}
}
if (p == n)
{
correct = false;
report.content += "发生错误！";
}
if (correct)
{
for (int i = 0; i < n; i++)
{
if (cpa[i] == cpa_ref[(p + i + n) % n])
{
m++;
}
}
int tm = 0;
for (int i = 0; i < n; i++)
{
if (cpa[i] == cpa_ref[(p - i + n) % n])
{
tm++;
}
}
if (tm > m)
{
m = tm;
}
report.content += "正确率：" + m.ToString() + "/" + n.ToString();
if (m == n)
{
report.content += "正确！";
}
else
{
report.content += "不正确！";
}
}
}
}
report.Show();
}
