Segfault for large problems when using Broyden

[jacobwilliams]

This line

fjac = prev_fjac + &
    matmul((delf-matmul(prev_fjac,delx)),&
    transpose(delx))/delxmag2

will segfault on intel for a large problem (e.g., n=1406, m=1200).

What can we do about it?

[jacobwilliams]

see this discussion: https://fortran-lang.discourse.group/t/automatic-arrays-and-intrinsic-array-operations-to-use-or-not-to-use/4070