[Pw_forum] A bug in the pp.x?

[Lorenzo Paulatto]

On Sun, June 8, 2008 13:46, zhaohscas wrote:
>  Hi forks,
Hi chopstick ;-),

>  In the abobe case, the e1 \bot e2, and e2 \bot e3, but e1 isn'
>  perpendicular to e3.  So, this will violate orthogonal condition.

You are right: the actual check performed by the code is:
  e_i \cdot \e_j > 1.d-6
which wrongly allows non-orthogonal vectors if their pairwise scalar
products are negative. Thank you for reporting, it will be fixed in the
next release.

bye

-- 
Lorenzo Paulatto
SISSA  &  DEMOCRITOS (Trieste)
+39 040 3787 511
http://people.sissa.it/~paulatto/


----------------------------------------------------------------
  SISSA Webmail https://webmail.sissa.it/
  Powered by SquirrelMail http://www.squirrelmail.org/