- [ ] Create a factory for solvers, and have that factory provide the `adept::Stack`. - [ ] Eliminate the serial version of mach (i.e. get rid of the relevant cpp define guards) - [x] Merge `dev` into `euler` fully - [ ] Consider using `res->GetEnergy()` for the change in entropy calculation. - [ ] We need a better way of updating the `dt` in the nonlinear mass integrator (does not generalize right now) - [ ] turn sandbox examples into regression tests before switching to parallel mfem. (See @tuckerbabcock comment below) - [ ] Remove steady option, and instead use a new ODE solver to implement Newton and pseudo-transient continuation. - [ ] The bilinear and nonlinear forms are initialized (i.e. pointers are not nullptr) regardless of whether or not they are used. Can we avoid this?
adept::Stack.devintoeulerfullyres->GetEnergy()for the change in entropy calculation.dtin the nonlinear mass integrator (does not generalize right now)