| /* |
| * Copyright (c) 1997, 2006, Oracle and/or its affiliates. All rights reserved. |
| * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. |
| * |
| * This code is free software; you can redistribute it and/or modify it |
| * under the terms of the GNU General Public License version 2 only, as |
| * published by the Free Software Foundation. Oracle designates this |
| * particular file as subject to the "Classpath" exception as provided |
| * by Oracle in the LICENSE file that accompanied this code. |
| * |
| * This code is distributed in the hope that it will be useful, but WITHOUT |
| * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| * version 2 for more details (a copy is included in the LICENSE file that |
| * accompanied this code). |
| * |
| * You should have received a copy of the GNU General Public License version |
| * 2 along with this work; if not, write to the Free Software Foundation, |
| * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. |
| * |
| * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA |
| * or visit www.oracle.com if you need additional information or have any |
| * questions. |
| */ |
| |
| package java.awt.geom; |
| |
| import java.awt.Shape; |
| import java.awt.Rectangle; |
| import java.util.Arrays; |
| import java.io.Serializable; |
| import sun.awt.geom.Curve; |
| |
| /** |
| * The <code>CubicCurve2D</code> class defines a cubic parametric curve |
| * segment in {@code (x,y)} coordinate space. |
| * <p> |
| * This class is only the abstract superclass for all objects which |
| * store a 2D cubic curve segment. |
| * The actual storage representation of the coordinates is left to |
| * the subclass. |
| * |
| * @author Jim Graham |
| * @since 1.2 |
| */ |
| public abstract class CubicCurve2D implements Shape, Cloneable { |
| |
| /** |
| * A cubic parametric curve segment specified with |
| * {@code float} coordinates. |
| * @since 1.2 |
| */ |
| public static class Float extends CubicCurve2D implements Serializable { |
| /** |
| * The X coordinate of the start point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float x1; |
| |
| /** |
| * The Y coordinate of the start point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float y1; |
| |
| /** |
| * The X coordinate of the first control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float ctrlx1; |
| |
| /** |
| * The Y coordinate of the first control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float ctrly1; |
| |
| /** |
| * The X coordinate of the second control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float ctrlx2; |
| |
| /** |
| * The Y coordinate of the second control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float ctrly2; |
| |
| /** |
| * The X coordinate of the end point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float x2; |
| |
| /** |
| * The Y coordinate of the end point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public float y2; |
| |
| /** |
| * Constructs and initializes a CubicCurve with coordinates |
| * (0, 0, 0, 0, 0, 0, 0, 0). |
| * @since 1.2 |
| */ |
| public Float() { |
| } |
| |
| /** |
| * Constructs and initializes a {@code CubicCurve2D} from |
| * the specified {@code float} coordinates. |
| * |
| * @param x1 the X coordinate for the start point |
| * of the resulting {@code CubicCurve2D} |
| * @param y1 the Y coordinate for the start point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate for the first control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate for the first control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate for the second control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate for the second control point |
| * of the resulting {@code CubicCurve2D} |
| * @param x2 the X coordinate for the end point |
| * of the resulting {@code CubicCurve2D} |
| * @param y2 the Y coordinate for the end point |
| * of the resulting {@code CubicCurve2D} |
| * @since 1.2 |
| */ |
| public Float(float x1, float y1, |
| float ctrlx1, float ctrly1, |
| float ctrlx2, float ctrly2, |
| float x2, float y2) |
| { |
| setCurve(x1, y1, ctrlx1, ctrly1, ctrlx2, ctrly2, x2, y2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getX1() { |
| return (double) x1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getY1() { |
| return (double) y1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getP1() { |
| return new Point2D.Float(x1, y1); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlX1() { |
| return (double) ctrlx1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlY1() { |
| return (double) ctrly1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getCtrlP1() { |
| return new Point2D.Float(ctrlx1, ctrly1); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlX2() { |
| return (double) ctrlx2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlY2() { |
| return (double) ctrly2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getCtrlP2() { |
| return new Point2D.Float(ctrlx2, ctrly2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getX2() { |
| return (double) x2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getY2() { |
| return (double) y2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getP2() { |
| return new Point2D.Float(x2, y2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public void setCurve(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2) |
| { |
| this.x1 = (float) x1; |
| this.y1 = (float) y1; |
| this.ctrlx1 = (float) ctrlx1; |
| this.ctrly1 = (float) ctrly1; |
| this.ctrlx2 = (float) ctrlx2; |
| this.ctrly2 = (float) ctrly2; |
| this.x2 = (float) x2; |
| this.y2 = (float) y2; |
| } |
| |
| /** |
| * Sets the location of the end points and control points |
| * of this curve to the specified {@code float} coordinates. |
| * |
| * @param x1 the X coordinate used to set the start point |
| * of this {@code CubicCurve2D} |
| * @param y1 the Y coordinate used to set the start point |
| * of this {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate used to set the first control point |
| * of this {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate used to set the first control point |
| * of this {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate used to set the second control point |
| * of this {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate used to set the second control point |
| * of this {@code CubicCurve2D} |
| * @param x2 the X coordinate used to set the end point |
| * of this {@code CubicCurve2D} |
| * @param y2 the Y coordinate used to set the end point |
| * of this {@code CubicCurve2D} |
| * @since 1.2 |
| */ |
| public void setCurve(float x1, float y1, |
| float ctrlx1, float ctrly1, |
| float ctrlx2, float ctrly2, |
| float x2, float y2) |
| { |
| this.x1 = x1; |
| this.y1 = y1; |
| this.ctrlx1 = ctrlx1; |
| this.ctrly1 = ctrly1; |
| this.ctrlx2 = ctrlx2; |
| this.ctrly2 = ctrly2; |
| this.x2 = x2; |
| this.y2 = y2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Rectangle2D getBounds2D() { |
| float left = Math.min(Math.min(x1, x2), |
| Math.min(ctrlx1, ctrlx2)); |
| float top = Math.min(Math.min(y1, y2), |
| Math.min(ctrly1, ctrly2)); |
| float right = Math.max(Math.max(x1, x2), |
| Math.max(ctrlx1, ctrlx2)); |
| float bottom = Math.max(Math.max(y1, y2), |
| Math.max(ctrly1, ctrly2)); |
| return new Rectangle2D.Float(left, top, |
| right - left, bottom - top); |
| } |
| |
| /* |
| * JDK 1.6 serialVersionUID |
| */ |
| private static final long serialVersionUID = -1272015596714244385L; |
| } |
| |
| /** |
| * A cubic parametric curve segment specified with |
| * {@code double} coordinates. |
| * @since 1.2 |
| */ |
| public static class Double extends CubicCurve2D implements Serializable { |
| /** |
| * The X coordinate of the start point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double x1; |
| |
| /** |
| * The Y coordinate of the start point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double y1; |
| |
| /** |
| * The X coordinate of the first control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double ctrlx1; |
| |
| /** |
| * The Y coordinate of the first control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double ctrly1; |
| |
| /** |
| * The X coordinate of the second control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double ctrlx2; |
| |
| /** |
| * The Y coordinate of the second control point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double ctrly2; |
| |
| /** |
| * The X coordinate of the end point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double x2; |
| |
| /** |
| * The Y coordinate of the end point |
| * of the cubic curve segment. |
| * @since 1.2 |
| * @serial |
| */ |
| public double y2; |
| |
| /** |
| * Constructs and initializes a CubicCurve with coordinates |
| * (0, 0, 0, 0, 0, 0, 0, 0). |
| * @since 1.2 |
| */ |
| public Double() { |
| } |
| |
| /** |
| * Constructs and initializes a {@code CubicCurve2D} from |
| * the specified {@code double} coordinates. |
| * |
| * @param x1 the X coordinate for the start point |
| * of the resulting {@code CubicCurve2D} |
| * @param y1 the Y coordinate for the start point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate for the first control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate for the first control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate for the second control point |
| * of the resulting {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate for the second control point |
| * of the resulting {@code CubicCurve2D} |
| * @param x2 the X coordinate for the end point |
| * of the resulting {@code CubicCurve2D} |
| * @param y2 the Y coordinate for the end point |
| * of the resulting {@code CubicCurve2D} |
| * @since 1.2 |
| */ |
| public Double(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2) |
| { |
| setCurve(x1, y1, ctrlx1, ctrly1, ctrlx2, ctrly2, x2, y2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getX1() { |
| return x1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getY1() { |
| return y1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getP1() { |
| return new Point2D.Double(x1, y1); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlX1() { |
| return ctrlx1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlY1() { |
| return ctrly1; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getCtrlP1() { |
| return new Point2D.Double(ctrlx1, ctrly1); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlX2() { |
| return ctrlx2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getCtrlY2() { |
| return ctrly2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getCtrlP2() { |
| return new Point2D.Double(ctrlx2, ctrly2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getX2() { |
| return x2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public double getY2() { |
| return y2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Point2D getP2() { |
| return new Point2D.Double(x2, y2); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public void setCurve(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2) |
| { |
| this.x1 = x1; |
| this.y1 = y1; |
| this.ctrlx1 = ctrlx1; |
| this.ctrly1 = ctrly1; |
| this.ctrlx2 = ctrlx2; |
| this.ctrly2 = ctrly2; |
| this.x2 = x2; |
| this.y2 = y2; |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Rectangle2D getBounds2D() { |
| double left = Math.min(Math.min(x1, x2), |
| Math.min(ctrlx1, ctrlx2)); |
| double top = Math.min(Math.min(y1, y2), |
| Math.min(ctrly1, ctrly2)); |
| double right = Math.max(Math.max(x1, x2), |
| Math.max(ctrlx1, ctrlx2)); |
| double bottom = Math.max(Math.max(y1, y2), |
| Math.max(ctrly1, ctrly2)); |
| return new Rectangle2D.Double(left, top, |
| right - left, bottom - top); |
| } |
| |
| /* |
| * JDK 1.6 serialVersionUID |
| */ |
| private static final long serialVersionUID = -4202960122839707295L; |
| } |
| |
| /** |
| * This is an abstract class that cannot be instantiated directly. |
| * Type-specific implementation subclasses are available for |
| * instantiation and provide a number of formats for storing |
| * the information necessary to satisfy the various accessor |
| * methods below. |
| * |
| * @see java.awt.geom.CubicCurve2D.Float |
| * @see java.awt.geom.CubicCurve2D.Double |
| * @since 1.2 |
| */ |
| protected CubicCurve2D() { |
| } |
| |
| /** |
| * Returns the X coordinate of the start point in double precision. |
| * @return the X coordinate of the start point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getX1(); |
| |
| /** |
| * Returns the Y coordinate of the start point in double precision. |
| * @return the Y coordinate of the start point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getY1(); |
| |
| /** |
| * Returns the start point. |
| * @return a {@code Point2D} that is the start point of |
| * the {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract Point2D getP1(); |
| |
| /** |
| * Returns the X coordinate of the first control point in double precision. |
| * @return the X coordinate of the first control point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getCtrlX1(); |
| |
| /** |
| * Returns the Y coordinate of the first control point in double precision. |
| * @return the Y coordinate of the first control point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getCtrlY1(); |
| |
| /** |
| * Returns the first control point. |
| * @return a {@code Point2D} that is the first control point of |
| * the {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract Point2D getCtrlP1(); |
| |
| /** |
| * Returns the X coordinate of the second control point |
| * in double precision. |
| * @return the X coordinate of the second control point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getCtrlX2(); |
| |
| /** |
| * Returns the Y coordinate of the second control point |
| * in double precision. |
| * @return the Y coordinate of the second control point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getCtrlY2(); |
| |
| /** |
| * Returns the second control point. |
| * @return a {@code Point2D} that is the second control point of |
| * the {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract Point2D getCtrlP2(); |
| |
| /** |
| * Returns the X coordinate of the end point in double precision. |
| * @return the X coordinate of the end point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getX2(); |
| |
| /** |
| * Returns the Y coordinate of the end point in double precision. |
| * @return the Y coordinate of the end point of the |
| * {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract double getY2(); |
| |
| /** |
| * Returns the end point. |
| * @return a {@code Point2D} that is the end point of |
| * the {@code CubicCurve2D}. |
| * @since 1.2 |
| */ |
| public abstract Point2D getP2(); |
| |
| /** |
| * Sets the location of the end points and control points of this curve |
| * to the specified double coordinates. |
| * |
| * @param x1 the X coordinate used to set the start point |
| * of this {@code CubicCurve2D} |
| * @param y1 the Y coordinate used to set the start point |
| * of this {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate used to set the first control point |
| * of this {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate used to set the first control point |
| * of this {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate used to set the second control point |
| * of this {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate used to set the second control point |
| * of this {@code CubicCurve2D} |
| * @param x2 the X coordinate used to set the end point |
| * of this {@code CubicCurve2D} |
| * @param y2 the Y coordinate used to set the end point |
| * of this {@code CubicCurve2D} |
| * @since 1.2 |
| */ |
| public abstract void setCurve(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2); |
| |
| /** |
| * Sets the location of the end points and control points of this curve |
| * to the double coordinates at the specified offset in the specified |
| * array. |
| * @param coords a double array containing coordinates |
| * @param offset the index of <code>coords</code> from which to begin |
| * setting the end points and control points of this curve |
| * to the coordinates contained in <code>coords</code> |
| * @since 1.2 |
| */ |
| public void setCurve(double[] coords, int offset) { |
| setCurve(coords[offset + 0], coords[offset + 1], |
| coords[offset + 2], coords[offset + 3], |
| coords[offset + 4], coords[offset + 5], |
| coords[offset + 6], coords[offset + 7]); |
| } |
| |
| /** |
| * Sets the location of the end points and control points of this curve |
| * to the specified <code>Point2D</code> coordinates. |
| * @param p1 the first specified <code>Point2D</code> used to set the |
| * start point of this curve |
| * @param cp1 the second specified <code>Point2D</code> used to set the |
| * first control point of this curve |
| * @param cp2 the third specified <code>Point2D</code> used to set the |
| * second control point of this curve |
| * @param p2 the fourth specified <code>Point2D</code> used to set the |
| * end point of this curve |
| * @since 1.2 |
| */ |
| public void setCurve(Point2D p1, Point2D cp1, Point2D cp2, Point2D p2) { |
| setCurve(p1.getX(), p1.getY(), cp1.getX(), cp1.getY(), |
| cp2.getX(), cp2.getY(), p2.getX(), p2.getY()); |
| } |
| |
| /** |
| * Sets the location of the end points and control points of this curve |
| * to the coordinates of the <code>Point2D</code> objects at the specified |
| * offset in the specified array. |
| * @param pts an array of <code>Point2D</code> objects |
| * @param offset the index of <code>pts</code> from which to begin setting |
| * the end points and control points of this curve to the |
| * points contained in <code>pts</code> |
| * @since 1.2 |
| */ |
| public void setCurve(Point2D[] pts, int offset) { |
| setCurve(pts[offset + 0].getX(), pts[offset + 0].getY(), |
| pts[offset + 1].getX(), pts[offset + 1].getY(), |
| pts[offset + 2].getX(), pts[offset + 2].getY(), |
| pts[offset + 3].getX(), pts[offset + 3].getY()); |
| } |
| |
| /** |
| * Sets the location of the end points and control points of this curve |
| * to the same as those in the specified <code>CubicCurve2D</code>. |
| * @param c the specified <code>CubicCurve2D</code> |
| * @since 1.2 |
| */ |
| public void setCurve(CubicCurve2D c) { |
| setCurve(c.getX1(), c.getY1(), c.getCtrlX1(), c.getCtrlY1(), |
| c.getCtrlX2(), c.getCtrlY2(), c.getX2(), c.getY2()); |
| } |
| |
| /** |
| * Returns the square of the flatness of the cubic curve specified |
| * by the indicated control points. The flatness is the maximum distance |
| * of a control point from the line connecting the end points. |
| * |
| * @param x1 the X coordinate that specifies the start point |
| * of a {@code CubicCurve2D} |
| * @param y1 the Y coordinate that specifies the start point |
| * of a {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate that specifies the first control point |
| * of a {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate that specifies the first control point |
| * of a {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate that specifies the second control point |
| * of a {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate that specifies the second control point |
| * of a {@code CubicCurve2D} |
| * @param x2 the X coordinate that specifies the end point |
| * of a {@code CubicCurve2D} |
| * @param y2 the Y coordinate that specifies the end point |
| * of a {@code CubicCurve2D} |
| * @return the square of the flatness of the {@code CubicCurve2D} |
| * represented by the specified coordinates. |
| * @since 1.2 |
| */ |
| public static double getFlatnessSq(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2) { |
| return Math.max(Line2D.ptSegDistSq(x1, y1, x2, y2, ctrlx1, ctrly1), |
| Line2D.ptSegDistSq(x1, y1, x2, y2, ctrlx2, ctrly2)); |
| |
| } |
| |
| /** |
| * Returns the flatness of the cubic curve specified |
| * by the indicated control points. The flatness is the maximum distance |
| * of a control point from the line connecting the end points. |
| * |
| * @param x1 the X coordinate that specifies the start point |
| * of a {@code CubicCurve2D} |
| * @param y1 the Y coordinate that specifies the start point |
| * of a {@code CubicCurve2D} |
| * @param ctrlx1 the X coordinate that specifies the first control point |
| * of a {@code CubicCurve2D} |
| * @param ctrly1 the Y coordinate that specifies the first control point |
| * of a {@code CubicCurve2D} |
| * @param ctrlx2 the X coordinate that specifies the second control point |
| * of a {@code CubicCurve2D} |
| * @param ctrly2 the Y coordinate that specifies the second control point |
| * of a {@code CubicCurve2D} |
| * @param x2 the X coordinate that specifies the end point |
| * of a {@code CubicCurve2D} |
| * @param y2 the Y coordinate that specifies the end point |
| * of a {@code CubicCurve2D} |
| * @return the flatness of the {@code CubicCurve2D} |
| * represented by the specified coordinates. |
| * @since 1.2 |
| */ |
| public static double getFlatness(double x1, double y1, |
| double ctrlx1, double ctrly1, |
| double ctrlx2, double ctrly2, |
| double x2, double y2) { |
| return Math.sqrt(getFlatnessSq(x1, y1, ctrlx1, ctrly1, |
| ctrlx2, ctrly2, x2, y2)); |
| } |
| |
| /** |
| * Returns the square of the flatness of the cubic curve specified |
| * by the control points stored in the indicated array at the |
| * indicated index. The flatness is the maximum distance |
| * of a control point from the line connecting the end points. |
| * @param coords an array containing coordinates |
| * @param offset the index of <code>coords</code> from which to begin |
| * getting the end points and control points of the curve |
| * @return the square of the flatness of the <code>CubicCurve2D</code> |
| * specified by the coordinates in <code>coords</code> at |
| * the specified offset. |
| * @since 1.2 |
| */ |
| public static double getFlatnessSq(double coords[], int offset) { |
| return getFlatnessSq(coords[offset + 0], coords[offset + 1], |
| coords[offset + 2], coords[offset + 3], |
| coords[offset + 4], coords[offset + 5], |
| coords[offset + 6], coords[offset + 7]); |
| } |
| |
| /** |
| * Returns the flatness of the cubic curve specified |
| * by the control points stored in the indicated array at the |
| * indicated index. The flatness is the maximum distance |
| * of a control point from the line connecting the end points. |
| * @param coords an array containing coordinates |
| * @param offset the index of <code>coords</code> from which to begin |
| * getting the end points and control points of the curve |
| * @return the flatness of the <code>CubicCurve2D</code> |
| * specified by the coordinates in <code>coords</code> at |
| * the specified offset. |
| * @since 1.2 |
| */ |
| public static double getFlatness(double coords[], int offset) { |
| return getFlatness(coords[offset + 0], coords[offset + 1], |
| coords[offset + 2], coords[offset + 3], |
| coords[offset + 4], coords[offset + 5], |
| coords[offset + 6], coords[offset + 7]); |
| } |
| |
| /** |
| * Returns the square of the flatness of this curve. The flatness is the |
| * maximum distance of a control point from the line connecting the |
| * end points. |
| * @return the square of the flatness of this curve. |
| * @since 1.2 |
| */ |
| public double getFlatnessSq() { |
| return getFlatnessSq(getX1(), getY1(), getCtrlX1(), getCtrlY1(), |
| getCtrlX2(), getCtrlY2(), getX2(), getY2()); |
| } |
| |
| /** |
| * Returns the flatness of this curve. The flatness is the |
| * maximum distance of a control point from the line connecting the |
| * end points. |
| * @return the flatness of this curve. |
| * @since 1.2 |
| */ |
| public double getFlatness() { |
| return getFlatness(getX1(), getY1(), getCtrlX1(), getCtrlY1(), |
| getCtrlX2(), getCtrlY2(), getX2(), getY2()); |
| } |
| |
| /** |
| * Subdivides this cubic curve and stores the resulting two |
| * subdivided curves into the left and right curve parameters. |
| * Either or both of the left and right objects may be the same |
| * as this object or null. |
| * @param left the cubic curve object for storing for the left or |
| * first half of the subdivided curve |
| * @param right the cubic curve object for storing for the right or |
| * second half of the subdivided curve |
| * @since 1.2 |
| */ |
| public void subdivide(CubicCurve2D left, CubicCurve2D right) { |
| subdivide(this, left, right); |
| } |
| |
| /** |
| * Subdivides the cubic curve specified by the <code>src</code> parameter |
| * and stores the resulting two subdivided curves into the |
| * <code>left</code> and <code>right</code> curve parameters. |
| * Either or both of the <code>left</code> and <code>right</code> objects |
| * may be the same as the <code>src</code> object or <code>null</code>. |
| * @param src the cubic curve to be subdivided |
| * @param left the cubic curve object for storing the left or |
| * first half of the subdivided curve |
| * @param right the cubic curve object for storing the right or |
| * second half of the subdivided curve |
| * @since 1.2 |
| */ |
| public static void subdivide(CubicCurve2D src, |
| CubicCurve2D left, |
| CubicCurve2D right) { |
| double x1 = src.getX1(); |
| double y1 = src.getY1(); |
| double ctrlx1 = src.getCtrlX1(); |
| double ctrly1 = src.getCtrlY1(); |
| double ctrlx2 = src.getCtrlX2(); |
| double ctrly2 = src.getCtrlY2(); |
| double x2 = src.getX2(); |
| double y2 = src.getY2(); |
| double centerx = (ctrlx1 + ctrlx2) / 2.0; |
| double centery = (ctrly1 + ctrly2) / 2.0; |
| ctrlx1 = (x1 + ctrlx1) / 2.0; |
| ctrly1 = (y1 + ctrly1) / 2.0; |
| ctrlx2 = (x2 + ctrlx2) / 2.0; |
| ctrly2 = (y2 + ctrly2) / 2.0; |
| double ctrlx12 = (ctrlx1 + centerx) / 2.0; |
| double ctrly12 = (ctrly1 + centery) / 2.0; |
| double ctrlx21 = (ctrlx2 + centerx) / 2.0; |
| double ctrly21 = (ctrly2 + centery) / 2.0; |
| centerx = (ctrlx12 + ctrlx21) / 2.0; |
| centery = (ctrly12 + ctrly21) / 2.0; |
| if (left != null) { |
| left.setCurve(x1, y1, ctrlx1, ctrly1, |
| ctrlx12, ctrly12, centerx, centery); |
| } |
| if (right != null) { |
| right.setCurve(centerx, centery, ctrlx21, ctrly21, |
| ctrlx2, ctrly2, x2, y2); |
| } |
| } |
| |
| /** |
| * Subdivides the cubic curve specified by the coordinates |
| * stored in the <code>src</code> array at indices <code>srcoff</code> |
| * through (<code>srcoff</code> + 7) and stores the |
| * resulting two subdivided curves into the two result arrays at the |
| * corresponding indices. |
| * Either or both of the <code>left</code> and <code>right</code> |
| * arrays may be <code>null</code> or a reference to the same array |
| * as the <code>src</code> array. |
| * Note that the last point in the first subdivided curve is the |
| * same as the first point in the second subdivided curve. Thus, |
| * it is possible to pass the same array for <code>left</code> |
| * and <code>right</code> and to use offsets, such as <code>rightoff</code> |
| * equals (<code>leftoff</code> + 6), in order |
| * to avoid allocating extra storage for this common point. |
| * @param src the array holding the coordinates for the source curve |
| * @param srcoff the offset into the array of the beginning of the |
| * the 6 source coordinates |
| * @param left the array for storing the coordinates for the first |
| * half of the subdivided curve |
| * @param leftoff the offset into the array of the beginning of the |
| * the 6 left coordinates |
| * @param right the array for storing the coordinates for the second |
| * half of the subdivided curve |
| * @param rightoff the offset into the array of the beginning of the |
| * the 6 right coordinates |
| * @since 1.2 |
| */ |
| public static void subdivide(double src[], int srcoff, |
| double left[], int leftoff, |
| double right[], int rightoff) { |
| double x1 = src[srcoff + 0]; |
| double y1 = src[srcoff + 1]; |
| double ctrlx1 = src[srcoff + 2]; |
| double ctrly1 = src[srcoff + 3]; |
| double ctrlx2 = src[srcoff + 4]; |
| double ctrly2 = src[srcoff + 5]; |
| double x2 = src[srcoff + 6]; |
| double y2 = src[srcoff + 7]; |
| if (left != null) { |
| left[leftoff + 0] = x1; |
| left[leftoff + 1] = y1; |
| } |
| if (right != null) { |
| right[rightoff + 6] = x2; |
| right[rightoff + 7] = y2; |
| } |
| x1 = (x1 + ctrlx1) / 2.0; |
| y1 = (y1 + ctrly1) / 2.0; |
| x2 = (x2 + ctrlx2) / 2.0; |
| y2 = (y2 + ctrly2) / 2.0; |
| double centerx = (ctrlx1 + ctrlx2) / 2.0; |
| double centery = (ctrly1 + ctrly2) / 2.0; |
| ctrlx1 = (x1 + centerx) / 2.0; |
| ctrly1 = (y1 + centery) / 2.0; |
| ctrlx2 = (x2 + centerx) / 2.0; |
| ctrly2 = (y2 + centery) / 2.0; |
| centerx = (ctrlx1 + ctrlx2) / 2.0; |
| centery = (ctrly1 + ctrly2) / 2.0; |
| if (left != null) { |
| left[leftoff + 2] = x1; |
| left[leftoff + 3] = y1; |
| left[leftoff + 4] = ctrlx1; |
| left[leftoff + 5] = ctrly1; |
| left[leftoff + 6] = centerx; |
| left[leftoff + 7] = centery; |
| } |
| if (right != null) { |
| right[rightoff + 0] = centerx; |
| right[rightoff + 1] = centery; |
| right[rightoff + 2] = ctrlx2; |
| right[rightoff + 3] = ctrly2; |
| right[rightoff + 4] = x2; |
| right[rightoff + 5] = y2; |
| } |
| } |
| |
| /** |
| * Solves the cubic whose coefficients are in the <code>eqn</code> |
| * array and places the non-complex roots back into the same array, |
| * returning the number of roots. The solved cubic is represented |
| * by the equation: |
| * <pre> |
| * eqn = {c, b, a, d} |
| * dx^3 + ax^2 + bx + c = 0 |
| * </pre> |
| * A return value of -1 is used to distinguish a constant equation |
| * that might be always 0 or never 0 from an equation that has no |
| * zeroes. |
| * @param eqn an array containing coefficients for a cubic |
| * @return the number of roots, or -1 if the equation is a constant. |
| * @since 1.2 |
| */ |
| public static int solveCubic(double eqn[]) { |
| return solveCubic(eqn, eqn); |
| } |
| |
| /** |
| * Solve the cubic whose coefficients are in the <code>eqn</code> |
| * array and place the non-complex roots into the <code>res</code> |
| * array, returning the number of roots. |
| * The cubic solved is represented by the equation: |
| * eqn = {c, b, a, d} |
| * dx^3 + ax^2 + bx + c = 0 |
| * A return value of -1 is used to distinguish a constant equation, |
| * which may be always 0 or never 0, from an equation which has no |
| * zeroes. |
| * @param eqn the specified array of coefficients to use to solve |
| * the cubic equation |
| * @param res the array that contains the non-complex roots |
| * resulting from the solution of the cubic equation |
| * @return the number of roots, or -1 if the equation is a constant |
| * @since 1.3 |
| */ |
| public static int solveCubic(double eqn[], double res[]) { |
| // From Numerical Recipes, 5.6, Quadratic and Cubic Equations |
| double d = eqn[3]; |
| if (d == 0.0) { |
| // The cubic has degenerated to quadratic (or line or ...). |
| return QuadCurve2D.solveQuadratic(eqn, res); |
| } |
| double a = eqn[2] / d; |
| double b = eqn[1] / d; |
| double c = eqn[0] / d; |
| int roots = 0; |
| double Q = (a * a - 3.0 * b) / 9.0; |
| double R = (2.0 * a * a * a - 9.0 * a * b + 27.0 * c) / 54.0; |
| double R2 = R * R; |
| double Q3 = Q * Q * Q; |
| a = a / 3.0; |
| if (R2 < Q3) { |
| double theta = Math.acos(R / Math.sqrt(Q3)); |
| Q = -2.0 * Math.sqrt(Q); |
| if (res == eqn) { |
| // Copy the eqn so that we don't clobber it with the |
| // roots. This is needed so that fixRoots can do its |
| // work with the original equation. |
| eqn = new double[4]; |
| System.arraycopy(res, 0, eqn, 0, 4); |
| } |
| res[roots++] = Q * Math.cos(theta / 3.0) - a; |
| res[roots++] = Q * Math.cos((theta + Math.PI * 2.0)/ 3.0) - a; |
| res[roots++] = Q * Math.cos((theta - Math.PI * 2.0)/ 3.0) - a; |
| fixRoots(res, eqn); |
| } else { |
| boolean neg = (R < 0.0); |
| double S = Math.sqrt(R2 - Q3); |
| if (neg) { |
| R = -R; |
| } |
| double A = Math.pow(R + S, 1.0 / 3.0); |
| if (!neg) { |
| A = -A; |
| } |
| double B = (A == 0.0) ? 0.0 : (Q / A); |
| res[roots++] = (A + B) - a; |
| } |
| return roots; |
| } |
| |
| /* |
| * This pruning step is necessary since solveCubic uses the |
| * cosine function to calculate the roots when there are 3 |
| * of them. Since the cosine method can have an error of |
| * +/- 1E-14 we need to make sure that we don't make any |
| * bad decisions due to an error. |
| * |
| * If the root is not near one of the endpoints, then we will |
| * only have a slight inaccuracy in calculating the x intercept |
| * which will only cause a slightly wrong answer for some |
| * points very close to the curve. While the results in that |
| * case are not as accurate as they could be, they are not |
| * disastrously inaccurate either. |
| * |
| * On the other hand, if the error happens near one end of |
| * the curve, then our processing to reject values outside |
| * of the t=[0,1] range will fail and the results of that |
| * failure will be disastrous since for an entire horizontal |
| * range of test points, we will either overcount or undercount |
| * the crossings and get a wrong answer for all of them, even |
| * when they are clearly and obviously inside or outside the |
| * curve. |
| * |
| * To work around this problem, we try a couple of Newton-Raphson |
| * iterations to see if the true root is closer to the endpoint |
| * or further away. If it is further away, then we can stop |
| * since we know we are on the right side of the endpoint. If |
| * we change direction, then either we are now being dragged away |
| * from the endpoint in which case the first condition will cause |
| * us to stop, or we have passed the endpoint and are headed back. |
| * In the second case, we simply evaluate the slope at the |
| * endpoint itself and place ourselves on the appropriate side |
| * of it or on it depending on that result. |
| */ |
| private static void fixRoots(double res[], double eqn[]) { |
| final double EPSILON = 1E-5; |
| for (int i = 0; i < 3; i++) { |
| double t = res[i]; |
| if (Math.abs(t) < EPSILON) { |
| res[i] = findZero(t, 0, eqn); |
| } else if (Math.abs(t - 1) < EPSILON) { |
| res[i] = findZero(t, 1, eqn); |
| } |
| } |
| } |
| |
| private static double solveEqn(double eqn[], int order, double t) { |
| double v = eqn[order]; |
| while (--order >= 0) { |
| v = v * t + eqn[order]; |
| } |
| return v; |
| } |
| |
| private static double findZero(double t, double target, double eqn[]) { |
| double slopeqn[] = {eqn[1], 2*eqn[2], 3*eqn[3]}; |
| double slope; |
| double origdelta = 0; |
| double origt = t; |
| while (true) { |
| slope = solveEqn(slopeqn, 2, t); |
| if (slope == 0) { |
| // At a local minima - must return |
| return t; |
| } |
| double y = solveEqn(eqn, 3, t); |
| if (y == 0) { |
| // Found it! - return it |
| return t; |
| } |
| // assert(slope != 0 && y != 0); |
| double delta = - (y / slope); |
| // assert(delta != 0); |
| if (origdelta == 0) { |
| origdelta = delta; |
| } |
| if (t < target) { |
| if (delta < 0) return t; |
| } else if (t > target) { |
| if (delta > 0) return t; |
| } else { /* t == target */ |
| return (delta > 0 |
| ? (target + java.lang.Double.MIN_VALUE) |
| : (target - java.lang.Double.MIN_VALUE)); |
| } |
| double newt = t + delta; |
| if (t == newt) { |
| // The deltas are so small that we aren't moving... |
| return t; |
| } |
| if (delta * origdelta < 0) { |
| // We have reversed our path. |
| int tag = (origt < t |
| ? getTag(target, origt, t) |
| : getTag(target, t, origt)); |
| if (tag != INSIDE) { |
| // Local minima found away from target - return the middle |
| return (origt + t) / 2; |
| } |
| // Local minima somewhere near target - move to target |
| // and let the slope determine the resulting t. |
| t = target; |
| } else { |
| t = newt; |
| } |
| } |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean contains(double x, double y) { |
| if (!(x * 0.0 + y * 0.0 == 0.0)) { |
| /* Either x or y was infinite or NaN. |
| * A NaN always produces a negative response to any test |
| * and Infinity values cannot be "inside" any path so |
| * they should return false as well. |
| */ |
| return false; |
| } |
| // We count the "Y" crossings to determine if the point is |
| // inside the curve bounded by its closing line. |
| double x1 = getX1(); |
| double y1 = getY1(); |
| double x2 = getX2(); |
| double y2 = getY2(); |
| int crossings = |
| (Curve.pointCrossingsForLine(x, y, x1, y1, x2, y2) + |
| Curve.pointCrossingsForCubic(x, y, |
| x1, y1, |
| getCtrlX1(), getCtrlY1(), |
| getCtrlX2(), getCtrlY2(), |
| x2, y2, 0)); |
| return ((crossings & 1) == 1); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean contains(Point2D p) { |
| return contains(p.getX(), p.getY()); |
| } |
| |
| /* |
| * Fill an array with the coefficients of the parametric equation |
| * in t, ready for solving against val with solveCubic. |
| * We currently have: |
| * <pre> |
| * val = P(t) = C1(1-t)^3 + 3CP1 t(1-t)^2 + 3CP2 t^2(1-t) + C2 t^3 |
| * = C1 - 3C1t + 3C1t^2 - C1t^3 + |
| * 3CP1t - 6CP1t^2 + 3CP1t^3 + |
| * 3CP2t^2 - 3CP2t^3 + |
| * C2t^3 |
| * 0 = (C1 - val) + |
| * (3CP1 - 3C1) t + |
| * (3C1 - 6CP1 + 3CP2) t^2 + |
| * (C2 - 3CP2 + 3CP1 - C1) t^3 |
| * 0 = C + Bt + At^2 + Dt^3 |
| * C = C1 - val |
| * B = 3*CP1 - 3*C1 |
| * A = 3*CP2 - 6*CP1 + 3*C1 |
| * D = C2 - 3*CP2 + 3*CP1 - C1 |
| * </pre> |
| */ |
| private static void fillEqn(double eqn[], double val, |
| double c1, double cp1, double cp2, double c2) { |
| eqn[0] = c1 - val; |
| eqn[1] = (cp1 - c1) * 3.0; |
| eqn[2] = (cp2 - cp1 - cp1 + c1) * 3.0; |
| eqn[3] = c2 + (cp1 - cp2) * 3.0 - c1; |
| return; |
| } |
| |
| /* |
| * Evaluate the t values in the first num slots of the vals[] array |
| * and place the evaluated values back into the same array. Only |
| * evaluate t values that are within the range <0, 1>, including |
| * the 0 and 1 ends of the range iff the include0 or include1 |
| * booleans are true. If an "inflection" equation is handed in, |
| * then any points which represent a point of inflection for that |
| * cubic equation are also ignored. |
| */ |
| private static int evalCubic(double vals[], int num, |
| boolean include0, |
| boolean include1, |
| double inflect[], |
| double c1, double cp1, |
| double cp2, double c2) { |
| int j = 0; |
| for (int i = 0; i < num; i++) { |
| double t = vals[i]; |
| if ((include0 ? t >= 0 : t > 0) && |
| (include1 ? t <= 1 : t < 1) && |
| (inflect == null || |
| inflect[1] + (2*inflect[2] + 3*inflect[3]*t)*t != 0)) |
| { |
| double u = 1 - t; |
| vals[j++] = c1*u*u*u + 3*cp1*t*u*u + 3*cp2*t*t*u + c2*t*t*t; |
| } |
| } |
| return j; |
| } |
| |
| private static final int BELOW = -2; |
| private static final int LOWEDGE = -1; |
| private static final int INSIDE = 0; |
| private static final int HIGHEDGE = 1; |
| private static final int ABOVE = 2; |
| |
| /* |
| * Determine where coord lies with respect to the range from |
| * low to high. It is assumed that low <= high. The return |
| * value is one of the 5 values BELOW, LOWEDGE, INSIDE, HIGHEDGE, |
| * or ABOVE. |
| */ |
| private static int getTag(double coord, double low, double high) { |
| if (coord <= low) { |
| return (coord < low ? BELOW : LOWEDGE); |
| } |
| if (coord >= high) { |
| return (coord > high ? ABOVE : HIGHEDGE); |
| } |
| return INSIDE; |
| } |
| |
| /* |
| * Determine if the pttag represents a coordinate that is already |
| * in its test range, or is on the border with either of the two |
| * opttags representing another coordinate that is "towards the |
| * inside" of that test range. In other words, are either of the |
| * two "opt" points "drawing the pt inward"? |
| */ |
| private static boolean inwards(int pttag, int opt1tag, int opt2tag) { |
| switch (pttag) { |
| case BELOW: |
| case ABOVE: |
| default: |
| return false; |
| case LOWEDGE: |
| return (opt1tag >= INSIDE || opt2tag >= INSIDE); |
| case INSIDE: |
| return true; |
| case HIGHEDGE: |
| return (opt1tag <= INSIDE || opt2tag <= INSIDE); |
| } |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean intersects(double x, double y, double w, double h) { |
| // Trivially reject non-existant rectangles |
| if (w <= 0 || h <= 0) { |
| return false; |
| } |
| |
| // Trivially accept if either endpoint is inside the rectangle |
| // (not on its border since it may end there and not go inside) |
| // Record where they lie with respect to the rectangle. |
| // -1 => left, 0 => inside, 1 => right |
| double x1 = getX1(); |
| double y1 = getY1(); |
| int x1tag = getTag(x1, x, x+w); |
| int y1tag = getTag(y1, y, y+h); |
| if (x1tag == INSIDE && y1tag == INSIDE) { |
| return true; |
| } |
| double x2 = getX2(); |
| double y2 = getY2(); |
| int x2tag = getTag(x2, x, x+w); |
| int y2tag = getTag(y2, y, y+h); |
| if (x2tag == INSIDE && y2tag == INSIDE) { |
| return true; |
| } |
| |
| double ctrlx1 = getCtrlX1(); |
| double ctrly1 = getCtrlY1(); |
| double ctrlx2 = getCtrlX2(); |
| double ctrly2 = getCtrlY2(); |
| int ctrlx1tag = getTag(ctrlx1, x, x+w); |
| int ctrly1tag = getTag(ctrly1, y, y+h); |
| int ctrlx2tag = getTag(ctrlx2, x, x+w); |
| int ctrly2tag = getTag(ctrly2, y, y+h); |
| |
| // Trivially reject if all points are entirely to one side of |
| // the rectangle. |
| if (x1tag < INSIDE && x2tag < INSIDE && |
| ctrlx1tag < INSIDE && ctrlx2tag < INSIDE) |
| { |
| return false; // All points left |
| } |
| if (y1tag < INSIDE && y2tag < INSIDE && |
| ctrly1tag < INSIDE && ctrly2tag < INSIDE) |
| { |
| return false; // All points above |
| } |
| if (x1tag > INSIDE && x2tag > INSIDE && |
| ctrlx1tag > INSIDE && ctrlx2tag > INSIDE) |
| { |
| return false; // All points right |
| } |
| if (y1tag > INSIDE && y2tag > INSIDE && |
| ctrly1tag > INSIDE && ctrly2tag > INSIDE) |
| { |
| return false; // All points below |
| } |
| |
| // Test for endpoints on the edge where either the segment |
| // or the curve is headed "inwards" from them |
| // Note: These tests are a superset of the fast endpoint tests |
| // above and thus repeat those tests, but take more time |
| // and cover more cases |
| if (inwards(x1tag, x2tag, ctrlx1tag) && |
| inwards(y1tag, y2tag, ctrly1tag)) |
| { |
| // First endpoint on border with either edge moving inside |
| return true; |
| } |
| if (inwards(x2tag, x1tag, ctrlx2tag) && |
| inwards(y2tag, y1tag, ctrly2tag)) |
| { |
| // Second endpoint on border with either edge moving inside |
| return true; |
| } |
| |
| // Trivially accept if endpoints span directly across the rectangle |
| boolean xoverlap = (x1tag * x2tag <= 0); |
| boolean yoverlap = (y1tag * y2tag <= 0); |
| if (x1tag == INSIDE && x2tag == INSIDE && yoverlap) { |
| return true; |
| } |
| if (y1tag == INSIDE && y2tag == INSIDE && xoverlap) { |
| return true; |
| } |
| |
| // We now know that both endpoints are outside the rectangle |
| // but the 4 points are not all on one side of the rectangle. |
| // Therefore the curve cannot be contained inside the rectangle, |
| // but the rectangle might be contained inside the curve, or |
| // the curve might intersect the boundary of the rectangle. |
| |
| double[] eqn = new double[4]; |
| double[] res = new double[4]; |
| if (!yoverlap) { |
| // Both y coordinates for the closing segment are above or |
| // below the rectangle which means that we can only intersect |
| // if the curve crosses the top (or bottom) of the rectangle |
| // in more than one place and if those crossing locations |
| // span the horizontal range of the rectangle. |
| fillEqn(eqn, (y1tag < INSIDE ? y : y+h), y1, ctrly1, ctrly2, y2); |
| int num = solveCubic(eqn, res); |
| num = evalCubic(res, num, true, true, null, |
| x1, ctrlx1, ctrlx2, x2); |
| // odd counts imply the crossing was out of [0,1] bounds |
| // otherwise there is no way for that part of the curve to |
| // "return" to meet its endpoint |
| return (num == 2 && |
| getTag(res[0], x, x+w) * getTag(res[1], x, x+w) <= 0); |
| } |
| |
| // Y ranges overlap. Now we examine the X ranges |
| if (!xoverlap) { |
| // Both x coordinates for the closing segment are left of |
| // or right of the rectangle which means that we can only |
| // intersect if the curve crosses the left (or right) edge |
| // of the rectangle in more than one place and if those |
| // crossing locations span the vertical range of the rectangle. |
| fillEqn(eqn, (x1tag < INSIDE ? x : x+w), x1, ctrlx1, ctrlx2, x2); |
| int num = solveCubic(eqn, res); |
| num = evalCubic(res, num, true, true, null, |
| y1, ctrly1, ctrly2, y2); |
| // odd counts imply the crossing was out of [0,1] bounds |
| // otherwise there is no way for that part of the curve to |
| // "return" to meet its endpoint |
| return (num == 2 && |
| getTag(res[0], y, y+h) * getTag(res[1], y, y+h) <= 0); |
| } |
| |
| // The X and Y ranges of the endpoints overlap the X and Y |
| // ranges of the rectangle, now find out how the endpoint |
| // line segment intersects the Y range of the rectangle |
| double dx = x2 - x1; |
| double dy = y2 - y1; |
| double k = y2 * x1 - x2 * y1; |
| int c1tag, c2tag; |
| if (y1tag == INSIDE) { |
| c1tag = x1tag; |
| } else { |
| c1tag = getTag((k + dx * (y1tag < INSIDE ? y : y+h)) / dy, x, x+w); |
| } |
| if (y2tag == INSIDE) { |
| c2tag = x2tag; |
| } else { |
| c2tag = getTag((k + dx * (y2tag < INSIDE ? y : y+h)) / dy, x, x+w); |
| } |
| // If the part of the line segment that intersects the Y range |
| // of the rectangle crosses it horizontally - trivially accept |
| if (c1tag * c2tag <= 0) { |
| return true; |
| } |
| |
| // Now we know that both the X and Y ranges intersect and that |
| // the endpoint line segment does not directly cross the rectangle. |
| // |
| // We can almost treat this case like one of the cases above |
| // where both endpoints are to one side, except that we may |
| // get one or three intersections of the curve with the vertical |
| // side of the rectangle. This is because the endpoint segment |
| // accounts for the other intersection in an even pairing. Thus, |
| // with the endpoint crossing we end up with 2 or 4 total crossings. |
| // |
| // (Remember there is overlap in both the X and Y ranges which |
| // means that the segment itself must cross at least one vertical |
| // edge of the rectangle - in particular, the "near vertical side" |
| // - leaving an odd number of intersections for the curve.) |
| // |
| // Now we calculate the y tags of all the intersections on the |
| // "near vertical side" of the rectangle. We will have one with |
| // the endpoint segment, and one or three with the curve. If |
| // any pair of those vertical intersections overlap the Y range |
| // of the rectangle, we have an intersection. Otherwise, we don't. |
| |
| // c1tag = vertical intersection class of the endpoint segment |
| // |
| // Choose the y tag of the endpoint that was not on the same |
| // side of the rectangle as the subsegment calculated above. |
| // Note that we can "steal" the existing Y tag of that endpoint |
| // since it will be provably the same as the vertical intersection. |
| c1tag = ((c1tag * x1tag <= 0) ? y1tag : y2tag); |
| |
| // Now we have to calculate an array of solutions of the curve |
| // with the "near vertical side" of the rectangle. Then we |
| // need to sort the tags and do a pairwise range test to see |
| // if either of the pairs of crossings spans the Y range of |
| // the rectangle. |
| // |
| // Note that the c2tag can still tell us which vertical edge |
| // to test against. |
| fillEqn(eqn, (c2tag < INSIDE ? x : x+w), x1, ctrlx1, ctrlx2, x2); |
| int num = solveCubic(eqn, res); |
| num = evalCubic(res, num, true, true, null, y1, ctrly1, ctrly2, y2); |
| |
| // Now put all of the tags into a bucket and sort them. There |
| // is an intersection iff one of the pairs of tags "spans" the |
| // Y range of the rectangle. |
| int tags[] = new int[num+1]; |
| for (int i = 0; i < num; i++) { |
| tags[i] = getTag(res[i], y, y+h); |
| } |
| tags[num] = c1tag; |
| Arrays.sort(tags); |
| return ((num >= 1 && tags[0] * tags[1] <= 0) || |
| (num >= 3 && tags[2] * tags[3] <= 0)); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean intersects(Rectangle2D r) { |
| return intersects(r.getX(), r.getY(), r.getWidth(), r.getHeight()); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean contains(double x, double y, double w, double h) { |
| if (w <= 0 || h <= 0) { |
| return false; |
| } |
| // Assertion: Cubic curves closed by connecting their |
| // endpoints form either one or two convex halves with |
| // the closing line segment as an edge of both sides. |
| if (!(contains(x, y) && |
| contains(x + w, y) && |
| contains(x + w, y + h) && |
| contains(x, y + h))) { |
| return false; |
| } |
| // Either the rectangle is entirely inside one of the convex |
| // halves or it crosses from one to the other, in which case |
| // it must intersect the closing line segment. |
| Rectangle2D rect = new Rectangle2D.Double(x, y, w, h); |
| return !rect.intersectsLine(getX1(), getY1(), getX2(), getY2()); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public boolean contains(Rectangle2D r) { |
| return contains(r.getX(), r.getY(), r.getWidth(), r.getHeight()); |
| } |
| |
| /** |
| * {@inheritDoc} |
| * @since 1.2 |
| */ |
| public Rectangle getBounds() { |
| return getBounds2D().getBounds(); |
| } |
| |
| /** |
| * Returns an iteration object that defines the boundary of the |
| * shape. |
| * The iterator for this class is not multi-threaded safe, |
| * which means that this <code>CubicCurve2D</code> class does not |
| * guarantee that modifications to the geometry of this |
| * <code>CubicCurve2D</code> object do not affect any iterations of |
| * that geometry that are already in process. |
| * @param at an optional <code>AffineTransform</code> to be applied to the |
| * coordinates as they are returned in the iteration, or <code>null</code> |
| * if untransformed coordinates are desired |
| * @return the <code>PathIterator</code> object that returns the |
| * geometry of the outline of this <code>CubicCurve2D</code>, one |
| * segment at a time. |
| * @since 1.2 |
| */ |
| public PathIterator getPathIterator(AffineTransform at) { |
| return new CubicIterator(this, at); |
| } |
| |
| /** |
| * Return an iteration object that defines the boundary of the |
| * flattened shape. |
| * The iterator for this class is not multi-threaded safe, |
| * which means that this <code>CubicCurve2D</code> class does not |
| * guarantee that modifications to the geometry of this |
| * <code>CubicCurve2D</code> object do not affect any iterations of |
| * that geometry that are already in process. |
| * @param at an optional <code>AffineTransform</code> to be applied to the |
| * coordinates as they are returned in the iteration, or <code>null</code> |
| * if untransformed coordinates are desired |
| * @param flatness the maximum amount that the control points |
| * for a given curve can vary from colinear before a subdivided |
| * curve is replaced by a straight line connecting the end points |
| * @return the <code>PathIterator</code> object that returns the |
| * geometry of the outline of this <code>CubicCurve2D</code>, |
| * one segment at a time. |
| * @since 1.2 |
| */ |
| public PathIterator getPathIterator(AffineTransform at, double flatness) { |
| return new FlatteningPathIterator(getPathIterator(at), flatness); |
| } |
| |
| /** |
| * Creates a new object of the same class as this object. |
| * |
| * @return a clone of this instance. |
| * @exception OutOfMemoryError if there is not enough memory. |
| * @see java.lang.Cloneable |
| * @since 1.2 |
| */ |
| public Object clone() { |
| try { |
| return super.clone(); |
| } catch (CloneNotSupportedException e) { |
| // this shouldn't happen, since we are Cloneable |
| throw new InternalError(); |
| } |
| } |
| } |