Hello, Is there any way of incorporating mass matrix for transient problems solved with built-in multistage schemes like:
scheme = BDF1(F, u, t, [bc]) solver = RKSolver(scheme)
Do i have to invert it manually and multiply by variational form F?
