All Downloads are FREE. Search and download functionalities are using the official Maven repository.

org.elasticsearch.geo.XShapeTestUtil Maven / Gradle / Ivy

There is a newer version: 8.16.0
Show newest version
/*
 * @notice
 * Licensed to the Apache Software Foundation (ASF) under one or more
 * contributor license agreements.  See the NOTICE file distributed with
 * this work for additional information regarding copyright ownership.
 * The ASF licenses this file to You under the Apache License, Version 2.0
 * (the "License"); you may not use this file except in compliance with
 * the License.  You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
package org.elasticsearch.geo;

import com.carrotsearch.randomizedtesting.RandomizedContext;
import com.carrotsearch.randomizedtesting.generators.BiasedNumbers;

import org.apache.lucene.geo.XYCircle;
import org.apache.lucene.geo.XYPolygon;
import org.apache.lucene.geo.XYRectangle;
import org.apache.lucene.tests.util.TestUtil;

import java.util.ArrayList;
import java.util.Random;

/** generates random cartesian geometry; heavy reuse of {@link org.apache.lucene.tests.geo.GeoTestUtil} */
public class XShapeTestUtil {

    // If we allow polygons centers to reach anywhere near Float.MAX, we will end up will all points rounding to the same value
    public static final float CENTER_SCALE_FACTOR = 0.001f;
    // When comparing two floats, consider very similar numbers as identical
    private static final float MIN_DIFF = 1e-10f;

    /** returns next pseudorandom polygon */
    public static XYPolygon nextPolygon() {
        if (random().nextBoolean()) {
            return surpriseMePolygon();
        } else if (random().nextInt(10) == 1) {
            // this poly is slow to create ... only do it 10% of the time:
            while (true) {
                int gons = TestUtil.nextInt(random(), 4, 500);
                try {
                    return createRegularPolygon(nextDouble() * CENTER_SCALE_FACTOR, nextDouble() * CENTER_SCALE_FACTOR, nextRadius(), gons);
                } catch (IllegalArgumentException iae) {
                    // we tried to cross dateline or pole ... try again
                }
            }
        }

        XYRectangle box = nextBoxInternal();
        if (random().nextBoolean()) {
            // box
            return boxPolygon(box);
        } else {
            // triangle
            return trianglePolygon(box);
        }
    }

    private static double nextRadius() {
        return random().nextDouble() * CENTER_SCALE_FACTOR * Float.MAX_VALUE + 1.0;
    }

    private static XYPolygon trianglePolygon(XYRectangle box) {
        final float[] polyX = new float[4];
        final float[] polyY = new float[4];
        polyX[0] = box.minX;
        polyY[0] = box.minY;
        polyX[1] = box.minX;
        polyY[1] = box.maxY;
        polyX[2] = box.maxX;
        polyY[2] = box.maxY;
        polyX[3] = box.minX;
        polyY[3] = box.minY;
        return validatePolygon(new XYPolygon(polyX, polyY));
    }

    public static XYRectangle nextBox() {
        return nextBoxInternal();
    }

    private static XYRectangle nextBoxInternal() {
        // prevent lines instead of boxes
        float x0 = nextFloat() * CENTER_SCALE_FACTOR;
        float x1 = nextFloat() * CENTER_SCALE_FACTOR;
        while (Math.abs(x1 - x0) < MIN_DIFF) {
            x1 = nextFloat() * CENTER_SCALE_FACTOR;
        }
        // prevent lines instead of boxes
        float y0 = nextFloat() * CENTER_SCALE_FACTOR;
        float y1 = nextFloat() * CENTER_SCALE_FACTOR;
        while (Math.abs(y1 - y0) < MIN_DIFF) {
            y1 = nextFloat() * CENTER_SCALE_FACTOR;
        }

        if (x1 < x0) {
            float x = x0;
            x0 = x1;
            x1 = x;
        }

        if (y1 < y0) {
            float y = y0;
            y0 = y1;
            y1 = y;
        }

        return new XYRectangle(x0, x1, y0, y1);
    }

    private static XYPolygon boxPolygon(XYRectangle box) {
        final float[] polyX = new float[5];
        final float[] polyY = new float[5];
        polyX[0] = box.minX;
        polyY[0] = box.minY;
        polyX[1] = box.minX;
        polyY[1] = box.maxY;
        polyX[2] = box.maxX;
        polyY[2] = box.maxY;
        polyX[3] = box.maxX;
        polyY[3] = box.minY;
        polyX[4] = box.minX;
        polyY[4] = box.minY;
        return validatePolygon(new XYPolygon(polyX, polyY));
    }

    private static XYPolygon validatePolygon(XYPolygon polygon) {
        if (validPolygon(polygon)) return polygon;
        throw new IllegalArgumentException("Invalid polygon: " + polygon);
    }

