Hi,
For large problems (around 100 million degrees of freedom in 3D) I get the error
File "Laplacian.py", line 40, in <module>
assemble_system(a, l, bc, A_tensor = A, b_tensor = b)
File "/ubc/cs/research/scl/people/mwathen/.hastdist/bld/profile/4nxzwfz7gakq/lib/python2.7/site-packages/dolfin/fem/assembling.py", line 261, in assemble_system
assembler.assemble(A_tensor, b_tensor)
RuntimeError:
*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
*** fenics@fenicsproject.org
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error: Unable to successfully call PETSc function 'MatMPIAIJSetPreallocation'.
*** Reason: PETSc error code is: 1.
*** Where: This error was encountered inside /ubc/cs/research/scl/people/mwathen/.hastdist/tmp/dolfin-2ckaxbyj3rot/dolfin/la/PETScMatrix.cpp.
*** Process: 0
***
*** DOLFIN version: 1.6.0
*** Git changeset:
*** -------
Is there a way of getting around this?
