Geometry.java

  1. /*
  2.  * Licensed to the Hipparchus project under one or more
  3.  * contributor license agreements.  See the NOTICE file distributed with
  4.  * this work for additional information regarding copyright ownership.
  5.  * The Hipparchus project licenses this file to You under the Apache License, Version 2.0
  6.  * (the "License"); you may not use this file except in compliance with
  7.  * the License.  You may obtain a copy of the License at
  8.  *
  9.  *      https://www.apache.org/licenses/LICENSE-2.0
  10.  *
  11.  * Unless required by applicable law or agreed to in writing, software
  12.  * distributed under the License is distributed on an "AS IS" BASIS,
  13.  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  14.  * See the License for the specific language governing permissions and
  15.  * limitations under the License.
  16.  */
  17. package org.hipparchus.geometry;

  18. import org.hipparchus.exception.LocalizedCoreFormats;
  19. import org.hipparchus.exception.MathIllegalArgumentException;

  20. import java.util.List;

  21. /** Utilities for geometry.
  22.  * @since 4.0
  23.  */
  24. public class Geometry {

  25.     /**
  26.      * Private constructor for a utility class.
  27.      */
  28.     private Geometry() {
  29.         // nothing to do
  30.     }

  31.     /**
  32.      * Compute the barycenter of n points.
  33.      * @param <S> Type of the space.
  34.      * @param <P> Type of the points in space.
  35.      * @param points points generating the barycenter
  36.      * @return barycenter of the points
  37.      */
  38.     public static <S extends Space, P extends Point<S, P>> P barycenter(final List<P> points) {

  39.         // safety check
  40.         if (points.isEmpty()) {
  41.             throw new MathIllegalArgumentException(LocalizedCoreFormats.NUMBER_OF_ELEMENTS_SHOULD_BE_POSITIVE, 0);
  42.         }

  43.         // compute barycenter by moving from point to point
  44.         P current = points.get(0);
  45.         for (int i = 1; i < points.size(); i++) {
  46.             current = points.get(i).moveTowards(current, ((double) i) / (i + 1));
  47.         }

  48.         return current;

  49.     }

  50. }