<div dir="ltr"><div>In my opinion, the postprocessing code "pp.x" is a mess beyond control that badly needs to be split, re-organized and cleaned up. This being said: your patch adds a negligible amount of entropy, so I do not personally see major reasons not to accept it. I would however avoid changes to the symmetry section: it is sufficient to call "sym_rho_init" just once, outside the loop over k-points and bands. I would also avoid adding new "plot_num" values since there are no new quantities to be plotted, just new options to plot them<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 7, 2016 at 6:05 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"><div style="word-wrap:break-word"><div>Dear QE developers,</div><div><br></div><div>this is a feature request, which I am sending to the mailing list in lieu of access to the QE bug tracker on <a href="http://qe-forge.org" target="_blank">qe-forge.org</a></div><div><br></div><div><br></div><div>I am sure many of you are familiar with the AiiDA project (<a href="http://www.aiida.net/" target="_blank">http://www.aiida.net/</a>), which provides not only a very convenient python infrastructure for running QE calculations, but also a database for a consistent storage of results.</div><div> </div><div>We have recently come across a "conflict of philosophy" between QE (more specifically, pp.x) and AiiDA that causes considerable computational/data overhead.</div><div><br></div><div>Basically, we need to plot lots of KS orbitals for different nanostructures. </div><div>pp.x is designed to produce 1 cube file per execution and AiiDA is designed to produce 1 working directory per calculation. </div><div>I.e. when computing 100 KS orbitals from a band structure calculation, using AiiDA the information contained in the band structure directory is replicated 100 times. </div><div> </div><div>Our solution so far has been to modify the source of pp.x in such a way that one run can produce N KS orbitals. This has the added benefit of saving a lot of "startup time" as compared to running pp.x N times (in a test run generating 21 KS orbitals, the speedup was ~3x but this does vary, of course).</div><div>We have already contacted AiiDA developers Nicolas Mounet and Giovanni Pizzi and they are in favor of this solution as well.</div><div><br></div><div>I attach a patch against QE 6.0 (git SHA 93ced686f67fe1f3250824358015e2<wbr>ffde364960) that adds plot_num 22 and 23 in pp.x.</div><div>plot_num 23 allows to plot N KS orbitals and plot_num 22 allows to plot the local density of states on a grid of energies between emin and emax (useful for scanning tunneling spectroscopy simulations). The patch starts with a commit message that contains sample input files.</div><div><br></div><div>I think the main question is, whether you would consider deviating from the "1 run = 1 output file" policy. </div><div>If yes, then one can think of the details regarding the implementation (the patch, including the input format for the KS orbitals, can certainly be improved).</div><div><br></div><div>Best wishes,</div><div>Leopold Talirz, Carlo Pignedoli and Sasha Yakutovich</div><div><br></div><div></div></div><br><div style="word-wrap:break-word"><div></div></div><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></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, 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>