Constraint.cc

Go to the documentation of this file.

/*
 * $Id: Constraint.cc 239 2014-11-30 14:11:47Z jdl3 $
 * Copyright (C) 2013, 2014 John D Lamb
 *
 * 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., 675 Mass Ave, Cambridge, MA 02139, USA.
 */

#ifdef HAVE_CONFIG_H
#  include<config.h>
#endif

#include"Constraint.hpp"

using namespace ipo;

Constraint::Constraint( detail::ModelBase& model, detail::SharedFunctionPtr function,
            std::string const& name )
  : Objective { model, function, name }{}

Constraint::Constraint( detail::ModelBase& model, ipo_function::Function* function,
            std::string const& name )
  : Constraint { model, detail::SharedFunctionPtr { function }, name }{}

Constraint::Constraint( detail::ModelBase& model, ipo_function::Function& function,
            std::string const& name )
  : Constraint { model, detail::SharedFunctionPtr { function }, name }{}

Constraint::Constraint( detail::ModelBase& model, detail::SharedFunctionPtr function,
            char const* const name )
  : Constraint { model, function, 0 == name ? std::string {} : name }{}

Constraint::Constraint( detail::ModelBase& model, ipo_function::Function* function,
            char const* const name )
  : Constraint { model, function, 0 == name ? std::string {} : name }{}

Constraint::Constraint( detail::ModelBase& model, ipo_function::Function& function,
            char const* const name )
  : Constraint { model, function, 0 == name ? std::string {} : name }{}

double
Constraint::getUpperBound() const {
  return data->upperBound;
}

double
Constraint::getLowerBound() const {
  return data->lowerBound;
}

void
Constraint::setUpperBound( double const upperBound ){
  if( upperBound <= data->lowerBound )
    throw IPOE( "ipo::Constraint::setUpperBound(): upper bound must exceed lower." );
  data->upperBound = upperBound;
}

void
Constraint::setLowerBound( double const lowerBound ){
  if( data->upperBound <= lowerBound )
    throw IPOE( "ipo::Constraint::setLowerBound(): upper bound must exceed lower." );
  data->lowerBound = lowerBound;
}

bool
ipo::operator==( Constraint const& lhs, Constraint const& rhs ){
  return lhs.data == rhs.data;
}


Constraint::Constraint( Constraint const& constraint )
  : Objective { constraint }{}

Constraint&
Constraint::operator=( Constraint const& constraint ){
  model = constraint.model;
  function = constraint.function;
  data = std::move( constraint.data );
  return *this;
}

void
Constraint::summary( std::ostream& ostream, std::string const& prefix ) const {
  auto f = [](double x)->std::string {
    std::ostringstream ost;
    double a { std::fabs( x ) };
    if( x < 0 ) ost << "−";
    if( a == std::numeric_limits<double>::infinity() )
      ost << "∞";
    else
      ost << x;
    return ost.str();
  };
  ostream << prefix << f( getLowerBound() ) << " ≤ " << getName() << " ( ";
  bool first { true };
  for( auto const& e : const_cast<ipo::Constraint&>( *this ).getVariables() ){
    if( first ) first = false; else ostream << ", ";
    ostream << e.getName();
  }
  ostream << " ) ≤ " <<  f( getUpperBound() ) << std::endl;
  std::cout << prefix << "  Constraint function value = " << value() << std::endl;
}