NAME

sample_from_gaussian_process_prior - Samples from the prior of a Gaussian process

SYNOPSIS

#include "gp/gp_gaussian_processes.h"

Example compile flags (system dependent):
  -DLINUX_X86_64 -DLINUX_X86_64_OPTERON  -DGNU_COMPILER 
   -I/home/kobus/include
   -L/home/kobus/misc/load/linux_x86_64_opteron -L/usr/lib/x86_64-linux-gnu
  -lKJB                               -lfftw3  -lgsl -lgslcblas -ljpeg  -lSVM -lstdc++                    -lpthread -lSLATEC -lg2c    -lacml -lacml_mv -lblas -lg2c      -lncursesw 


int sample_from_gaussian_process_prior
(
	Vector_vector **sample,
	const Vector_vector *indices,
	int (*mean_func)(Vector **,const Vector *,int),
	int (*cov_func)(Matrix **,const Vector *,const Vector *,const void *,int),
	const void *hyper_params,
	int d
);

DESCRIPTION

This routine samples from the prior distribution of a Gaussian process, whose mean function and covariance function are given by mean_func and cov_func, respectively. The sample will be in the indices given by the vectors of indices, and will each be of dimension d. Naturally, the vectors of indices must all be of equal length, and mean_func and cov_func must return a d-vector and a square matrix of dimension d, respectively. The library provides a few widely-used mean and covariance functions. Also, sample will have indices->length vectors, each of length d. If the vector pointed to by sample is NULL, then a vector of the appropriate size is created. If it exists, but is the wrong size, then it is recreated. Otherwise, the storage is recycled.

RETURNS

If the routine fails (due to storage allocation, an error in the mean or covariance functions, or a mismatch in the sizes of the indices), then ERROR is returned with and error message being set. Otherwise NO_ERROR is returned.

RELATED

squared_exponential_covariance_function, zero_mean_function, sample_from_gaussian_process_prior_i

DISCLAIMER

This software is not adequatedly tested. It is recomended that results are checked independantly where appropriate.

AUTHOR

Ernesto Brau

DOCUMENTER

Ernesto Brau

SEE ALSO

fill_covariance_matrix , fill_mean_vector , sample_from_gaussian_process_prior_i , sample_from_gaussian_process_predictive , sample_from_gaussian_process_predictive_i , get_gaussian_process_predictive_distribution , get_gaussian_process_predictive_distribution_i , get_gaussian_process_posterior_distribution , get_gaussian_process_posterior_distribution_i , compute_gaussian_process_likelihood , compute_gaussian_process_likelihood_i , compute_gaussian_process_marginal_likelihood , compute_gaussian_process_marginal_likelihood_i , compute_gaussian_process_marginal_log_likelihood , compute_gaussian_process_marginal_log_likelihood_i