[ase-users] how to get the pressure

[陈冠宇]

hi ,I have some  questions about ASE.how can i get the pressure of system ,such as h2.
my scribe is following:
 
from ase import Atoms
from gpaw import GPAW
a = Atoms('H2', positions=[(0, 0, 0), (0, 0, 1.1)],
           calculator=GPAW())
stress=a.get_stress()
pressure=get_isotropic_pressure()
print stress ,pressure
but it does not work.
what should i do?
print h2.get_forces()
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listserv.fysik.dtu.dk/pipermail/ase-users/attachments/20110420/a00e61e4/attachment.html>