Dear Prof. Hutter,<br><br>Many thanks,<br><br>Stas<br><br><div class="gmail_quote">On 16 June 2010 15:31, Juerg Hutter <span dir="ltr"><<a href="mailto:hut...@pci.uzh.ch">hut...@pci.uzh.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi<br>
<br>
the problem was related to the calculation of electrons<br>
(= occupied states). We used a piece of code like<br>
 N = INT(A-C)<br>
where A = B + C. This is numerically not stable. I replaced<br>
it with<br>
 N = NINT(A-C).<br>
<br>
regards<br>
<br>
Juerg Hutter<br>
<br>
----------------------------------------------------------<br>
Juerg Hutter                   Phone : ++41 44 635 4491<br>
Physical Chemistry Institute   FAX   : ++41 44 635 6838<br>
University of Zurich           E-mail: <a href="mailto:hut...@pci.uzh.ch" target="_blank">hut...@pci.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zurich, Switzerland<br>
----------------------------------------------------------<div><div></div><div class="h5"><br>
<br>
<br>
On Sun, 13 Jun 2010, Stanislav Avdoshenko wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div></div><div class="h5">
Dear all,<br>
<br>
I have faced to a problem of some "magic numbers" in the core_correction(CC,<br>
- PL1 variables ) function. In attached you can find an example. In this<br>
run, as you can see, is three different values of PL variable:  0.71, 0.72<br>
and 0.73.  For 0.71 and 0.73 according to my  outputs I had 5 electrons on<br>
net, but for 0.72 I had 4.  what kind of the re-normalization it might be,<br>
if for 0.7200001 I have 5 electrons as well? I have tried to find an answer<br>
in Google.net, but  without much results, I hope that community might be<br>
able to help me.<br>
<br>
thanks in advance, Stas<br>
<br></div></div>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
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>.<br>
<br>
<br>
</blockquote>
<br>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
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>.<br>
<br>
</blockquote></div><br>