@@ -9,7 +9,13 @@ if (knitr::is_html_output()) {
99}
1010```
1111
12- ## integrate_ode_rk45, integrate_ode_adams, integrate_ode_bdf ODE integrators {#functions-old-ode-solver}
12+ ``` {r results='asis', echo=FALSE}
13+ if (knitr::is_html_output()) {
14+ cat(' * <a href="cov_exp_quad.html">cov_exp_quad</a>\n')
15+ }
16+ ```
17+
18+ ## integrate_ode_rk45, integrate_ode_adams, integrate_ode_bdf ODE Integrators {#functions-old-ode-solver}
1319
1420These ODE integrator functions have been replaced by those described in:
1521
@@ -24,7 +30,7 @@ if (knitr::is_html_output()) {
2430A system of ODEs is specified as an ordinary function in Stan within
2531the functions block. The ODE system function must have this function
2632signature:
27-
33+
2834``` stan
2935array[] real ode(real time, array[] real state, array[] real theta,
3036 array[] real x_r, array[] int x_i);
@@ -117,7 +123,7 @@ are as follows.
117123
118124* * ` ode ` * : function literal referring to a function specifying the
119125system of differential equations with signature:
120-
126+
121127```
122128(real, array[] real, array[] real, data array[] real, data array[] int):array[] real
123129```
@@ -140,7 +146,7 @@ derivatives with respect to time of the state,
140146
141147For more fine-grained control of the ODE solvers, these parameters can
142148also be provided:
143-
149+
144150 * ` data ` * ` rel_tol ` * : relative tolerance for the ODE solver, type
145151` real ` , data only,
146152
@@ -159,7 +165,7 @@ with values consisting of solutions at the specified times.
159165
160166The sizes must match, and in particular, the following groups are of
161167the same size:
162-
168+
163169 * state variables passed into the system function, derivatives
164170returned by the system function, initial state passed into the
165171solver, and rows of the return value of the solver,
@@ -169,3 +175,65 @@ solver,
169175
170176* parameters, real data and integer data passed to the solver will
171177be passed to the system function
178+
179+
180+ ## Exponentiated quadratic covariance functions {#cov_exp_quad}
181+
182+ These covariance functions have been replaced by those described in:
183+
184+ ``` {r results='asis', echo=FALSE}
185+ if (knitr::is_html_output()) {
186+ cat(' * <a href="gaussian-process-covariance-functions.html">Gaussian Process Covariance Functions</a>\n')
187+ }
188+ ```
189+
190+ With magnitude $\alpha$ and length scale $l$, the exponentiated quadratic kernel is:
191+
192+ $$
193+ k(x_i, x_j) = \alpha^2 \exp \left(-\dfrac{1}{2\rho^2} \sum_{d=1}^D (x_{i,d} - x_{j,d})^2 \right)
194+ $$
195+
196+ <!-- matrix; cov_exp_quad; (row_vectors x, real alpha, real rho); -->
197+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (row\_ vectors x, real alpha, real rho): matrix}|hyperpage}
198+
199+ ` matrix ` ** ` cov_exp_quad ` ** ` (row_vectors x, real alpha, real rho) ` <br >\newline
200+ The covariance matrix with an exponentiated quadratic kernel of x.
201+ ` r since("2.16") `
202+
203+ <!-- matrix; cov_exp_quad; (vectors x, real alpha, real rho); -->
204+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (vectors x, real alpha, real rho): matrix}|hyperpage}
205+
206+ ` matrix ` ** ` cov_exp_quad ` ** ` (vectors x, real alpha, real rho) ` <br >\newline
207+ The covariance matrix with an exponentiated quadratic kernel of x.
208+ ` r since("2.16") `
209+
210+ <!-- matrix; cov_exp_quad; (array[] real x, real alpha, real rho); -->
211+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (array[ ] real x, real alpha, real rho): matrix}|hyperpage}
212+
213+ ` matrix ` ** ` cov_exp_quad ` ** ` (array[] real x, real alpha, real rho) ` <br >\newline
214+ The covariance matrix with an exponentiated quadratic kernel of x.
215+ ` r since("2.16") `
216+
217+ <!-- matrix; cov_exp_quad; (row_vectors x1, row_vectors x2, real alpha, real rho); -->
218+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (row\_ vectors x1, row\_ vectors x2, real alpha, real rho): matrix}|hyperpage}
219+
220+ ` matrix ` ** ` cov_exp_quad ` ** ` (row_vectors x1, row_vectors x2, real alpha, real rho) ` <br >\newline
221+ The covariance matrix with an exponentiated quadratic kernel of x1 and
222+ x2.
223+ ` r since("2.18") `
224+
225+ <!-- matrix; cov_exp_quad; (vectors x1, vectors x2, real alpha, real rho); -->
226+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (vectors x1, vectors x2, real alpha, real rho): matrix}|hyperpage}
227+
228+ ` matrix ` ** ` cov_exp_quad ` ** ` (vectors x1, vectors x2, real alpha, real rho) ` <br >\newline
229+ The covariance matrix with an exponentiated quadratic kernel of x1 and
230+ x2.
231+ ` r since("2.18") `
232+
233+ <!-- matrix; cov_exp_quad; (array[] real x1, array[] real x2, real alpha, real rho); -->
234+ \index{{\tt \bfseries cov\_ exp\_ quad }!{\tt (array[ ] real x1, array[ ] real x2, real alpha, real rho): matrix}|hyperpage}
235+
236+ ` matrix ` ** ` cov_exp_quad ` ** ` (array[] real x1, array[] real x2, real alpha, real rho) ` <br >\newline
237+ The covariance matrix with an exponentiated quadratic kernel of x1 and
238+ x2.
239+ ` r since("2.18") `
0 commit comments