Hi, Flo Yes. "cp_fm_invert" is only called by electronic coupling subroutines. One concern I have is on the value of the projection matrix (W) elements, which certainly are influenced by the sign of molecular overlap SAB, the determinant of orbital overlap matrix. Thanks. Hanning <div class="gmail_quote">On Wed, Nov 24, 2010 at 2:55 AM, flo <<a href="mailto:fsch...@pci.uzh.ch">fsch...@pci.uzh.ch</a>> wrote: <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"> Hi Hanning, I just noticed, when I wrote the subroutine, I committed the wrong description for the subroutine. I will update this soon. Indeed the sign can be wrong, as it depends on the entries in ipivot. As far as I know CP2K uses the determinant only for computing the electronic coupling matrix elements and here it is done by calling the pure serial routine in mathlib.f, where it is corrected. In all other cases the sign of the determinant does not matter. In principle it is possible to correct this bug according to the way it is done in mathlib.f. The point why it hasn't been done yet is, that it wasn't needed and involves quite a bit of thinking about the transformation between local positions in ipivot and the global positions. Cheers Flo -- You received this message because you are subscribed to the Google Groups "cp2k" group. To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>. To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>. For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>. </blockquote></div>