root / ase / test / rotate.py @ 4
Historique | Voir | Annoter | Télécharger (336 octet)
| 1 | 1 | tkerber | from ase.utils import rotate, irotate |
|---|---|---|---|
| 2 | 1 | tkerber | |
| 3 | 1 | tkerber | def test(xyz): |
| 4 | 1 | tkerber | a = rotate(xyz) |
| 5 | 1 | tkerber | ixyz = '%sx,%sy,%sz' % irotate(a)
|
| 6 | 1 | tkerber | a2 = rotate(ixyz) |
| 7 | 1 | tkerber | print xyz
|
| 8 | 1 | tkerber | print ixyz
|
| 9 | 1 | tkerber | #print np.around(a-a2, 5)
|
| 10 | 1 | tkerber | assert abs(a-a2).max() < 1e-10 |
| 11 | 1 | tkerber | |
| 12 | 1 | tkerber | test('10z')
|
| 13 | 1 | tkerber | test('155x,43y,190z')
|
| 14 | 1 | tkerber | test('55x,90y,190z')
|
| 15 | 1 | tkerber | test('180x,-90y,45z')
|
| 16 | 1 | tkerber | test('-180y')
|
| 17 | 1 | tkerber | test('40z,50x') |