erdc / proteus-mprans

Module for solving 3D, non-hydrostatic, turbulent, incompressible air/water flow (multi-phase RANS)
7 stars 6 forks source link

suboff2D example hangs in initialization if discretization is too coarse #27

Open mfarthin opened 10 years ago

mfarthin commented 10 years ago

When I run with the refinement flag set to 4, the code hangs in the calculateElementQuadrature function in RANS2P.py during the
femSpace.elementMaps.getJacobianValues call. Running in the debugger, the q['J'] value looks fine before the call to

femSpace.elementMaps.getValues(self.elementQuadraturePoints,self.q['x'])

After this call, pdb hangs if I try to print out anything on q['J']