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
*** -------------------------------------------------------------------------