------- Comment #2 from kargl at gcc dot gnu dot org 2010-09-09 22:25 ------- There is no way to fix this problem unless you would like +inf along the diagonal.
gfortran will constant fold 1./alpha if alpha has the parameter attribute. After all, this attribute tells the compiler that alpha is a named constant, ie., it will not change. PS: Fortran bugs are never given a severity high than normal unless the bug breaks bootstrapping the compiler. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=45624