private void CD_po_together(fault2d f1, fault2d f2)
{
util u1 = new util();
PO f1s, f1d, f2s, f2d, ptm, p1t, p2t, p3t;
f1s = f1.l0.p1;
f1d = f1.l0.p2;
f2s = f2.l0.p1;
f2d = f2.l0.p2;
//判断f1、f2的状态
#region judge topo
if (f1s.Me_to_po_length(f2d) == 0)
{
ptm = f2d;
f2d = f2s;
f2s = ptm;
}
//u1.InFile(u1.infopath, "ihere11");
if (f1d.Me_to_po_length(f2s) == 0)
{
ptm = f1d;
f1d = f1s;
f1s = ptm;
}
//u1.InFile(u1.infopath, "ihere12");
if (f1d.Me_to_po_length(f2d) == 0)
{
ptm = f1d;
f1d = f1s;
f1s = ptm;
ptm = f2d;
f2d = f2s;
f2s = ptm;
}
//u1.InFile(u1.infopath, "ihere13");
#endregion judge topo
p1t = u1.Po_vec_f(f1s, u1.cal_vec(f1s, f1d).unitVector(), 4);
p2t = u1.Po_vec_f(f2s, u1.cal_vec(f2s, f2d).unitVector(), 4);
p3t = new PO((p1t.x + p2t.x) / 2, (p1t.y + p2t.y) / 2);
p1t.key = 1;
p2t.key = 1;
p3t.key = 1;
PAINT.ll_tag.Insert(new LINE(f1s, p3t));
PAINT.pl_tag.Insert(p1t);
PAINT.pl_tag.Insert(p2t);
PAINT.pl_tag.Insert(p3t);
Cal_new_edge_po_insert(f1, f1s, p3t);
Cal_new_edge_po_insert(f2, f2s, p3t);
}
public override void GenerateNetPo()
{
int sitan;
float PI = 3.1415926f;
util u1 = new util();
sitan = 10;
PO vec, vec2, pc;//沿井身
vec = new PO(0, 0, 1);
vec2 = new PO();
pi = new PO();
pj = new PO();
for (int i = 5; i < 20; i = i + 4)
{
for (int j = 0; j < sitan; j++)
{
vec2.x = (float)System.Math.Cos(2 * PI * j / sitan);
vec2.y = (float)System.Math.Sin(2 * PI * j / sitan);
vec2.z = 0;
pc = u1.Po_vec_f(u1.vec_f(vec, 0), vec2, i);
if (i == 17 && j == 1)
{
pi = pc;
}
if (i == 17 && j == 2)
{
pj = pc;
}
PAINT.GetPL(pc);
ge_p_l.Insert(pc);
}
}
//PAINT.GetLL(new LINE(pi,pj));
}
public override void GenerateNetPo()
{
int sitan;
float PI = 3.1415926f;
util u1 = new util();
sitan = 10;
PO vec, vec2;//沿井身
vec = new PO(0, 0, 1);
vec2 = new PO();
pi = new PO();
pj = new PO();
pk = new PO();
for (int k = 0; k < 8; k = k + 2)
{
for (int i = -10; i < 20; i = i + 5)
{
for (int j = 0; j < sitan; j++)
{
vec2.x = (float)System.Math.Cos(2 * PI * j / sitan);
vec2.y = (float)System.Math.Sin(2 * PI * j / sitan);
vec2.z = 0;
if (k == 6 && i == -10 && j == 1)
{
pi = u1.Po_vec_f(u1.vec_f(vec, i), vec2, k);
}
if (k == 6 && i == -10 && j == 2)
{
pj = u1.Po_vec_f(u1.vec_f(vec, i), vec2, k);
}
if (k == 6 && i == -5 && j == 2)
{
pk = u1.Po_vec_f(u1.vec_f(vec, i), vec2, k);
}
ge_p_l.Insert(u1.Po_vec_f(u1.vec_f(vec, i), vec2, k));
}
}
}
}