* `gp`: Implement sparse gaussian process methods (cf. `SparseGaussianProcess`)
* Python binding: `SparseGpMix`, see doc/tutorial
* GP/SGP API
* hyperparameter tuning : initial theta guess and bounds can be specified (`theta_init`, `theta_bounds`)
* `n_start` controls the number of optimization multistart
* In GP/SGP `rayon` is used to make parallel optimization multistart