|
Java example source code file (WelzlEncloser3DTest.java)
This example Java source code file (WelzlEncloser3DTest.java) is included in the alvinalexander.com
"Java Source Code
Warehouse" project. The intent of this project is to help you "Learn
Java by Example" TM.
Learn more about this Java project at its project page.
The WelzlEncloser3DTest.java Java example source code
/*
* 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.apache.commons.math3.geometry.enclosing;
import java.io.IOException;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.List;
import org.apache.commons.math3.geometry.euclidean.threed.Euclidean3D;
import org.apache.commons.math3.geometry.euclidean.threed.SphereGenerator;
import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
import org.apache.commons.math3.random.RandomGenerator;
import org.apache.commons.math3.random.UnitSphereRandomVectorGenerator;
import org.apache.commons.math3.random.Well1024a;
import org.junit.Assert;
import org.junit.Test;
public class WelzlEncloser3DTest {
@Test
public void testNullList() {
SphereGenerator generator = new SphereGenerator();
WelzlEncloser<Euclidean3D, Vector3D> encloser =
new WelzlEncloser<Euclidean3D, Vector3D>(1.0e-10, generator);
EnclosingBall<Euclidean3D, Vector3D> ball = encloser.enclose(null);
Assert.assertTrue(ball.getRadius() < 0);
}
@Test
public void testNoPoints() {
SphereGenerator generator = new SphereGenerator();
WelzlEncloser<Euclidean3D, Vector3D> encloser =
new WelzlEncloser<Euclidean3D, Vector3D>(1.0e-10, generator);
EnclosingBall<Euclidean3D, Vector3D> ball = encloser.enclose(new ArrayList());
Assert.assertTrue(ball.getRadius() < 0);
}
@Test
public void testReducingBall() {
List<Vector3D> list =
Arrays.asList(new Vector3D(-7.140397329568118, -16.571661242582177, 11.714458961735405),
new Vector3D(-7.137986707455888, -16.570767323375720, 11.708602108715928),
new Vector3D(-7.139185068549035, -16.570891204702250, 11.715554057357394),
new Vector3D(-7.142682716997507, -16.571609818234290, 11.710787934580328),
new Vector3D(-7.139018392423351, -16.574405614157020, 11.710518716711425),
new Vector3D(-7.140870659936730, -16.567993074240455, 11.710914678204503),
new Vector3D(-7.136350173659562, -16.570498228820930, 11.713965225900928),
new Vector3D(-7.141675762759172, -16.572852471407028, 11.714033471449508),
new Vector3D(-7.140453077221105, -16.570212820780647, 11.708624578004980),
new Vector3D(-7.140322188726825, -16.574152894557717, 11.710305611121410),
new Vector3D(-7.141116131477088, -16.574061164624560, 11.712938509321699));
WelzlEncloser<Euclidean3D, Vector3D> encloser =
new WelzlEncloser<Euclidean3D, Vector3D>(1.0e-10, new SphereGenerator());
EnclosingBall<Euclidean3D, Vector3D> ball = encloser.enclose(list);
Assert.assertTrue(ball.getRadius() > 0);
}
@Test
public void testInfiniteLoop() {
// this test used to generate an infinite loop
List<Vector3D> list =
Arrays.asList(new Vector3D( -0.89227075512164380, -2.89317694645713900, 14.84572323743355500),
new Vector3D( -0.92099498940693580, -2.31086108263908940, 12.92071026467688300),
new Vector3D( -0.85227999411005200, -3.06314731441320730, 15.40163831651287000),
new Vector3D( -1.77399413020785970, -3.65630391378114260, 14.13190097751873400),
new Vector3D( 0.33157833272465354, -2.22813591757792160, 14.21225234159008200),
new Vector3D( -1.53065579165484400, -1.65692084770139570, 14.61483055714788500),
new Vector3D( -1.08457093941217140, -1.96100325935602980, 13.09265170575555000),
new Vector3D( 0.30029469589708850, -3.05470831395667370, 14.56352400426342600),
new Vector3D( -0.95007443938638460, -1.86810946486118360, 15.14491234340057000),
new Vector3D( -1.89661503804130830, -2.17004080885185860, 14.81235128513927000),
new Vector3D( -0.72193328761607530, -1.44513142833618270, 14.52355724218561800),
new Vector3D( -0.26895980939606550, -3.69512371522084140, 14.72272846327652000),
new Vector3D( -1.53501693431786170, -3.25055166611021900, 15.15509062584274800),
new Vector3D( -0.71727553535519410, -3.62284279460799100, 13.26256700929380700),
new Vector3D( -0.30220950676137365, -3.25410412500779070, 13.13682612771606000),
new Vector3D( -0.04543996608267075, -1.93081853923797750, 14.79497997883171400),
new Vector3D( -1.53348892951571640, -3.66688919703524900, 14.73095600812074200),
new Vector3D( -0.98034899533935820, -3.34004481162763960, 13.03245014017556800));
WelzlEncloser<Euclidean3D, Vector3D> encloser =
new WelzlEncloser<Euclidean3D, Vector3D>(1.0e-10, new SphereGenerator());
EnclosingBall<Euclidean3D, Vector3D> ball = encloser.enclose(list);
Assert.assertTrue(ball.getRadius() > 0);
}
@Test
public void testLargeSamples() throws IOException {
RandomGenerator random = new Well1024a(0x35ddecfc78131e1dl);
final UnitSphereRandomVectorGenerator sr = new UnitSphereRandomVectorGenerator(3, random);
for (int k = 0; k < 50; ++k) {
// define the reference sphere we want to compute
double d = 25 * random.nextDouble();
double refRadius = 10 * random.nextDouble();
Vector3D refCenter = new Vector3D(d, new Vector3D(sr.nextVector()));
// set up a large sample inside the reference sphere
int nbPoints = random.nextInt(1000);
List<Vector3D> points = new ArrayList();
for (int i = 0; i < nbPoints; ++i) {
double r = refRadius * random.nextDouble();
points.add(new Vector3D(1.0, refCenter, r, new Vector3D(sr.nextVector())));
}
// test we find a sphere at most as large as the one used for random drawings
checkSphere(points, refRadius);
}
}
private void checkSphere(List<Vector3D> points, double refRadius) {
EnclosingBall<Euclidean3D, Vector3D> sphere = checkSphere(points);
// compare computed sphere with bounding sphere
Assert.assertTrue(sphere.getRadius() <= refRadius);
// check removing any point of the support Sphere fails to enclose the point
for (int i = 0; i < sphere.getSupportSize(); ++i) {
List<Vector3D> reducedSupport = new ArrayList();
int count = 0;
for (Vector3D s : sphere.getSupport()) {
if (count++ != i) {
reducedSupport.add(s);
}
}
EnclosingBall<Euclidean3D, Vector3D> reducedSphere =
new SphereGenerator().ballOnSupport(reducedSupport);
boolean foundOutside = false;
for (int j = 0; j < points.size() && !foundOutside; ++j) {
if (!reducedSphere.contains(points.get(j), 1.0e-10)) {
foundOutside = true;
}
}
Assert.assertTrue(foundOutside);
}
}
private EnclosingBall<Euclidean3D, Vector3D> checkSphere(List points) {
WelzlEncloser<Euclidean3D, Vector3D> encloser =
new WelzlEncloser<Euclidean3D, Vector3D>(1.0e-10, new SphereGenerator());
EnclosingBall<Euclidean3D, Vector3D> Sphere = encloser.enclose(points);
// all points are enclosed
for (Vector3D v : points) {
Assert.assertTrue(Sphere.contains(v, 1.0e-10));
}
for (Vector3D v : points) {
boolean inSupport = false;
for (Vector3D s : Sphere.getSupport()) {
if (v == s) {
inSupport = true;
}
}
if (inSupport) {
// points on the support should be outside of reduced ball
Assert.assertFalse(Sphere.contains(v, -0.001));
}
}
return Sphere;
}
}
Other Java examples (source code examples)
Here is a short list of links related to this Java WelzlEncloser3DTest.java source code file:
|