/// <summary>
///
/// </summary>
private static Orient3d GetHalfedgeTransform(E hedge, Func <V, Vector3d> getPosition, double unrollFactor)
{
if (unrollFactor < 0.0)
{
return(Orient3d.Identity);
}
var he0 = hedge;
var he1 = he0.Twin;
Vector3d p0 = getPosition(he0.Start);
Vector3d p1 = getPosition(he1.Start);
Vector3d p2 = (getPosition(he0.Previous.Start) + getPosition(he0.Next.End)) * 0.5;
Vector3d p3 = (getPosition(he1.Previous.Start) + getPosition(he1.Next.End)) * 0.5;
Vector3d x = p1 - p0;
Vector3d y0 = p2 - p0;
Vector3d y1 = (unrollFactor < 1.0) ? y0.SlerpTo(p1 - p3, unrollFactor) : p1 - p3; // TODO handle anti-parallel case
var t0 = new Orient3d(OrthoBasis3d.CreateFromXY(x, y0), p0);
var t1 = new Orient3d(OrthoBasis3d.CreateFromXY(x, y1), p0);
return(Orient3d.CreateFromTo(ref t0, ref t1));
}
/// <summary>
///
/// </summary>
private static Orient3d GetHalfedgeTransform(E hedge, Func <V, Vector3d> getPosition)
{
var he0 = hedge;
var he1 = he0.Twin;
Vector3d p0 = getPosition(he0.Start);
Vector3d p1 = getPosition(he1.Start);
Vector3d p2 = (getPosition(he0.Previous.Start) + getPosition(he0.Next.End)) * 0.5;
Vector3d p3 = (getPosition(he1.Previous.Start) + getPosition(he1.Next.End)) * 0.5;
Vector3d x = p1 - p0;
var t0 = new Orient3d(OrthoBasis3d.CreateFromXY(x, p2 - p0), p0);
var t1 = new Orient3d(OrthoBasis3d.CreateFromXY(x, p1 - p3), p0);
return(Orient3d.CreateFromTo(ref t0, ref t1));
}
/// <summary>
///
/// </summary>
/// <param name="plane"></param>
/// <returns></returns>
public static Orient3d ToOrient3d(this Plane plane)
{
return(new Orient3d(OrthoBasis3d.CreateFromXY(plane.XAxis, plane.YAxis), plane.Origin));
}
