`get_discrete`

This method is deprecated in favour of the newer solution interface, and has much worse performance. These methods will now only work with `sol.original_sol`

.

`MethodOfLines.jl`

exports a helper function `get_discrete`

, which returns a `Dict`

with the keys being the independent and dependent variables, and the values their corresponding discrete grid, and discretized variables used in the discretization. It is used as following:

```
grid = get_discrete(pdesys, discretization)
discrete_x = grid[x]
# Retrieve shaped solution
u_sol = [map(d -> sol[d][i], grid[u(t, x)]) for i in 1:length(sol[t])]
```