Solver() create problem in 3d mesh

Hope so now you people get a good idea and able to help me .

code:

``````from dolfin import *

# Optimization options:
parameters["form_compiler"]["cpp_optimize"]       = True
parameters['form_compiler']['cpp_optimize_flags'] = '-O3 -ffast-math -march=native'
parameters["form_compiler"]["optimize"]           = True
parameters["refinement_algorithm"] = "plaza_with_parent_facets"
parameters["ghost_mode"]           = "shared_facet"

mesh = Mesh()
.
.
.
.
solve(F == 0, up, bc, solver_parameters={"newton_solver":{"linear_solver":'mumps'},"newton_solver":{"relative_tolerance":1e-7}})
u, p, lamda = up.split()
``````

system space:

``````              total        used        free      shared  buff/cache   available
Mem:            62Gi       3.5Gi        49Gi        77Mi       9.5Gi        58Gi
Swap:          6.5Gi          0B       6.5Gi
``````

error;

``````<meshio mesh object>
Number of points: 26272
Number of cells:
triangle: 4740
triangle: 3250
triangle: 556
triangle: 3222
triangle: 566
triangle: 704
triangle: 680
triangle: 702
triangle: 690
triangle: 4732
tetra: 132339
Cell sets: gmsh:bounding_entities
Point data: gmsh:dim_tags
Cell data: gmsh:physical, gmsh:geometrical

Calling FFC just-in-time (JIT) compiler, this may take some time.
WARNING: The number of integration points for each cell will be: 125
Consider using the option 'quadrature_degree' to reduce the number of points
WARNING: The number of integration points for each cell will be: 125
Consider using the option 'quadrature_degree' to reduce the number of points
Calling FFC just-in-time (JIT) compiler, this may take some time.
WARNING: The number of integration points for each cell will be: 125
Consider using the option 'quadrature_degree' to reduce the number of points
WARNING: The number of integration points for each cell will be: 125
Consider using the option 'quadrature_degree' to reduce the number of points

UMFPACK V5.7.9 (Oct 20, 2019): ERROR: out of memory

Traceback (most recent call last):
File "/home/ayush/flow.py", line 116, in <module>
solve(F == 0, up, bc, solver_parameters={"newton_solver":{"linear_solver":'mumps'},"newton_solver":{"relative_tolerance":1e-7}})
File "/usr/lib/petsc/lib/python3/dist-packages/dolfin/fem/solving.py", line 233, in solve
_solve_varproblem(*args, **kwargs)
File "/usr/lib/petsc/lib/python3/dist-packages/dolfin/fem/solving.py", line 314, in _solve_varproblem
solver.solve()
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
***
***
*** 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 'KSPSolve'.
*** Reason:  PETSc error code is: 76 (Error in external library).
*** Where:   This error was encountered inside ./dolfin/la/PETScKrylovSolver.cpp.
*** Process: 0
***
*** DOLFIN version: 2019.2.0.dev0
*** Git changeset:  ubuntu
*** -------------------------------------------------------------------------
``````

Your code is not sufficient to reproduce the error.