# solve DGP with one of the models {dgpsystem, dgp, dgpmax}.mod
#   data file in eiter gph/ or protein/

## choose a model
param formulation symbolic;
model dgpsystem.mod;
let formulation := "dgpsystem";
#
#model dgp.mod;
#let formulation := "dgp";
#
#model dgpmax.mod;
#let formulation := "dgpmax";

## choose embedding dimension (usually 2 or 3)
let Kdim := 2;

## choose a data file
#data gph/dmdgp-3_10.dat;
#data gph/dmdgp-3_20.dat;
#data gph/dmdgp-3_30.dat;
#data gph/dmdgp-3_50.dat;
#data gph/euclid-10_0.2.dat;
#data gph/euclid-10_0.5.dat;
#data gph/euclid-10_0.8.dat;
#data gph/euclid-20_0.2.dat;
#data gph/euclid-20_0.5.dat;
#data gph/euclid-20_0.8.dat;
#data gph/euclid-30_0.2.dat;
#data gph/euclid-30_0.5.dat;
#data gph/euclid-30_0.8.dat;
#data gph/euclid-50_0.2.dat;
#data gph/euclid-50_0.5.dat;
#data gph/euclid-50_0.8.dat;
#data gph/hypercube-4.dat;
#data gph/hypercube-5.dat;
#data gph/hypercube-6.dat;
#data gph/norm-10_0.5_1.dat;
#data gph/norm-20_0.5_1.dat;
#data gph/norm-30_0.5_1.dat;
#data gph/norm-50_0.5_1.dat;
#data protein/tiny_gph.dat;
data protein/pept_gph.dat;
#data protein/1guu_gph.dat;

## choose a solver
# name of the solver
param solver symbolic;
let solver := "baron";
option solver baron;
option baron_options "outlev=1 epsr=1e-3 epsa=1e-4 prfreq=50";
solve;
#
#let solver := "knitro";
#option solver knitro;
#solve;
#
#let solver := "snopt";
#option solver snopt;
#let {i in V, k in K} x[i,k] := Uniform(-10,10); #starting point
#solve;
#
#let solver := "multistart";
#include "dgp_ms.run";

## output
include "dgp-postprocess.run";