Hi everyone,
I'm applying ISSM on some Patagonian glaciers, but I'm running into an error when I use the solve function.
When I run md = solve(md,'Stressbalance');, I get this error message:
solver residue: norm(KU-F)/norm(F)=-nan
??? Error in ==> ./modules/Solverx/Solverx.cpp:63
checkconvergence error message: Solver residue is NaN
??? Error in ==> ./modules/Solverx/Solverx.cpp:63
checkconvergence error message: Solver residue is NaN
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpiexec detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[24149,1],0]
Exit code: 1
--------------------------------------------------------------------------
loading results from cluster
[Warning:
============================================================
Binary file HPEGlaciers.outbin not found!
This typically happens when the run crashed.
Please check for error messages above or in the outlog
============================================================
This error, however, only occurs when I run my code on the HPC of my university, not when I run it on my Windows laptop. I'm using exactly the same input data, Matlab version (R2019b), ISSM version (4.20) and scripts on both the HPC and my Windows PC. md = solve(md,'Stressbalance') and md = solve(md, 'transient') work perfectly fine on my Windows laptop without an NaN value for the solver residue, but not on our HPC.
Do you know where 'solver residue: norm(KU-F)/norm(F)=-nan' might be relate to?
Best,
Matthias
