<div dir="ltr"><div>Hi Leopold, thank you for the patch. I just committed it, removed the undocumented and presumably unused case plot_num=13,14,15 and related variable "epsilon", fixed example05. The problem in example06 shows up only in parallel execution and is an infamous "S matirx not positive definite", for which no easy solution is available, I am afraid<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 6, 2017 at 11:11 PM, Leopold Talirz <span dir="ltr"><<a href="mailto:leopold.talirz@gmail.com" target="_blank">leopold.talirz@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Paolo,<br>
<br>
thanks for moving forward so quickly with the patch.<br>
I have now verified that the specific namelist behavior required for the patch is provided by the recent gfortran, ifort, pgf90 and nag compilers (although I have compiled full QE only with gfortran and ifort).<br>
<br>
While testing the examples provided with PP I stumbled upon a problem in my patch (namely, when pp.x is used with an empty INPUTPP section). The attached patch resolves the problem and also includes some fixes to the documentation (which I have finally been able to compile).<br>
<br>
While running the examples, I noticed two further things that are not related to PP:<br>
<br>
* example05:  run_example (line 190) wants to copy a file feo_af.occup that doesn't seem to be produced by pw.x. Perhaps this output was removed at some point?<br>
* example06:  run_example (line 177) stops with<br>
  running the band calculation for Fe with PAW and spin-orbit...[cli_0]: aborting job:<br>
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<br>
<br>
Apart from this, the examples run through fine.<br>
<br>
<br>
Best,<br>
Leopold<br>
<br>
<br><br>
<br>
<br>
> On 11 Nov 2016, at 13:10 , Nicola Marzari <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:<br>
><br>
><br>
><br>
> Dear All,<br>
><br>
><br>
> slightly tangential to this discussion, but just wanted to remind/alert everyone<br>
> that Giovanni Pizzi (and especially his Japanese colleagues) have made a major<br>
> effort in devising standard paths in k-space for any structure, in ways that are<br>
> fully compliant with current crystallography standards (this was not the case<br>
> up to now).<br>
><br>
> Great if we all could adhere to this framework - see page below for a link to<br>
> the paper, an online tool that provides all points, labels and paths given a<br>
> pwscf input, and the link to the github repository if you want to use directly<br>
> the python source:<br>
><br>
> <a href="http://www.materialscloud.org/tools/seekpath/input_structure/" rel="noreferrer" target="_blank">http://www.materialscloud.org/<wbr>tools/seekpath/input_<wbr>structure/</a><br>
><br>
><br>
>                               nicola<br>
><br>
><br>
><br>
> On 11/11/2016 12:41, Pignedoli, Carlo wrote:<br>
>> Dear all,<br>
>> I understand it could be more complex, but since<br>
>> you are going to consider this useful modification<br>
>> I would spend an effort in trying to allow flexibility in the selection of k-points:<br>
>><br>
>> normally I will plot k-points after a bandstructure calculation<br>
>> performed with several k-points (could easily be 100 in one directions)<br>
>><br>
>> I iwl probably need to plot soem selected bands at selected k-points<br>
>><br>
>> now a "stride option" like plot every 3 kpoint would not be enough<br>
>> because most likely it would exclude some high simmetry points...<br>
>><br>
>> thus a format that allows to specify<br>
>><br>
>> 1 5 10 15 20 25 34<br>
>> in a calculation where e.g. 1 is G and 34 is X<br>
>> would be extremely useful.<br>
>><br>
>> Kind regards<br>
>><br>
>> Carlo<br>
>><br>
>><br>
>><br>
>><br>
>> Dr. Carlo Antonio Pignedoli<br>
>> Empa<br>
>> Swiss Federal Laboratories for Materials Science and Technology<br>
>> Überlandstrasse 129<br>
>> 8600 Dübendorf<br>
>> Switzerland<br>
>> Tel +41 58 765 4206<br>
>> Fax +41 58 765 4031<br>
>> e-mail <a href="mailto:carlo.pignedoli@empa.ch">carlo.pignedoli@empa.ch</a><br>
>> <a href="http://www.surfaces.ch" rel="noreferrer" target="_blank">www.surfaces.ch</a><br>
>><br>
>><br>
>>> -----Original Message-----<br>
>>> From: Leopold Talirz [mailto:<a href="mailto:leopold.talirz@gmail.com">leopold.talirz@gmail.<wbr>com</a>]<br>
>>> Sent: Thursday, November 10, 2016 12:55<br>
>>> To: <a href="mailto:tone.kokalj@ijs.si">tone.kokalj@ijs.si</a><br>
>>> Cc: General discussion list for Quantum ESPRESSO developers; Yakutovich,<br>
>>> Aliaksandr; Pignedoli, Carlo<br>
>>> Subject: Re: [Q-e-developers] Plotting many KS orbitals in one pp.x run (and<br>
>>> more)<br>
>>><br>
>>> Dear Tone,<br>
>>><br>
>>> I agree that Guido's suggestion is the best compromise.<br>
>>><br>
>>> Regarding your example, I think one should keep requiring the user to use<br>
>>> the kpoint keyword as well (otherwise one might end up by lots of orbitals by<br>
>>> accident).<br>
>>> But I guess this was not the point you were making. Regarding input format,<br>
>>> the same code will allow both<br>
>>><br>
>>> kband(1) = 21<br>
>>> kband(2) = 38<br>
>>> spin_component(1) = 0<br>
>>> spin_component(2) = 3<br>
>>><br>
>>> and<br>
>>><br>
>>> kband = 21 38<br>
>>> spin_component = 0 3<br>
>>><br>
>>> and<br>
>>><br>
>>> kband = 21<br>
>>> spin_component = 0<br>
>>><br>
>>> (at least with gfortran and intel, I haven't tested other compilers so far).<br>
>>><br>
>>> Best wishes,<br>
>>> Leopold<br>
>>><br>
>>><br>
>>>> On 10 Nov 2016, at 11:38 , Tone Kokalj <<a href="mailto:tone.kokalj@ijs.si">tone.kokalj@ijs.si</a>> wrote:<br>
>>>><br>
>>>> On Thu, 2016-11-10 at 11:20 +0000, Leopold Talirz wrote:<br>
>>>>> Dear Tone,<br>
>>>>><br>
>>>>>> In addition to lsign, I suggest also to make the spin_component<br>
>>>>>> constant: if more components are needed it can be done with few pp.x<br>
>>>>>> calculations.<br>
>>>>><br>
>>>>> From just the input file perspective I don't see a particular reason<br>
>>>>> to handle spin_component differently than kpoint and kband. Or am I<br>
>>>>> missing something?<br>
>>>><br>
>>>> Indeed you are right; if somebody is interested in more components it<br>
>>>> would make sense to make the spin_component follow the Guido<br>
>>>> suggestion. To illustrate: I found this useful (Guido scheme)<br>
>>>><br>
>>>> kband(1) = 21<br>
>>>> kband(2) = 38<br>
>>>> spin_component(1) = 0<br>
>>>> spin_component(2) = 3<br>
>>>><br>
>>>> according to Guide scheme this should make the total and x,y,z<br>
>>>> components for bands from 21 to 38. But I found this too fine grained:<br>
>>>><br>
>>>> kband = 21 22 23 24<br>
>>>> spin_component = 1 1 2 1<br>
>>>><br>
>>>><br>
>>>> Best regards, Tone<br>
>>>> --<br>
>>>> J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia<br>
>>>> (tel: +386-1-477-3523 // fax: +386-1-251-9385)<br>
>><br>
>><br>
>> ______________________________<wbr>_________________<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/<wbr>listinfo/q-e-developers</a><br>
>><br>
><br>
> --<br>
> ------------------------------<wbr>------------------------------<wbr>----------<br>
> Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
> Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br>
> <a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/<wbr>Contact</a> <a href="http://nccr-marvel.ch/en/project" rel="noreferrer" target="_blank">http://nccr-marvel.ch/en/<wbr>project</a><br>
<br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, <span>Dip. Scienze Matematiche Informatiche e Fisiche</span>,<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>