PhaseIObjectiveFunctionAndDerivatives.cc

Go to the documentation of this file.

/*
 * $Id: PhaseIObjectiveFunctionAndDerivatives.cc 137 2013-06-29 15:10:31Z jdl3 $
 * Copyright (C) 2013 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"PhaseIObjectiveFunctionAndDerivatives.hpp"

using namespace ipo_function::detail;

PhaseIObjectiveFunctionAndDerivatives::
PhaseIObjectiveFunctionAndDerivatives( Function const& function )
  : FunctionBase { 0 == function.getSize() ? 0 : function.getSize() + 1 }
{
  if( size != 0 ){
    functionGradient = gsl::vector( size );
    functionGradient.set_all( 1.0 );
    functionHessian = gsl::matrix { size, size };
    functionHessian.set_all( 0.0 );
  }
}

double
PhaseIObjectiveFunctionAndDerivatives::operator()( gsl::vector const& vector ){
  size_t const VSIZE { vector.size() };
  if( 0 == VSIZE ) return 0; // do something sane
  return vector[VSIZE - 1];
}

gsl::vector
PhaseIObjectiveFunctionAndDerivatives::gradient( gsl::vector const& vector ){
  size_t const SIZE { 0 == getSize() ? vector.size() : getSize() };
  if( 0 == getSize() )
    if( functionGradient.size() != SIZE ) functionGradient = gsl::vector( SIZE );

  for( size_t i { 0 }; i < SIZE; ++i )
    functionGradient[i] = (i + 1 == SIZE) ? 1.0 : 0.0;
  return functionGradient;
}

gsl::matrix
PhaseIObjectiveFunctionAndDerivatives::hessian( gsl::vector const& vector ){
  size_t const SIZE { 0 == getSize() ? vector.size() : getSize() };
  if( 0 == getSize() ){
    if( functionHessian.size1() != SIZE or functionHessian.size2() != SIZE ){
      functionHessian = gsl::matrix( SIZE, SIZE );
      functionHessian.set_all( 0.0 );
    }
  }
  return functionHessian;
}

void
PhaseIObjectiveFunctionAndDerivatives::setVector( gsl::vector const& vector ){
  size_t const SIZE { 0 == getSize() ? vector.size() : getSize() };
  if( 0 == getSize() ){
    if( functionGradient.size() != SIZE ) functionGradient = gsl::vector( SIZE );
    if( functionHessian.size1() != SIZE or functionHessian.size2() != SIZE ){
      functionHessian = gsl::matrix( SIZE, SIZE );
      functionHessian.set_all( 0 );
    }
  }

  double const VSIZE = vector.size();
  functionValue = 0 == VSIZE ? 0 : vector[VSIZE - 1];

  for( size_t i = 0; i < SIZE; ++i )
    functionGradient[i] = (i + 1 == SIZE) ? 1 : 0;
}