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"
parameters["form_compiler"]["quadrature_degree"] = 8

# reading the mesh made using gmsh :
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
***
***     fenics-support@googlegroups.com
***
*** 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.
Please see:

https://fenicsproject.org/qa/4177/reason-petsc-error-code-is-76/

and make a minimal reproducible example.