File: test_h2.py

package info (click to toggle)
python-ase 3.26.0-2
  • links: PTS, VCS
  • area: main
  • in suites:
  • size: 15,484 kB
  • sloc: python: 148,112; xml: 2,728; makefile: 110; javascript: 47
file content (13 lines) | stat: -rw-r--r-- 345 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
# fmt: off
from ase import Atoms
from ase.calculators.emt import EMT
from ase.calculators.fd import calculate_numerical_forces


def test_h2():

    h2 = Atoms('H2', positions=[(0, 0, 0), (0, 0, 1.1)],
               calculator=EMT())
    f1 = calculate_numerical_forces(h2, 0.0001)
    f2 = h2.get_forces()
    assert abs(f1 - f2).max() < 1e-6