Coverage for gpaw/test/core/test_interpolate.py: 100%

1import pytest 

2from gpaw.core import UGDesc 

3 

4 

5@pytest.mark.ci 

6def test_fft_interpolation(): 

7 a = UGDesc(cell=[1, 1, 1], size=(4, 4, 4)).zeros() 

8 b = UGDesc(cell=[1, 1, 1], size=(8, 8, 8)).zeros() 

9 a.data[2, 2, 2] = 1.0 

10 a.interpolate(out=b) 

11 assert (b.data[::2, ::2, ::2] == a.data).all() 

12 

13 b.fft_restrict(out=a) 

14 assert a.integrate() == pytest.approx(b.integrate(), abs=1e-12)