Function _stop

python/simulation.py:5300–5321

(sim)

    }

    def _stop(sim):
        fabs = abs(sim.get_field_point(c, pt)) * abs(sim.get_field_point(c, pt))
        closure["cur_max"] = max(closure["cur_max"], fabs)

        if sim.round_time() <= dt + closure["t0"]:
            return False
        else:
            old_cur = closure["cur_max"]
            closure["cur_max"] = 0
            closure["t0"] = sim.round_time()
            closure["max_abs"] = max(closure["max_abs"], old_cur)
            if closure["max_abs"] != 0 and verbosity.meep > 0:
                fmt = "field decay(t = {}): {} / {} = {}"
                print(
                    fmt.format(
                        sim.meep_time(),
                        old_cur,
                        closure["max_abs"],
                        old_cur / closure["max_abs"],
                    )
                )
            return old_cur <= closure["max_abs"] * decay_by

    return _stop