virtual void cycle(int e) final { n[e] = (n[e] + 1) % n_tlev - n_tlev; // -n_tlev so that n+1 does not give out of bounds if (e == n_eqns - 1) mem->cycle(rank); }