Method test_constant_diffusion

tests/commit/physics/test_diffuse.py:30–35

(self)

            # diffuse.implicit(grid, diffusivity, 1, order=2)  # not yet supported

    def test_constant_diffusion(self):
        grid = CenteredGrid(1, extrapolation.PERIODIC, x=5, y=5)
        explicit = diffuse.explicit(grid, 1, 1, substeps=10)
        implicit = diffuse.implicit(grid, 1, 1)
        fourier = diffuse.fourier(grid, 1, 1)
        math.assert_close(grid.values, explicit.values, implicit.values, fourier.values)

    def test_equality_1d_periodic(self):
        DIFFUSIVITY = 0.5