<div dir="ltr">It seems like that you have exceeded stack space with your local array. Did you run under valgrind?<div><br></div><div>  Matt</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Feb 7, 2016 at 3:51 AM, Adrian Jackson <span dir="ltr"><<a href="mailto:adrianj@epcc.ed.ac.uk" target="_blank">adrianj@epcc.ed.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
We're using PETSc in a fortran code and get an error if we compile the<br>
code using the -fstack-arrays flag when using the GNU fortran compilers.<br>
 This is the error we are getting:<br>
<br>
[22]PETSC ERROR: --------------------- Error Message<br>
--------------------------------------------------------------<br>
[22]PETSC ERROR: Arguments are incompatible<br>
[22]PETSC ERROR: Zero diagonal on row 23003<br>
....<br>
[22]PETSC ERROR: #1 MatInvertDiagonal_SeqAIJ() line 1841 in<br>
/b/cray-petsc/.cray-build/GNU/49/sandybridge/cray-petsc-base-dynamic/petsc-3.5.2/src/mat/impls/aij/seq/aij.c<br>
[22]PETSC ERROR: #2 MatSOR_SeqAIJ() line 1873 in<br>
/b/cray-petsc/.cray-build/GNU/49/sandybridge/cray-petsc-base-dynamic/petsc-3.5.2/src/mat/impls/aij/seq/aij.c<br>
....<br>
[22]PETSC ERROR: #14 KSPSolve() line 459 in<br>
/b/cray-petsc/.cray-build/GNU/49/sandybridge/cray-petsc-base-dynamic/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c<br>
 WARNING: Failed to converge.<br>
 PETSc did not converge for matrix solve of: DeltaU<br>
 Reason for non-convergence is undefined:            0<br>
 Number of iterations:            0<br>
 Sending signal to dump and finish<br>
<br>
<br>
I've checked the right hand side vector passed to KSPSolve in both cases<br>
and it's the same.  I've not fully verified the ksp context is exactly<br>
the same but it doesn't seem to be corrupted when using -fstack-arrays.<br>
 Is this a problem people have encountered before?<br>
<br>
thanks<br>
<br>
adrianj<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
The University of Edinburgh is a charitable body, registered in<br>
Scotland, with registration number SC005336.<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div>