sghyperplane.hpp

// This file is part of the SGSolve library for stochastic games
// Copyright (C) 2019 Benjamin A. Brooks
// 
// SGSolve 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 3 of the License, or
// (at your option) any later version.
// 
// SGSolve 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, see
// <http://www.gnu.org/licenses/>.
// 
// Benjamin A. Brooks
// ben@benjaminbrooks.net
// Chicago, IL
 
#ifndef _SGHYPERPLANE_HPP
#define _SGHYPERPLANE_HPP
 
#include "sgcommon.hpp"
#include "sgpoint.hpp"
#include "sgexception.hpp"
 
 
class SGHyperplane
{
private:
  SGPoint normal; 
  vector<double> levels; 
public:
  SGHyperplane(int numStates): levels(numStates,0)
  {}
 
  SGHyperplane(const SGPoint & _normal,
           const vector<double> & _levels):
    normal(_normal),
    levels(_levels)
  {}
 
  SGHyperplane(): normal(0,0), levels(1,0) {}
  ~SGHyperplane() {}
 
 
  double expectation(const vector<double> & prob) const;
 
  static double distance(const SGHyperplane& t0,
             const SGHyperplane& t1);
  
  void push_back(double level);
  void clear();
 
  int size() const { return levels.size(); }
 
  void setNormal(const SGPoint & newNorm)
  { normal = newNorm; }
  const SGPoint & getNormal() const
  { return normal; }
  const vector<double> & getLevels() const
  { return levels; }
  
  // Operators
  double& operator[](int state);
  const double& operator[](int state) const;
  SGHyperplane& operator=(const SGHyperplane & rhs);
  SGHyperplane& operator+=(const SGHyperplane & rhs);
  SGHyperplane& operator-=(const SGHyperplane & rhs);
 
  SGHyperplane& operator+=(double rhs);
  SGHyperplane& operator*=(double d);
  friend SGHyperplane operator*(double d,const SGHyperplane & hp);
  friend SGHyperplane operator*(const SGHyperplane & hp,double d);
  friend ostream& operator<<(ostream& out, const SGHyperplane& rhs);
 
  void round(double tol);
 
  template<class Archive>
  void serialize(Archive &ar, const unsigned int version)
  {
    ar & normal & levels;
  }
 
  friend class boost::serialization::access;
}; // SGHyperplane
 
#endif