Hello All,
Symmetric structures may give non-symmetric interaction
matrices, something that is unexpected and suggests a
bug somewhere. The effect as I have seen it is small, it
would disappear in the rounding errors if computations
were made with 4 byte floating point.
The asymmetry is caused by the normalization of
QJRATI(I) in the routine QJRATV.
QJRATI(1) is set to (1.0, 0.0) and the higher members
are scaled by the original complex value of QJRATI(1).
This way both amplitude and phase are scaled. When
QJRATI is used in JNQSET, only the real part is used
and that causes the symmetry error.
Symmetry is restored if QJRATI is replaced by its
magnitude but I do not quite understand what the code
is intended to do so I do not know if that would
be the correct fix.
Is this a known bug?
Leif Asbrink
-- The NEC-List mailing list NEC-List_at_robomod.net http://www.robomod.net/mailman/listinfo/nec-listReceived on Wed Mar 28 2007 - 14:33:32 EDT
This archive was generated by hypermail 2.2.0 : Sat Oct 02 2010 - 00:10:46 EDT