PhaseIBoundedVariableFunctionAndDerivatives.cc

Go to the documentation of this file.

/*
 * $Id: PhaseIBoundedVariableFunctionAndDerivatives.cc 236 2014-11-30 14:09:06Z jdl3 $
 * Copyright (C) 2013–14 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"PhaseIBoundedVariableFunctionAndDerivatives.hpp"

using namespace ipo_function::detail;

PhaseIBoundedVariableFunctionAndDerivatives::
PhaseIBoundedVariableFunctionAndDerivatives( double const b,
                         bool const upperBound )
  : FunctionBase { 2 }, function( function ), b { b },
  upperBound { upperBound }
{
  functionGradient = gsl::vector( 2 );
  if( upperBound ) functionGradient[0] = 1.0;
  else functionGradient[0] = -1.0;
  functionGradient[1] = -1.0;
  functionHessian = gsl::matrix { 2, 2 };
  functionHessian.set_all( 0.0 );
}

double
PhaseIBoundedVariableFunctionAndDerivatives::operator()( gsl::vector const& vector ){
  if( vector.size() < 2  ) return 0; // do something sane

  double const& x { vector[0] };
  double const& s { vector[1] };
  if( upperBound )
    return x - b - s;
  else
    return b - x - s;
}

gsl::vector
PhaseIBoundedVariableFunctionAndDerivatives::gradient( gsl::vector const& vector ){
  return functionGradient;
}

gsl::matrix
PhaseIBoundedVariableFunctionAndDerivatives::hessian( gsl::vector const& vector ){
  return functionHessian;
}

void
PhaseIBoundedVariableFunctionAndDerivatives::setVector( gsl::vector const& vector ){
  if( vector.size() < 2  ){ functionValue = 0; return; } // do something sane

  double const& x { vector[0] };
  double const& s { vector[1] };
  if( upperBound )
    functionValue = x - b - s;
  else
    functionValue = b - x - s;
}