<br><font size=2 face="sans-serif">Thank you for you quick answer.</font>
<br>
<br><font size=2 face="sans-serif">I enclose the file, that is the exact
file provided in the downloaded package to run the example (14 in versin
4.2.1), with only a small modification.</font>
<br><font size=2 face="sans-serif">Indeed, as you will see, I just substitute
the X-points coordinates by some other coordinates, for the d3.x (line
199) and the phonon calculation (line 181).</font>
<br><font size=2 face="sans-serif">Or perhaps this is only why I think
I am doing...</font>
<br>
<br><font size=2 face="sans-serif">Once again, I have no trumble at all
to run it on k-points belonging to high symmetry axis (1,0,0   0.5,0.5,0.5,
  0.2,0.2,0.2   0.3,0.3,0    ....).</font>
<br><font size=2 face="sans-serif">The k-point here in the file is 0.3,0.2,0.1.</font>
<br>
<br><font size=2 face="sans-serif">Regards<br>
<br>
<br>
</font>
<table>
<tr valign=top>
<td rowspan=2><img src=cid:_2_0785B0100785AC28003124B6C1257B48>
<td><font size=1 face="sans-serif"> </font>
<tr valign=top>
<td><font size=1 color=#ff8100 face="Arial"><b>Hugo LEVARD</b></font><font size=1 color=#ff8100 face="Arial"><b><br>
Doctorant</b></font><font size=1 color=#0062e1 face="Arial"><br>
EDF – R&D<br>
IRDEP - EFESE<br>
6 Quai Watier - BP 49<br>
78401 Chatou</font><font size=1 color=#0062e1 face="sans-serif"><br>
 </font>
<br><font size=1 color=#0062e1 face="Arial"><b>hugo.levard@edf.fr</b></font>
<br><font size=1 color=#0062e1 face="Arial">Tél. : 01 30 87 90 28<br>
Fax : 01 30 87 85 65</font>
<tr>
<td>
<div align=right><img src=cid:_1_0785BECC0785BB10003124B6C1257B48></div>
<td><font size=1 color=#0062e1 face="Arial">Un geste simple pour l'environnement,
n'imprimez ce message que si vous en avez l'utilité.</font></table>
<br>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>lorenzo.paulatto@impmc.upmc.fr</b>
</font>
<br><font size=1 face="sans-serif">Envoyé par : pw_forum-bounces@pwscf.org</font>
<p><font size=1 face="sans-serif">09/04/2013 10:33</font>
<table border>
<tr valign=top>
<td bgcolor=white>
<div align=center><font size=1 face="sans-serif">Veuillez répondre à<br>
lorenzo.paulatto@impmc.upmc.fr; Veuillez répondre à<br>
pw_forum@pwscf.org</font></div></table>
<br>
<td width=59%>
<table width=100%>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">A</font></div>
<td><font size=1 face="sans-serif">pw_forum@pwscf.org</font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">Objet</font></div>
<td><font size=1 face="sans-serif">Re: [Pw_forum] d3.x on various q-points</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br><font size=3>On 04/09/2013 09:07 AM, Hugo LEVARD wrote:</font>
<br><font size=2 face="sans-serif">I am trying to understand an error occuring
on third order anharmonic calculation with d3.x.</font><font size=3> </font><font size=2 face="sans-serif"><br>
The default example works perfectly (q=(1 0 0)), but when I just change
the q-point on which the calculation is done, I do not always get a result.</font><font size=3>
</font><font size=2 face="sans-serif"><br>
With, for instance, q=(0.3 0.2 0.1)</font><font size=3> </font><font size=2 face="sans-serif"><br>
+ In version 4.3.2, the tensor components are simply "NaN".</font><font size=3>
</font><font size=2 face="sans-serif"><br>
+ In version 5.0.2, the calculation fail with an I/O error ("error
in routine davcio").</font><font size=3> </font><font size=2 face="sans-serif"><br>
From the panel of q -points I have tried, it seems to me that this problem
arises when leaving the high symmetry line in the BZ.</font><font size=3>
</font>
<br><font size=3><br>
Dear Hugo,<br>
the d3.x calculation has to be done following a strict order. Did you precede
the d3.x calculation with a phonon calculation at the same point? With
the same prefix and outdir?<br>
<br>
Anyway, please send me the three inputs (pw, ph, d3) as well as any script
you are using to run them and I'll try to reproduce your problem: without
it is impossible to say anything definitive.<br>
<br>
cheers<br>
</font>
<br><tt><font size=3>-- <br>
Dr. Lorenzo Paulatto <br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
phone: +33 (0)1 44275 084 / skype: paulatz<br>
www:   </font></tt><a href="http://www-int.impmc.upmc.fr/~paulatto/"><tt><font size=3 color=blue><u>http://www-int.impmc.upmc.fr/~paulatto/</u></font></tt></a><tt><font size=3><br>
mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris
Cédex 05<br>
</font></tt><tt><font size=2>_______________________________________________<br>
Pw_forum mailing list<br>
Pw_forum@pwscf.org<br>
http://pwscf.org/mailman/listinfo/pw_forum</font></tt>
<br>