[Pw_forum] mistake in impose_deviatoric_strain_2d

[Niels Walet]

Hope this helps someone:
in pw.x (espresso 5.1.2), with
                  cell_dofree = '2Dshape' ,
the area (and volume) of the cell are actually not conserved.

Digging around a bit, it looks like the power is wrong in flib impose_deviatoric_strain_2d:
In line 59, the power of 1/3 should be 1/2 (i.e., the final 3.d0 should be 2.d0), otherwise area is not conserved!!!
Niels

---
Prof. Niels R. Walet                             Phone:  +44(0)1613063693
School of Physics and Astronomy         Fax:    +44(0)1613064303
The University of Manchester               Mobile: +44(0)7516622121
Manchester, M13 9PL,  UK                  room 7.7, Schuster Building
email: Niels.Walet at manchester.ac.uk
web:   http://www.theory.physics.manchester.ac.uk/~mccsnrw<../../owa/redir.aspx?C=8f9098c1943846dcba17536b8618a38a&URL=http%3a%2f%2fwww.theory.physics.manchester.ac.uk%2f%7emccsnrw>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20150814/3f5dfe41/attachment.html>