Coverage for gpaw/test/radial/test_integral4.py: 100%

1import pytest 

2import numpy.random as ra 

3from gpaw.setup import create_setup 

4from gpaw.xc import XC 

5 

6 

7@pytest.mark.ci 

8def test_radial_integral4(): 

9 rng = ra.default_rng(8) 

10 xc = XC('LDA') 

11 s = create_setup('H', xc) 

12 ni = s.ni 

13 nii = ni * (ni + 1) // 2 

14 D_p = 0.1 * rng.random((1, nii)) + 0.2 

15 

16 def f(x): 

17 return x 

18 

19 s.xc_correction.four_phi_integrals(D_p, f) 

20 

21 # Check integrals using two_phi_integrals function and finite differences: 

22 pass