Use of save_idxs solver option with MethodOfLines/ModelingToolkit

Open an issue. I don’t think that works with the PDE indexing right now.