gfortran matrix multiplication result size bug

foo.f90

program foo 

real, parameter :: A(2, 3) = reshape((/ 0, 0, 0, 0, 0, 0 /), (/ 2, 3 /)) 
real, parameter :: B(2) = (/ 0, 0 /)
real, parameter :: C(3) = matmul(B, A)
!                      1
!Error: Different shape for array assignment at (1) on dimension 1 (3 and 2)
!f951: internal compiler error: in gfc_simplify_reshape, at fortran/simplify.c:5104

! This works in some versions of gfortran:
real, parameter :: D(3) = reshape(matmul(B, A), (/ 3 /))
! Others give:
!f951: internal compiler error: in gfc_simplify_reshape, at fortran/simplify.c:5104

! This outputs 3 values:
print *, matmul(B, A)

end program foo