<div dir="ltr"><div>It is a known problem that was raised in the pw_forum mailing list. I have no strong opinions about this (well, I have one on some compilers, but this is of no help) but I tend to agree that the best (or the only) solution is to put the routine dv_of_drho into a module. <br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 11, 2016 at 5:29 PM, Timrov Iurii <span dir="ltr"><<a href="mailto:iurii.timrov@epfl.ch" target="_blank">iurii.timrov@epfl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Lorenzo,<br>
<br>
You are totally right! In fact, one of the QE users also noticed such a problem with the routine dv_of_drho when using xlf.<br>
<br>
It is my mistake to introduce such a bug in dv_of_drho. Sorry about that!<br>
<br>
In my opinion, the first option which you suggest is better. Personally I do not like the second option, because in several places of PHonon and in TDDFPT the routine dv_of_drho is called without the need to pass the response core charge density "drhoc". So, in the case of the second option, one would need to pass a zero array drhoc=(0.0d0, 0.0d0) to dv_of_drho, which has no sense and moreover it consumes some memory without any need (add_nlcc=.false. in that case, so such a zero array will not be used).<br>
<br>
Best regards,<br>
Iurii<br>
<br>
<br>
--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
Swiss Federal Institute of Technology Lausanne (EPFL)<br>
Laboratory of Theory and Simulation of Materials (THEOS)<br>
CH-1015 Lausanne, Switzerland<br>
<br>
<br>
________________________________________<br>
From: <a href="mailto:q-e-developers-bounces@qe-forge.org">q-e-developers-bounces@qe-forge.org</a> <<a href="mailto:q-e-developers-bounces@qe-forge.org">q-e-developers-bounces@qe-forge.org</a>> on behalf of Lorenzo Paulatto <<a href="mailto:lorenzo.paulatto@impmc.upmc.fr">lorenzo.paulatto@impmc.upmc.fr</a>><br>
Sent: Wednesday, May 11, 2016 5:00 PM<br>
To: discussion developers<br>
Subject: [Q-e-developers] dv_of_drho optional arguments<br>
<div class="HOEnZb"><div class="h5"><br>
Hello,<br>
I have noticed that the subroutine dv_of_drho in LR_Modules has an optional<br>
argument. However, this subroutine does not have an explicit interface (i.e.<br>
it is not in a module), which can cause problems. It does indeed cause<br>
problems with xlf.<br>
<br>
There are two possible solutions:<br>
1. put the subroutine in a module<br>
2. make the argument non optional<br>
<br>
As the subroutine is used in several codes (phonon, tddft) I want to ask your<br>
opinion before proceeding<br>
<br>
cheers<br>
<br>
--<br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
<a href="tel:%2B33%20%280%291%2044%20275%20084" value="+33144275084">+33 (0)1 44 275 084</a> / skype: paulatz<br>
<a href="http://www.impmc.upmc.fr/~paulatto/" rel="noreferrer" target="_blank">http://www.impmc.upmc.fr/~paulatto/</a><br>
23-24/4é16 Boîte courrier 115,<br>
4 place Jussieu 75252 Paris Cédex 05<br>
<br>
_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
<br>
_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>