gaussian_process_normal.hpp

Go to the documentation of this file.

/*
-------------------------------------------------------------------------
   This file is part of BayesOpt, an efficient C++ library for 
   Bayesian optimization.

   Copyright (C) 2011-2015 Ruben Martinez-Cantin <rmcantin@unizar.es>
 
   BayesOpt is free software: you can redistribute it and/or modify it 
   under the terms of the GNU Affero General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   BayesOpt 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 Affero General Public License for more details.

   You should have received a copy of the GNU Affero General Public License
   along with BayesOpt.  If not, see <http://www.gnu.org/licenses/>.
------------------------------------------------------------------------
*/


#ifndef __GAUSSIAN_PROCESS_NORMAL_HPP__
#define __GAUSSIAN_PROCESS_NORMAL_HPP__

#include "gauss_distribution.hpp"
#include "gaussian_process_hierarchical.hpp"


namespace bayesopt
{
  
  class GaussianProcessNormal: public HierarchicalGaussianProcess
  {
  public:
    GaussianProcessNormal(size_t dim, Parameters params, const Dataset& data, 
              MeanModel& mean, randEngine& eng);
    virtual ~GaussianProcessNormal();

    ProbabilityDistribution* prediction(const vectord &query);

  private:

    double negativeLogLikelihood();

    void precomputePrediction();

  private:
    vectord mWMap;                      
    vectord mW0;
    vectord mInvVarW;
    vectord mVf;
    matrixd mKF, mD;     

    GaussianDistribution* d_;      
  };

} //namespace bayesopt

#endif