|
30 | 30 | "\n", |
31 | 31 | "import matplotlib.pyplot as plt \n", |
32 | 32 | "import copy\n", |
33 | | - "import GPflow\n", |
34 | | - "import GPflowOpt\n", |
| 33 | + "import gpflow\n", |
| 34 | + "import gpflowopt\n", |
35 | 35 | "import tensorflow as tf" |
36 | 36 | ] |
37 | 37 | }, |
|
49 | 49 | " return f[:,None] + rng.rand(X.shape[0], 1) * 0.25\n", |
50 | 50 | "\n", |
51 | 51 | "# Setup input domain\n", |
52 | | - "domain = GPflowOpt.domain.ContinuousParameter('x1', -3, 3) + \\\n", |
53 | | - " GPflowOpt.domain.ContinuousParameter('x2', -2, 2)" |
| 52 | + "domain = gpflowopt.domain.ContinuousParameter('x1', -3, 3) + \\\n", |
| 53 | + " gpflowopt.domain.ContinuousParameter('x2', -2, 2)" |
54 | 54 | ] |
55 | 55 | }, |
56 | 56 | { |
|
77 | 77 | }, |
78 | 78 | "outputs": [], |
79 | 79 | "source": [ |
80 | | - "class AugmentedEI(GPflowOpt.acquisition.ExpectedImprovement):\n", |
| 80 | + "class AugmentedEI(gpflowopt.acquisition.ExpectedImprovement):\n", |
81 | 81 | " def __init__(self, model):\n", |
82 | 82 | " super(AugmentedEI, self).__init__(model)\n", |
83 | 83 | "\n", |
|
115 | 115 | } |
116 | 116 | ], |
117 | 117 | "source": [ |
118 | | - "design = GPflowOpt.design.LatinHyperCube(9, domain)\n", |
| 118 | + "design = gpflowopt.design.LatinHyperCube(9, domain)\n", |
119 | 119 | "X = design.generate()\n", |
120 | 120 | "Y = camelback(X)\n", |
121 | | - "m = GPflow.gpr.GPR(X, Y, GPflow.kernels.Matern52(2, ARD=True, lengthscales=[10,10], variance=10000))\n", |
| 121 | + "m = gpflow.gpr.GPR(X, Y, gpflow.kernels.Matern52(2, ARD=True, lengthscales=[10,10], variance=10000))\n", |
122 | 122 | "m.likelihood.variance = 1\n", |
123 | 123 | "m.likelihood.variance.fixed = True\n", |
124 | | - "ei = GPflowOpt.acquisition.ExpectedImprovement(m)\n", |
125 | | - "m = GPflow.gpr.GPR(X, Y, GPflow.kernels.Matern52(2, ARD=True, lengthscales=[10,10], variance=10000))\n", |
| 124 | + "ei = gpflowopt.acquisition.ExpectedImprovement(m)\n", |
| 125 | + "m = gpflow.gpr.GPR(X, Y, gpflow.kernels.Matern52(2, ARD=True, lengthscales=[10,10], variance=10000))\n", |
126 | 126 | "m.likelihood.variance = 1\n", |
127 | 127 | "m.likelihood.variance.fixed = False\n", |
128 | 128 | "aei = AugmentedEI(m)\n", |
129 | 129 | "\n", |
130 | | - "opt = GPflowOpt.optim.StagedOptimizer([GPflowOpt.optim.MCOptimizer(domain, 200), \n", |
131 | | - " GPflowOpt.optim.SciPyOptimizer(domain)])\n", |
| 130 | + "opt = gpflowopt.optim.StagedOptimizer([gpflowopt.optim.MCOptimizer(domain, 200), \n", |
| 131 | + " gpflowopt.optim.SciPyOptimizer(domain)])\n", |
132 | 132 | "\n", |
133 | | - "bopt1 = GPflowOpt.BayesianOptimizer(domain, ei, optimizer=opt)\n", |
| 133 | + "bopt1 = gpflowopt.BayesianOptimizer(domain, ei, optimizer=opt)\n", |
134 | 134 | "with bopt1.silent():\n", |
135 | 135 | " bopt1.optimize(camelback, n_iter=50)\n", |
136 | 136 | "\n", |
137 | | - "bopt2 = GPflowOpt.BayesianOptimizer(domain, aei, optimizer=opt)\n", |
| 137 | + "bopt2 = gpflowopt.BayesianOptimizer(domain, aei, optimizer=opt)\n", |
138 | 138 | "with bopt2.silent():\n", |
139 | 139 | " bopt2.optimize(camelback, n_iter=50)\n", |
140 | 140 | "\n", |
141 | 141 | "f, axes = plt.subplots(1,2, figsize=(14,7))\n", |
142 | 142 | "\n", |
143 | | - "Xeval = GPflowOpt.design.FactorialDesign(101, domain).generate()\n", |
| 143 | + "Xeval = gpflowopt.design.FactorialDesign(101, domain).generate()\n", |
144 | 144 | "Yeval = camelback(Xeval)\n", |
145 | 145 | "titles = ['EI', 'AEI']\n", |
146 | 146 | "shape = (101, 101)\n", |
|
172 | 172 | "name": "python", |
173 | 173 | "nbconvert_exporter": "python", |
174 | 174 | "pygments_lexer": "ipython3", |
175 | | - "version": "3.6.1" |
| 175 | + "version": "3.5.2" |
176 | 176 | } |
177 | 177 | }, |
178 | 178 | "nbformat": 4, |
|
0 commit comments