The function returned by jit_compile(), jit_compile_linear(), functional_graident(), custom_gradient() now dynamically choose the backend when invoked.
Overhaul solve API
Add TensorLike, implemented by Fields, Extrapolations and more.
Remove simple field math functions, directly call math functions instead.
Update demos
Changes to linear solves
Move implementation to Backend
Add 'CG-adaptive' method
Initial guesses for solves are now part of the Solve object
Add record_trajectory and suppress errors arguments to Solve
Solves now return a SolveResult object which is no longer referenced by Solve
Changes to linear solves