/// <summary>
/// Fruchterman-Reingold-Layoutを適用します。
/// </summary>
/// <param name="width">描画領域幅</param>
/// <param name="height">描画領域高</param>
public void Fruchterman_Reingold_Layout(float width, float height)
{
if (vertexNum == 1)
return;
Vector2D[] Vpos = new Vector2D[vertexNum];
Vector2D[] Vdisp = new Vector2D[vertexNum];
Vector2D delta = new Vector2D();
float rate = 0.9f;
float k = (float)Math.Sqrt((double)(width * rate * height * rate) / vertexNum);
float t = 0.1f;
float d;
for (int i = 0; i < vertexNum; i++)
{
Vpos[i] = new Vector2D(((VertexView)vertex[i]).GetPosition().X, ((VertexView)vertex[i]).GetPosition().Y);
Vdisp[i] = new Vector2D();
}
for (int count = 0; count < 50; count++)
{
for (int i = 0; i < Vpos.Length; i++)
{
Vdisp[i].Init();
for (int j = 0; j < Vpos.Length; j++)
{
if (i != j)
{
delta = Vpos[i] - Vpos[j];
d = delta.Getnorm();
delta.Normalization();
Vdisp[i] += delta * f_r(d, k);
}
}
}
for (int i = 0; i < edgeNum; i++)
{
int n1 = vertex.IndexOf((VertexView)edge[i].GetEdgeVertex(Edge.VERTEX_FIRST));
int n2 = vertex.IndexOf((VertexView)edge[i].GetEdgeVertex(Edge.VERTEX_SECOND));
if (n1 != -1 && n2 != -1)
{
delta = Vpos[n1] - Vpos[n2];
d = delta.Getnorm();
delta.Normalization();
Vdisp[n1] -= delta * f_a(d, k);
Vdisp[n2] += delta * f_a(d, k);
}
}
for (int i = 0; i < Vpos.Length; i++)
{
d = Vdisp[i].Getnorm();
Vpos[i] += (Vdisp[i] / d) * Math.Min(d, t);
Vpos[i].X = (float)Math.Min(width * rate, Math.Max(width * (1 - rate), Vpos[i].X));
Vpos[i].Y = (float)Math.Min(height * rate, Math.Max(height * (1 - rate), Vpos[i].Y));
((VertexView)vertex[i]).SetPosition(Vpos[i].X, Vpos[i].Y);
}
}
}
/// <summary>
/// 指定された座標が Edge 上にあるかどうかを判断します。
/// </summary>
/// <param name="x">X座標</param>
/// <param name="y">Y座標</param>
public bool CheckMouseOn(float x, float y)
{
Vector2D vector1 = new Vector2D(((VertexView)vertex2).GetPosition().X - ((VertexView)vertex1).GetPosition().X, ((VertexView)vertex2).GetPosition().Y - ((VertexView)vertex1).GetPosition().Y);
Vector2D vector2 = new Vector2D(x - ((VertexView)vertex1).GetPosition().X, y - ((VertexView)vertex1).GetPosition().Y);
if (vector1.Getnorm() < vector2.Getnorm())
return false;
vector1.Normalization();
vector2.Normalization();
return (Vector2D.Inner(vector1, vector2) > 0.9999) ? true : false;
}
/// <summary>
/// Spring-Electrical-Layoutを適用します。
/// </summary>
/// <param name="width">描画領域幅</param>
/// <param name="height">描画領域高</param>
public void Spring_Electrical(float width, float height)
{
Vector2D[] velocity = new Vector2D[vertexNum]; //頂点ごとの速度
Vector2D force = new Vector2D(); //頂点に作用する力
Vector2D delta = new Vector2D(); //差分
float rate = 0.9f; //比率
for (int i = 0; i < vertexNum; i++)
{
force.Init();
for (int j = 0; j < vertexNum; j++)
{
if (i != j)
{
delta.Set(((VertexView)vertex[i]).GetPosition().X - ((VertexView)vertex[j]).GetPosition().X, ((VertexView)vertex[i]).GetPosition().Y - ((VertexView)vertex[j]).GetPosition().Y);
force += 500 * delta / (delta.Getnorm() * delta.Getnorm());
}
}
for (int j = 0; j < edgeNum; j++)
{
int index;
if ((index = edge[j].FindContainIndex(vertex[i])) != -1)
{
float X1, Y1, X2, Y2;
X1 = ((VertexView)vertex[i]).GetPosition().X;
Y1 = ((VertexView)vertex[i]).GetPosition().Y;
X2 = ((VertexView)edge[j].GetEdgeVertex((index + 1) % 2)).GetPosition().X;
Y2 = ((VertexView)edge[j].GetEdgeVertex((index + 1) % 2)).GetPosition().Y;
delta.Set(X2 - X1, Y2 - Y1);
force += 0.06f * delta;
}
}
velocity[i] = new Vector2D(force.X * 0.85f, force.Y * 0.85f);
}
for (int i = 0; i < vertexNum; i++)
{
delta.Set(((VertexView)vertex[i]).GetPosition().X, ((VertexView)vertex[i]).GetPosition().Y);
delta += velocity[i];
delta.X = (float)Math.Min(width * rate, Math.Max(width * (1 - rate), delta.X));
delta.Y = (float)Math.Min(height * rate, Math.Max(height * (1 - rate), delta.Y));
((VertexView)vertex[i]).SetPosition(delta.X, delta.Y);
}
}
