Laurence Marks
2018-10-01 18:03:19 UTC
Line 38 of ScaleDiag should be:
call dgemm('N', 'N', MAXMIX, MEMORY, MEMORY, 1.D0, Y, MAXMIX, W,
MEMORY, 0.0D0, W2, MAXMIX)
It was "0.0, W2" without the D0, which probably does not matter but is not
right.
Thanks to Pavel OndraÄka for finding this one.
call dgemm('N', 'N', MAXMIX, MEMORY, MEMORY, 1.D0, Y, MAXMIX, W,
MEMORY, 0.0D0, W2, MAXMIX)
It was "0.0, W2" without the D0, which probably does not matter but is not
right.
Thanks to Pavel OndraÄka for finding this one.
--
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu ; Corrosion in 4D: MURI4D.numis.northwestern.edu
Partner of the CFW 100% program for gender equity, www.cfw.org/100-percent
Co-Editor, Acta Cryst A
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu ; Corrosion in 4D: MURI4D.numis.northwestern.edu
Partner of the CFW 100% program for gender equity, www.cfw.org/100-percent
Co-Editor, Acta Cryst A