Method test_implicit_stability

tests/commit/physics/test_diffuse.py:57–66

(self)

        field.assert_close(grid, back_explicit, back_implicit, back_fourier, rel_tolerance=0, abs_tolerance=0.1)

    def test_implicit_stability(self):
        DIFFUSIVITY = 10
        grid = CenteredGrid((1,) * 3 + (0,) * 3, extrapolation.PERIODIC, x=21)
        try:
            implicit = diffuse.implicit(grid, DIFFUSIVITY, 1)
            print(implicit.values)
            field.assert_close(0 <= implicit <= 1.0001, True)
        except NotConverged as err:
            print(err)
            pass  # solve_linear did not converge

    def test_explicit_centered_non_isotropic(self):
        grid = CenteredGrid(wrap([[0, 0, 0], [0, 1, 0], [0, 0, 0]], spatial('x,y')), 0)