nlesolver-fortran
nlesolver-fortran copied to clipboard
Segfault for large problems when using Broyden
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?