root / ase / test / exciting.py @ 1
Historique | Voir | Annoter | Télécharger (733 octet)
1 | 1 | tkerber | import os |
---|---|---|---|
2 | 1 | tkerber | from ase import Atoms |
3 | 1 | tkerber | from ase.io import read, write |
4 | 1 | tkerber | from ase.calculators import Exciting |
5 | 1 | tkerber | from ase.units import Bohr, Hartree |
6 | 1 | tkerber | from ase.test import NotAvailable |
7 | 1 | tkerber | |
8 | 1 | tkerber | try:
|
9 | 1 | tkerber | import lxml |
10 | 1 | tkerber | except ImportError: |
11 | 1 | tkerber | raise NotAvailable('This test need lxml module.') |
12 | 1 | tkerber | |
13 | 1 | tkerber | a = Atoms('N3O',
|
14 | 1 | tkerber | [(0, 0, 0), (1, 0, 0), (0, 0, 1), (0.5, 0.5, 0.5)], |
15 | 1 | tkerber | pbc=True)
|
16 | 1 | tkerber | |
17 | 1 | tkerber | raise NotAvailable('Problem with lxml module.') |
18 | 1 | tkerber | |
19 | 1 | tkerber | write('geo.exi', a)
|
20 | 1 | tkerber | b = read('geo.exi')
|
21 | 1 | tkerber | |
22 | 1 | tkerber | print a
|
23 | 1 | tkerber | print a.get_positions()
|
24 | 1 | tkerber | print b
|
25 | 1 | tkerber | print b.get_positions()
|
26 | 1 | tkerber | |
27 | 1 | tkerber | calculator = Exciting(dir='excitingtestfiles',
|
28 | 1 | tkerber | kpts=(4, 4, 3), |
29 | 1 | tkerber | maxscl=3,
|
30 | 1 | tkerber | #bin='/fshome/chm/git/exciting/bin/excitingser'
|
31 | 1 | tkerber | ) |