├── .gitignore ├── Doc ├── ExamplesMap.html ├── Manuscrit_durantin_final.pdf └── Soutenance_Durantin_v2.pdf ├── Example ├── CORS_Branin_unconstrained.m ├── CoRBF_ForresterTest.m ├── Cokriging_ForresterTest.m ├── EGO_cheap_constraint.m ├── EGO_constrained.m ├── EGO_constrained_qualitative.m ├── EGO_unconstrained.m ├── EGO_unconstrained_qualitative.m ├── EI_Cokriging_Branin.m ├── EI_Cokriging_ForresterTest.m ├── Error_prediction_Branin.m ├── Gutmann_Branin.m ├── Gutmann_cheap_constraint.m ├── Kriging_Branin.m ├── Kriging_Qualitative.m ├── MGDA_prob_convex.m ├── Multiobjective_EHVI.m ├── Multiobjective_EI_MGDA.m ├── Multiobjective_EI_MGDA_qualitative.m ├── NSGA_2_constrained.m ├── RBF_Branin.m ├── Robust_EGO_1D.m ├── Robust_EGO_2D_cons.m ├── Robust_EGO_multiobj_2D.m ├── Sobol_Branin.m └── Thesis_examples │ ├── Chapitre_2 │ ├── B.1.Plan_experience │ │ ├── Full_fractionnal.m │ │ ├── LHS_bad.m │ │ ├── OLHS.m │ │ ├── Random_DOE.m │ │ └── Sobol_kmeans.m │ ├── B.2.Kriging │ │ ├── Bruit_blanc_GP.m │ │ ├── Exp_GP.m │ │ ├── Gauss_GP.m │ │ ├── Interpolation_property.m │ │ ├── Intervalle_confiance.m │ │ ├── Matern32_GP.m │ │ ├── Matern52_GP.m │ │ ├── Maximum_vraisemblance_gauss.m │ │ ├── Maximum_vraisemblance_matern52.m │ │ └── Variance_property.m │ ├── B.3.RBF │ │ ├── Cubique_RBF.m │ │ ├── Inv_multiquadrique_RBF.m │ │ ├── LOO_RBF_gauss.m │ │ ├── LOO_RBF_matern52.m │ │ ├── Lineaire_RBF.m │ │ ├── Multiquadrique_RBF.m │ │ ├── RBF_exemple_xsinx.m │ │ └── ThinPlateSpline_RBF.m │ ├── C.Cokriging │ │ └── Cokriging_exemple_forrester.m │ └── D.Benchmark │ │ ├── Cellule_PA │ │ ├── DOE_FLNS_20_1.mat │ │ ├── DOE_FLNS_80.mat │ │ ├── DOE_Kreuzer_20_1.mat │ │ ├── DOE_Kreuzer_300.mat │ │ ├── Plot_error_FLNS.m │ │ ├── Plot_error_kreuzer.m │ │ ├── Plot_validation_FLNS_RBF.m │ │ ├── Plot_validation_FLNS_kriging.m │ │ ├── Plot_validation_kreuzer_RBF.m │ │ ├── Plot_validation_kreuzer_kriging.m │ │ ├── Results_FLNS_20_RBF.mat │ │ ├── Results_FLNS_20_kriging.mat │ │ ├── Results_Kreuzer_80_RBF.mat │ │ └── Results_Kreuzer_80_kriging.mat │ │ ├── Coupleur │ │ ├── DOE_coupler_lambda_1.mat │ │ ├── DOE_coupler_lambda_10.mat │ │ ├── Plot_error_lambda.m │ │ ├── Plot_validation_RBF.m │ │ ├── Plot_validation_kriging.m │ │ ├── Results_lambda_RBF.mat │ │ └── Results_lambda_kriging.mat │ │ ├── Exemple_non_fonctionnel │ │ ├── Bencharmk_RBF_hotplate.m │ │ ├── Benchmark_kriging_cellulePA_FLNS.m │ │ └── Benchmark_kriging_coupleur.m │ │ ├── Hotplate │ │ ├── DOE_setBias_100_CVT_1.mat │ │ ├── DOE_setBias_300_CVT_test.mat │ │ ├── Hyperparametres_analyse.m │ │ ├── Plot_error_prediction.m │ │ ├── Plot_validation_RBF.m │ │ ├── Plot_validation_kriging.m │ │ ├── Results_setBias_100_CVT_RBF.mat │ │ └── Results_setBias_100_CVT_kriging.mat │ │ └── figtitle.m │ ├── Chapitre_3 │ ├── A │ │ ├── Branin_IMSE_krig.m │ │ └── Branin_MSE_krig.m │ ├── B │ │ ├── Bench_optim_constrained │ │ │ ├── EI_times_PF_opt.m │ │ │ ├── EI_times_PF_result.mat │ │ │ ├── EI_versus_PF_opt.m │ │ │ ├── EI_versus_PF_result.mat │ │ │ ├── GUT_opt.m │ │ │ ├── Gutmann_result.mat │ │ │ ├── Opt_Camel_real.mat │ │ │ ├── Plot_camel.m │ │ │ ├── Postprocess_results.m │ │ │ └── Test_points_Camel.mat │ │ ├── Bench_optim_unconstrained │ │ │ ├── CORS_bench.m │ │ │ ├── CORS_result.mat │ │ │ ├── EI_bench.m │ │ │ ├── EI_result.mat │ │ │ ├── Gutmann_bench.m │ │ │ ├── Gutmann_result.mat │ │ │ ├── Opt_Branin_real.m │ │ │ ├── Opt_Branin_real.mat │ │ │ ├── PI_bench.m │ │ │ ├── PI_result.mat │ │ │ ├── Plot_Branin.m │ │ │ ├── Postprocess_results.m │ │ │ └── Test_points_Branin.mat │ │ ├── Gradient │ │ │ ├── Plot_gradient_result.m │ │ │ ├── pas_opt.m │ │ │ ├── section_dor.m │ │ │ ├── square_2D.m │ │ │ └── steepest_slope_step_opt.m │ │ ├── Plot_EI_PI.m │ │ └── Plot_Gutmann.m │ ├── C │ │ ├── bench_MO │ │ │ ├── EHVI_opt.m │ │ │ ├── EHVI_results.mat │ │ │ ├── EI_euclid_opt.m │ │ │ ├── EI_euclid_results.mat │ │ │ ├── NSGA_opt.m │ │ │ ├── NSGA_results.mat │ │ │ └── Post_traitement_resultats.m │ │ └── kmeans_exemple.m │ └── D │ │ └── Exemple_robust_measure.m │ ├── Chapitre_4 │ ├── Bench_coRBF │ │ ├── DOE_borehole_500.mat │ │ ├── DOE_curretal_200.mat │ │ ├── borehole_HF.m │ │ ├── borehole_LF.m │ │ ├── borehole_cokrig.m │ │ ├── borehole_cokrig_LOO.m │ │ ├── borehole_corbf.m │ │ ├── coRBF_borehole_result.mat │ │ ├── coRBF_curretal_result.mat │ │ ├── coRBF_exemple.m │ │ ├── cokrig_borehole_result.mat │ │ ├── cokrig_borehole_result_LOO.mat │ │ ├── cokrig_curretal_result.mat │ │ ├── curretal88exp.m │ │ ├── curretal_HF.m │ │ ├── curretal_LF.m │ │ ├── curretal_cokrig.m │ │ ├── curretal_corbf.m │ │ ├── plot_borehole.m │ │ ├── plot_borehole_cokriging_LOO.m │ │ └── plot_curretal.m │ └── bench_EI_MGDA │ │ ├── EHVI_opt.m │ │ ├── EHVI_results.mat │ │ ├── EI_MGDA_opt.m │ │ ├── EI_MGDA_results.mat │ │ └── Plot_results.m │ └── Read_me_(french).txt ├── IMPORTANT_LICENSE_NOTICE ├── LICENSE ├── Lib ├── +ooDACE │ ├── +basisfunctions │ │ ├── corrcubic.m │ │ ├── correxp.m │ │ ├── corrgauss.m │ │ ├── corrgaussp.m │ │ ├── corrlin.m │ │ ├── corrmatern32.m │ │ ├── corrmatern32_iso.m │ │ ├── corrmatern32_variance.m │ │ ├── corrmatern32iso.m │ │ ├── corrmatern52.m │ │ ├── corrspherical.m │ │ ├── corrspline.m │ │ └── covmatern32.m │ ├── +tools │ │ ├── averageEuclideanError.m │ │ ├── buildVandermondeMatrix.m │ │ ├── cfix.m │ │ ├── makeEvalGrid.m │ │ ├── makeGrid.m │ │ ├── mergeStruct.m │ │ ├── plotScatteredData.m │ │ └── powerBase.m │ ├── @BasicGaussianProcess │ │ ├── BasicGaussianProcess.m │ │ ├── correlationFunction.m │ │ ├── cvpe.m │ │ ├── extrinsicCorrelationMatrix.m │ │ ├── fit.m │ │ ├── generateDegrees.m │ │ ├── getExpression.m │ │ ├── imse.m │ │ ├── intrinsicCovarianceMatrix.m │ │ ├── likelihood.m │ │ ├── marginalLikelihood.m │ │ ├── mseTestset.m │ │ ├── plotLikelihood.m │ │ ├── plotVariogram.m │ │ ├── predict.m │ │ ├── predict_derivatives.m │ │ ├── pseudoLikelihood.m │ │ ├── rcValues.m │ │ ├── regressionFunction.m │ │ ├── regressionMatrix.m │ │ ├── tuneParameters.m │ │ ├── updateModel.m │ │ ├── updateRegression.m │ │ └── updateStochasticProcess.m │ ├── @CmaesOptimizer │ │ ├── CmaesOptimizer.m │ │ └── optimize.m │ ├── @CoKriging │ │ ├── CoKriging.m │ │ ├── correlationFunction.m │ │ ├── extrinsicCorrelationMatrix.m │ │ ├── fit.m │ │ ├── intrinsicCovarianceMatrix.m │ │ ├── regressionFunction.m │ │ ├── regressionMatrix.m │ │ └── setData.m │ ├── @Kriging │ │ ├── Kriging.m │ │ ├── cvpe.m │ │ ├── getExpression.m │ │ ├── predict.m │ │ ├── predict_derivatives.m │ │ └── setData.m │ ├── @MatlabOptimizer │ │ ├── MatlabOptimizer.m │ │ ├── optimize.m │ │ └── setInputConstraints.m │ ├── @Optimizer │ │ ├── Optimizer.m │ │ ├── getBounds.m │ │ ├── getHint.m │ │ ├── getInitialPopulation.m │ │ ├── getInputDimension.m │ │ ├── getOutputDimension.m │ │ ├── getPopulationSize.m │ │ ├── optimize.m │ │ ├── setBounds.m │ │ ├── setDimensions.m │ │ ├── setHint.m │ │ ├── setInitialPopulation.m │ │ └── setInputConstraints.m │ ├── @SQPLabOptimizer │ │ ├── SQPLabOptimizer.m │ │ ├── optimize.m │ │ └── simulator.m │ ├── agpl-3.0.txt │ ├── license.txt │ ├── oodacefit.m │ └── plotKrigingModel.m ├── cholincsp.m ├── cmaes.m ├── ipdm.m ├── license_ipdm.txt ├── sqplab-0.4.5-distrib │ ├── COPYING │ ├── VERSIONS.txt │ ├── example │ │ ├── README │ │ ├── hanging_chain.m │ │ ├── hanging_chain_data.m │ │ ├── hanging_chain_dimopt.m │ │ ├── hanging_chain_initopt.m │ │ ├── hanging_chain_plot.m │ │ ├── hanging_chain_plot_params.m │ │ ├── hanging_chain_postopt.m │ │ └── hanging_chain_simul.m │ └── src │ │ ├── sqplab.m │ │ ├── sqplab_armijo.m │ │ ├── sqplab_badsimul.m │ │ ├── sqplab_bfgs.m │ │ ├── sqplab_bfgs_inv.m │ │ ├── sqplab_checkoptions.m │ │ ├── sqplab_dnorm.m │ │ ├── sqplab_finish.m │ │ ├── sqplab_givens.m │ │ ├── sqplab_lgivens.m │ │ ├── sqplab_loop.m │ │ ├── sqplab_lsmult.m │ │ ├── sqplab_merit.m │ │ ├── sqplab_optimality.m │ │ ├── sqplab_options.m │ │ ├── sqplab_pnorm.m │ │ ├── sqplab_prelim.m │ │ ├── sqplab_qrg.m │ │ ├── sqplab_rgivens.m │ │ ├── sqplab_rqn.m │ │ ├── sqplab_sigma.m │ │ ├── sqplab_step.m │ │ ├── sqplab_tcg.m │ │ ├── sqplab_tr.m │ │ └── sqplab_wolfe.m └── stk │ ├── AUTHORS │ ├── CITATION │ ├── COPYING │ ├── ChangeLog │ ├── NEWS │ ├── README.md │ ├── arrays │ ├── @stk_dataframe │ │ ├── abs.m │ │ ├── acos.m │ │ ├── acosd.m │ │ ├── acosh.m │ │ ├── and.m │ │ ├── apply.m │ │ ├── asin.m │ │ ├── asind.m │ │ ├── asinh.m │ │ ├── atan.m │ │ ├── atand.m │ │ ├── atanh.m │ │ ├── bsxfun.m │ │ ├── cat.m │ │ ├── cos.m │ │ ├── cosd.m │ │ ├── cosh.m │ │ ├── ctranspose.m │ │ ├── diff.m │ │ ├── disp.m │ │ ├── display.m │ │ ├── double.m │ │ ├── end.m │ │ ├── eq.m │ │ ├── exp.m │ │ ├── expm1.m │ │ ├── fieldnames.m │ │ ├── ge.m │ │ ├── get.m │ │ ├── gt.m │ │ ├── horzcat.m │ │ ├── isempty.m │ │ ├── isequal.m │ │ ├── ismember.m │ │ ├── ldivide.m │ │ ├── le.m │ │ ├── length.m │ │ ├── log.m │ │ ├── log10.m │ │ ├── log1p.m │ │ ├── log2.m │ │ ├── logical.m │ │ ├── lt.m │ │ ├── max.m │ │ ├── mean.m │ │ ├── median.m │ │ ├── min.m │ │ ├── minus.m │ │ ├── mode.m │ │ ├── mtimes.m │ │ ├── ne.m │ │ ├── openvar.m │ │ ├── or.m │ │ ├── plot.m │ │ ├── plus.m │ │ ├── power.m │ │ ├── private │ │ │ ├── get_column_number.c │ │ │ ├── get_column_number.m │ │ │ ├── get_column_number.mexa64 │ │ │ ├── get_column_number.mexa64.info │ │ │ ├── get_column_number.mexw64 │ │ │ ├── get_column_number.mexw64.info │ │ │ ├── names_errmsg.m │ │ │ ├── process_cell_indices.m │ │ │ ├── process_char_indices.m │ │ │ └── reserved_field_names.m │ │ ├── prod.m │ │ ├── quantile.m │ │ ├── rdivide.m │ │ ├── realpow.m │ │ ├── reshape.m │ │ ├── set.m │ │ ├── sin.m │ │ ├── sind.m │ │ ├── sinh.m │ │ ├── size.m │ │ ├── sort.m │ │ ├── sqrt.m │ │ ├── std.m │ │ ├── stk_boundingbox.m │ │ ├── stk_conditioning.m │ │ ├── stk_dataframe.m │ │ ├── stk_generate_samplepaths.m │ │ ├── stk_length.m │ │ ├── stk_normalize.m │ │ ├── stk_rescale.m │ │ ├── stk_sprintf.m │ │ ├── stk_sprintf_colnames.m │ │ ├── stk_sprintf_rownames.m │ │ ├── stk_sprintf_sizetype.m │ │ ├── subsasgn.m │ │ ├── subsref.m │ │ ├── sum.m │ │ ├── tan.m │ │ ├── tand.m │ │ ├── tanh.m │ │ ├── times.m │ │ ├── transpose.m │ │ ├── uminus.m │ │ ├── uplus.m │ │ ├── var.m │ │ ├── vertcat.m │ │ ├── xlim.m │ │ ├── xor.m │ │ └── ylim.m │ ├── @stk_factorialdesign │ │ ├── contour.m │ │ ├── fieldnames.m │ │ ├── get.m │ │ ├── isequal.m │ │ ├── mesh.m │ │ ├── meshc.m │ │ ├── meshz.m │ │ ├── ndgrid.m │ │ ├── pcolor.m │ │ ├── private │ │ │ └── plot_surfmeshcontour.m │ │ ├── set.m │ │ ├── stk_boundingbox.m │ │ ├── stk_factorialdesign.m │ │ ├── stk_normalize.m │ │ ├── stk_rescale.m │ │ ├── stk_sprintf.m │ │ ├── stk_sprintf_levels.m │ │ ├── subsasgn.m │ │ ├── surf.m │ │ └── uminus.m │ ├── @stk_hrect │ │ ├── axis.m │ │ ├── disp.m │ │ ├── display.m │ │ ├── get.m │ │ ├── horzcat.m │ │ ├── ismember.m │ │ ├── stk_boundingbox.m │ │ ├── stk_hrect.m │ │ ├── stk_normalize.m │ │ ├── stk_rescale.m │ │ ├── subsasgn.m │ │ ├── subsref.m │ │ └── vertcat.m │ ├── @struct │ │ ├── double.m │ │ ├── stk_boundingbox.m │ │ └── stk_length.m │ └── generic │ │ ├── stk_boundingbox.m │ │ ├── stk_commonsize.m │ │ ├── stk_feval.m │ │ ├── stk_length.m │ │ ├── stk_normalize.m │ │ ├── stk_rescale.m │ │ ├── stk_sprintf.m │ │ ├── stk_sprintf_colnames.m │ │ ├── stk_sprintf_rownames.m │ │ └── stk_sprintf_sizetype.m │ ├── core │ ├── @stk_kreq_qr │ │ ├── disp.m │ │ ├── display.m │ │ ├── get.m │ │ ├── isequal.m │ │ ├── linsolve.m │ │ ├── private │ │ │ └── compute_P_scaling.m │ │ ├── stk_kreq_qr.m │ │ ├── stk_set_righthandside.m │ │ ├── stk_squared_seminorm.m │ │ ├── stk_update.m │ │ └── subsref.m │ ├── @stk_model_gpposterior │ │ ├── disp.m │ │ ├── display.m │ │ ├── fieldnames.m │ │ ├── get.m │ │ ├── get_input_data.m │ │ ├── get_output_data.m │ │ ├── get_prior_model.m │ │ ├── horzcat.m │ │ ├── isequal.m │ │ ├── set.m │ │ ├── stk_isnoisy.m │ │ ├── stk_model_gpposterior.m │ │ ├── stk_model_update.m │ │ ├── stk_predict.m │ │ ├── stk_predict_leaveoneout.m │ │ ├── subsasgn.m │ │ ├── subsref.m │ │ └── vertcat.m │ ├── stk_cholcov.m │ ├── stk_isnoisy.m │ ├── stk_make_matcov.m │ ├── stk_model.m │ ├── stk_model_fixlm.m │ ├── stk_model_update.m │ ├── stk_ortho_func.m │ ├── stk_predict.m │ └── stk_predict_leaveoneout.m │ ├── covfcs │ ├── rbf │ │ ├── stk_rbf_exponential.m │ │ ├── stk_rbf_gauss.m │ │ ├── stk_rbf_matern.m │ │ ├── stk_rbf_matern32.m │ │ ├── stk_rbf_matern52.m │ │ ├── stk_rbf_spherical.m │ │ ├── stk_sf_gausscorr.m │ │ ├── stk_sf_matern.m │ │ ├── stk_sf_matern32.m │ │ └── stk_sf_matern52.m │ ├── stk_discretecov.m │ ├── stk_expcov_aniso.m │ ├── stk_expcov_iso.m │ ├── stk_gausscov_aniso.m │ ├── stk_gausscov_iso.m │ ├── stk_materncov32_aniso.m │ ├── stk_materncov32_iso.m │ ├── stk_materncov52_aniso.m │ ├── stk_materncov52_iso.m │ ├── stk_materncov_aniso.m │ ├── stk_materncov_iso.m │ ├── stk_noisecov.m │ ├── stk_sphcov_aniso.m │ └── stk_sphcov_iso.m │ ├── doc │ ├── dev │ │ └── model.texi │ └── html │ │ ├── COPYING.html │ │ ├── NEWS.html │ │ ├── css │ │ └── stk.css │ │ ├── function │ │ ├── @stk_dataframe │ │ │ └── stk_dataframe.html │ │ ├── @stk_factorialdesign │ │ │ ├── ndgrid.html │ │ │ └── stk_factorialdesign.html │ │ ├── @stk_hrect │ │ │ └── stk_hrect.html │ │ ├── @stk_lm_affine │ │ │ └── stk_lm_affine.html │ │ ├── @stk_lm_constant │ │ │ └── stk_lm_constant.html │ │ ├── @stk_lm_cubic │ │ │ └── stk_lm_cubic.html │ │ ├── @stk_lm_matrix │ │ │ └── stk_lm_matrix.html │ │ ├── @stk_lm_null │ │ │ └── stk_lm_null.html │ │ ├── @stk_lm_quadratic │ │ │ └── stk_lm_quadratic.html │ │ ├── @stk_model_gpposterior │ │ │ ├── get_input_data.html │ │ │ ├── get_output_data.html │ │ │ ├── get_prior_model.html │ │ │ ├── stk_model_gpposterior.html │ │ │ ├── stk_model_update.html │ │ │ ├── stk_predict.html │ │ │ └── stk_predict_leaveoneout.html │ │ ├── @stk_optim_fmincon │ │ │ └── stk_optim_fmincon.html │ │ ├── @stk_optim_fminsearch │ │ │ └── stk_optim_fminsearch.html │ │ ├── @stk_optim_octavesqp │ │ │ └── stk_optim_octavesqp.html │ │ ├── stk_boundingbox.html │ │ ├── stk_conditioning.html │ │ ├── stk_dataset_twobumps.html │ │ ├── stk_discretecov.html │ │ ├── stk_dist.html │ │ ├── stk_distrib_normal_ei.html │ │ ├── stk_distrib_student_ei.html │ │ ├── stk_dominatedhv.html │ │ ├── stk_example_doe01.html │ │ ├── stk_example_doe02.html │ │ ├── stk_example_doe03.html │ │ ├── stk_example_doe04.html │ │ ├── stk_example_kb01.html │ │ ├── stk_example_kb01n.html │ │ ├── stk_example_kb02.html │ │ ├── stk_example_kb02n.html │ │ ├── stk_example_kb03.html │ │ ├── stk_example_kb05.html │ │ ├── stk_example_kb06.html │ │ ├── stk_example_kb07.html │ │ ├── stk_example_kb08.html │ │ ├── stk_example_kb09.html │ │ ├── stk_example_kb10.html │ │ ├── stk_example_misc01.html │ │ ├── stk_example_misc02.html │ │ ├── stk_example_misc03.html │ │ ├── stk_example_misc04.html │ │ ├── stk_expcov_aniso.html │ │ ├── stk_expcov_iso.html │ │ ├── stk_feval.html │ │ ├── stk_filldist.html │ │ ├── stk_gausscov_aniso.html │ │ ├── stk_gausscov_iso.html │ │ ├── stk_generate_samplepaths.html │ │ ├── stk_init.html │ │ ├── stk_isdominated.html │ │ ├── stk_isnoisy.html │ │ ├── stk_length.html │ │ ├── stk_lm_polynomial.html │ │ ├── stk_make_matcov.html │ │ ├── stk_materncov32_aniso.html │ │ ├── stk_materncov32_iso.html │ │ ├── stk_materncov52_aniso.html │ │ ├── stk_materncov52_iso.html │ │ ├── stk_materncov_aniso.html │ │ ├── stk_materncov_iso.html │ │ ├── stk_maxabscorr.html │ │ ├── stk_mindist.html │ │ ├── stk_minimize_boxconstrained.html │ │ ├── stk_minimize_unconstrained.html │ │ ├── stk_model.html │ │ ├── stk_model_update.html │ │ ├── stk_noisecov.html │ │ ├── stk_normalize.html │ │ ├── stk_optim_isavailable.html │ │ ├── stk_options_get.html │ │ ├── stk_options_set.html │ │ ├── stk_param_estim.html │ │ ├── stk_param_gls.html │ │ ├── stk_param_init.html │ │ ├── stk_param_init_lnv.html │ │ ├── stk_param_relik.html │ │ ├── stk_paretofind.html │ │ ├── stk_phipcrit.html │ │ ├── stk_plot1d.html │ │ ├── stk_plot_histnormres.html │ │ ├── stk_plot_predvsobs.html │ │ ├── stk_pmisclass.html │ │ ├── stk_predict.html │ │ ├── stk_predict_leaveoneout.html │ │ ├── stk_rbf_exponential.html │ │ ├── stk_rbf_gauss.html │ │ ├── stk_rbf_matern.html │ │ ├── stk_rbf_matern32.html │ │ ├── stk_rbf_matern52.html │ │ ├── stk_rbf_spherical.html │ │ ├── stk_rescale.html │ │ ├── stk_runtests.html │ │ ├── stk_sampcrit_ehvi_eval.html │ │ ├── stk_sampcrit_ei_eval.html │ │ ├── stk_sampcrit_emmi_eval.html │ │ ├── stk_sampling_halton_rr2.html │ │ ├── stk_sampling_maximinlhs.html │ │ ├── stk_sampling_olhs.html │ │ ├── stk_sampling_randomlhs.html │ │ ├── stk_sampling_randunif.html │ │ ├── stk_sampling_regulargrid.html │ │ ├── stk_sampling_sobol.html │ │ ├── stk_sampling_vdc_rr2.html │ │ ├── stk_sphcov_aniso.html │ │ ├── stk_sphcov_iso.html │ │ ├── stk_testfun_borehole.html │ │ ├── stk_testfun_braninhoo.html │ │ ├── stk_testfun_goldsteinprice.html │ │ ├── stk_testfun_twobumps.html │ │ └── stk_version.html │ │ ├── images │ │ ├── download.png │ │ ├── news.png │ │ └── stk_logo.png │ │ └── index.html │ ├── examples │ ├── 01_kriging_basics │ │ ├── stk_example_kb01.m │ │ ├── stk_example_kb01n.m │ │ ├── stk_example_kb02.m │ │ ├── stk_example_kb02n.m │ │ ├── stk_example_kb03.m │ │ ├── stk_example_kb04.m │ │ ├── stk_example_kb05.m │ │ ├── stk_example_kb06.m │ │ ├── stk_example_kb07.m │ │ ├── stk_example_kb08.m │ │ ├── stk_example_kb09.m │ │ └── stk_example_kb10.m │ ├── 02_design_of_experiments │ │ ├── stk_example_doe01.m │ │ ├── stk_example_doe02.m │ │ ├── stk_example_doe03.m │ │ └── stk_example_doe04.m │ ├── 03_miscellaneous │ │ ├── stk_example_misc01.m │ │ ├── stk_example_misc02.m │ │ ├── stk_example_misc03.m │ │ └── stk_example_misc04.m │ ├── datasets │ │ └── stk_dataset_twobumps.m │ └── test_functions │ │ ├── stk_testfun_borehole.m │ │ ├── stk_testfun_braninhoo.m │ │ ├── stk_testfun_goldsteinprice.m │ │ └── stk_testfun_twobumps.m │ ├── lm │ ├── @stk_lm_affine │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_affine.m │ ├── @stk_lm_constant │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_constant.m │ ├── @stk_lm_cubic │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_cubic.m │ ├── @stk_lm_matrix │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_matrix.m │ ├── @stk_lm_null │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_null.m │ ├── @stk_lm_quadratic │ │ ├── feval.m │ │ ├── isequal.m │ │ └── stk_lm_quadratic.m │ └── stk_lm_polynomial.m │ ├── misc │ ├── benchmarks │ │ ├── stk_benchmark_bivnorm.m │ │ ├── stk_benchmark_cov1.m │ │ ├── stk_benchmark_cov2.m │ │ ├── stk_benchmark_examples.m │ │ ├── stk_benchmark_init.m │ │ ├── stk_benchmark_logdet.m │ │ └── stk_benchmark_paramestim.m │ ├── design │ │ ├── stk_maxabscorr.m │ │ └── stk_phipcrit.m │ ├── dist │ │ ├── private │ │ │ ├── stk_dist_matrixx.c │ │ │ ├── stk_dist_matrixx.mexa64 │ │ │ ├── stk_dist_matrixx.mexa64.info │ │ │ ├── stk_dist_matrixx.mexw64 │ │ │ ├── stk_dist_matrixx.mexw64.info │ │ │ ├── stk_dist_matrixy.c │ │ │ ├── stk_dist_matrixy.mexa64 │ │ │ ├── stk_dist_matrixy.mexa64.info │ │ │ ├── stk_dist_matrixy.mexw64 │ │ │ ├── stk_dist_matrixy.mexw64.info │ │ │ ├── stk_dist_pairwise.c │ │ │ ├── stk_dist_pairwise.mexa64 │ │ │ ├── stk_dist_pairwise.mexa64.info │ │ │ ├── stk_dist_pairwise.mexw64 │ │ │ ├── stk_dist_pairwise.mexw64.info │ │ │ ├── stk_filldist_discr_mex.c │ │ │ ├── stk_filldist_discr_mex.mexa64 │ │ │ ├── stk_filldist_discr_mex.mexa64.info │ │ │ ├── stk_filldist_discr_mex.mexw64 │ │ │ ├── stk_filldist_discr_mex.mexw64.info │ │ │ ├── stk_gpquadform_matrixx.c │ │ │ ├── stk_gpquadform_matrixx.mexa64 │ │ │ ├── stk_gpquadform_matrixx.mexa64.info │ │ │ ├── stk_gpquadform_matrixx.mexw64 │ │ │ ├── stk_gpquadform_matrixx.mexw64.info │ │ │ ├── stk_gpquadform_matrixy.c │ │ │ ├── stk_gpquadform_matrixy.mexa64 │ │ │ ├── stk_gpquadform_matrixy.mexa64.info │ │ │ ├── stk_gpquadform_matrixy.mexw64 │ │ │ ├── stk_gpquadform_matrixy.mexw64.info │ │ │ ├── stk_gpquadform_pairwise.c │ │ │ ├── stk_gpquadform_pairwise.mexa64 │ │ │ ├── stk_gpquadform_pairwise.mexa64.info │ │ │ ├── stk_gpquadform_pairwise.mexw64 │ │ │ ├── stk_gpquadform_pairwise.mexw64.info │ │ │ ├── stk_mindist_mex.c │ │ │ ├── stk_mindist_mex.mexa64 │ │ │ ├── stk_mindist_mex.mexa64.info │ │ │ ├── stk_mindist_mex.mexw64 │ │ │ └── stk_mindist_mex.mexw64.info │ │ ├── stk_dist.m │ │ ├── stk_filldist.m │ │ ├── stk_filldist_discretized.m │ │ ├── stk_filldist_exact.m │ │ ├── stk_gpquadform.m │ │ └── stk_mindist.m │ ├── distrib │ │ ├── README │ │ ├── private │ │ │ ├── stk_distrib_bivnorm0_cdf.c │ │ │ ├── stk_distrib_bivnorm0_cdf.mexa64 │ │ │ ├── stk_distrib_bivnorm0_cdf.mexa64.info │ │ │ ├── stk_distrib_bivnorm0_cdf.mexw64 │ │ │ └── stk_distrib_bivnorm0_cdf.mexw64.info │ │ ├── stk_distrib_bivnorm_cdf.m │ │ ├── stk_distrib_normal_cdf.m │ │ ├── stk_distrib_normal_ei.m │ │ ├── stk_distrib_normal_pdf.m │ │ ├── stk_distrib_student_cdf.m │ │ ├── stk_distrib_student_ei.m │ │ └── stk_distrib_student_pdf.m │ ├── error │ │ ├── stk_assert_box.m │ │ └── stk_error.m │ ├── include │ │ ├── pareto.h │ │ ├── primes.h │ │ └── stk_mex.h │ ├── mole │ │ ├── README │ │ ├── graphics_toolkit │ │ │ └── graphics_toolkit.m │ │ ├── isoctave │ │ │ └── isoctave.m │ │ ├── isrow │ │ │ └── isrow.m │ │ ├── linsolve │ │ │ └── linsolve.m │ │ ├── matlab │ │ │ ├── fflush.m │ │ │ ├── file_in_loadpath.m │ │ │ ├── file_in_path.m │ │ │ ├── index.m │ │ │ └── stdout.m │ │ └── quantile │ │ │ └── quantile.m │ ├── optim │ │ ├── @stk_optim_baseclass │ │ │ ├── isequal.m │ │ │ ├── stk_optim_baseclass.m │ │ │ └── stk_optim_isavailable.m │ │ ├── @stk_optim_fmincon │ │ │ ├── isequal.m │ │ │ ├── stk_minimize_boxconstrained.m │ │ │ └── stk_optim_fmincon.m │ │ ├── @stk_optim_fminsearch │ │ │ ├── isequal.m │ │ │ ├── stk_minimize_boxconstrained.m │ │ │ ├── stk_minimize_unconstrained.m │ │ │ └── stk_optim_fminsearch.m │ │ ├── @stk_optim_octavesqp │ │ │ ├── isequal.m │ │ │ ├── private │ │ │ │ ├── sqp_quadprog.m │ │ │ │ ├── sqp_quadprog_testf.m │ │ │ │ └── sqp_quadprog_testg.m │ │ │ ├── stk_minimize_boxconstrained.m │ │ │ ├── stk_minimize_unconstrained.m │ │ │ └── stk_optim_octavesqp.m │ │ ├── stk_minimize_boxconstrained.m │ │ ├── stk_minimize_unconstrained.m │ │ ├── stk_optim_hasfmincon.m │ │ ├── stk_optim_isavailable.m │ │ ├── stk_optim_testmin_box.m │ │ ├── stk_optim_testmin_unc.m │ │ └── stk_select_optimizer.m │ ├── options │ │ ├── stk_options_get.m │ │ └── stk_options_set.m │ ├── parallel │ │ ├── @stk_parallel_engine_none │ │ │ ├── isequal.m │ │ │ ├── stk_parallel_engine_none.m │ │ │ ├── stk_parallel_feval.m │ │ │ └── stk_parallel_stop.m │ │ ├── @stk_parallel_engine_parfor │ │ │ ├── isequal.m │ │ │ ├── stk_parallel_engine_parfor.m │ │ │ ├── stk_parallel_feval.m │ │ │ └── stk_parallel_stop.m │ │ ├── stk_parallel_cutblocks.m │ │ ├── stk_parallel_engine_get.m │ │ ├── stk_parallel_engine_set.m │ │ ├── stk_parallel_feval.m │ │ ├── stk_parallel_start.m │ │ └── stk_parallel_stop.m │ ├── pareto │ │ ├── private │ │ │ ├── stk_dominatedhv_mex.c │ │ │ ├── stk_dominatedhv_mex.mexa64 │ │ │ ├── stk_dominatedhv_mex.mexa64.info │ │ │ ├── stk_dominatedhv_mex.mexw64 │ │ │ ├── stk_dominatedhv_mex.mexw64.info │ │ │ ├── stk_isdominated_mex.c │ │ │ ├── stk_isdominated_mex.mexa64 │ │ │ ├── stk_isdominated_mex.mexa64.info │ │ │ ├── stk_isdominated_mex.mexw64 │ │ │ ├── stk_isdominated_mex.mexw64.info │ │ │ ├── stk_paretofind_mex.c │ │ │ ├── stk_paretofind_mex.mexa64 │ │ │ ├── stk_paretofind_mex.mexa64.info │ │ │ ├── stk_paretofind_mex.mexw64 │ │ │ ├── stk_paretofind_mex.mexw64.info │ │ │ ├── wfg.README │ │ │ ├── wfg.c │ │ │ └── wfg.h │ │ ├── stk_dominatedhv.m │ │ ├── stk_isdominated.m │ │ └── stk_paretofind.m │ ├── plot │ │ ├── stk_axes.m │ │ ├── stk_figure.m │ │ ├── stk_labels.m │ │ ├── stk_plot1d.m │ │ ├── stk_plot_getaxesarg.m │ │ ├── stk_plot_histnormres.m │ │ ├── stk_plot_predvsobs.m │ │ ├── stk_plot_shadedci.m │ │ ├── stk_subplot.m │ │ ├── stk_title.m │ │ ├── stk_xlabel.m │ │ ├── stk_ylabel.m │ │ └── stk_zlabel.m │ ├── test │ │ ├── stk_is_lhs.m │ │ ├── stk_isequal_tolabs.m │ │ ├── stk_isequal_tolrel.m │ │ ├── stk_runtests.m │ │ ├── stk_test.m │ │ ├── stk_test_class.m │ │ └── stk_test_dfbinaryop.m │ └── text │ │ ├── stk_disp_examplewelcome.m │ │ ├── stk_disp_isloose.m │ │ ├── stk_disp_progress.m │ │ ├── stk_sprintf_colvect.m │ │ ├── stk_sprintf_colvect_fixedpoint.m │ │ ├── stk_sprintf_colvect_scientific.m │ │ └── stk_sprintf_framed.m │ ├── paramestim │ ├── private │ │ ├── stk_get_optimizable_parameters.m │ │ └── stk_set_optimizable_parameters.m │ ├── stk_param_estim.m │ ├── stk_param_getdefaultbounds.m │ ├── stk_param_gls.m │ ├── stk_param_init.m │ ├── stk_param_init_lnv.m │ └── stk_param_relik.m │ ├── sampling │ ├── private │ │ ├── stk_sampling_sobol_mex.c │ │ ├── stk_sampling_sobol_mex.mexa64 │ │ ├── stk_sampling_sobol_mex.mexa64.info │ │ ├── stk_sampling_sobol_mex.mexw64 │ │ └── stk_sampling_sobol_mex.mexw64.info │ ├── stk_sampcrit_ehvi_eval.m │ ├── stk_sampcrit_ei_eval.m │ ├── stk_sampcrit_emmi_eval.m │ ├── stk_sampling_halton_rr2.m │ ├── stk_sampling_maximinlhs.m │ ├── stk_sampling_olhs.m │ ├── stk_sampling_randomlhs.m │ ├── stk_sampling_randunif.m │ ├── stk_sampling_regulargrid.m │ ├── stk_sampling_sobol.m │ ├── stk_sampling_vdc_rr2.c │ ├── stk_sampling_vdc_rr2.m │ ├── stk_sampling_vdc_rr2.mexa64 │ ├── stk_sampling_vdc_rr2.mexa64.info │ ├── stk_sampling_vdc_rr2.mexw64 │ └── stk_sampling_vdc_rr2.mexw64.info │ ├── stk_init.m │ ├── stk_version.m │ └── utils │ ├── stk_conditioning.m │ ├── stk_generate_samplepaths.m │ └── stk_pmisclass.m ├── README.md └── Src ├── Adaptive_Sampling ├── Derived_objects │ ├── @EI_cheap_constraint │ │ ├── EI_cheap_constraint.m │ │ ├── EI_constrained.m │ │ ├── EI_unconstrained.m │ │ └── Opt_crit.m │ ├── @EI_multifi │ │ ├── EI_multifi.m │ │ └── Find_min_value.m │ ├── @EI_times_PF │ │ ├── EI_constrained.m │ │ ├── EI_times_PF.m │ │ ├── Opt_crit.m │ │ └── Set_options_optim.m │ ├── @Gutmann_RBF_cheap_constraint │ │ ├── Gutmann_RBF_cheap_constraint.m │ │ ├── Gutmann_crit.m │ │ ├── Obj_predict.m │ │ └── Opt_crit.m │ ├── @Gutmann_multifi │ │ ├── Find_min_value.m │ │ └── Gutmann_multifi.m │ ├── @Q_expected_improvement │ │ ├── Find_min_value.m │ │ ├── Opt_crit.m │ │ └── Q_expected_improvement.m │ ├── @Q_multi_obj_EI_MGDA │ │ ├── Find_min_value.m │ │ ├── Obj_func_EI.m │ │ ├── Obj_func_nsga.m │ │ ├── Opt_crit.m │ │ └── Q_multi_obj_EI_MGDA.m │ └── @Robust_EGO_multiobj │ │ ├── Eval_rob_meas.m │ │ └── Robust_EGO_multiobj.m └── Main_objects │ ├── @Adaptive_sampling │ ├── Adaptive_sampling.m │ ├── Conv_check.m │ ├── Eval.m │ ├── K_filtering.m │ ├── Opt.m │ └── Restart.m │ ├── @CORS_RBF │ ├── CORS_RBF.m │ ├── Conv_check_crit.m │ ├── Find_min_value.m │ ├── Opt_crit.m │ ├── Prediction.m │ └── Set_options_optim.m │ ├── @Error_prediction │ ├── Conv_check_crit.m │ ├── Error_crit.m │ ├── Error_prediction.m │ ├── Meta_int1.m │ ├── Meta_int2.m │ ├── Meta_int3.m │ ├── Opt_crit.m │ └── Set_options_optim.m │ ├── @Expected_improvement │ ├── Conv_check_crit.m │ ├── EI_constrained.m │ ├── EI_unconstrained.m │ ├── Expected_improvement.m │ ├── Find_min_value.m │ ├── Opt_crit.m │ └── Set_options_optim.m │ ├── @Gutmann_RBF │ ├── Conv_check_crit.m │ ├── Find_min_value.m │ ├── Gutmann_RBF.m │ ├── Gutmann_crit.m │ ├── Opt_crit.m │ ├── Prediction.m │ └── Set_options_optim.m │ ├── @Multi_obj_EHVI │ ├── Conv_check_crit.m │ ├── EI_euclid.m │ ├── Find_min_value.m │ ├── Multi_obj_EHVI.m │ ├── Obj_func_EHVI.m │ ├── Obj_func_EI_euclid.m │ ├── Obj_func_nsga.m │ ├── Opt_crit.m │ └── Set_options_optim.m │ ├── @Multi_obj_EI_MGDA │ ├── Conv_check_crit.m │ ├── Find_min_value.m │ ├── Multi_obj_EI_MGDA.m │ ├── Obj_func_EI.m │ ├── Obj_func_nsga.m │ ├── Opt_crit.m │ └── Set_options_optim.m │ └── @Robust_EGO │ ├── Compute_CRN.m │ ├── Conv_check_crit.m │ ├── EI_constrained.m │ ├── EI_unconstrained.m │ ├── Eval_rob_meas.m │ ├── Find_min_value_opt.m │ ├── Find_min_value_prob.m │ ├── Mean_meas.m │ ├── Mixed_meas.m │ ├── Opt_crit.m │ ├── Optim_env.m │ ├── Optim_meas.m │ ├── Robust_EGO.m │ ├── Set_options_optim.m │ ├── Update_CRN.m │ ├── Var_meas.m │ ├── Worstcase_meas.m │ └── X_to_rob.m ├── Analytical_test_function ├── BQQV │ ├── Q_sin_cos.m │ ├── emse_1.m │ ├── emse_2.m │ ├── test_ego_constrained.m │ ├── test_ego_unconstrained.m │ └── test_multi_obj.m ├── Branin.m ├── Branin_LF.m ├── Camel_constrained.m ├── Cheap_cons_branin.m ├── MO_constrained.m ├── MO_convex.m ├── MO_convex_nsga.m ├── Multifi_1D_HF.m ├── Multifi_1D_LF.m ├── Robust_1D.m └── test_multi_obj_ind.m ├── Metamodel ├── @CoRBF │ ├── Clean.m │ ├── CoRBF.m │ ├── Def_hyp_corr.m │ ├── Obj_diff_opt.m │ ├── Predict.m │ └── Train.m ├── @Cokriging │ ├── Clean.m │ ├── Cokriging.m │ ├── Def_hyp_corr.m │ ├── Predict.m │ └── Train.m ├── @Kriging │ ├── Clean.m │ ├── Kriging.m │ ├── Predict.m │ └── Train.m ├── @Metamodel │ ├── Mape.m │ ├── Metamodel.m │ ├── Nmse.m │ ├── Plot.m │ ├── Plot_error.m │ ├── Predict.m │ ├── R2.m │ ├── Raae.m │ ├── Rmae.m │ ├── Rmse.m │ ├── Train.m │ └── Update_training_data.m ├── @Q_kriging │ ├── Build_tau.m │ ├── Clean.m │ ├── Ext_corr.m │ ├── Fit.m │ ├── Generate_degrees.m │ ├── Info_display.m │ ├── Init_tau_id.m │ ├── Int_cov.m │ ├── Likelihood.m │ ├── Marginal_likelihood.m │ ├── Nrmse.m │ ├── Plot.m │ ├── Predict.m │ ├── Predict_comp.m │ ├── Predict_deriv.m │ ├── Preproc_tau.m │ ├── Q_kriging.m │ ├── Regression_matrix.m │ ├── Set_data.m │ ├── Set_optim.m │ ├── Tau_wrap.m │ ├── Train.m │ ├── Tune_parameters.m │ ├── Update_model.m │ ├── Update_reg.m │ └── Update_sto.m └── @RBF │ ├── Clean.m │ ├── Corrcubic.m │ ├── Corrgauss.m │ ├── Corrinvmultiquadric.m │ ├── Corrlinear.m │ ├── Corrmatern32.m │ ├── Corrmatern52.m │ ├── Corrmultiquadric.m │ ├── Corrthinplatespline.m │ ├── Loo_error.m │ ├── Norm_theta.m │ ├── Predict.m │ ├── RBF.m │ └── Train.m ├── Optimization ├── @MGDA │ ├── Eval_function.m │ ├── Get_grad_shared.m │ ├── Get_x_finite_diff.m │ ├── MGDA.m │ ├── Opt.m │ └── Permute_info.m └── @NSGA_2 │ ├── Croisement.m │ ├── Crowding_distance_assignement.m │ ├── Crowding_distance_sorting.m │ ├── Domination_sorting.m │ ├── Fast_non_dominated_sort.m │ ├── Mutation.m │ ├── NSGA_2.m │ ├── Opt.m │ ├── Restart.m │ └── Selection.m ├── Problem ├── @Problem │ ├── Add_data.m │ ├── Eval.m │ ├── Eval_error_handling.m │ ├── Get_design.m │ ├── Input_assert.m │ ├── Output_assert.m │ ├── Problem.m │ └── Sampling.m ├── @Problem_multifi │ ├── Check_pb.m │ ├── Eval.m │ ├── Get_design.m │ └── Problem_multifi.m └── @Q_problem │ ├── Add_data.m │ ├── Eval.m │ ├── Eval_error_handling.m │ ├── Get_design.m │ ├── Input_assert.m │ ├── Q_problem.m │ ├── SLHS.m │ └── Sampling.m ├── Sensitivity_analysis └── @Sobol │ ├── Logical_matrix.m │ ├── Partial_comp.m │ ├── Sampling.m │ ├── Sobol.m │ └── Total_comp.m ├── Tools ├── BQQV │ ├── Q_basis_functions │ │ ├── Q_Exp.m │ │ ├── Q_Gauss.m │ │ ├── Q_Matern3_2.m │ │ └── Q_Matern5_2.m │ └── Tools │ │ ├── Q_pareto_points.m │ │ ├── hypersph2cart.m │ │ ├── ind2subVect.m │ │ ├── inter_min.m │ │ ├── q2qval.m │ │ ├── qval2qsub.m │ │ ├── reg_polytopes.m │ │ ├── sub2indVect.m │ │ ├── sub2tril.m │ │ ├── vect2full.m │ │ └── vect2tril.m ├── Kmeans_SBDO.m ├── Kmeans_analysis_SBDO.m ├── Pareto_points.m ├── Rand_gauss.m ├── Round_data.m ├── Scale_data.m ├── Spacing_PF.m ├── Spread_PF.m ├── Theta_bound.m ├── Theta_bound_normal.m ├── Unscale_data.m └── copy.m └── Unit_test ├── Test_metamodel.m ├── Test_problem.m ├── Test_problem_multifi.m ├── Test_q_kriging.m └── Test_q_problem.m /.gitignore: -------------------------------------------------------------------------------- 1 | # Temporary files MATLAB 2 | *~ 3 | *.asv 4 | 5 | # Data files MATLAB 6 | *.mat 7 | 8 | # PDF 9 | *.pdf 10 | 11 | # Image 12 | *.eps 13 | *.png 14 | *.jpg 15 | *.jpeg 16 | -------------------------------------------------------------------------------- /Doc/Manuscrit_durantin_final.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Doc/Manuscrit_durantin_final.pdf -------------------------------------------------------------------------------- /Doc/Soutenance_Durantin_v2.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Doc/Soutenance_Durantin_v2.pdf -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/Full_fractionnal.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | % Grille de points 8 | [x1 , x2] = meshgrid(linspace(1.25,8.75,5)); 9 | 10 | x1 = reshape(x1,25,1); 11 | x2 = reshape(x2,25,1); 12 | 13 | figure 14 | hold on 15 | plot(x1,x2,'bo','MarkerFaceColor','b') 16 | axis([0 10 0 10]) 17 | axis square 18 | box on 19 | xlabel('$x_1$','interpreter','latex') 20 | ylabel('$x_2$','interpreter','latex') 21 | hold off -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/LHS_bad.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | figure 8 | hold on 9 | plot(0.5:9.5,0.5:9.5,'bo','MarkerFaceColor','b') 10 | axis([0 10 0 10]) 11 | ax = gca; 12 | ax.XTick = 0:10; 13 | ax.YTick = 0:10; 14 | axis square 15 | box on 16 | grid on 17 | xlabel('$x_1$','interpreter','latex') 18 | ylabel('$x_2$','interpreter','latex') 19 | hold off -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/OLHS.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | % OLHS 8 | x_temp = stk_sampling_maximinlhs( 25, 2, [0 0; 10 10], 500 ); 9 | x = x_temp.data; 10 | 11 | figure 12 | hold on 13 | plot(x(:,1),x(:,2),'bo','MarkerFaceColor','b') 14 | axis([0 10 0 10]) 15 | axis square 16 | box on 17 | xlabel('$x_1$','interpreter','latex') 18 | ylabel('$x_2$','interpreter','latex') 19 | hold off -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/Random_DOE.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/Random_DOE.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.1.Plan_experience/Sobol_kmeans.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | % Sobol 8 | x_temp = stk_sampling_sobol( 10^5, 2, [0 0; 10 10], true ); 9 | x_sampling = x_temp.data; 10 | 11 | % k means 12 | [~,x] = kmeans( x_sampling, 25, 'MaxIter', 500); 13 | 14 | % Avec le kmeans de SBDOT : 15 | % x = Kmeans_SBDO( x_sampling, 25); 16 | 17 | figure 18 | hold on 19 | plot(x(:,1),x(:,2),'bo','MarkerFaceColor','b') 20 | axis([0 10 0 10]) 21 | axis square 22 | box on 23 | xlabel('$x_1$','interpreter','latex') 24 | ylabel('$x_2$','interpreter','latex') 25 | hold off -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Bruit_blanc_GP.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.2.Kriging/Bruit_blanc_GP.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Exp_GP.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.2.Kriging/Exp_GP.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Gauss_GP.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.2.Kriging/Gauss_GP.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Intervalle_confiance.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Set random seed for reproductibility 6 | rng(1) 7 | % Define problem structure 8 | n_x = 1; 9 | m_y = 1; 10 | m_g = 0; 11 | lb = 0; 12 | ub = 10; 13 | 14 | % Create Problem object 15 | prob = Problem( @(x)x.*sin(x), n_x, m_y, m_g, lb, ub , 'parallel', true); 16 | 17 | % Evaluate the model 18 | prob.Get_design( 6 ,'LHS' ) 19 | 20 | % Kriging 21 | krig1 = Kriging ( prob , 1 , [] ); 22 | 23 | % Variance prediction 24 | x_pred = linspace(0,10,200)'; 25 | 26 | [y_pred1, s_pred1] = krig1.Predict( x_pred ); 27 | 28 | % Confidence interval 29 | confidence_interval = [y_pred1 - 1.96*sqrt(s_pred1) ; flipud(y_pred1 + 1.96*sqrt(s_pred1))]; 30 | 31 | % Figure 32 | figure 33 | hold on 34 | plot( prob.x, prob.y, 'bo') 35 | plot( x_pred, y_pred1, 'k-') 36 | fill( [x_pred ; flipud(x_pred)], confidence_interval, 'g','FaceAlpha',0.3,'EdgeColor','none') 37 | box on 38 | xlabel('$x$','interpreter','latex') 39 | ylabel('$\mathcal{M}$','interpreter','latex') 40 | legend({'Donn\''{e}es d''entrainement $\mathcal{D}$','$\hat y$','Intervalle de confiance \`{a} 95\%'},... 41 | 'Interpreter','latex','Location','northwest') 42 | axis([0 10 -8 10]) 43 | hold off 44 | 45 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Matern32_GP.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.2.Kriging/Matern32_GP.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.2.Kriging/Matern52_GP.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.2.Kriging/Matern52_GP.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.3.RBF/Cubique_RBF.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | theta1 = 0.1 ; % Hyperparametre 1 8 | theta2 = 0.2 ; % Hyperparametre 2 9 | theta3 = 0.3 ; % Hyperparametre 3 10 | 11 | % Fonction radiale cubique 12 | x_R = linspace(0,2,400)'; 13 | R1 = (x_R.*theta1).^3; 14 | R2 = (x_R.*theta2).^3; 15 | R3 = (x_R.*theta3).^3; 16 | 17 | figure 18 | hold on 19 | plot( x_R, R1, 'b-') 20 | plot( x_R, R2, 'g-') 21 | plot( x_R, R3, 'r-') 22 | box on 23 | xlabel('$(x -x'')$','interpreter','latex') 24 | ylabel('$\varphi_R({x}-{x}'')$','interpreter','latex') 25 | legend({'$\theta = 0.1$','$\theta = 0.2$','$\theta = 0.3$'},... 26 | 'interpreter','latex','Location','northwest') 27 | title('Cubique','interpreter','latex') 28 | hold off 29 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.3.RBF/Inv_multiquadrique_RBF.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | theta1 = 1 ; % Hyperparametre 1 8 | theta2 = 1.5 ; % Hyperparametre 2 9 | theta3 = 2 ; % Hyperparametre 3 10 | 11 | % Fonction radiale inverse multiquadrique 12 | x_R = linspace(0,1,400)'; 13 | R1 = 1./sqrt((1/theta1).^2+x_R.^2); 14 | R2 = 1./sqrt((1/theta2).^2+x_R.^2); 15 | R3 = 1./sqrt((1/theta3).^2+x_R.^2); 16 | 17 | figure 18 | hold on 19 | plot( x_R, R1, 'b-') 20 | plot( x_R, R2, 'g-') 21 | plot( x_R, R3, 'r-') 22 | box on 23 | xlabel('$(x -x'')$','interpreter','latex') 24 | ylabel('$\varphi_R({x}-{x}'')$','interpreter','latex') 25 | legend({'$\theta = 1$','$\theta = 1.5$','$\theta = 2$'},... 26 | 'interpreter','latex','Location','northeast') 27 | title('Inverse multiquadrique','interpreter','latex') 28 | hold off 29 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.3.RBF/Lineaire_RBF.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/B.3.RBF/Lineaire_RBF.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.3.RBF/Multiquadrique_RBF.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | theta1 = 1 ; % Hyperparametre 1 8 | theta2 = 1.5 ; % Hyperparametre 2 9 | theta3 = 2 ; % Hyperparametre 3 10 | 11 | % Fonction radiale multiquadrique 12 | x_R = linspace(0,1,400)'; 13 | R1 = sqrt((1/theta1).^2+x_R.^2); 14 | R2 = sqrt((1/theta2).^2+x_R.^2); 15 | R3 = sqrt((1/theta3).^2+x_R.^2); 16 | 17 | figure 18 | hold on 19 | plot( x_R, R1, 'b-') 20 | plot( x_R, R2, 'g-') 21 | plot( x_R, R3, 'r-') 22 | box on 23 | xlabel('$(x -x'')$','interpreter','latex') 24 | ylabel('$\varphi_R({x}-{x}'')$','interpreter','latex') 25 | legend({'$\theta = 1$','$\theta = 1.5$','$\theta = 2$'},... 26 | 'interpreter','latex','Location','northwest') 27 | title('Multiquadrique','interpreter','latex') 28 | hold off 29 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/B.3.RBF/ThinPlateSpline_RBF.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | rng(1) 6 | 7 | theta1 = 1 ; % Hyperparametre 1 8 | theta2 = 1.5 ; % Hyperparametre 2 9 | theta3 = 2 ; % Hyperparametre 3 10 | 11 | % Fonction radiale spline en plaque mince 12 | x_R = linspace(0,1,400)'; 13 | R1 = ((x_R.*theta1).^2).*log(x_R.*theta1); 14 | R2 = ((x_R.*theta2).^2).*log(x_R.*theta2); 15 | R3 = ((x_R.*theta3).^2).*log(x_R.*theta3); 16 | 17 | figure 18 | hold on 19 | plot( x_R, R1, 'b-') 20 | plot( x_R, R2, 'g-') 21 | plot( x_R, R3, 'r-') 22 | box on 23 | xlabel('$(x -x'')$','interpreter','latex') 24 | ylabel('$\varphi_R({x}-{x}'')$','interpreter','latex') 25 | legend({'$\theta = 1$','$\theta = 1.5$','$\theta = 2$'},... 26 | 'interpreter','latex','Location','northwest') 27 | hold off 28 | title('Spline en plaque mince','interpreter','latex') 29 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_FLNS_20_1.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_FLNS_20_1.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_FLNS_80.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_FLNS_80.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_Kreuzer_20_1.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_Kreuzer_20_1.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_Kreuzer_300.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/DOE_Kreuzer_300.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Plot_error_FLNS.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Plot_error_FLNS.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Plot_error_kreuzer.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Plot_error_kreuzer.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_FLNS_20_RBF.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_FLNS_20_RBF.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_FLNS_20_kriging.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_FLNS_20_kriging.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_Kreuzer_80_RBF.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_Kreuzer_80_RBF.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_Kreuzer_80_kriging.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Cellule_PA/Results_Kreuzer_80_kriging.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/DOE_coupler_lambda_1.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/DOE_coupler_lambda_1.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/DOE_coupler_lambda_10.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/DOE_coupler_lambda_10.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Plot_error_lambda.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Load DOE training 6 | load('DOE_coupler_lambda_1.mat') 7 | 8 | % Create Problem 9 | problem=Problem(@(x)Coupleur(x),5,1,0,lb,ub); 10 | problem.Add_data(x,y,[]); 11 | 12 | % Build Kriging / RBF 13 | krig=Kriging(problem,1,[],... 14 | 'estim_hyp',@pseudoLikelihood,'corr','correxp','reg',false); 15 | 16 | rbf=RBF(problem,1,[],... 17 | 'optimizer','CMAES','corr','Corrlinear'); 18 | 19 | % Load DOE test points 20 | load('DOE_coupler_lambda_10.mat') 21 | 22 | krig.Plot_error(x,y(:,1)) 23 | title('Ecart du taux de couplage') 24 | 25 | rbf.Plot_error(x,y(:,1)) 26 | title('Ecart du taux de couplage') 27 | 28 | % Cut view param 1&2 29 | krig.Plot([1 2],[1.25 0.6 1.31]) 30 | hold on 31 | xlabel('$W_1$','interpreter','latex') 32 | ylabel('$W_2$','interpreter','latex') 33 | zlabel('$T_c$','interpreter','latex') 34 | hold off 35 | 36 | % Cut view param 3&4 37 | krig.Plot([3 4],[0.4 0.4 1.31]) 38 | hold on 39 | xlabel('$g_1$','interpreter','latex') 40 | ylabel('$g_2$','interpreter','latex') 41 | zlabel('$T_c$','interpreter','latex') 42 | hold off 43 | 44 | % Cut view param 4&5 45 | krig.Plot([4 5],[0.4 0.4 1.25]) 46 | hold on 47 | xlabel('$g_2$','interpreter','latex') 48 | ylabel('$\lambda_c$','interpreter','latex') 49 | zlabel('$T_c$','interpreter','latex') 50 | hold off 51 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Plot_validation_RBF.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Load result file 6 | load Results_lambda_RBF.mat 7 | 8 | % Figure options 9 | name_legende = {'Matern 52', 'Matern 32', 'Gaussienne', 'Lin\''eaire' , 'Spline plaque mince' , 'Multiquadrique' , 'Inverse multiquadrique', 'Cubique'}; 10 | positions = [1 1.25 1.5 1.75 2 2.25 2.5 2.75]; 11 | group = [1,2,3,4,5,6,7,8]; 12 | 13 | %% 14 | figure 15 | 16 | subplot(1,2,1) 17 | hold on 18 | boxplot(RAAE,group,'positions',positions,'Colors','r') 19 | set(gca,'xticklabel', name_legende ,'TickLabelInterpreter','latex') 20 | set(gca,'XTickLabelRotation', 45) 21 | ylabel('RAAE','interpreter','latex') 22 | hold off 23 | 24 | subplot(1,2,2) 25 | hold on 26 | boxplot(RMAE,group,'positions',positions,'Colors','r') 27 | set(gca,'xticklabel', name_legende ,'TickLabelInterpreter','latex') 28 | set(gca,'XTickLabelRotation', 45) 29 | ylabel('RMAE','interpreter','latex') 30 | hold off 31 | 32 | figtitle('Ecart du taux de couplage','interpreter','latex') 33 | set(gcf, 'Position', [0 0 1000 600]) -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Plot_validation_kriging.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Load result file 6 | load Results_lambda_kriging.mat 7 | 8 | % Figure options 9 | name_legende = {'Matern 52 ','Exponentielle','Gaussienne','Matern 32'}; 10 | positions = [1 1.25 1.5 1.75 2 2.25 2.5 2.75]; 11 | group = [1,2,3,4,5,6,7,8]; 12 | 13 | %% 14 | figure 15 | 16 | subplot(1,2,1) 17 | hold on 18 | boxplot(RAAE,group,'positions',positions,'Colors','rb') 19 | set(gca,'xtick',[1.125 1.625 2.125 2.625]) 20 | set(gca,'xticklabel', name_legende ,'TickLabelInterpreter','latex') 21 | set(gca,'XTickLabelRotation', 45) 22 | ylabel('RAAE','interpreter','latex') 23 | box_vars = findall(gca,'Tag','Box'); 24 | legend(box_vars([2,1]), {'MLE','LOO'},'Location','Northwest' ,'Interpreter','latex'); 25 | hold off 26 | 27 | subplot(1,2,2) 28 | hold on 29 | boxplot(RMAE,group,'positions',positions,'Colors','rb') 30 | set(gca,'xtick',[1.125 1.625 2.125 2.625]) 31 | set(gca,'xticklabel', name_legende ,'TickLabelInterpreter','latex') 32 | set(gca,'XTickLabelRotation', 45) 33 | ylabel('RMAE','interpreter','latex') 34 | hold off 35 | 36 | figtitle('Ecart du taux de couplage','interpreter','latex') 37 | set(gcf, 'Position', [0 0 1000 600]) -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Results_lambda_RBF.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Results_lambda_RBF.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Results_lambda_kriging.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Coupleur/Results_lambda_kriging.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Bencharmk_RBF_hotplate.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Bencharmk_RBF_hotplate.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Benchmark_kriging_cellulePA_FLNS.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Benchmark_kriging_cellulePA_FLNS.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Benchmark_kriging_coupleur.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Exemple_non_fonctionnel/Benchmark_kriging_coupleur.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/DOE_setBias_100_CVT_1.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/DOE_setBias_100_CVT_1.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/DOE_setBias_300_CVT_test.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/DOE_setBias_300_CVT_test.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Plot_error_prediction.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Plot_error_prediction.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Results_setBias_100_CVT_RBF.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Results_setBias_100_CVT_RBF.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Results_setBias_100_CVT_kriging.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_2/D.Benchmark/Hotplate/Results_setBias_100_CVT_kriging.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/A/Branin_IMSE_krig.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Set random seed for reproductibility 6 | rng(1) 7 | 8 | % Define problem structure 9 | m_x = 2; % Number of parameters 10 | m_y = 1; % Number of objectives 11 | m_g = 1; % Number of constraints 12 | lb = [-5 0]; % Lower bound 13 | ub = [10 15]; % Upper bound 14 | 15 | % Create Problem object with optionnal parallel input as true 16 | prob_krig = Problem( 'Branin', m_x, m_y, m_g, lb, ub , 'parallel', true); 17 | 18 | % Evaluate the model on 20 points created with LHS 19 | prob_krig.Get_design( 20 ,'LHS' ); 20 | 21 | %% 22 | options_optim.DispModulo = 0; 23 | options_optim.Restarts = 0; 24 | 25 | % Adaptive sampling 26 | obj1 = Error_prediction(prob_krig, 1, [], @Kriging ,'iter_max',10,'crit_type','IMSE','options_optim',options_optim); 27 | 28 | figure 29 | hold on 30 | plot(prob_krig.x(1:20,1),prob_krig.x(1:20,2),'bo','MarkerFaceColor','b','MarkerSize',8) 31 | plot(prob_krig.x(21:end,1),prob_krig.x(21:end,2),'ro','MarkerFaceColor','r','MarkerSize',8) 32 | xlabel('$x_1$','interpreter','latex') 33 | ylabel('$x_2$','interpreter','latex') 34 | axis([-5 10 0 15]) 35 | box on 36 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/A/Branin_MSE_krig.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Set random seed for reproductibility 6 | rng(1) 7 | 8 | % Define problem structure 9 | m_x = 2; % Number of parameters 10 | m_y = 1; % Number of objectives 11 | m_g = 1; % Number of constraints 12 | lb = [-5 0]; % Lower bound 13 | ub = [10 15]; % Upper bound 14 | 15 | % Create Problem object with optionnal parallel input as true 16 | prob_krig = Problem( 'Branin', m_x, m_y, m_g, lb, ub , 'parallel', true); 17 | 18 | % Evaluate the model on 20 points created with LHS 19 | prob_krig.Get_design( 20 ,'LHS' ); 20 | 21 | %% 22 | options_optim.DispModulo = 0; 23 | options_optim.Restarts = 0; 24 | 25 | % Adaptive sampling 26 | obj1 = Error_prediction(prob_krig, 1, [], @Kriging ,'iter_max',10,'crit_type','MSE','options_optim',options_optim); 27 | 28 | figure 29 | hold on 30 | plot(prob_krig.x(1:20,1),prob_krig.x(1:20,2),'bo','MarkerFaceColor','b','MarkerSize',8) 31 | plot(prob_krig.x(21:end,1),prob_krig.x(21:end,2),'ro','MarkerFaceColor','r','MarkerSize',8) 32 | xlabel('$x_1$','interpreter','latex') 33 | ylabel('$x_2$','interpreter','latex') 34 | axis([-5 10 0 15]) 35 | box on -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/EI_times_PF_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/EI_times_PF_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/EI_versus_PF_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/EI_versus_PF_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Gutmann_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Gutmann_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Opt_Camel_real.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Opt_Camel_real.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Plot_camel.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Plot Camel contour lines 6 | 7 | [ X_plot1, X_plot2 ] = meshgrid( ... 8 | linspace( -2, 2, 500), ... 9 | linspace( -2, 2, 500) ); 10 | 11 | X_plot = [ reshape( X_plot1, size( X_plot1, 1)^2, 1)... 12 | reshape( X_plot2, size( X_plot2, 1 )^2, 1) ]; 13 | 14 | reshape_size = size(X_plot1); 15 | 16 | [Y_plot , G_plot] = Camel_constrained( X_plot ); 17 | 18 | Y_plot(G_plot >0 ,:) = NaN; 19 | 20 | figure 21 | hold on 22 | contourf( X_plot1, X_plot2, reshape( Y_plot, reshape_size ), 21 ); 23 | contour( X_plot1, X_plot2, reshape( G_plot, reshape_size ), [0,0] ,'r' ); 24 | plot(0.02,0.7,'ro','MarkerFaceColor','r') 25 | xlabel('$x_1$','interpreter','latex') 26 | ylabel('$x_2$','interpreter','latex') 27 | colorbar 28 | colormap('jet') 29 | box on 30 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Test_points_Camel.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_constrained/Test_points_Camel.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/CORS_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/CORS_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/EI_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/EI_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Gutmann_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Gutmann_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Opt_Branin_real.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | options_optim.TolX = 1e-8; %Tolerance on final optimized input 6 | options_optim.TolFun = 1e-8; %Tolerance on final optimized output 7 | options_optim.DispModulo = 0; % Display message or not during optimization 8 | options_optim.Restarts = 4; % Number of restar of the algorithm 9 | options_optim.LBounds = [-5 0]; 10 | options_optim.UBounds = [10 15]; 11 | options_optim.IncPopSize = 4; 12 | 13 | [ x_min_real , y_min_real ] = cmaes(@Branin, [-3 12], [7 7]',options_optim); 14 | 15 | save('Opt_Branin_real','x_min_real','y_min_real') 16 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Opt_Branin_real.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Opt_Branin_real.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/PI_bench.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | for i = 1 : 20 6 | 7 | % Set random seed for reproductibility 8 | rng(i) 9 | 10 | % Define problem structure 11 | m_x = 2; % Number of parameters 12 | m_y = 1; % Number of objectives 13 | m_g = 1; % Number of constraint 14 | lb = [-5 0]; % Lower bound 15 | ub = [10 15]; % Upper bound 16 | 17 | % Create Problem object with optionnal parallel input as true 18 | prob = Problem( 'Branin', m_x, m_y, m_g, lb, ub , 'parallel', true); 19 | 20 | % Get design 21 | prob.Get_design( 20 ,'LHS' ) 22 | 23 | % Instantiate optimization object 24 | obj = Expected_improvement(prob, 1, [], @Kriging, 'fcall_max', 40 ,'criterion', 'PI'); 25 | 26 | % Launch optimization 27 | obj.Opt(); 28 | 29 | load Test_points_Branin.mat 30 | 31 | RMSE_result(i,:) = obj.meta_y.Rmse(x_test,y_test); 32 | 33 | failure(i,:) = obj.failed; 34 | x_min_result(i,:) = obj.x_min; 35 | y_min_result(i,:) = obj.y_min; 36 | fcall_result(i,:) = obj.fcall_num; 37 | 38 | end 39 | 40 | save('PI_result','RMSE_result','failure','x_min_result','y_min_result','fcall_result') 41 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/PI_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/PI_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Plot_Branin.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | save_plot = true; 6 | s = hgexport('readstyle','manuscrit'); 7 | s.Format = 'png'; 8 | 9 | [ X_plot1, X_plot2 ] = meshgrid( ... 10 | linspace( -5,10, 200), ... 11 | linspace( 0, 15, 200) ); 12 | 13 | X_plot = [ reshape( X_plot1, size( X_plot1, 1)^2, 1)... 14 | reshape( X_plot2, size( X_plot2, 1 )^2, 1) ]; 15 | 16 | reshape_size = size(X_plot1); 17 | 18 | Y_plot = Branin( X_plot ); 19 | 20 | X = [pi 2.275 ; 9.42478 2.475]; 21 | 22 | figure 23 | hold on 24 | contourf( X_plot1, X_plot2, reshape( Y_plot, reshape_size ), 21 ); 25 | plot(X(:,1), X(:,2), 'wx', 'MarkerSize', 14); 26 | plot(-3.18,12.36,'ro','MarkerFaceColor','r') 27 | xlabel('$x_1$','interpreter','latex') 28 | ylabel('$x_2$','interpreter','latex') 29 | colorbar 30 | colormap('jet') 31 | box on 32 | 33 | if save_plot 34 | hgexport(gcf,'Branin_contour',s); 35 | end 36 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Test_points_Branin.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/B/Bench_optim_unconstrained/Test_points_Branin.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Gradient/Plot_gradient_result.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Set random seed for reproductibility 6 | rng(1) 7 | 8 | % Optimization 9 | [ f_min1, x_min1 , iter1, f_call1,x_hist1]=steepest_slope_step_opt('square_2D',[-3 2],[-5 -5],[5 5]); 10 | 11 | [X1_plot,X2_plot]=meshgrid(-5:0.01:5); 12 | X_plot=[reshape(X1_plot,size(X1_plot,1)^2,1),reshape(X2_plot,size(X1_plot,1)^2,1)]; 13 | y_plot=square_2D(X_plot); 14 | y_plot=reshape(y_plot,size(X1_plot,1),size(X1_plot,1)); 15 | 16 | figure 17 | hold on 18 | plot(x_hist1(:,1),x_hist1(:,2),'ko','MarkerFaceColor','k','MarkerSize',8) 19 | contour(X1_plot,X2_plot,y_plot,[600 500 400 300 200 100 50 20 15 5 1 0.1]) 20 | plot(-3,2,'ko','MarkerFaceColor','k','MarkerSize',8') 21 | xlabel('$x_1$','interpreter','latex') 22 | ylabel('$x_2$','interpreter','latex') 23 | plot([-3;x_hist1(:,1)],[2;x_hist1(:,2)],'k-') 24 | legend({'$\bf{x}^{k}$'},... 25 | 'Interpreter','latex','Location','Northeast') 26 | box on 27 | colorbar 28 | hold off 29 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Gradient/pas_opt.m: -------------------------------------------------------------------------------- 1 | function objectif=pas_opt(x,functname,x0,pas) 2 | 3 | 4 | objectif=feval(functname,x0+x.*pas); 5 | 6 | end -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/B/Gradient/square_2D.m: -------------------------------------------------------------------------------- 1 | function y=square_2D(x) 2 | 3 | y=(x(:,1).^4)/4+x(:,2).^4; 4 | 5 | end -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/EHVI_opt.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | save_plot = true; 6 | s = hgexport('readstyle','manuscrit'); 7 | s.Format = 'png'; 8 | 9 | for i = 1 : 20 10 | 11 | % Set random seed for reproductibility 12 | rng(i) 13 | 14 | % Define problem structure 15 | m_x = 2; % Number of parameters 16 | m_y = 2; % Number of objectives 17 | m_g = 0; % Number of constraint 18 | lb = [0 0]; % Lower bound 19 | ub = [2 2]; % Upper bound 20 | 21 | % Create Problem object with optionnal parallel input as true 22 | prob = Problem( 'MO_convex', m_x, m_y, m_g, lb, ub , 'parallel', true); 23 | 24 | % Evaluate the model on 20 points created with LHS 25 | prob.Get_design( 20 ,'LHS' ) 26 | 27 | % Instantiate optimization object 28 | EGO = Multi_obj_EHVI( prob, [1 2], [], @Kriging, 5, ... 29 | 'fcall_max',40 ,'criterion', 'EHVI' ); 30 | 31 | pareto_final = Pareto_points( prob.y(:,[1 2]) ); 32 | 33 | spread_meas(i,:) = Spread_PF(pareto_final); 34 | spacing_meas(i,:) = Spacing_PF(pareto_final); 35 | hyp_meas(i,:) = stk_dominatedhv (pareto_final, [10 10]); 36 | nbr_pareto(i,:) = size(pareto_final,1); 37 | 38 | end 39 | 40 | save('EHVI_results','spread_meas','spacing_meas','hyp_meas','nbr_pareto') 41 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/EHVI_results.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/C/bench_MO/EHVI_results.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/EI_euclid_opt.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | save_plot = true; 6 | s = hgexport('readstyle','manuscrit'); 7 | s.Format = 'png'; 8 | 9 | for i = 1 : 20 10 | 11 | % Set random seed for reproductibility 12 | rng(i) 13 | 14 | % Define problem structure 15 | m_x = 2; % Number of parameters 16 | m_y = 2; % Number of objectives 17 | m_g = 0; % Number of constraint 18 | lb = [0 0]; % Lower bound 19 | ub = [2 2]; % Upper bound 20 | 21 | % Create Problem object with optionnal parallel input as true 22 | prob = Problem( 'MO_convex', m_x, m_y, m_g, lb, ub , 'parallel', true); 23 | 24 | % Evaluate the model on 20 points created with LHS 25 | prob.Get_design( 20 ,'LHS' ) 26 | 27 | % Instantiate optimization object 28 | EGO = Multi_obj_EHVI( prob, [1 2], [], @Kriging, 5, ... 29 | 'fcall_max',40 ,'criterion', 'EI_euclid' ); 30 | 31 | pareto_final = Pareto_points( prob.y(:,[1 2]) ); 32 | 33 | spread_meas(i,:) = Spread_PF(pareto_final); 34 | spacing_meas(i,:) = Spacing_PF(pareto_final); 35 | hyp_meas(i,:) = stk_dominatedhv (pareto_final, [10 10]); 36 | nbr_pareto(i,:) = size(pareto_final,1); 37 | 38 | end 39 | 40 | save('EI_euclid_results','spread_meas','spacing_meas','hyp_meas','nbr_pareto') 41 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/EI_euclid_results.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/C/bench_MO/EI_euclid_results.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/NSGA_opt.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | for i = 1 : 20 6 | 7 | % Set random seed for reproductibility 8 | rng(i) 9 | 10 | % Define problem structure 11 | m_x = 2; % Number of parameters 12 | m_y = 2; % Number of objectives 13 | m_g = 0; % Number of constraint 14 | lb = [0 0]; % Lower bound 15 | ub = [2 2]; % Upper bound 16 | 17 | % Instantiate optimization object 18 | nsga = NSGA_2('MO_convex_nsga',lb,ub,2,'n_pop',50,'max_gen',100,'display',true); 19 | 20 | pareto_final = Pareto_points( nsga.y(:,[1 2]) ); 21 | 22 | spread_meas(i,:) = Spread_PF(pareto_final); 23 | spacing_meas(i,:) = Spacing_PF(pareto_final); 24 | hyp_meas(i,:) = stk_dominatedhv (pareto_final, [10 10]); 25 | nbr_pareto(i,:) = size(pareto_final,1); 26 | 27 | end 28 | 29 | save('NSGA_results','spread_meas','spacing_meas','hyp_meas','nbr_pareto') 30 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/C/bench_MO/NSGA_results.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/C/bench_MO/NSGA_results.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_3/D/Exemple_robust_measure.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_3/D/Exemple_robust_measure.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/DOE_borehole_500.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/DOE_borehole_500.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/DOE_curretal_200.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/DOE_curretal_200.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_HF.m: -------------------------------------------------------------------------------- 1 | function ye = borehole_HF(xe) 2 | 3 | rw = xe(:,1); 4 | r = xe(:,2); 5 | Tu = xe(:,3); 6 | Hu = xe(:,4); 7 | Tl = xe(:,5); 8 | Hl = xe(:,6); 9 | L = xe(:,7); 10 | Kw = xe(:,8); 11 | 12 | frac1 = 2 * pi .* Tu .* (Hu-Hl); 13 | 14 | frac2a = 2.*L.*Tu ./ (log(r./rw).*rw.^2.*Kw); 15 | frac2b = Tu ./ Tl; 16 | frac2 = log(r./rw) .* (1+frac2a+frac2b); 17 | 18 | ye = frac1 ./ frac2; 19 | 20 | end 21 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_LF.m: -------------------------------------------------------------------------------- 1 | function yc = borehole_LF(xc) 2 | 3 | rw = xc(:,1); 4 | r = xc(:,2); 5 | Tu = xc(:,3); 6 | Hu = xc(:,4); 7 | Tl = xc(:,5); 8 | Hl = xc(:,6); 9 | L = xc(:,7); 10 | Kw = xc(:,8); 11 | 12 | 13 | frac1 = 5 * pi .* Tu .* (Hu-Hl); 14 | 15 | frac2a = 2.*L.*Tu ./ (log(r./rw).*rw.^2.*Kw); 16 | frac2b = Tu ./ Tl; 17 | frac2 = log(r./rw) .* (1.5+frac2a+frac2b); 18 | 19 | yc = frac1 ./ frac2; 20 | 21 | end 22 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_cokrig.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_cokrig.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_cokrig_LOO.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_cokrig_LOO.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_corbf.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/borehole_corbf.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/coRBF_borehole_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/coRBF_borehole_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/coRBF_curretal_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/coRBF_curretal_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_borehole_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_borehole_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_borehole_result_LOO.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_borehole_result_LOO.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_curretal_result.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/cokrig_curretal_result.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal88exp.m: -------------------------------------------------------------------------------- 1 | 2 | 3 | function [y] = curretal88exp(xx) 4 | 5 | x1 = xx(:,1); 6 | x2 = xx(:,2); 7 | 8 | fact1 = 1 - exp(-1./(2.*x2)); 9 | fact2 = 2300*x1.^3 + 1900*x1.^2 + 2092*x1 + 60; 10 | fact3 = 100*x1.^3 + 500*x1.^2 + 4*x1 + 20; 11 | 12 | y = fact1 .* fact2 ./fact3; 13 | 14 | end 15 | 16 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_HF.m: -------------------------------------------------------------------------------- 1 | function ye = curretal_HF(xe) 2 | 3 | 4 | x1 = xe(:,1); 5 | x2 = xe(:,2); 6 | 7 | fact1 = 1 - exp(-1./(2*x2)); 8 | fact2 = 2300*x1.^3 + 1900*x1.^2 + 2092*x1 + 60; 9 | fact3 = 100*x1.^3 + 500*x1.^2 + 4*x1 + 20; 10 | 11 | ye = fact1 .* fact2./fact3; 12 | 13 | 14 | end 15 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_LF.m: -------------------------------------------------------------------------------- 1 | function yc = curretal_LF(xc) 2 | 3 | x1 = xc(:,1); 4 | x2 = xc(:,2); 5 | 6 | maxarg = max([zeros(length(x1),1), x2-1./20],2); 7 | 8 | yh1 = curretal88exp([x1+1/20, x2+1/20]); 9 | yh2 = curretal88exp([x1+1/20, maxarg]); 10 | yh3 = curretal88exp([x1-1/20, x2+1/20]); 11 | yh4 = curretal88exp([x1-1/20, maxarg]); 12 | 13 | 14 | yc = (yh1 + yh2 + yh3 + yh4) / 4; 15 | 16 | end 17 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_cokrig.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_cokrig.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_corbf.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/Bench_coRBF/curretal_corbf.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/Bench_coRBF/plot_borehole_cokriging_LOO.m: -------------------------------------------------------------------------------- 1 | clear all 2 | close all 3 | clc 4 | 5 | % Load results 6 | load cokrig_borehole_result_LOO.mat 7 | mean_raae_cokrig(:,:) = mean(RAAE_borehole,3); 8 | std_raae_cokrig(:,:) = std(RAAE_borehole,[],3); 9 | 10 | % Error figures 11 | figure 12 | imagesc([40 240],[10 130],mean_raae_cokrig) 13 | axis xy 14 | set(gca,'XTick',[40 80 120 160 200 240],'YTick',[10 30 50 70 90 110 130]); 15 | colorbar 16 | %colorbar('Limits',[0.002 1]); 17 | %caxis([0.002 1]) 18 | colormap gray 19 | xlabel('Nombre de points LF','interpreter','latex') 20 | ylabel('Nombre de points HF','interpreter','latex') 21 | title('Moyenne','interpreter','latex') 22 | 23 | figure 24 | imagesc([40 240],[10 130],std_raae_cokrig) 25 | axis xy 26 | set(gca,'XTick',[40 80 120 160 200 240],'YTick',[10 30 50 70 90 110 130]); 27 | colorbar 28 | %colorbar('Limits',[5e-4 0.30]); 29 | %caxis([5e-4 0.30]) 30 | colormap gray 31 | xlabel('Nombre de points LF','interpreter','latex') 32 | ylabel('Nombre de points HF','interpreter','latex') 33 | title('Ecart type','interpreter','latex') 34 | 35 | 36 | -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EHVI_opt.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EHVI_opt.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EHVI_results.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EHVI_results.mat -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EI_MGDA_opt.m: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EI_MGDA_opt.m -------------------------------------------------------------------------------- /Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EI_MGDA_results.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Example/Thesis_examples/Chapitre_4/bench_EI_MGDA/EI_MGDA_results.mat -------------------------------------------------------------------------------- /IMPORTANT_LICENSE_NOTICE: -------------------------------------------------------------------------------- 1 | IMPORTANT LICENSE NOTICE 2 | ------------------------ 3 | 4 | SBDOT relies on the following packages or pieces of software : 5 | - cholincsp http://theoval.cmp.uea.ac.uk/matlab/ 6 | - cmaes https://www.lri.fr/~hansen/cmaes.m 7 | - ipdm https://fr.mathworks.com/matlabcentral/fileexchange/18937-ipdm--inter-point-distance-matrix?focused=5243203&tab=function 8 | - ooDACE http://www.sumo.intec.ugent.be/ooDACE 9 | - SQP Lab https://who.rocq.inria.fr/Jean-Charles.Gilbert/modulopt/optimization-routines/sqplab/sqplab.html 10 | - STK http://kriging.sourceforge.net/htmldoc/ 11 | 12 | cmaes and ooDACE have been adapted to fit our needs. 13 | 14 | Please note that, due to the License Terms of the ooDACE toolbox, RESTRICTIONS APPLY to commercial use or funded academic research. 15 | Conditions can be found on the ooDACE toolbox website http://sumo.intec.ugent.be/ooDACE. 16 | -------------------------------------------------------------------------------- /Lib/+ooDACE/@CmaesOptimizer/CmaesOptimizer.m: -------------------------------------------------------------------------------- 1 | classdef CmaesOptimizer < ooDACE.Optimizer 2 | 3 | % private members 4 | properties (SetAccess = 'private', GetAccess = 'private') 5 | opt; 6 | end 7 | 8 | methods 9 | 10 | function this = CmaesOptimizer(varargin) 11 | % call superclass 12 | this = this@ooDACE.Optimizer(varargin{:}); 13 | 14 | % default case 15 | if(nargin == 1) 16 | config = varargin{1}; 17 | 18 | % Create custom options structure 19 | this.opt.TolX = 1e-4; 20 | this.opt.Restarts = 2; 21 | this.opt.IncPopSize = 2; 22 | this.opt.TolFun = 1e-4; 23 | this.opt.DispModulo = 0; 24 | 25 | this.debug = config.self.getBooleanOption('debug', false); 26 | elseif(nargin == 3) 27 | % First 2 are parsed by base class 28 | %nvar = varargin{1}; 29 | %nobj = varargin{2}; 30 | this.opt = varargin{3}; 31 | else 32 | error('Invalid number of arguments given'); 33 | end 34 | 35 | end % constructor 36 | 37 | [this, xmin, fmin] = optimize(this, arg ); 38 | 39 | end % methods 40 | end % classdef 41 | -------------------------------------------------------------------------------- /Lib/+ooDACE/@CmaesOptimizer/optimize.m: -------------------------------------------------------------------------------- 1 | function [this, xmin, fmin] = optimize(this, arg ) 2 | 3 | if isa( arg, 'Model' ) 4 | func = @(x) evaluate(arg,x); 5 | else % assume function handle 6 | func = arg; 7 | end 8 | 9 | [LB, UB] = this.getBounds(); 10 | pop = this.getInitialPopulation(); 11 | this.opt.LBounds = LB; 12 | this.opt.UBounds = UB; 13 | 14 | % Actually run the the optimization routine 15 | [xmin,fmin] = cmaes( func, pop, [], this.opt ); 16 | 17 | end 18 | -------------------------------------------------------------------------------- /Lib/license_ipdm.txt: -------------------------------------------------------------------------------- 1 | Copyright (c) 2009, John D'Errico 2 | All rights reserved. 3 | 4 | Redistribution and use in source and binary forms, with or without 5 | modification, are permitted provided that the following conditions are 6 | met: 7 | 8 | * Redistributions of source code must retain the above copyright 9 | notice, this list of conditions and the following disclaimer. 10 | * Redistributions in binary form must reproduce the above copyright 11 | notice, this list of conditions and the following disclaimer in 12 | the documentation and/or other materials provided with the distribution 13 | 14 | THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 15 | AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 16 | IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 17 | ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE 18 | LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 19 | CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 20 | SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 21 | INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 22 | CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 23 | ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 24 | POSSIBILITY OF SUCH DAMAGE. -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/example/README: -------------------------------------------------------------------------------- 1 | Hanging chain example 2 | ^^^^^^^^^^^^^^^^^^^^^ 3 | The example 'hanging_chain' is provided to show how 'sqplab' can be 4 | used to solve an equality and inequality constrained problem. The 5 | problem consists in finding the static equilibrium position of a chain 6 | made of rigid bars. This problem is further explained and discussed in 7 | the book whose reference is given below. 8 | 9 | To run the example: 10 | - launch matlab 11 | - enter 'hanging_chain' in the execution window. 12 | 13 | The output of SQPlab should then be close to that in 14 | 'hanging_chain.res' and the generated figure should be close to the one 15 | in 'hanging_chain.eps'. 16 | 17 | Reference 18 | ^^^^^^^^^ 19 | @book{bonnans-gilbert-lemarechal-sagastizabal-2006, 20 | author={J.F. Bonnans and J.Ch. Gilbert and C. Lemar\'echal and C. 21 | Sagastiz\'abal}, 22 | title={Numerical Optimization -- Theoretical and Practical Aspects} # 23 | seconde_edition, 24 | series={Universitext}, 25 | publisher={Springer Verlag, Berlin}, 26 | year=2006 27 | } 28 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/example/hanging_chain_data.m: -------------------------------------------------------------------------------- 1 | function [n,nb,mi,me,xy0,a,b,lb,r,s,paxis] = hanging_chain_data (); 2 | 3 | % 4 | % [n,nb,mi,me,xy0,a,b,lb,r,s,paxis] = hanging_chain_data (); 5 | % 6 | % Returns data for the hanging chain problem. 7 | % 8 | % Same as data_1a, but with triple floor. 9 | 10 | xy0 = []; 11 | a = []; 12 | b = []; 13 | lb = []; 14 | r = []; 15 | s = []; 16 | 17 | % Initial position of the chain 18 | 19 | xy0 = [ 0.1; 0.3; 0.5; 0.7; 20 | -0.2; -0.4; -0.4; -0.5]; 21 | 22 | % Coordinates of the second hook 23 | 24 | a = 1; 25 | b = -0.1; 26 | 27 | % Desired length of the bars 28 | 29 | lb = [0.4 0.3 0.25 0.2 0.4]'; 30 | 31 | % Floor 32 | 33 | r = [-0.20; -0.40; -0.60]; 34 | s = [-0.60; -0.10; 0.20]; 35 | 36 | % Plot axis 37 | 38 | paxis = zeros(1,4); 39 | paxis(1) = -0.05; % xmin 40 | paxis(2) = 1.05; % xmax 41 | paxis(3) = -0.60; % ymin 42 | paxis(4) = 0.30; % ymax 43 | 44 | % Dimensions 45 | 46 | n = length(xy0); 47 | nb = 0; 48 | mi = length(r)*n/2; 49 | me = length(lb); 50 | 51 | end 52 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/example/hanging_chain_dimopt.m: -------------------------------------------------------------------------------- 1 | function [n,nb,mi,me] = hanging_chain_dimopt (); 2 | 3 | % [n,nb,mi,me] = hanging_chain_dimopt (); 4 | % 5 | % This function returns the dimensions of the problem. More specifically 6 | % . n = number of variables, 7 | % . nb = number of variables with bounds, 8 | % . mi = number of inequality constraints, 9 | % . me = number of equality constraints. 10 | 11 | % Read the data file 12 | 13 | [n,nb,mi,me,xy0,a,b,lb,r,s] = hanging_chain_data(); 14 | 15 | % Check dimensiions 16 | 17 | if mi ~= (length(lb)-1)*length(r) 18 | fprintf('(hanging_chain_dimopt) >>> corrupted data file ''%s''\n',pdat); 19 | fprintf('(hanging_chain_dimopt) >>> mi (=%0i) is not equal to (length(lb)-1)*length(r) (=%0i)\n', ... 20 | mi,(length(lb)-1)*length(r)); 21 | exit; 22 | end 23 | 24 | if me ~= length(lb) 25 | fprintf('(hanging_chain_dimopt) >>> corrupted data file ''%s''\n',pdat); 26 | fprintf('(hanging_chain_dimopt) >>> me (=%0i) is not equal to length(lb) (=%0i)\n',me,length(lb)); 27 | exit; 28 | end 29 | 30 | return 31 | 32 | end 33 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/example/hanging_chain_plot_params.m: -------------------------------------------------------------------------------- 1 | function [coord,ground,init,inter,final] = hanging_chain_plot_params (); 2 | 3 | % 4 | % [coord,ground,init,inter,final] = hanging_chain_plot_params (); 5 | % 6 | % This function can be used to set a few plot parameters. To modify the 7 | % parameters, copy this file in the working directory and change the 8 | % settings (files in the Libopt environment should not be modified). 9 | 10 | % Specifications for axis 11 | 12 | coord.present = 1; % with (1) or without (0) axis 13 | 14 | % Specifications for the floor/ground 15 | 16 | ground.present = 1; % with (1) or without (0) the floor/ground (in case there is indeed a floor) 17 | ground.color = 0.81*[1 1 1]; % higher for lighter (0.81 is the Matlab gray) 18 | 19 | % Specifications for the initial hanging chain 20 | 21 | init.present = 1; % with (1) or without (0) initial hanging chain 22 | init.type = '-'; 23 | init.color = 'r'; 24 | init.width = 1; 25 | 26 | % Specifications for the intermediate hanging chains 27 | 28 | inter.present = 1; % with (1) or without (0) intermediate hanging chains 29 | inter.type = '--'; 30 | inter.color = 'm'; 31 | inter.width = 1; 32 | 33 | % Specifications for the final hanging chain 34 | 35 | final.present = 1; % with (1) or without (0) final hanging chain 36 | final.type = '-'; 37 | final.color = 'b'; 38 | final.width = 2; 39 | 40 | end 41 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_badsimul.m: -------------------------------------------------------------------------------- 1 | function [info] = sqplab_badsimul (outdic,info,options,values) 2 | 3 | % 4 | % [info] = sqplab_badsimul (outdic,info,options,values) 5 | % 6 | % Print a message, modify info, and return 7 | 8 | %----------------------------------------------------------------------- 9 | % 10 | % Author: Jean Charles Gilbert, INRIA. 11 | % 12 | % Copyright 2008, 2009, INRIA. 13 | % 14 | % SQPlab is distributed under the terms of the Q Public License version 15 | % 1.0. 16 | % 17 | % This program is distributed in the hope that it will be useful, but 18 | % WITHOUT ANY WARRANTY; without even the implied warranty of 19 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 20 | % License version 1.0 for more details. 21 | % 22 | % You should have received a copy of the Q Public License version 1.0 23 | % along with this program. If not, see 24 | % . 25 | % 26 | %----------------------------------------------------------------------- 27 | 28 | % Let's go 29 | 30 | if outdic == 2 31 | if options.verbose; fprintf(options.fout,'\n\n### sqplab: the simulator wants to stop\n'); end; 32 | info.flag = values.stop_on_simul; 33 | elseif outdic > 2 34 | if options.verbose; fprintf(options.fout,'\n\n### sqplab: error with the simulator (outdic = %0i)\n',outdic); end; 35 | info.flag = values.fail_on_simul; 36 | end 37 | 38 | return 39 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_dnorm.m: -------------------------------------------------------------------------------- 1 | function [vn] = sqplab_dnorm (v); 2 | 3 | % 4 | % [vn] = sqplab_dnorm (v); 5 | % 6 | % Computes the dual norm of the vector v. This norm must be the dual 7 | % norm of the one computed by sqplab_pnorm (). 8 | % 9 | % In the present case, this is the infinity norm. 10 | 11 | %----------------------------------------------------------------------- 12 | % 13 | % Author: Jean Charles Gilbert, INRIA. 14 | % 15 | % Copyright 2008, 2009, INRIA. 16 | % 17 | % SQPlab is distributed under the terms of the Q Public License version 18 | % 1.0. 19 | % 20 | % This program is distributed in the hope that it will be useful, but 21 | % WITHOUT ANY WARRANTY; without even the implied warranty of 22 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 23 | % License version 1.0 for more details. 24 | % 25 | % You should have received a copy of the Q Public License version 1.0 26 | % along with this program. If not, see 27 | % . 28 | % 29 | %----------------------------------------------------------------------- 30 | 31 | vn = norm(v,inf); 32 | 33 | return 34 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_givens.m: -------------------------------------------------------------------------------- 1 | function [c,s,r] = sqplab_givens (a,b) 2 | 3 | % 4 | % [c,s,r] = sqplab_givens (a,b) 5 | % 6 | % Being given a vector [a;b], sqplab_givens returns the elements 'c' 7 | % and 's' of the orthogonal Givens matrix G = [c, -s; s, c], as well as 8 | % the scalar 'r' such that G*[a;b] = [r;0]. 9 | 10 | %----------------------------------------------------------------------- 11 | % 12 | % Author: Jean Charles Gilbert, INRIA. 13 | % 14 | % Copyright 2008, 2009, INRIA. 15 | % 16 | % SQPlab is distributed under the terms of the Q Public License version 17 | % 1.0. 18 | % 19 | % This program is distributed in the hope that it will be useful, but 20 | % WITHOUT ANY WARRANTY; without even the implied warranty of 21 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 22 | % License version 1.0 for more details. 23 | % 24 | % You should have received a copy of the Q Public License version 1.0 25 | % along with this program. If not, see 26 | % . 27 | % 28 | %----------------------------------------------------------------------- 29 | 30 | if b == 0 31 | c = 1; s = 0; r = a; 32 | else 33 | if abs(b) > abs(a) 34 | t = -a/b; s = -1/sqrt(1+t^2); c = t*s; 35 | else 36 | t = -b/a; c = 1/sqrt(1+t^2); s = t*c; 37 | end 38 | end 39 | r = c*a-s*b; 40 | 41 | return 42 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_lgivens.m: -------------------------------------------------------------------------------- 1 | function A = sqplab_lgivens (A,n,c,s) 2 | 3 | % 4 | % Left multiplication by a Givens rotation matrix. 5 | % 6 | % A = sqplab_lgivens (A,n,c,s) 7 | % 8 | % The matrix A(2,n) is replaced by G*A, where G is the Givens rotation 9 | % matrix G = [c, -s; s, c], where c^2 + s^2 = 1 (note the change of 10 | % sign with respect to rgivens). 11 | 12 | %----------------------------------------------------------------------- 13 | % 14 | % Author: Jean Charles Gilbert, INRIA. 15 | % 16 | % Copyright 2008, 2009, INRIA. 17 | % 18 | % SQPlab is distributed under the terms of the Q Public License version 19 | % 1.0. 20 | % 21 | % This program is distributed in the hope that it will be useful, but 22 | % WITHOUT ANY WARRANTY; without even the implied warranty of 23 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 24 | % License version 1.0 for more details. 25 | % 26 | % You should have received a copy of the Q Public License version 1.0 27 | % along with this program. If not, see 28 | % . 29 | % 30 | %----------------------------------------------------------------------- 31 | 32 | for j = 1:n 33 | t1 = A(1,j); 34 | t2 = A(2,j); 35 | A(1,j) = c*t1-s*t2; 36 | A(2,j) = s*t1+c*t2; 37 | end 38 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_pnorm.m: -------------------------------------------------------------------------------- 1 | function [vn] = sqplab_pnorm (v); 2 | 3 | % 4 | % [vn] = sqplab_pnorm (v); 5 | % 6 | % Computes the (primal) norm of the vector v. This norm is the one in 7 | % the merit function used by SQPlab. 8 | % 9 | % In the present case, this is the 1-norm. 10 | 11 | %----------------------------------------------------------------------- 12 | % 13 | % Author: Jean Charles Gilbert, INRIA. 14 | % 15 | % Copyright 2008, 2009, INRIA. 16 | % 17 | % SQPlab is distributed under the terms of the Q Public License version 18 | % 1.0. 19 | % 20 | % This program is distributed in the hope that it will be useful, but 21 | % WITHOUT ANY WARRANTY; without even the implied warranty of 22 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 23 | % License version 1.0 for more details. 24 | % 25 | % You should have received a copy of the Q Public License version 1.0 26 | % along with this program. If not, see 27 | % . 28 | % 29 | %----------------------------------------------------------------------- 30 | 31 | vn = norm(v,1); 32 | 33 | return 34 | -------------------------------------------------------------------------------- /Lib/sqplab-0.4.5-distrib/src/sqplab_rgivens.m: -------------------------------------------------------------------------------- 1 | function A = sqplab_rgivens (A,m,c,s) 2 | 3 | % 4 | % Right multiplication by a Givens rotation matrix. 5 | % 6 | % A = sqplab_rgivens (A,m,c,s) 7 | % 8 | % The matrix A(m,2) is replaced by A*G, where G is the Givens rotation 9 | % matrix G = [c, s; -s, c], where c^2 + s^2 = 1 (note the change of 10 | % sign with respect to lgivens). 11 | 12 | %----------------------------------------------------------------------- 13 | % 14 | % Author: Jean Charles Gilbert, INRIA. 15 | % 16 | % Copyright 2008, 2009, INRIA. 17 | % 18 | % SQPlab is distributed under the terms of the Q Public License version 19 | % 1.0. 20 | % 21 | % This program is distributed in the hope that it will be useful, but 22 | % WITHOUT ANY WARRANTY; without even the implied warranty of 23 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Q Public 24 | % License version 1.0 for more details. 25 | % 26 | % You should have received a copy of the Q Public License version 1.0 27 | % along with this program. If not, see 28 | % . 29 | % 30 | %----------------------------------------------------------------------- 31 | 32 | for i = 1:m, 33 | t1 = A(i,1); 34 | t2 = A(i,2); 35 | A(i,1) = c*t1-s*t2; 36 | A(i,2) = s*t1+c*t2; 37 | end 38 | -------------------------------------------------------------------------------- /Lib/stk/AUTHORS: -------------------------------------------------------------------------------- 1 | STK is mainly developped and maintained by (alphabetical order): 2 | 3 | BECT Julien 4 | VAZQUEZ Emmanuel 5 | 6 | The following people have helped develop the toolbox in various ways 7 | (from providing ideas and bug reports to actually contributing code): 8 | 9 | ALEKSOVSKA Ivana 10 | ASSOULINE Tom 11 | AUTRET Florent 12 | BENASSI Romain 13 | DABOUSSI Elias 14 | DRAUG Carnë 15 | DUHAMEL Stephano 16 | DUTRIEUX Héloïse 17 | FELIOT Paul 18 | FRASNEDO Sophie 19 | JAN Benoit 20 | KETTANI Othmane 21 | KRAUTH Alexandra 22 | LI Ling 23 | PIERA-MARTINEZ Miguel 24 | RAHALI Elham 25 | RAVISANKAR Ashwin 26 | RESSEGUIER Valentin 27 | STROH Rémi 28 | VILLEMONTEIX Julien 29 | 30 | Individual copyright notices are provided at the beginning of each 31 | file in the toolbox. 32 | -------------------------------------------------------------------------------- /Lib/stk/CITATION: -------------------------------------------------------------------------------- 1 | -------------------------------------------------------------------- 2 | 3 | To cite this release of STK in publications use: 4 | 5 | Julien Bect, Emmanuel Vazquez and others (2017). 6 | STK: a Small (Matlab/Octave) Toolbox for Kriging. Release 2.4. 7 | URL http://kriging.sourceforge.net 8 | 9 | A BibTeX entry for LaTeX users is: 10 | 11 | @misc{, 12 | author = {Bect, Julien and Vazquez, Emmanuel and others}, 13 | title = {{STK}: a {S}mall ({M}atlab/{O}ctave) {T}oolbox 14 | for {K}riging. {R}elease 2.4}, 15 | year = {2014}, 16 | url = {http://kriging.sourceforge.net} 17 | } 18 | 19 | We have invested a lot of time and effort in the development of STK. 20 | Please cite it if you use it. 21 | 22 | -------------------------------------------------------------------- 23 | 24 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/abs.m: -------------------------------------------------------------------------------- 1 | % ABS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = abs (x) 30 | 31 | ydata = abs (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = abs (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, abs (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/acos.m: -------------------------------------------------------------------------------- 1 | % ACOS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = acos (x) 30 | 31 | ydata = acos (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = acos (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, acos (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/acosd.m: -------------------------------------------------------------------------------- 1 | % ACOSD [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = acosd (x) 30 | 31 | ydata = acosd (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = acosd (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, acosd (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/acosh.m: -------------------------------------------------------------------------------- 1 | % ACOSH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = acosh (x) 30 | 31 | ydata = acosh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = acosh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, acosh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/and.m: -------------------------------------------------------------------------------- 1 | % AND [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = and(x1, x2) 30 | 31 | y = bsxfun(@and, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@and, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@and, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@and, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/asin.m: -------------------------------------------------------------------------------- 1 | % ASIN [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = asin (x) 30 | 31 | ydata = asin (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = asin (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, asin (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/asind.m: -------------------------------------------------------------------------------- 1 | % ASIND [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = asind (x) 30 | 31 | ydata = asind (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = asind (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, asind (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/asinh.m: -------------------------------------------------------------------------------- 1 | % ASINH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = asinh (x) 30 | 31 | ydata = asinh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = asinh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, asinh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/atan.m: -------------------------------------------------------------------------------- 1 | % ATAN [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = atan (x) 30 | 31 | ydata = atan (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = atan (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, atan (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/atand.m: -------------------------------------------------------------------------------- 1 | % ATAND [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = atand (x) 30 | 31 | ydata = atand (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = atand (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, atand (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/atanh.m: -------------------------------------------------------------------------------- 1 | % ATANH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = atanh (x) 30 | 31 | ydata = atanh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = atanh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, atanh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/cos.m: -------------------------------------------------------------------------------- 1 | % COS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = cos (x) 30 | 31 | ydata = cos (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = cos (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, cos (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/cosd.m: -------------------------------------------------------------------------------- 1 | % COSD [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = cosd (x) 30 | 31 | ydata = cosd (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = cosd (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, cosd (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/cosh.m: -------------------------------------------------------------------------------- 1 | % ABS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = cosh (x) 30 | 31 | ydata = cosh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = cosh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, cosh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/diff.m: -------------------------------------------------------------------------------- 1 | % DIFF [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function z = diff (x, varargin) 30 | 31 | z = diff (x.data, varargin{:}); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/eq.m: -------------------------------------------------------------------------------- 1 | % EQ [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = eq(x1, x2) 30 | 31 | y = bsxfun(@eq, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@eq, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@eq, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@eq, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/exp.m: -------------------------------------------------------------------------------- 1 | % EXP [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = exp (x) 30 | 31 | ydata = exp (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = exp (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, exp (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/expm1.m: -------------------------------------------------------------------------------- 1 | % EXPM1 [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = expm1 (x) 30 | 31 | ydata = expm1 (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = expm1 (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, expm1 (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/fieldnames.m: -------------------------------------------------------------------------------- 1 | % FIELDNAMES [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function fn = fieldnames (x) 30 | 31 | fn = [x.colnames reserved_field_names()]; 32 | 33 | end % function 34 | 35 | %!test 36 | %! x = stk_dataframe (rand (3, 2), {'u' 'v'}); 37 | %! s1 = sort (fieldnames (x)); 38 | %! s2 = {'colnames' 'data' 'info' 'rownames' 'u' 'v'}; 39 | %! assert (all (strcmp (s1, s2))); 40 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/ge.m: -------------------------------------------------------------------------------- 1 | % GE [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = ge(x1, x2) 30 | 31 | y = bsxfun(@ge, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@ge, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@ge, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@ge, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/gt.m: -------------------------------------------------------------------------------- 1 | % GT [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = gt(x1, x2) 30 | 31 | y = bsxfun(@gt, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@gt, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@gt, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@gt, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/isempty.m: -------------------------------------------------------------------------------- 1 | % ISEMPTY [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function b = isempty(x) 30 | 31 | b = isempty(x.data); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/ldivide.m: -------------------------------------------------------------------------------- 1 | % LDIVIDE [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = ldivide(x1, x2) 30 | 31 | y = bsxfun(@ldivide, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@ldivide, 1 + rand(7, 2), rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@ldivide, 1 + rand(7, 2), pi); 37 | %!error stk_test_dfbinaryop(@ldivide, 1 + rand(7, 2), rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/le.m: -------------------------------------------------------------------------------- 1 | % LE [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = le(x1, x2) 30 | 31 | y = bsxfun(@le, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@le, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@le, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@le, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/length.m: -------------------------------------------------------------------------------- 1 | % LENGTH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | 30 | function N = length (x) %#ok 31 | 32 | errmsg = ['length() is not defined for objects of class ' class(x)]; 33 | stk_error (errmsg, 'MethodNotDefined'); 34 | 35 | end % function 36 | 37 | %!error length (stk_dataframe ([1 2; 3 4; 5 6])) 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/log.m: -------------------------------------------------------------------------------- 1 | % LOG [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = log (x) 30 | 31 | ydata = log (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = log (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, log (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/log10.m: -------------------------------------------------------------------------------- 1 | % LOG10 [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = log10 (x) 30 | 31 | ydata = log10 (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = log10 (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, log10 (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/log1p.m: -------------------------------------------------------------------------------- 1 | % LOG1P [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = log1p (x) 30 | 31 | ydata = log1p (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = log1p (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, log1p (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/log2.m: -------------------------------------------------------------------------------- 1 | % LOG2 [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = log2 (x) 30 | 31 | ydata = log2 (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = log2 (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, log2 (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/lt.m: -------------------------------------------------------------------------------- 1 | % LT [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = lt(x1, x2) 30 | 31 | y = bsxfun(@lt, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@lt, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@lt, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@lt, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/minus.m: -------------------------------------------------------------------------------- 1 | % MINUS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = minus(x1, x2) 30 | 31 | y = bsxfun(@minus, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@minus, rand(7, 2), rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@minus, rand(7, 2), pi); 37 | %!error stk_test_dfbinaryop(@minus, rand(7, 2), rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/ne.m: -------------------------------------------------------------------------------- 1 | % NE [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = ne(x1, x2) 30 | 31 | y = bsxfun(@ne, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@ne, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@ne, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@ne, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/or.m: -------------------------------------------------------------------------------- 1 | % OR [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = or(x1, x2) 30 | 31 | y = bsxfun(@or, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@or, floor(3*rand(7, 2)), floor(3*rand(7, 2))); 36 | %!test stk_test_dfbinaryop(@or, floor(3*rand(7, 2)), 1.0); 37 | %!error stk_test_dfbinaryop(@or, floor(3*rand(7, 2)), floor(3*rand(7, 3))); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/plus.m: -------------------------------------------------------------------------------- 1 | % PLUS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = plus(x1, x2) 30 | 31 | y = bsxfun(@plus, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@plus, rand(7, 2), rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@plus, rand(7, 2), pi); 37 | %!error stk_test_dfbinaryop(@plus, rand(7, 2), rand(7, 3)); 38 | 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/power.m: -------------------------------------------------------------------------------- 1 | % POWER [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = power(x1, x2) 30 | 31 | y = bsxfun(@power, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@power, rand(7, 2), .1 + rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@power, rand(7, 2), .1); 37 | %!error stk_test_dfbinaryop(@power, rand(7, 2), .1 + rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/get_column_number.m: -------------------------------------------------------------------------------- 1 | % GET_COLUMN_NUMBER [STK internal] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | % Do NOT remove this seemingly useless M-file, otherwise some not-so-old 30 | % versions of Octave would faild to detect the corresponding mex-file. 31 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexa64 -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexw64 -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/get_column_number.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/private/reserved_field_names.m: -------------------------------------------------------------------------------- 1 | % RESERVED_FIELD_NAMES [STK internal] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013, 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function c = reserved_field_names () 30 | 31 | c = {'data', 'info', 'rownames', 'colnames'}; 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/rdivide.m: -------------------------------------------------------------------------------- 1 | % RDIVIDE [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = rdivide(x1, x2) 30 | 31 | y = bsxfun(@rdivide, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@rdivide, rand(7, 2), 1 + rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@rdivide, rand(7, 2), pi); 37 | %!error stk_test_dfbinaryop(@rdivide, rand(7, 2), 1 + rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/realpow.m: -------------------------------------------------------------------------------- 1 | % REALPOW [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = realpow(x1, x2) 30 | 31 | y = bsxfun(@realpow, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@realpow, rand(7, 2), .1 + rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@realpow, rand(7, 2), .1); 37 | %!error stk_test_dfbinaryop(@realpow, rand(7, 2), .1 + rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/sin.m: -------------------------------------------------------------------------------- 1 | % SIN [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = sin (x) 30 | 31 | ydata = sin (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = sin (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, sin (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/sind.m: -------------------------------------------------------------------------------- 1 | % SIND [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = sind (x) 30 | 31 | ydata = sind (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = sind (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, sind (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/sinh.m: -------------------------------------------------------------------------------- 1 | % SINH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = sinh (x) 30 | 31 | ydata = sinh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = sinh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, sinh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/sqrt.m: -------------------------------------------------------------------------------- 1 | % SQRT [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = sqrt (x) 30 | 31 | ydata = sqrt (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = sqrt (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, sqrt (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/stk_length.m: -------------------------------------------------------------------------------- 1 | % STK_LENGTH [overload STK function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Authors: Julien Bect 8 | % Emmanuel Vazquez 9 | 10 | % Copying Permission Statement 11 | % 12 | % This file is part of 13 | % 14 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 15 | % (http://sourceforge.net/projects/kriging) 16 | % 17 | % STK is free software: you can redistribute it and/or modify it under 18 | % the terms of the GNU General Public License as published by the Free 19 | % Software Foundation, either version 3 of the License, or (at your 20 | % option) any later version. 21 | % 22 | % STK is distributed in the hope that it will be useful, but WITHOUT 23 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 25 | % License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | % along with STK. If not, see . 29 | 30 | function l = stk_length (x) 31 | 32 | l = size (x.data, 1); 33 | 34 | end % function 35 | 36 | %!test 37 | %! x = stk_dataframe ([1 2; 3 4; 5 6]); 38 | %! assert (isequal (stk_length (x), 3)); 39 | 40 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/stk_sprintf_sizetype.m: -------------------------------------------------------------------------------- 1 | % STK_SPRINTF_SIZETYPE prints the size and type into a string 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013, 2014 SUPELEC 6 | % 7 | % Authors: Julien Bect 8 | % Emmanuel Vazquez 9 | 10 | % Copying Permission Statement 11 | % 12 | % This file is part of 13 | % 14 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 15 | % (http://sourceforge.net/projects/kriging) 16 | % 17 | % STK is free software: you can redistribute it and/or modify it under 18 | % the terms of the GNU General Public License as published by the Free 19 | % Software Foundation, either version 3 of the License, or (at your 20 | % option) any later version. 21 | % 22 | % STK is distributed in the hope that it will be useful, but WITHOUT 23 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 25 | % License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | % along with STK. If not, see . 29 | 30 | function s = stk_sprintf_sizetype (x) 31 | 32 | [n, d] = size (x); 33 | 34 | s = sprintf ('%dx%d %s array', n, d, class (x)); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/tan.m: -------------------------------------------------------------------------------- 1 | % TAN [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = tan (x) 30 | 31 | ydata = tan (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = tan (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, tan (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/tand.m: -------------------------------------------------------------------------------- 1 | % TAND [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = tand (x) 30 | 31 | ydata = tand (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = tand (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, tand (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/tanh.m: -------------------------------------------------------------------------------- 1 | % TANH [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function ydata = tanh (x) 30 | 31 | ydata = tanh (x.data); 32 | 33 | end % function 34 | 35 | 36 | %!test 37 | %! u = rand (4, 3); x = stk_dataframe (u); v = tanh (x); 38 | %! assert (strcmp (class (v), class (u)) && isequal (v, tanh (u))) 39 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/times.m: -------------------------------------------------------------------------------- 1 | % TIMES [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function y = times(x1, x2) 30 | 31 | y = bsxfun(@times, x1, x2); 32 | 33 | end % function 34 | 35 | %!test stk_test_dfbinaryop(@times, rand(7, 2), rand(7, 2)); 36 | %!test stk_test_dfbinaryop(@times, rand(7, 2), pi); 37 | %!error stk_test_dfbinaryop(@times, rand(7, 2), rand(7, 3)); 38 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/uminus.m: -------------------------------------------------------------------------------- 1 | % UMINUS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function x = uminus (x) 30 | 31 | x.data = - x.data; 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_dataframe/uplus.m: -------------------------------------------------------------------------------- 1 | % UPLUS [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function x = uplus (x) 30 | 31 | end % function 32 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/contour.m: -------------------------------------------------------------------------------- 1 | % CONTOUR [overload base function] 2 | % 3 | % See also: stk_factorialdesign/surf, stk_factorialdesign/mesh, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = contour (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@contour, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/mesh.m: -------------------------------------------------------------------------------- 1 | % MESH [overload base function] 2 | % 3 | % See also: stk_factorialdesign/contour, stk_factorialdesign/surf, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = mesh (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@mesh, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/meshc.m: -------------------------------------------------------------------------------- 1 | % MESHC [overload base function] 2 | % 3 | % See also: stk_factorialdesign/contour, stk_factorialdesign/surf, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = meshc (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@meshc, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/meshz.m: -------------------------------------------------------------------------------- 1 | % MESHZ [overload base function] 2 | % 3 | % See also: stk_factorialdesign/contour, stk_factorialdesign/surf, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = meshz (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@meshz, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/pcolor.m: -------------------------------------------------------------------------------- 1 | % PCOLOR [overload base function] 2 | % 3 | % See also: stk_factorialdesign/contour, stk_factorialdesign/mesh, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = pcolor (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@pcolor, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_factorialdesign/surf.m: -------------------------------------------------------------------------------- 1 | % SURF [overload base function] 2 | % 3 | % See also: stk_factorialdesign/contour, stk_factorialdesign/mesh, ... 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % Copyright (C) 2014 SUPELEC 9 | % 10 | % Author: Julien Bect 11 | 12 | % Copying Permission Statement 13 | % 14 | % This file is part of 15 | % 16 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 17 | % (http://sourceforge.net/projects/kriging) 18 | % 19 | % STK is free software: you can redistribute it and/or modify it under 20 | % the terms of the GNU General Public License as published by the Free 21 | % Software Foundation, either version 3 of the License, or (at your 22 | % option) any later version. 23 | % 24 | % STK is distributed in the hope that it will be useful, but WITHOUT 25 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 26 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 27 | % License for more details. 28 | % 29 | % You should have received a copy of the GNU General Public License 30 | % along with STK. If not, see . 31 | 32 | function h_plot = surf (varargin) 33 | 34 | h_plot = plot_surfmeshcontour (@surf, varargin{:}); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_hrect/disp.m: -------------------------------------------------------------------------------- 1 | % DISP [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function disp (hr) 30 | 31 | disp (hr.stk_dataframe); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@stk_hrect/display.m: -------------------------------------------------------------------------------- 1 | % DISPLAY [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function display (x) 30 | 31 | fprintf ( ... 32 | '\n%s = %d-dimensional hyper-rectangle (stk_hrect object):\n\n', ... 33 | inputname (1), size (x.stk_dataframe, 2)); 34 | 35 | disp (x); 36 | 37 | fprintf ('\n'); 38 | 39 | end % function 40 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@struct/double.m: -------------------------------------------------------------------------------- 1 | % DOUBLE [overload base function, deprecated] 2 | % 3 | % DEPRECATION WARNING: The use of .a structures is deprecated. 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2013, 2016 SUPELEC 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function xdata = double (x) 32 | 33 | xdata = x.a; 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/arrays/@struct/stk_length.m: -------------------------------------------------------------------------------- 1 | % STK_LENGTH [overload STK function, deprecated] 2 | % 3 | % DEPRECATION WARNING: The use of .a structures is deprecated. 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2015, 2016 CentraleSupelec 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function l = stk_length (x) 32 | 33 | l = size (x.a, 1); 34 | 35 | end % function 36 | 37 | 38 | %!test 39 | %! x = struct ('a', [5; 2]); 40 | %! assert (isequal (stk_length (x), 2)); 41 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_kreq_qr/disp.m: -------------------------------------------------------------------------------- 1 | % DISP [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function disp (x) 30 | 31 | disp (struct (x)); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_kreq_qr/display.m: -------------------------------------------------------------------------------- 1 | % DISPLAY [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function display (x) 30 | 31 | fprintf ('\n%s = <%s>\n\n', inputname (1), stk_sprintf_sizetype (x)); 32 | 33 | disp (x); 34 | 35 | fprintf ('\n'); 36 | 37 | end % function 38 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_kreq_qr/private/compute_P_scaling.m: -------------------------------------------------------------------------------- 1 | % COMPUTE_P_SCALING [STK internal] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Authors: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function s = compute_P_scaling (Kii, Pi) 30 | 31 | t = sum (Pi .^ 2); 32 | 33 | s = ones (1, size (Pi, 2)); 34 | s(t > 0) = sqrt ((max (sum (Kii .^ 2))) ./ t(t > 0)); 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/display.m: -------------------------------------------------------------------------------- 1 | % DISPLAY [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015, 2017 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function display (x) 30 | 31 | name = inputname (1); 32 | if isempty (name) 33 | name = 'ans'; 34 | end 35 | 36 | if stk_disp_isloose 37 | fprintf ('\n'); 38 | end 39 | 40 | fprintf ('%s = ', name); 41 | 42 | disp (x); 43 | 44 | end % function 45 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/fieldnames.m: -------------------------------------------------------------------------------- 1 | % FIELDNAMES [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2016 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function fn = fieldnames (model) 30 | 31 | fn = {'prior_model'; 'input_data'; 'output_data'}; 32 | 33 | % Note: kreq is a 'hidden' field that might change in future versions of STK, 34 | % we don't want ordinary users to see it 35 | 36 | end % function 37 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/get_input_data.m: -------------------------------------------------------------------------------- 1 | % GET_INTPUT_DATA returns the input data of the model 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015, 2016 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function input_data = get_input_data (model) 30 | 31 | input_data = model.input_data; 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/get_output_data.m: -------------------------------------------------------------------------------- 1 | % GET_OUTPUT_DATA returns the output data of the model 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015, 2016 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function output_data = get_output_data (model) 30 | 31 | output_data = model.output_data; 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/get_prior_model.m: -------------------------------------------------------------------------------- 1 | % GET_PRIOR_MODEL returns the prior_model structure 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2016 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function prior_model = get_prior_model (model) 30 | 31 | prior_model = model.prior_model; 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/horzcat.m: -------------------------------------------------------------------------------- 1 | % HORZCAT [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function varargout = horzcat (varargin) 30 | 31 | stk_error (['Arrays of stk_model_gpposterior objects are not supported. ', ... 32 | 'Use cell arrays instead.'], 'IllegalOperation'); 33 | 34 | end % function 35 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/stk_isnoisy.m: -------------------------------------------------------------------------------- 1 | % STK_ISNOISY [overload STK function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2016 CentraleSupelec & LNE 6 | % 7 | % Authors: Julien Bect 8 | % Rémi Stroh 9 | 10 | % Copying Permission Statement 11 | % 12 | % This file is part of 13 | % 14 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 15 | % (http://sourceforge.net/projects/kriging) 16 | % 17 | % STK is free software: you can redistribute it and/or modify it under 18 | % the terms of the GNU General Public License as published by the Free 19 | % Software Foundation, either version 3 of the License, or (at your 20 | % option) any later version. 21 | % 22 | % STK is distributed in the hope that it will be useful, but WITHOUT 23 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 25 | % License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | % along with STK. If not, see . 29 | 30 | function b = stk_isnoisy (model) 31 | 32 | b = stk_isnoisy (model.prior_model); 33 | 34 | end % function 35 | -------------------------------------------------------------------------------- /Lib/stk/core/@stk_model_gpposterior/vertcat.m: -------------------------------------------------------------------------------- 1 | % VERTCAT [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015 CentraleSupelec 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function varargout = vertcat (varargin) 30 | 31 | stk_error (['Arrays of stk_model_gpposterior objects are not supported. ', ... 32 | 'Use cell arrays instead.'], 'IllegalOperation'); 33 | 34 | end % function 35 | -------------------------------------------------------------------------------- /Lib/stk/covfcs/rbf/stk_sf_gausscorr.m: -------------------------------------------------------------------------------- 1 | % STK_SF_GAUSSCORR is deprecated, use stk_rbf_gauss instead 2 | % 3 | % See also: stk_rbf_gauss 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function k = stk_sf_gausscorr (h, diff) 32 | 33 | k = stk_rbf_gauss (h, diff); 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/covfcs/rbf/stk_sf_matern.m: -------------------------------------------------------------------------------- 1 | % STK_SF_MATERN is deprecated, use stk_rbf_matern instead 2 | % 3 | % See also: stk_rbf_matern 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function k = stk_sf_matern (nu, h, diff) 32 | 33 | k = stk_rbf_matern (h, diff); 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/covfcs/rbf/stk_sf_matern32.m: -------------------------------------------------------------------------------- 1 | % STK_SF_MATERN32 is deprecated, use stk_rbf_matern32 instead 2 | % 3 | % See also: stk_rbf_matern32 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function k = stk_sf_matern32 (h, diff) 32 | 33 | k = stk_rbf_matern32 (h, diff); 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/covfcs/rbf/stk_sf_matern52.m: -------------------------------------------------------------------------------- 1 | % STK_SF_MATERN52 is deprecated, use stk_rbf_matern52 instead 2 | % 3 | % See also: stk_rbf_matern52 4 | 5 | % Copyright Notice 6 | % 7 | % Copyright (C) 2016 CentraleSupelec 8 | % 9 | % Author: Julien Bect 10 | 11 | % Copying Permission Statement 12 | % 13 | % This file is part of 14 | % 15 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 16 | % (http://sourceforge.net/projects/kriging) 17 | % 18 | % STK is free software: you can redistribute it and/or modify it under 19 | % the terms of the GNU General Public License as published by the Free 20 | % Software Foundation, either version 3 of the License, or (at your 21 | % option) any later version. 22 | % 23 | % STK is distributed in the hope that it will be useful, but WITHOUT 24 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 25 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 26 | % License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | % along with STK. If not, see . 30 | 31 | function k = stk_sf_matern52 (h, diff) 32 | 33 | k = stk_rbf_matern52 (h, diff); 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/doc/html/images/download.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/doc/html/images/download.png -------------------------------------------------------------------------------- /Lib/stk/doc/html/images/news.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/doc/html/images/news.png -------------------------------------------------------------------------------- /Lib/stk/doc/html/images/stk_logo.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/doc/html/images/stk_logo.png -------------------------------------------------------------------------------- /Lib/stk/examples/01_kriging_basics/stk_example_kb04.m: -------------------------------------------------------------------------------- 1 | % STK_EXAMPLE_KB04 Estimating the variance of the noise 2 | % 3 | % This example no longer exists. See stk_example_kb02n instead. 4 | % 5 | % See also: stk_example_kb01n, stk_example_kb02n 6 | 7 | % Copyright Notice 8 | % 9 | % Copyright (C) 2016 CentraleSupelec 10 | % 11 | % Author: Julien Bect 12 | 13 | % Copying Permission Statement 14 | % 15 | % This file is part of 16 | % 17 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 18 | % (http://sourceforge.net/projects/kriging) 19 | % 20 | % STK is free software: you can redistribute it and/or modify it under 21 | % the terms of the GNU General Public License as published by the Free 22 | % Software Foundation, either version 3 of the License, or (at your 23 | % option) any later version. 24 | % 25 | % STK is distributed in the hope that it will be useful, but WITHOUT 26 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 27 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 28 | % License for more details. 29 | % 30 | % You should have received a copy of the GNU General Public License 31 | % along with STK. If not, see . 32 | 33 | stk_disp_examplewelcome (); 34 | -------------------------------------------------------------------------------- /Lib/stk/lm/@stk_lm_affine/feval.m: -------------------------------------------------------------------------------- 1 | % FEVAL [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2012-2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function z = feval (lm, x) %#ok 30 | 31 | n = size (x, 1); 32 | z = horzcat (ones (n, 1), x); 33 | 34 | end % function 35 | 36 | 37 | %!test 38 | %! n = 15; d = 4; 39 | %! x = stk_sampling_randunif (n, d); 40 | %! P = feval (stk_lm_affine (), x); 41 | %! assert (isequal (size (P), [n, d + 1])); 42 | -------------------------------------------------------------------------------- /Lib/stk/lm/@stk_lm_constant/feval.m: -------------------------------------------------------------------------------- 1 | % FEVAL [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2012-2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function z = feval (lm, x) %#ok 30 | 31 | n = size (x, 1); 32 | z = ones (n, 1); 33 | 34 | end % function 35 | 36 | 37 | %!test 38 | %! n = 15; d = 4; 39 | %! x = stk_sampling_randunif (n, d); 40 | %! P = feval (stk_lm_constant (), x); 41 | %! assert (isequal (size (P), [n, 1])); 42 | -------------------------------------------------------------------------------- /Lib/stk/lm/@stk_lm_matrix/feval.m: -------------------------------------------------------------------------------- 1 | % FEVAL [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2012-2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function z = feval (lm, x) 30 | 31 | x = double (x); 32 | 33 | if isempty (lm.data), 34 | n = size (x, 1); 35 | z = zeros (n, 0); 36 | else 37 | z = lm.data(x, :); 38 | end 39 | 40 | end % function 41 | -------------------------------------------------------------------------------- /Lib/stk/lm/@stk_lm_null/feval.m: -------------------------------------------------------------------------------- 1 | % FEVAL [overload base function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2012-2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function z = feval (lm, x) %#ok 30 | 31 | n = size (x, 1); 32 | z = zeros (n, 0); 33 | 34 | end % function 35 | 36 | %!test 37 | %! n = 15; d = 4; 38 | %! x = stk_sampling_randunif (n, d); 39 | %! P = feval (stk_lm_null (), x); 40 | %! assert (isequal (size (P), [n, 0])); 41 | -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixx.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_matrixx.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixx.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixx.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_matrixx.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixx.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixy.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_matrixy.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixy.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixy.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_matrixy.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_matrixy.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_pairwise.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_pairwise.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_pairwise.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_pairwise.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_dist_pairwise.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_dist_pairwise.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_filldist_discr_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixx.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_matrixy.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_gpquadform_pairwise.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_mindist_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_mindist_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_mindist_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_mindist_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/dist/private/stk_mindist_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/dist/private/stk_mindist_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/distrib/private/stk_distrib_bivnorm0_cdf.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/mole/README: -------------------------------------------------------------------------------- 1 | What is the MOLE ? 2 | ================== 3 | 4 | MOLE stands for "Matlab/Octave Langage Extensions" 5 | 6 | This directory is meant to contain general purpose functions that either extend 7 | the Matlab/Octave langage or provided a replacement for a function that is 8 | missing or buggy in one of the two langages. 9 | 10 | 11 | Long-term goal 12 | ============== 13 | 14 | The MOLE might one day be released as a package on its own, independently of 15 | STK. Therefore, here, functions DO NOT have their name prefixed by stk_. 16 | 17 | 18 | Directories 19 | =========== 20 | 21 | ./matlab contains functions that are missing in Matlab: 22 | 23 | o fflush 24 | o file_in_loadpath 25 | o file_in_path 26 | o index 27 | o page_screen_output 28 | o stdout 29 | 30 | Other functions are provided in single-function subdirectories: 31 | 32 | o isoctave 33 | o corr 34 | o linsolve 35 | o quantile 36 | -------------------------------------------------------------------------------- /Lib/stk/misc/mole/isrow/isrow.m: -------------------------------------------------------------------------------- 1 | % ISROW returns true for row vectors and false otherwise 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2014 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function b = isrow (x) 30 | 31 | n = length (x); 32 | b = isequal (size (x), [1 n]); 33 | 34 | end % function 35 | -------------------------------------------------------------------------------- /Lib/stk/misc/optim/@stk_optim_octavesqp/private/sqp_quadprog_testf.m: -------------------------------------------------------------------------------- 1 | function obj = sqp_quadprog_testf (x) 2 | 3 | obj = exp (prod (x)) - 0.5*(x(1)^3 + x(2)^3 + 1)^2; 4 | 5 | end % function 6 | -------------------------------------------------------------------------------- /Lib/stk/misc/optim/@stk_optim_octavesqp/private/sqp_quadprog_testg.m: -------------------------------------------------------------------------------- 1 | function r = sqp_quadprog_testg (x) 2 | 3 | r = [sum(abs(x).^2)-10; ... 4 | x(2)*x(3)-5*x(4)*x(5); ... 5 | x(1)^3+x(2)^3+1 ]; 6 | 7 | end % function 8 | -------------------------------------------------------------------------------- /Lib/stk/misc/parallel/@stk_parallel_engine_none/stk_parallel_engine_none.m: -------------------------------------------------------------------------------- 1 | % STK_PARALLEL_ENGINE_NONE [overload STK function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2017 CentraleSupelec 6 | % Copyright (C) 2013 SUPELEC 7 | % 8 | % Author: Julien Bect 9 | 10 | % Copying Permission Statement 11 | % 12 | % This file is part of 13 | % 14 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 15 | % (http://sourceforge.net/projects/kriging) 16 | % 17 | % STK is free software: you can redistribute it and/or modify it under 18 | % the terms of the GNU General Public License as published by the Free 19 | % Software Foundation, either version 3 of the License, or (at your 20 | % option) any later version. 21 | % 22 | % STK is distributed in the hope that it will be useful, but WITHOUT 23 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 25 | % License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | % along with STK. If not, see . 29 | 30 | function pareng = stk_parallel_engine_none() 31 | 32 | pareng = class(struct(), 'stk_parallel_engine_none'); 33 | 34 | end % function 35 | 36 | 37 | %!test stk_test_class ('stk_parallel_engine_none') 38 | -------------------------------------------------------------------------------- /Lib/stk/misc/parallel/@stk_parallel_engine_none/stk_parallel_stop.m: -------------------------------------------------------------------------------- 1 | % STK_PARALLEL_STOP [overload STK function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function stk_parallel_stop(eng) %#ok 30 | 31 | % nothing to do ! 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/misc/parallel/@stk_parallel_engine_parfor/stk_parallel_stop.m: -------------------------------------------------------------------------------- 1 | % STK_PARALLEL_STOP [overload STK function] 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function stk_parallel_stop(eng) %#ok 30 | 31 | if matlabpool('size') == 0, 32 | warning('There was no worker pool open.'); 33 | else 34 | matlabpool close; 35 | end 36 | 37 | stk_parallel_engine_set(stk_parallel_engine_none()); 38 | 39 | end % function 40 | -------------------------------------------------------------------------------- /Lib/stk/misc/parallel/stk_parallel_engine_get.m: -------------------------------------------------------------------------------- 1 | % STK_PARALLEL_ENGINE_GET returns the current parallelization engine. 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function pareng = stk_parallel_engine_get() 30 | 31 | pareng = stk_parallel_engine_set(); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/misc/parallel/stk_parallel_stop.m: -------------------------------------------------------------------------------- 1 | % STK_PARALLEL_STOP stops the parallelization engine. 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function stk_parallel_stop () 30 | 31 | stk_parallel_stop (stk_parallel_engine_get ()); 32 | 33 | end % function 34 | -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_dominatedhv_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_isdominated_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_isdominated_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_isdominated_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_isdominated_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_isdominated_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_isdominated_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_paretofind_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_paretofind_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_paretofind_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_paretofind_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/misc/pareto/private/stk_paretofind_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/misc/pareto/private/stk_paretofind_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/misc/text/stk_sprintf_framed.m: -------------------------------------------------------------------------------- 1 | % STK_SPRINTF_FRAMED writes formatted text into a frame. 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2012, 2013 SUPELEC 6 | % 7 | % Author: Julien Bect 8 | 9 | % Copying Permission Statement 10 | % 11 | % This file is part of 12 | % 13 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 14 | % (http://sourceforge.net/projects/kriging) 15 | % 16 | % STK is free software: you can redistribute it and/or modify it under 17 | % the terms of the GNU General Public License as published by the Free 18 | % Software Foundation, either version 3 of the License, or (at your 19 | % option) any later version. 20 | % 21 | % STK is distributed in the hope that it will be useful, but WITHOUT 22 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 23 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 24 | % License for more details. 25 | % 26 | % You should have received a copy of the GNU General Public License 27 | % along with STK. If not, see . 28 | 29 | function s = stk_sprintf_framed(fmt, varargin) 30 | 31 | a = sprintf([' ' fmt ' '], varargin{:}); 32 | b = repmat('=', 1, length(a)); 33 | s = sprintf('#%s#\n#%s#\n#%s#\n', b, a, b); 34 | 35 | end % function 36 | -------------------------------------------------------------------------------- /Lib/stk/sampling/private/stk_sampling_sobol_mex.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/sampling/private/stk_sampling_sobol_mex.mexa64 -------------------------------------------------------------------------------- /Lib/stk/sampling/private/stk_sampling_sobol_mex.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/sampling/private/stk_sampling_sobol_mex.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/sampling/private/stk_sampling_sobol_mex.mexw64 -------------------------------------------------------------------------------- /Lib/stk/sampling/private/stk_sampling_sobol_mex.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/sampling/stk_sampling_vdc_rr2.mexa64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/sampling/stk_sampling_vdc_rr2.mexa64 -------------------------------------------------------------------------------- /Lib/stk/sampling/stk_sampling_vdc_rr2.mexa64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.0.0.341360 (R2016a) -------------------------------------------------------------------------------- /Lib/stk/sampling/stk_sampling_vdc_rr2.mexw64: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/freeSBDO/SBDOT/f249949604a57b8dae73906c4282a5971e7974ba/Lib/stk/sampling/stk_sampling_vdc_rr2.mexw64 -------------------------------------------------------------------------------- /Lib/stk/sampling/stk_sampling_vdc_rr2.mexw64.info: -------------------------------------------------------------------------------- 1 | MATLAB 9.2.0.556344 (R2017a) -------------------------------------------------------------------------------- /Lib/stk/stk_version.m: -------------------------------------------------------------------------------- 1 | % STK_VERSION returns STK's version number 2 | 3 | % Copyright Notice 4 | % 5 | % Copyright (C) 2015-2017 CentraleSupelec 6 | % Copyright (C) 2013, 2014 SUPELEC 7 | % 8 | % Author: Julien Bect 9 | 10 | % Copying Permission Statement 11 | % 12 | % This file is part of 13 | % 14 | % STK: a Small (Matlab/Octave) Toolbox for Kriging 15 | % (http://sourceforge.net/projects/kriging) 16 | % 17 | % STK is free software: you can redistribute it and/or modify it under 18 | % the terms of the GNU General Public License as published by the Free 19 | % Software Foundation, either version 3 of the License, or (at your 20 | % option) any later version. 21 | % 22 | % STK is distributed in the hope that it will be useful, but WITHOUT 23 | % ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 | % or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public 25 | % License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | % along with STK. If not, see . 29 | 30 | function v = stk_version () 31 | 32 | v = '2.4.2'; 33 | 34 | end % function 35 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # Surrogate Based Design Optimization Toolbox 2 | 3 | --- 4 | 5 | MATLAB Toolbox for metamodeling and solving optimization problem. 6 | The function to surrogate and optimize can be evaluated from an external numerical simulator. 7 | 8 | > SBDOT Copyright (C) 2017 CEA - LETI, DOPT. Author: C. Durantin 9 | 10 | General informations 11 | -- 12 | 13 | Metamodeling features 14 | 15 | * Kriging (gaussian process based surrogate model) 16 | * Radial Basis Function 17 | * Co-kriging (multifidelity) 18 | * CoRBF (multifidelity) 19 | * BQQV Kriging (for qualitative variables) 20 | 21 | Optimization features 22 | 23 | * CMAES (constrained single objective) 24 | * NSGA-II (constrained multiobjective) 25 | * MGDA (unconstrained multiobjective) 26 | 27 | Adaptive Sampling strategy 28 | 29 | * Expected improvement (constrained global optimization) 30 | * Gutmann criterion (constrained global optimization) 31 | * Expected hypervolume improvement (multiobjective) 32 | * Robust efficient global optimization 33 | 34 | MATLAB toolboxes requirements 35 | -- 36 | 37 | * Optimization toolbox 38 | * Statistics toolbox 39 | -------------------------------------------------------------------------------- /Src/Adaptive_Sampling/Main_objects/@Robust_EGO/Mean_meas.m: -------------------------------------------------------------------------------- 1 | function [ meas ] = Mean_meas( obj, data, nb_points ) 2 | % MEAN_MEAS 3 | % Compute the mean robustness measure 4 | 5 | meas = mean( reshape( data, obj.CRN_samples, nb_points ), 1 )'; 6 | 7 | end 8 | 9 | 10 | 11 | 12 | 13 | % ========================================================================== 14 | % 15 | % This file is part of SBDOT. 16 | % 17 | % SBDOT is free software: you can redistribute it and/or modify 18 | % it under the terms of the GNU General Public License as published by 19 | % the Free Software Foundation, either version 3 of the License, or 20 | % (at your option) any later version. 21 | % 22 | % SBDOT is distributed in the hope that it will be useful, 23 | % but WITHOUT ANY WARRANTY; without even the implied warranty of 24 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 25 | % GNU General Public License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | %    along with SBDOT.  If not, see . 29 | % 30 | % Use of SBDOT is free for personal, non-profit, pure academic research 31 | % and educational purposes. Restrictions apply on commercial or funded 32 | % research use. Please read the IMPORTANT_LICENCE_NOTICE file. 33 | % 34 | % ========================================================================== 35 | 36 | 37 | -------------------------------------------------------------------------------- /Src/Adaptive_Sampling/Main_objects/@Robust_EGO/Var_meas.m: -------------------------------------------------------------------------------- 1 | function [ meas ] = Var_meas( obj, data, nb_points ) 2 | % VAR_MEAS 3 | % Compute the variance robustness measure 4 | 5 | meas = var( reshape( data, obj.CRN_samples, nb_points ), 1 )'; 6 | 7 | end 8 | 9 | 10 | 11 | 12 | 13 | % ========================================================================== 14 | % 15 | % This file is part of SBDOT. 16 | % 17 | % SBDOT is free software: you can redistribute it and/or modify 18 | % it under the terms of the GNU General Public License as published by 19 | % the Free Software Foundation, either version 3 of the License, or 20 | % (at your option) any later version. 21 | % 22 | % SBDOT is distributed in the hope that it will be useful, 23 | % but WITHOUT ANY WARRANTY; without even the implied warranty of 24 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 25 | % GNU General Public License for more details. 26 | % 27 | % You should have received a copy of the GNU General Public License 28 | %    along with SBDOT.  If not, see . 29 | % 30 | % Use of SBDOT is free for personal, non-profit, pure academic research 31 | % and educational purposes. Restrictions apply on commercial or funded 32 | % research use. Please read the IMPORTANT_LICENCE_NOTICE file. 33 | % 34 | % ========================================================================== 35 | 36 | 37 | -------------------------------------------------------------------------------- /Src/Analytical_test_function/MO_convex_nsga.m: -------------------------------------------------------------------------------- 1 | function [ y, cons ] = MO_convex_nsga( x ) 2 | % MO_CONVEX 3 | % For example script 4 | 5 | 6 | y(:,1)=x(:,1).^2+x(:,2).^2; 7 | y(:,2)=(x(:,1)-2).^2+(x(:,2)-2).^2; 8 | 9 | cons=[]; 10 | 11 | end 12 | 13 | 14 | 15 | 16 | 17 | % ========================================================================== 18 | % 19 | % This file is part of SBDOT. 20 | % 21 | % SBDOT is free software: you can redistribute it and/or modify 22 | % it under the terms of the GNU General Public License as published by 23 | % the Free Software Foundation, either version 3 of the License, or 24 | % (at your option) any later version. 25 | % 26 | % SBDOT is distributed in the hope that it will be useful, 27 | % but WITHOUT ANY WARRANTY; without even the implied warranty of 28 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 29 | % GNU General Public License for more details. 30 | % 31 | % You should have received a copy of the GNU General Public License 32 | %    along with SBDOT.  If not, see . 33 | % 34 | % Use of SBDOT is free for personal, non-profit, pure academic research 35 | % and educational purposes. Restrictions apply on commercial or funded 36 | % research use. Please read the IMPORTANT_LICENCE_NOTICE file. 37 | % 38 | % ========================================================================== 39 | 40 | 41 | -------------------------------------------------------------------------------- /Src/Analytical_test_function/Robust_1D.m: -------------------------------------------------------------------------------- 1 | function [ y ] = Robust_1D( x ) 2 | % ROBUST_1D 3 | % For example script 4 | 5 | y = 5 - ... 6 | 0.7 * exp( -(x-0.2).^2 ./ 0.04 ) - ... 7 | 2*exp( -(x-0.5).^2./1 ) - ... 8 | exp( -(x-0.8).^2 ./ 0.01); 9 | 10 | end 11 | 12 | 13 | 14 | 15 | 16 | % ========================================================================== 17 | % 18 | % This file is part of SBDOT. 19 | % 20 | % SBDOT is free software: you can redistribute it and/or modify 21 | % it under the terms of the GNU General Public License as published by 22 | % the Free Software Foundation, either version 3 of the License, or 23 | % (at your option) any later version. 24 | % 25 | % SBDOT is distributed in the hope that it will be useful, 26 | % but WITHOUT ANY WARRANTY; without even the implied warranty of 27 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 28 | % GNU General Public License for more details. 29 | % 30 | % You should have received a copy of the GNU General Public License 31 | %    along with SBDOT.  If not, see . 32 | % 33 | % Use of SBDOT is free for personal, non-profit, pure academic research 34 | % and educational purposes. Restrictions apply on commercial or funded 35 | % research use. Please read the IMPORTANT_LICENCE_NOTICE file. 36 | % 37 | % ========================================================================== 38 | 39 | 40 | -------------------------------------------------------------------------------- /Src/Tools/Rand_gauss.m: -------------------------------------------------------------------------------- 1 | function [ x_norm ] = Rand_gauss( n, m, mu, sigma ) 2 | %RAND_GAUSS Summary of this function goes here 3 | % Detailed explanation goes here 4 | 5 | x = randn( n, m ); 6 | x_norm = bsxfun( @plus, bsxfun( @times, x, sigma ), mu ); 7 | 8 | end 9 | 10 | 11 | 12 | 13 | 14 | % ========================================================================== 15 | % 16 | % This file is part of SBDOT. 17 | % 18 | % SBDOT is free software: you can redistribute it and/or modify 19 | % it under the terms of the GNU General Public License as published by 20 | % the Free Software Foundation, either version 3 of the License, or 21 | % (at your option) any later version. 22 | % 23 | % SBDOT is distributed in the hope that it will be useful, 24 | % but WITHOUT ANY WARRANTY; without even the implied warranty of 25 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 26 | % GNU General Public License for more details. 27 | % 28 | % You should have received a copy of the GNU General Public License 29 | %    along with SBDOT.  If not, see . 30 | % 31 | % Use of SBDOT is free for personal, non-profit, pure academic research 32 | % and educational purposes. Restrictions apply on commercial or funded 33 | % research use. Please read the IMPORTANT_LICENCE_NOTICE file. 34 | % 35 | % ========================================================================== 36 | 37 | 38 | --------------------------------------------------------------------------------