<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>I added two lines in<i> atom_write.f </i>in <i>SRC_dstart</i>
and commented one line (changes in bold):</p>
<p><br>
</p>
<p><br>
if (ia.eq.1) then<br>
! if(myid .eq. 0)then<br>
write(51,1970) iscf
<br>
write(51,78) title,lattic,aa,bb,cc,jrj(1)
<br>
<b># write(51,77)<br>
write(51,2022) lmmax(ia+1)</b><br>
! end if<br>
end if<br>
! if(myid .eq. 0)then<br>
write(51,1990) ia <br>
write(51,2771) lmmax(ia) <br>
<br>
do lm1=1,lmmax(ia) <br>
<br>
write(51,2011) lm(1,lm1,ia),lm(2,lm1,ia) <br>
write(51,2021) ( rholm(j,lm1), j=1,jrj(ia) ) <br>
write(51,2031)
<br>
!enddo<br>
enddo<br>
write(51,2033) <br>
! end if <br>
! <br>
77 format(1x,9(f9.7))
<br>
78 format(1x,a20,a4,3f10.6,i5)<br>
1970 format(3x,' TOTAL CHARGE DENSITY GENERATED ',
&<br>
'BY',i3,'. ITERATION (NORM: CLM=CLM*R*R)') <br>
1980 format(3x)
<br>
2771 format(3x,'NUMBER OF LM',I3//)<br>
<b> 2022 FORMAT(1x,'LM-max:'1x,I3)</b></p>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Maybe the experts can quickly check, if
I did all correctly .. Thanks in advance!</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Am 10.07.2025 um 12:44 schrieb Fecher,
Gerhard:<br>
</div>
<blockquote type="cite"
cite="mid:89b8cc8d7c5b49549f7b5019f593df11@uni-mainz.de">
<pre wrap="" class="moz-quote-pre">Fine that it works also for other cases
you have to check in atom_write.f of dstart
Ciao
Gerhard
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:
"I think the problem, to be quite honest with you,
is that you have never actually known what the question is."
====================================
Dr. Gerhard H. Fecher
Institut of Physics
Johannes Gutenberg - University
55099 Mainz
________________________________________
Von: Wien [<a class="moz-txt-link-abbreviated" href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] im Auftrag von Michael Fechtelkord [<a class="moz-txt-link-abbreviated" href="mailto:Michael.Fechtelkord@ruhr-uni-bochum.de">Michael.Fechtelkord@ruhr-uni-bochum.de</a>]
Gesendet: Donnerstag, 10. Juli 2025 12:29
An: A Mailing list for WIEN2k users
Betreff: Re: [Wien] LAPW0 -- again
Hello Gerhard,
I just tried your solution on my MgF2 case. lapw0 runs here without any
problems when I insert LM-MAX into the initial case.clmsum
TOTAL CHARGE DENSITY GENERATED BY 0. ITERATION (NORM:
CLM=CLM*R*R)
MgF2 P 8.736207 8.736207 5.767446 781
ATOMNUMBER = 1
NUMBER OF LM 10
CLM(R) FOR L 0 M= 0
...
adding
LM-max: 16
in the third line makes lapw0 run, even with omp_num_threads larger than 1
TOTAL CHARGE DENSITY GENERATED BY 0. ITERATION (NORM:
CLM=CLM*R*R)
MgF2 P 8.736207 8.736207 5.767446 781
LM-max: 16
ATOMNUMBER = 1
NUMBER OF LM 10
CLM(R) FOR L 0 M= 0
....
It should be easy to add that line by dstart in the initial *.clmsum.
Best regards,
Michael
Am 10.07.2025 um 11:54 schrieb Fecher, Gerhard:
</pre>
<blockquote type="cite">
<pre wrap="" class="moz-quote-pre">Dear all,
I am curious how lapw0 should know in which cycle it is,
there must be a difference between the first and all other cycles,
that means something is not or not completely initialised in the 1st cycle.
There are not too many files to be inspected that are used by lapw0 as input,
case.clsum contains in the first line the number of the iteration, which is "0." in the first cycle and "1." after the first cycle, and so on .....
However, there is another difference:
after dstart, the third line of case.clmsum is empty
TOTAL CHARGE DENSITY GENERATED BY 0. ITERATION
Li a_exp B 6.632941 6.632941 6.632941 781
ATOMNUMBER = 1
NUMBER OF LM 5
...
whereas after the first cycle it is, in my simple test case
TOTAL CHARGE DENSITY GENERATED BY 1. ITERATION
Li a_exp B 6.632941 6.632941 6.632941 781
LM-max: 5
ATOMNUMBER = 1
NUMBER OF LM 5
...
Now, when I insert ' LM-max: 5' into the empty line of a fresh case.clmsum after dstart, then lapw0 runs without problem altready in the first cycle.
The line with LM-max is also missing after dstart in the case.clmup/dn files of spin polarised cakculations, but there it seems that it doesn't matter .
Solution, if the missing line is the cause for the segmentation faults (maybe through wrong array size by allocation):
1) dstart writes the LM-max line with the correct value into case.clmsum
or
2) lapw0 needs to be inspected in more detail:
lmmax ix used for various allocations, maybe that is the reason why something goes wrong if the empty input line is interpreted wrong
check in lapw0 line 336 ff (and the format lines)
! READ SCFDATA UNTIL BOTTOM AND INCREASE CYCLENUMBER BY ONE
if(myid.eq.0) then
READ(8,2044) ISCF
read(8,2043) lmmaxx
endif
#ifdef Parallel
call MPI_Bcast(iscf,1,MPI_INTEGER,0,MPI_COMM_WORLD,ierr)
call MPI_Bcast(lmmaxx,1,MPI_INTEGER,0,MPI_COMM_WORLD,ierr)
#endif
ISCF=ISCF+1
if(lmmaxx.eq.0) lmmaxx=ncom
IF(ISCF.GT.999) ISCF=1
...
(note: ncom=261 from params.inc)
It might also be that trying to read an integer from the empty third line causes also later some other problems with the input, which was not a problem with ifort or gfortran but appears in ifx
(It is not known that the old comments about the OMP problems in lapw0 whether they appeared only in the 1st or also in other cycles,
however from the comments in lapw0.F and the changes file in SRC_lapw0 from 2021, it seems to be not just a problem of ifx)
Ciao
Gerhard
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:
"I think the problem, to be quite honest with you,
is that you have never actually known what the question is."
====================================
Dr. Gerhard H. Fecher
Institut of Physics
Johannes Gutenberg - University
55099 Mainz
_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
SEARCH the MAILING-LIST at: <a class="moz-txt-link-freetext" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>
</pre>
</blockquote>
<pre wrap="" class="moz-quote-pre">
--
Dr. Michael Fechtelkord
Institut für Geologie, Mineralogie und Geophysik
Ruhr-Universität Bochum
Universitätsstr. 150
D-44780 Bochum
Phone: +49 (234) 32-24380
Fax: +49 (234) 32-04380
Email: <a class="moz-txt-link-abbreviated" href="mailto:Michael.Fechtelkord@ruhr-uni-bochum.de">Michael.Fechtelkord@ruhr-uni-bochum.de</a>
Web Page: <a class="moz-txt-link-freetext" href="https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/">https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/</a>
_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
SEARCH the MAILING-LIST at: <a class="moz-txt-link-freetext" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>
_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
SEARCH the MAILING-LIST at: <a class="moz-txt-link-freetext" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
Dr. Michael Fechtelkord
Institut für Geologie, Mineralogie und Geophysik
Ruhr-Universität Bochum
Universitätsstr. 150
D-44780 Bochum
Phone: +49 (234) 32-24380
Fax: +49 (234) 32-04380
Email: <a class="moz-txt-link-abbreviated" href="mailto:Michael.Fechtelkord@ruhr-uni-bochum.de">Michael.Fechtelkord@ruhr-uni-bochum.de</a>
Web Page: <a class="moz-txt-link-freetext" href="https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/">https://www.ruhr-uni-bochum.de/kristallographie/kc/mitarbeiter/fechtelkord/</a>
</pre>
</body>
</html>