get_discrete
This method is deprecated in favor 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 of their corresponding discrete grid, and discretized variables used in the discretization. It is used as follows:
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])]