MFEM + PyMFEM (finite element method library)
Project description
MFEM + PyMFEM (FEM library)
This package provides MFEM and its Python wrapper (PyMFEM). MFEM is a high performance parallel finite element method (FEM) library (http://mfem.org/).
Installer downloads a couple of external libraries and build them. By default, "pip install mfem" downloads and builds the serial version of MFEM and PyMFEM. See more detail below for other configurations
Install
pip install mfem # binary install is available only on linux platforms (Py36-39)
pip install mfem --no-binary mfem # install serial MFEM + wrapper
The setup script accept various options. TO use it, please download the package and run the script manually. For example, this below download and build parallel version of MFEM library (linked with Metis and Hypre) and install under /mfem
$ pip3 download mfem
(expand tar.gz file and move to the downloaded directory)
$ python setup.py install --with-parallel # it download and build metis/hypre/mfem
Choosing compiler
$ python setup.py install --with-parallel --CC=icc --CXX=icpc --MPICC=mpiicc --MPICXX=mpiicpc
For other configurations, see docs/install.txt or help
$ python setup.py install --help
Usage
Here is an example to solve div(grad(f)) = 1 in a square and to plot the result with matplotlib (modified from ex1.cpp). Use the badge above to open this in Binder.
import mfem.ser as mfem
# create sample mesh for square shape
mesh = mfem.Mesh(10, 10, "TRIANGLE")
# create finite element function space
fec = mfem.H1_FECollection(1, mesh.Dimension()) # H1 order=1
fespace = mfem.FiniteElementSpace(mesh, fec)
#
ess_tdof_list = mfem.intArray()
ess_bdr = mfem.intArray([1]*mesh.bdr_attributes.Size())
fespace.GetEssentialTrueDofs(ess_bdr, ess_tdof_list)
# constant coefficient
one = mfem.ConstantCoefficient(1.0)
# define Bilinear and Linear operator
a = mfem.BilinearForm(fespace)
a.AddDomainIntegrator(mfem.DiffusionIntegrator(one))
a.Assemble()
b = mfem.LinearForm(fespace)
b.AddDomainIntegrator(mfem.DomainLFIntegrator(one))
b.Assemble()
# create gridfunction, which is where the solution vector is stored
x = mfem.GridFunction(fespace);
x.Assign(0.0)
# form linear equation (AX=B)
A = mfem.OperatorPtr()
B = mfem.Vector()
X = mfem.Vector()
a.FormLinearSystem(ess_tdof_list, x, b, A, X, B);
print("Size of linear system: " + str(A.Height()))
# solve it using PCG solver and store the solution to x
AA = mfem.OperatorHandle2SparseMatrix(A)
M = mfem.GSSmoother(AA)
mfem.PCG(AA, M, B, X, 1, 200, 1e-12, 0.0)
a.RecoverFEMSolution(X, b, x)
# extract vertices and solution as numpy array
verts = mesh.GetVertexArray()
sol = x.GetDataArray()
# plot solution using Matplotlib
import matplotlib.pyplot as plt
import matplotlib.tri as tri
triang = tri.Triangulation(verts[:,0], verts[:,1])
fig1, ax1 = plt.subplots()
ax1.set_aspect('equal')
tpc = ax1.tripcolor(triang, sol, shading='gouraud')
fig1.colorbar(tpc)
plt.show()
License
PyMFEM is licensed under BSD-3. Please refer the developers' web sites for the external libraries
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distributions
Hashes for mfem-4.3.0.1-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 9a6a38829cebec8b5f2010638171e3e7826cb128ffea5acaf86deb8530c671b9 |
|
MD5 | 44cac70182c31a9c57f6ad5afe6d0bc5 |
|
BLAKE2b-256 | b8518b878f7a34b566f7ab4f72c06a99a67bded6f17db736fbc15f3f38ce935d |
Hashes for mfem-4.3.0.1-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | dc1e24e8fa3b174f7e9010d118e59ffd07203554020955115ac250b0d5b6d26d |
|
MD5 | 9731b83c10358dea94c20cfe9bc04665 |
|
BLAKE2b-256 | 9d12e20ad3a73835c16fd2a55db5bf608c75f45ac5720ada81c289ab97f3883b |
Hashes for mfem-4.3.0.1-cp37-cp37m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 71d5d3654b4cf269c796c7003a865ec48de6cb7fdad48092dd14e2f4bd96d255 |
|
MD5 | 2a59348e379e252c67b37137226ad2ad |
|
BLAKE2b-256 | 5f8e87eb809a68bb619af9aa0bf6016a2642ebc9e312c72bd11d620c83202074 |
Hashes for mfem-4.3.0.1-cp36-cp36m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 87212af97ddc58b3b6562c1ba851fe3085dd60efeaa17d933a389e20c78059d0 |
|
MD5 | 6c13c88eb7561b166f30fc2ead6ee56c |
|
BLAKE2b-256 | 7fc08e6d8b62ecff6b366b3ca8f206a0a29b0f1cdfeccb636ef51a0a955a6bee |