<div dir="ltr">Many tests (also with smearing) are compiled with different compilers (including gfortran 10.2) and run during the night (CET time): <a href="http://bionano2.dimi.uniud.it:8010/#/console">http://bionano2.dimi.uniud.it:8010/#/console</a>. Unless you have more evidence, what you report seems to indicate that the problem is in the compiler, not in the code. <div><br></div><div>Paolo</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 23, 2021 at 8:11 AM Bainglass, Edan <<a href="mailto:Edan.Bainglass@unt.edu">Edan.Bainglass@unt.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-US" style="overflow-wrap: break-word;">
<div class="gmail-m_-6212740695167873857WordSection1">
<p class="MsoNormal"><b>Hello,<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>For gfortran(10.3.0)-compiled QE 6.8, I’m getting the following error for jobs that include smearing:<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>Program received signal SIGSEGV: Segmentation fault - invalid memory reference.<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Backtrace for this error:<u></u><u></u></p>
<p class="MsoNormal">#0 0x146941220d21 in ???<u></u><u></u></p>
<p class="MsoNormal">#1 0x14694121fef5 in ???<u></u><u></u></p>
<p class="MsoNormal">#2 0x1469414dd3bf in ???<u></u><u></u></p>
<p class="MsoNormal">#3 0x7ffd48ddd9cc in ???<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>Job terminates right after this line:<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> iteration # 1 ecut= 35.00 Ry beta= 0.70<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>This does not occur with the latest Intel compiler. It also does not occur with either compiler on QE 6.7.<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>For completeness, here are a few jobs that fail:<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>1.<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">&CONTROL<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> calculation = 'scf'<u></u><u></u></p>
<p class="MsoNormal"> restart_mode = 'from_scratch'<u></u><u></u></p>
<p class="MsoNormal"> tprnfor = .TRUE.<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&SYSTEM<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> ibrav = 8<u></u><u></u></p>
<p class="MsoNormal"> celldm(1) = 10.6881<u></u><u></u></p>
<p class="MsoNormal"> celldm(2) = 0.866025<u></u><u></u></p>
<p class="MsoNormal"> celldm(3) = 3.95422<u></u><u></u></p>
<p class="MsoNormal"> nat = 10<u></u><u></u></p>
<p class="MsoNormal"> ntyp = 3<u></u><u></u></p>
<p class="MsoNormal"> ecutwfc = 35<u></u><u></u></p>
<p class="MsoNormal"> ecutrho = 300<u></u><u></u></p>
<p class="MsoNormal"> occupations = 'smearing'<u></u><u></u></p>
<p class="MsoNormal"> degauss = 0.03<u></u><u></u></p>
<p class="MsoNormal"> smearing = 'mv'<u></u><u></u></p>
<p class="MsoNormal"> nbnd = 80<u></u><u></u></p>
<p class="MsoNormal"> tot_charge = 0<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&ELECTRONS<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> conv_thr = 5.D-6<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">K_POINTS gamma<u></u><u></u></p>
<p class="MsoNormal">ATOMIC_SPECIES<u></u><u></u></p>
<p class="MsoNormal">C 1 C.pbe-rrkjus.UPF<u></u><u></u></p>
<p class="MsoNormal">O 1 O.pbe-rrkjus.UPF<u></u><u></u></p>
<p class="MsoNormal">Pt 1 Pt.pbe-n-rrkjus_psl.0.1.UPF<u></u><u></u></p>
<p class="MsoNormal">ATOMIC_POSITIONS (bohr)<u></u><u></u></p>
<p class="MsoNormal">C 5.335084148 4.646723426 12.901029877<u></u><u></u></p>
<p class="MsoNormal">O 5.335009643 4.619623254 15.079854269<u></u><u></u></p>
<p class="MsoNormal">Pt 8.061327071 0.098057998 8.992142901<u></u><u></u></p>
<p class="MsoNormal">Pt 2.608989366 0.098058283 8.992140585<u></u><u></u></p>
<p class="MsoNormal">Pt 0.000036609 4.720846294 8.968756935<u></u><u></u></p>
<p class="MsoNormal">Pt 5.335159557 4.721612729 9.380196435<u></u><u></u></p>
<p class="MsoNormal">Pt 0.000041121 7.802951963 4.604626508<u></u><u></u></p>
<p class="MsoNormal">Pt 5.335161233 7.697749113 4.753489408<u></u><u></u></p>
<p class="MsoNormal">Pt 2.697860636 3.152173889 4.688412329<u></u><u></u></p>
<p class="MsoNormal">Pt 7.972463687 3.152174491 4.688415209<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>2.<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">&CONTROL<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> calculation = 'scf'<u></u><u></u></p>
<p class="MsoNormal"> restart_mode = 'from_scratch'<u></u><u></u></p>
<p class="MsoNormal"> tprnfor = .TRUE.<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&SYSTEM<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> ecutrho = 300<u></u><u></u></p>
<p class="MsoNormal"> ecutwfc = 30<u></u><u></u></p>
<p class="MsoNormal"> ibrav = 1<u></u><u></u></p>
<p class="MsoNormal"> celldm( 1 ) = 20<u></u><u></u></p>
<p class="MsoNormal"> nat = 2<u></u><u></u></p>
<p class="MsoNormal"> ntyp = 2<u></u><u></u></p>
<p class="MsoNormal"> nbnd = 6<u></u><u></u></p>
<p class="MsoNormal"> tot_charge = 0<u></u><u></u></p>
<p class="MsoNormal"> nspin = 2<u></u><u></u></p>
<p class="MsoNormal"> starting_magnetization(1) = 1<u></u><u></u></p>
<p class="MsoNormal"> occupations = 'smearing'<u></u><u></u></p>
<p class="MsoNormal"> smearing = 'gaussian'<u></u><u></u></p>
<p class="MsoNormal"> degauss = 0.01<u></u><u></u></p>
<p class="MsoNormal"> nosym = .true.<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&ELECTRONS<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal"> conv_thr = 5.D-3<u></u><u></u></p>
<p class="MsoNormal"> !<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">K_POINTS (gamma)<u></u><u></u></p>
<p class="MsoNormal">ATOMIC_SPECIES <u></u><u></u></p>
<p class="MsoNormal">C 12 C.pbe-rrkjus.UPF <u></u><u></u></p>
<p class="MsoNormal">H 1 H.pbe-rrkjus.UPF <u></u><u></u></p>
<p class="MsoNormal">ATOMIC_POSITIONS (angstrom)<u></u><u></u></p>
<p class="MsoNormal">C 0.00000000 0.00000000 0.16007400 <u></u><u></u></p>
<p class="MsoNormal">H 0.00000000 0.00000000 -0.96044600 <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b>3.<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">&CONTROL<u></u><u></u></p>
<p class="MsoNormal"> calculation = 'scf',<u></u><u></u></p>
<p class="MsoNormal"> tprnfor = .true.,<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&SYSTEM<u></u><u></u></p>
<p class="MsoNormal"> nspin = 1,<u></u><u></u></p>
<p class="MsoNormal"> ecutwfc = 35,<u></u><u></u></p>
<p class="MsoNormal"> nat = 2,<u></u><u></u></p>
<p class="MsoNormal"> occupations = 'smearing',<u></u><u></u></p>
<p class="MsoNormal"> ibrav = 0,<u></u><u></u></p>
<p class="MsoNormal"> degauss = 0.01,<u></u><u></u></p>
<p class="MsoNormal"> smearing = 'marzari-vanderbilt',<u></u><u></u></p>
<p class="MsoNormal"> tot_charge = 0,<u></u><u></u></p>
<p class="MsoNormal"> ntyp = 1,<u></u><u></u></p>
<p class="MsoNormal"> ecutrho = 350,<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&ELECTRONS<u></u><u></u></p>
<p class="MsoNormal"> mixing_beta = 0.3<u></u><u></u></p>
<p class="MsoNormal"> conv_thr = 1e-06,<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">&IONS<u></u><u></u></p>
<p class="MsoNormal">/<u></u><u></u></p>
<p class="MsoNormal">ATOMIC_SPECIES<u></u><u></u></p>
<p class="MsoNormal">Ag 107.868000 Ag.pbe-n-rrkjus_psl.1.0.0.UPF<u></u><u></u></p>
<p class="MsoNormal">ATOMIC_POSITIONS {bohr}<u></u><u></u></p>
<p class="MsoNormal">Ag 2.771763 2.771763 18.897260 <u></u><u></u></p>
<p class="MsoNormal"> Ag 0.000000 0.000000 22.740340 <u></u><u></u></p>
<p class="MsoNormal">K_POINTS {automatic}<u></u><u></u></p>
<p class="MsoNormal"> 1 1 1 0 0 0<u></u><u></u></p>
<p class="MsoNormal">CELL_PARAMETERS {bohr}<u></u><u></u></p>
<p class="MsoNormal">5.543526 0.000000 0.000000<u></u><u></u></p>
<p class="MsoNormal">0.000000 5.543526 0.000000<u></u><u></u></p>
<p class="MsoNormal">0.000000 0.000000 65.08548<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b><u></u> <u></u></b></p>
<p class="MsoNormal"><b>Thanks in advance for any info,<u></u><u></u></b></p>
<p class="MsoNormal"><b><u></u> <u></u></b></p>
<p class="MsoNormal"><b>Dr. Edan Bainglass<u></u><u></u></b></p>
<p class="MsoNormal"><b>Postdoctoral Research Assistant<u></u><u></u></b></p>
<p class="MsoNormal"><b>Environ Developer<u></u><u></u></b></p>
<p class="MsoNormal"><b>MATERIALab Group<u></u><u></u></b></p>
<p class="MsoNormal"><b>Physics Department<u></u><u></u></b></p>
<p class="MsoNormal"><b>University of North Texas<u></u><u></u></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" 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 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>