Only used/save cached bem matrices explicitly
Add output of Mallinson switching time error for ODE LLG