    private static XYPolygon surpriseMePolygon() {
        // repeat until we get a poly that doesn't cross dateline:
        while (true) {
            // System.out.println("\nPOLY ITER");
            double centerX = nextDouble() * CENTER_SCALE_FACTOR;
            double centerY = nextDouble() * CENTER_SCALE_FACTOR;
            double radius = nextRadius();
            double radiusDelta = random().nextDouble();

            ArrayList xList = new ArrayList<>();
            ArrayList yList = new ArrayList<>();
            double angle = 0.0;
            while (true) {
                angle += random().nextDouble() * 40.0;
                // System.out.println(" angle " + angle);
                if (angle > 360) {
                    break;
                }
                double len = radius * (1.0 - radiusDelta + radiusDelta * random().nextDouble());
                double maxX = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerX), StrictMath.abs(-Float.MAX_VALUE - centerX));
                double maxY = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerY), StrictMath.abs(-Float.MAX_VALUE - centerY));

                len = StrictMath.min(len, StrictMath.min(maxX, maxY));

                // System.out.println(" len=" + len);
                float x = (float) (centerX + len * Math.cos(Math.toRadians(angle)));
                float y = (float) (centerY + len * Math.sin(Math.toRadians(angle)));

                xList.add(x);
                yList.add(y);

                // System.out.println(" lat=" + lats.get(lats.size()-1) + " lon=" + lons.get(lons.size()-1));
            }

            // close it
            xList.add(xList.get(0));
            yList.add(yList.get(0));

            float[] xArray = new float[xList.size()];
            float[] yArray = new float[yList.size()];
            for (int i = 0; i < xList.size(); i++) {
                xArray[i] = xList.get(i);
                yArray[i] = yList.get(i);
            }
            XYPolygon polygon = new XYPolygon(xArray, yArray);
            if (validPolygon(polygon)) return polygon;
        }
    }

    /** Makes an n-gon, centered at the provided x/y, and each vertex approximately
     *  distanceMeters away from the center.
     *
     * Do not invoke me across the dateline or a pole!! */
    public static XYPolygon createRegularPolygon(double centerX, double centerY, double radius, int gons) {

        double maxX = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerX), StrictMath.abs(-Float.MAX_VALUE - centerX));
        double maxY = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerY), StrictMath.abs(-Float.MAX_VALUE - centerY));

        radius = StrictMath.min(radius, StrictMath.min(maxX, maxY));

        float[][] result = new float[2][];
        result[0] = new float[gons + 1];
        result[1] = new float[gons + 1];
        // System.out.println("make gon=" + gons);
        for (int i = 0; i < gons; i++) {
            double angle = 360.0 - i * (360.0 / gons);
            // System.out.println(" angle " + angle);
            double x = Math.cos(StrictMath.toRadians(angle));
            double y = Math.sin(StrictMath.toRadians(angle));
            result[0][i] = (float) (centerY + y * radius);
            result[1][i] = (float) (centerX + x * radius);
        }

        // close poly
        result[0][gons] = result[0][0];
        result[1][gons] = result[1][0];

        return validatePolygon(new XYPolygon(result[0], result[1]));
    }

    private static boolean validPolygon(XYPolygon polygon) {
        if (sameValues(polygon.getPolyX()) || sameValues(polygon.getPolyY())) {
            return false;
        }
        return true;
    }

    private static boolean sameValues(float[] values) {
        int consecutiveSame = 0;
        for (int i = 1; i < values.length; i++) {
            if (Math.abs(values[i] - values[i - 1]) <= MIN_DIFF) consecutiveSame++;
        }
        return consecutiveSame > values.length / 2;
    }

    public static double nextDouble() {
        return BiasedNumbers.randomDoubleBetween(random(), -Float.MAX_VALUE, Float.MAX_VALUE);
    }

    public static float nextFloat() {
        return BiasedNumbers.randomFloatBetween(random(), -Float.MAX_VALUE, Float.MAX_VALUE);
    }

    /** Keep it simple, we don't need to take arbitrary Random for geo tests */
    private static Random random() {
        return RandomizedContext.current().getRandom();
    }

    public static XYCircle nextCircle() {
        float centerX = nextFloat() * CENTER_SCALE_FACTOR;
        float centerY = nextFloat() * CENTER_SCALE_FACTOR;
        float radius = (float) nextRadius();
        double maxX = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerX), StrictMath.abs(-Float.MAX_VALUE - centerX));
        double maxY = StrictMath.min(StrictMath.abs(Float.MAX_VALUE - centerY), StrictMath.abs(-Float.MAX_VALUE - centerY));
        radius = Math.min(radius, (float) maxX);
        radius = Math.min(radius, (float) maxY);
        return new XYCircle(centerX, centerY, radius);
    }
}




© 2015 - 2025 Weber Informatics LLC | Privacy Policy