[Wien] Patch to fix ROTDEF error in lapw0 for Fd3m et al.

Martin Kroeker martin at ruby.chemie.uni-freiburg.de
Sat Jan 23 20:02:38 CET 2021


We have recently encountered an unexpected problem with lapw0 returning
"ROTDEF error", claiming "symmetry operation not found" for equivalent
positions of an X,X,X site in spacegroup Fd-3m. (One simplified example
is attached to this message to allow reproducing the problem.)
As it appears, SRC_lapw0/rotdef.f takes turns applying the centering
operations to the calculated distance vector looking for a match - and
then undoes that operation before trying the next, instead of simply
copying back the original value. This demonstrably leads to rounding
errors that can exceed the hardcoded tolerance threshold of 1.5e-7

The attached patch avoids the problem by trivially saving x1 before
the cascade of tests for centered lattices, and restoring it after
each that did not lead to a successful exit from the loop.
-- 
Dr. Martin Kroeker            martin at ruby.chemie.uni-freiburg.de
c/o Prof.Dr. Caroline Roehr
Institut fuer Anorganische und Analytische Chemie der Universitaet Freiburg

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: rotdef.struct
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20210123/6189ebd3/attachment.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: lapw0_rotdef.patch
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20210123/6189ebd3/attachment-0001.ksh>


More information about the Wien mailing list