Mohit,
can't help you there, but I have exactly the same problem with a
struct file (Mo3Sb7), which I used 2 years ago with the pre-compiled
linux executables.
Now with ifort/mkl, x symmetry (compiled fine thanks to Peter's and
Gerhard's hints) does the same to me, to be precise, the copy
struct_st option empties the struct file except for its first four
lines. This doesn't happen all the time, for example the TiC example
is ok. Kind of, because x symmetry can't DETERMINE LOCAL ROTATION
MATRIX, which it should for this NaCl structure.
On the same machine, the version compiled with gfortran/goto runs
fine. Including Mo3Sb7, TiC, and the local rotation matrices.
If somebody wants to have a look, I am copying the files below.
Holger
[All with wien2k_9.2]
I) TiC.outputs
1) ifort
For Ti:
pointgroup is m3m (pos. iatnr!!)
axes should be: any
z-rotation vector: 0. 0. 1.
y-rotation vector: 0. 0. 0.0
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
For C:
pointgroup is m3m (pos. iatnr!!)
axes should be: any
z-rotation vector: 0. 0. 1.
y-rotation vector: 1. 0. 0.0
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
2) gfortran
For Ti:
pointgroup is m3m (pos. iatnr!!)
axes should be: any
z-rotation vector: 0. 0. 1.
y-rotation vector: 0. 0. 0.0(same vectors as above)
LOCAL ROT MATRIX: NEWOLD
1.000 0.000 0.000 1.000 0.000 0.000
0.000 1.000 0.000 0.000 1.000 0.000
0.000 0.000 1.000 0.000 0.000 1.000
** IATNR IN STRUCT_ST CHANGED TO A POSITIVE NUMBER **
For C:
pointgroup is m3m (pos. iatnr!!)
axes should be: any
z-rotation vector: 0. 0. 1.
y-rotation vector: 0. 0. 0.0 (different y vector)
WARNING: LOCAL ROTATION MATRIX CHANGED
LOCAL ROT MATRIX: NEWOLD
1.000 0.000 0.000 0.000 0.000 0.000
0.000 1.000 0.000 0.000 0.000 0.000
0.000 0.000 1.000 0.000 0.000 0.000
II) Mo3Sb7.outputs
1) ifort
pointgroup is 4mm (neg. iatnr!!)
axes should be: 4 || z, m n y
z-rotation vector: 1. 0. 0.
y-rotation vector: 0. 1. 0.2
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
z-rotation vector: 1. 0. 0.
y-rotation vector: 0. 1. 1.2
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
z-rotation vector: 1. 0. 0.
y-rotation vector: 0. -1. 1.2
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
z-rotation vector: 1. 0. 0.
y-rotation vector: 0. 0. 1.2
WARNING !!! DET(CB)=0
COULD NOT DETERMINE LOCAL ROTATION MATRIX!!!
and end of file
- struct:
Mo3Sb7
B LATTICE,NONEQUIV.ATOMS: 3229_Im-3m
MODE OF CALC=RELA unit=ang
18.063899 18.063899 18.063899 90.00 90.00 90.00
(that's it)
2) gfortran
pointgroup is 4mm (neg. iatnr!!)
axes should be: 4 || z, m n y
z-rotation vector: 1. 0. 0.
y-rotation vector: 0. 1. 0.2
WARNING: LOCAL ROTATION MATRIX CHANGED
LOCAL ROT MATRIX: NEWOLD
0.000 0.000 1.000 1.000 0.000 0.000
0.000 1.000 0.000 0.000 1.000 0.000
-1.000 0.000 0.000 0.000 0.000 1.000
- struct:
Mo3Sb7
B LATTICE,NONEQUIV.ATOMS: 3229_Im-3m
MODE OF CALC=RELA unit=ang
18.063899 18.063899 18.063899 90.00 90.00
90.00
ATOM -1: X=0.3432 Y=0. Z=0.
MULT= 6 ISPLIT=-2
-1: X=0.6568 Y=0. Z=0.
-1: X=0. Y=0.3432 Z=0.
-1: X=0. Y=0.6568 Z=0.
-1: X=0. Y=0. Z=0.3432
-1: X=0. Y=0. Z=0.6568
Mo1NPT= 781 R0=0.0001 RMT=2.5000 Z: 42.0
LOCAL ROT MATRIX:0.000 0.000 1.000
0.000 1.000 0.000
-1.000 0.000 0.000
ATOM -2: X=0.2500 Y=0. Z=0.5000
MULT= 6 ISPLIT=-2
-2: X=0.7500 Y=0. Z=0.5000
-2: X=0. Y=0.2500 Z=0.5000
-2: X=0. Y=0.7500 Z=0.5000
-2: X=0.5000 Y=0. Z=0.2500
-2: X=0.5000 Y=0. Z=0.7500
Sb1NPT= 781 R0=0.0001 RMT=2.5000 Z: 51.0
LOCAL ROT MATRIX:0.000 0.000 1.000