Sum.cc

Go to the documentation of this file.

/*
 * $Id: Sum.cc 138 2013-06-29 15:10:53Z 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"Sum.hpp"

using namespace ipo_function::concrete;

Sum::Sum( size_t const size )
  : FunctionBase { size }
{
  if( size != 0 ){
    functionGradient = gsl::vector( size );
    functionGradient.set_all( 1.0 );
    functionHessian = gsl::matrix { size, size };
    functionHessian.set_all( 0.0 );
  }
}

double
Sum::operator()( gsl::vector const& vector ){
  // Attempt sane behaviour even if vector has wrong size
  return std::accumulate( vector.begin(), vector.end(), 0.0 );
}

gsl::vector
Sum::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 );
      functionGradient.set_all( 1.0 );
    }
  }
  return functionGradient;
}

gsl::matrix
Sum::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
Sum::setVector( gsl::vector const& vector ){
  if( 0 == getSize() ){
    size_t const SIZE { vector.size() };
    if( functionGradient.size() != SIZE ) functionGradient = gsl::vector( SIZE );
    if( functionHessian.size1() != SIZE or functionHessian.size2() != SIZE )
      functionHessian = gsl::matrix { SIZE, SIZE };
    functionGradient.set_all( 1.0 );
    functionHessian.set_all( 0.0 );
  }
  functionValue = std::accumulate( vector.begin(), vector.end(), 0.0 );
}