hoy.h

/*****************************************************************************
 *   Copyright (C) 2004-2015 The PaGMO development team,                     *
 *   Advanced Concepts Team (ACT), European Space Agency (ESA)               *
 *                                                                           *
 *   https://github.com/esa/pagmo                                            *
 *                                                                           *
 *   act@esa.int                                                             *
 *                                                                           *
 *   This program is free software; you can redistribute it and/or modify    *
 *   it under the terms of the GNU General Public License as published by    *
 *   the Free Software Foundation; either version 2 of the License, or       *
 *   (at your option) any later version.                                     *
 *                                                                           *
 *   This program 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 for more details.                            *
 *                                                                           *
 *   You should have received a copy of the GNU General Public License       *
 *   along with this program; if not, write to the                           *
 *   Free Software Foundation, Inc.,                                         *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.               *
 *****************************************************************************/

#ifndef PAGMO_UTIL_HV_ALGORITHM_HOY_H
#define PAGMO_UTIL_HV_ALGORITHM_HOY_H

#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <iterator>

#include "base.h"
#include "../hypervolume.h"

namespace pagmo { namespace util { namespace hv_algorithm {


class __PAGMO_VISIBLE hoy : public base {
public:
        hoy();
        double compute(std::vector<fitness_vector> &, const fitness_vector &) const;
        void verify_before_compute(const std::vector<fitness_vector> &, const fitness_vector &) const;
        base_ptr clone() const;
        std::string get_name() const;

private:
        inline bool covers(const double cub[], const double reg_low[]) const;
        inline bool part_covers(const double cub[], const double reg_up[]) const;
        inline int contains_boundary(const double cub[], const double reg_low[], const int split) const;
        inline double get_measure(const double reg_low[], const double reg_up[]) const;
        inline int is_pile(const double cub[], const double reg_low[]) const;
        inline double compute_trellis(const double reg_low[], const double reg_up[], const double trellis[]) const;
        inline double get_median(double* bounds, unsigned int n) const;
        inline void stream(double m_region_low[], double m_region_up[], double** points, const unsigned int n_points, int split, double cover, unsigned int rec_level) const;

        mutable int             m_dimension;
        mutable int             m_total_size;
        mutable double  m_sqrt_size;
        mutable double  m_volume;
        mutable double  *m_region_up;
        mutable double  *m_region_low;
        mutable int             *m_piles;
        mutable double  *m_trellis;
        mutable double  *m_boundaries;
        mutable double  *m_no_boundaries;
        mutable std::vector<double**> m_child_points;

        friend class boost::serialization::access;
        template <class Archive>
        void serialize(Archive &ar, const unsigned int)
        {
                ar & boost::serialization::base_object<base>(*this);
        }
};

} } }

BOOST_CLASS_EXPORT_KEY(pagmo::util::hv_algorithm::hoy)

#endif