<div dir="ltr"><br><div class="gmail_extra">Thank you very much Prof Gavin Abo for nice explanation.<br></div><div class="gmail_extra">Now, I got it.<br></div><div class="gmail_extra"><br><br></div><div class="gmail_extra"><div><div><div dir="ltr"><div><div dir="ltr"><div><div>Bhamu<br></div><div dir="ltr"><b><span style="color:rgb(0,0,255)"></span></b></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Fri, Sep 4, 2015 at 1:53 AM, Gavin Abo <span dir="ltr"><<a href="mailto:gsabo@crimson.ua.edu" target="_blank">gsabo@crimson.ua.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<font color="#000099">See comments below:<br>
<br>
</font>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div><b>For TiC</b> (I already posted this
query few month ago and got a quick reply
from Prof. Marks)<br>
<br>
<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg12624.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg12624.html</a><br>
<br>
</div>
<div dir="ltr">Reply from Prom Marks was<br>
<br>
"You do not have a problem! The charge
convergence is very sensitive to
numerical issues, so using -cc 0.0001
demands a lot and is probably not a
good idea"<br>
<pre></pre>
<pre>I ran calculation via w2web and also with terminal with -cc 0.0001 (suggested in UG) or -ec 0.0001.
</pre>
<pre>It conversed again in 14 cycles.
</pre>
<pre>Any comment why it is not conversing at 11 iterations? </pre>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
<font color="#000099">This should be simply because your calculation
takes 14 cycles to satisfy your convergence criteria (-cc 0.0001
or -ec 0.0001).</font><br>
<br>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div dir="ltr">
<pre>I want it in with 11 iterations.
</pre>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
<font color="#000099">You could force it to stop in the 11th
iteration by using the -i option (e.g., run_lapw -i 11 </font><font color="#000099">-cc 0.0001), but then it would not reach </font><font color="#000099">the convergence criteria (-cc 0.0001).<br>
<br>
Not being able to reproduce </font><font color="#000099"><font color="#000099">'exactly' the results of</font></font><font color="#000099"><font color="#000099"> </font>the examples like
TiC in the userguide is normal.<br>
<br>
The TiC calculation in the usersguide should be from one of the
old WIEN2k versions. It was also likely done using an old version
of ifort/mkl.<br>
<br>
The reason it converges in 14 instead of 11 cycles is likely
because a combination of:<br>
<br>
a. You might have slightly different input parameters somewhere.<br>
<br>
b. Changes to the WIEN2k code (fixes and improvements have been
made since then to programs like lapw1, lapw2, mixer [
<a href="http://www.wien2k.at/reg_user/updates/" target="_blank">http://www.wien2k.at/reg_user/updates/</a> ]; the examples might have
been produced without the -in1ef switch behavior, but it is now
coded to always use the -in1ef behavior [
<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg08632.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg08632.html</a>
]).<br>
<br>
c. Numerical result</font><font color="#000099"><font color="#000099"> </font>differences </font><font color="#000099"><font color="#000099"><font color="#000099">(from
rounding errors) </font></font>due to hardware and software
(compiler and math libraries) [
<a href="https://software.intel.com/en-us/articles/intel-composer-xe-2015-compilers-fixes-list" target="_blank">https://software.intel.com/en-us/articles/intel-composer-xe-2015-compilers-fixes-list</a>
, <a href="https://software.intel.com/en-us/forums/topic/299617" target="_blank">https://software.intel.com/en-us/forums/topic/299617</a> ,
<a href="http://science.slashdot.org/story/13/07/28/137209/same-programs--different-computers--different-weather-forecasts" target="_blank">http://science.slashdot.org/story/13/07/28/137209/same-programs--different-computers--different-weather-forecasts</a>
].</font><br>
</div>
<br>_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br></blockquote></div><br></div></div>