#773: McLachlan dissipation should be disabled by default --------------------------+------------------------------------------------- Reporter: hinder | Owner: Type: enhancement | Status: new Priority: major | Milestone: Component: Other | Version: Resolution: | Keywords: --------------------------+-------------------------------------------------
Comment (by eschnett):
What about modifying the way in which McLachlan iterates over variables? One calculation per variable should be easy to realise: {{{ dissCalcs = Table[ { Name -> BSSN <> "_Dissipation_" <> ToString[var /. {Tensor[n_,__] -> n}], Schedule -> {"IN " <> BSSN <> "_evolCalcGroup " <> "AFTER (" <> BSSN <> "_RHS1 " <> BSSN <> "_RHS2)"}, ConditionalOnKeyword -> {"apply_dissipation", "always"}, Where -> InteriorNoSync, Shorthands -> {epsdiss[ua]}, Equations -> { epsdiss[ua] -> EpsDiss, dot[var] -> dot[var] + epsdiss[ux] PDdiss[var,lx] } }, {var, {phi, gt[la,lb], Xt[ui], trK, At[la,lb], alpha, A, beta[ua], B[ua]}} ]; }}}
Otherwise, this is okay to change.
Do you then also want to change the default epsdiss value to be the same as for thorn Dissipation?