IDualSurfaceCurve[] IDualSurfaceCurve.Split(double v) { ICurve[] splitted = curve3D.Split(v); if (splitted == null || splitted.Length != 2) { return(null); } IDualSurfaceCurve dsc1 = null, dsc2 = null; if (curve3D is InterpolatedDualSurfaceCurve) { dsc1 = new DualSurfaceCurve(splitted[0], surface1, (splitted[0] as InterpolatedDualSurfaceCurve).CurveOnSurface1, surface2, (splitted[0] as InterpolatedDualSurfaceCurve).CurveOnSurface2); dsc2 = new DualSurfaceCurve(splitted[1], surface1, (splitted[1] as InterpolatedDualSurfaceCurve).CurveOnSurface1, surface2, (splitted[1] as InterpolatedDualSurfaceCurve).CurveOnSurface2); } else { dsc1 = new DualSurfaceCurve(splitted[0], surface1, surface1.GetProjectedCurve(splitted[0], 0.0), surface2, surface2.GetProjectedCurve(splitted[0], 0.0)); dsc2 = new DualSurfaceCurve(splitted[1], surface1, surface1.GetProjectedCurve(splitted[1], 0.0), surface2, surface2.GetProjectedCurve(splitted[1], 0.0)); } return(new IDualSurfaceCurve[] { dsc1, dsc2 }); //GeoPoint2D uv1 = surface1.PositionOf(splitted[0].EndPoint); //GeoPoint2D uv2 = surface2.PositionOf(splitted[0].EndPoint); //double u1 = curve2D1.PositionOf(uv1); //double u2 = curve2D2.PositionOf(uv2); //GeoVector2D dir12d = curve2D1.DirectionAt(u1); //GeoVector2D dir22d = curve2D2.DirectionAt(u2); //GeoVector dir1 = dir12d.x*surface1.UDirection(uv1)+ dir12d.y * surface1.VDirection(uv1); //GeoVector dir2 = dir22d.x*surface2.UDirection(uv2)+ dir22d.y * surface2.VDirection(uv2); //GeoVector dir = curve3D.DirectionAt(v); //// dir1 should be in direction of curve3d, dir2 in the opposite dierection //if (dir1 * dir < 0) curve2D1.Reverse(); //if (dir2 * dir > 0) curve2D2.Reverse(); //ICurve2D[] splitted1 = curve2D1.Split(u1); //ICurve2D[] splitted2 = curve2D2.Split(u2); //if (splitted1 == null || splitted1.Length != 2) return null; //if (splitted2 == null || splitted2.Length != 2) return null; //IDualSurfaceCurve dsc1 = null, dsc2 = null; //dsc1 = new DualSurfaceCurve(splitted[0], surface1, splitted1[0], surface2, splitted2[1]); //dsc2 = new DualSurfaceCurve(splitted[1], surface1, splitted1[1], surface2, splitted2[0]); //return new IDualSurfaceCurve[] { dsc1, dsc2 }; }
/// <summary> /// Overrides <see cref="CADability.GeoObject.ISurfaceImpl.GetPlaneIntersection (PlaneSurface, double, double, double, double, double)"/> /// </summary> /// <param name="pl"></param> /// <param name="umin"></param> /// <param name="umax"></param> /// <param name="vmin"></param> /// <param name="vmax"></param> /// <param name="precision"></param> /// <returns></returns> public override IDualSurfaceCurve[] GetPlaneIntersection(PlaneSurface pl, double umin, double umax, double vmin, double vmax, double precision) { if (firstCurve is Line && secondCurve is Line) { if (Precision.IsPerpendicular(firstCurve.StartDirection, pl.Normal, false) && Precision.IsPerpendicular(secondCurve.StartDirection, pl.Normal, false)) { // eine Ebene, die zu beiden Linien parallel ist GeoPoint sp1, ep1, sp2, ep2; sp1 = firstCurve.PointAt(umin); ep1 = secondCurve.PointAt(umin); sp2 = firstCurve.PointAt(umax); ep2 = secondCurve.PointAt(umax); GeoPoint2D[] ip1 = pl.GetLineIntersection(sp1, ep1 - sp1); GeoPoint2D[] ip2 = pl.GetLineIntersection(sp2, ep2 - sp2); if (ip1.Length == 1 && ip2.Length == 1) { GeoPoint sp = pl.PointAt(ip1[0]); GeoPoint ep = pl.PointAt(ip2[0]); double v = Geometry.LinePar(sp1, ep1, sp); Line line = Line.Construct(); line.SetTwoPoints(sp, ep); DualSurfaceCurve dsc = new DualSurfaceCurve(line, this, new Line2D(new GeoPoint2D(umin, v), new GeoPoint2D(umax, v)), pl, new Line2D(ip1[0], ip2[0])); return(new IDualSurfaceCurve[] { dsc }); } } } if (firstCurve is Ellipse && secondCurve is Ellipse) { Ellipse e1 = firstCurve as Ellipse; Ellipse e2 = secondCurve as Ellipse; if (Precision.SameDirection(pl.Normal, e1.Plane.Normal, false) && Precision.SameDirection(pl.Normal, e2.Plane.Normal, false)) { if (e1.IsCircle && e2.IsCircle) { GeoPoint sp1, ep1, sp2, ep2, spm, epm; sp1 = firstCurve.PointAt(umin); ep1 = secondCurve.PointAt(umin); sp2 = firstCurve.PointAt(umax); ep2 = secondCurve.PointAt(umax); spm = firstCurve.PointAt((umin + umax) / 2.0); epm = secondCurve.PointAt((umin + umax) / 2.0); GeoPoint2D[] ip1 = pl.GetLineIntersection(sp1, ep1 - sp1); GeoPoint2D[] ip2 = pl.GetLineIntersection(sp2, ep2 - sp2); GeoPoint2D[] ipm = pl.GetLineIntersection(spm, epm - spm); if (ip1.Length == 1 && ip2.Length == 1 && ipm.Length == 1) { Ellipse e3 = Ellipse.Construct(); e3.SetArc3Points(pl.PointAt(ip1[0]), pl.PointAt(ipm[0]), pl.PointAt(ip2[0]), pl.Plane); double v = Geometry.LinePar(sp1, ep1, pl.PointAt(ip1[0])); DualSurfaceCurve dsc = new DualSurfaceCurve(e3, this, new Line2D(new GeoPoint2D(umin, v), new GeoPoint2D(umax, v)), pl, e3.GetProjectedCurve(pl.Plane)); return(new IDualSurfaceCurve[] { dsc }); } } } } PlanarState ps1 = firstCurve.GetPlanarState(); PlanarState ps2 = secondCurve.GetPlanarState(); if ((ps1 == PlanarState.UnderDetermined || ps1 == PlanarState.Planar) && (ps2 == PlanarState.UnderDetermined || ps2 == PlanarState.Planar)) { if (Precision.IsPerpendicular(firstCurve.StartDirection, pl.Normal, false) && Precision.IsPerpendicular(secondCurve.StartDirection, pl.Normal, false)) { // beide Kurven sind eben und parallel zur Schnittebene, wir haben also ein festes v und somit eine Zwischenkurve GeoPoint ip = pl.Plane.Intersect(firstCurve.StartPoint, secondCurve.StartPoint - firstCurve.StartPoint); double v = Geometry.LinePar(firstCurve.StartPoint, secondCurve.StartPoint, ip); ICurve cv = FixedV(v, umin, umax); DualSurfaceCurve dsc = new DualSurfaceCurve(cv, this, new Line2D(new GeoPoint2D(umin, v), new GeoPoint2D(umax, v)), pl, cv.GetProjectedCurve(pl.Plane)); return(new IDualSurfaceCurve[] { dsc }); } } return(base.GetPlaneIntersection(pl, umin, umax, vmin, vmax, precision)); }