This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

mshr causing PETSc error after upgrade to 2016.2

0 votes

Hello,

Up to now, I was using FEniCS through the automated build script, and my last build was 1.7.0dev. I just updated to work with FEniCS under docker, so I followed instructions to work with the stable image 2016.2, and everything seemed to work fine (checked some demos and they worked fine). However, I found that in the container, when I try to run one of my old scripts using mshr, I get the following error (minimal code is reproduced at the end). I am stuck cause I am unable to find where this comes from, so I would be grateful if anyone could help me.

Thank you very much for your help.

[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run 
[0]PETSC ERROR: to get more information on the crash.
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=59
:
system msg for write_line failure : Bad file descriptor

Minimal test code:

from dolfin import *
import mshr 
if __name__ == "__main__":
    Center = Point(0.0, 0.0); Radius = 6.35E-3 
    circle1 = mshr.Circle(Center,Radius)
    circle2 = mshr.Circle(Center,2.5*Radius)
    circle2.set_subdomain(1,circle1)
    circle2.set_subdomain(2,circle2-circle1)
    mesh=mshr.generate_mesh(circle2,20)
asked Mar 12, 2017 by Juan Córcoles FEniCS Novice (250 points)

I can confirm that this works with quay.io/fenicsproject/stable:2016.1.0 but not with quay.io/fenicsproject/stable:current or quay.io/fenicsproject/dev. It works when increasing the Radius to 0.024. You might consider reporting an issue for mshr.

...