esys.finley Package¶
Our most general domain representation. Imports submodules into its namespace
Classes¶
- class esys.finley.FinleyDomain¶
A concrete class representing a domain. For more details, please consult the C++ documentation.
- __init__((object)arg1, (FinleyDomain)arg2) None ¶
- MPIBarrier((FinleyDomain)arg1) None : ¶
Wait until all processes have reached this point
- addPDEToLumpedSystem((FinleyDomain)arg1, (Data)arg2, (Data)arg3, (Data)mat, (Data)D, (object)d) None : ¶
adds a PDE onto the lumped stiffness matrix
- Parameters
mat (
Data
) –D (
Data
) –d (
Data
) –useHRZ (bool) –
- addPDEToRHS((FinleyDomain)arg1, (Data)arg2, (Data)rhs, (Data)X, (Data)Y, (Data)y, (Data)y_contact) None : ¶
adds a PDE onto the stiffness matrix mat and a rhs
- Parameters
rhs (
Data
) –X (
Data
) –Y (
Data
) –y (
Data
) –y_contact (
Data
) –
- addPDEToSystem((FinleyDomain)arg1, (Operator)arg2, (Data)arg3, (Data)mat, (Data)rhs, (Data)A, (Data)B, (Data)C, (Data)D, (Data)X, (Data)Y, (Data)d, (Data)y, (Data)d_contact, (Data)y_contact) None : ¶
adds a PDE onto the stiffness matrix mat and a rhs
- Parameters
mat (
OperatorAdapter
) –rhs (
Data
) –A (
Data
) –B (
Data
) –C (
Data
) –D (
Data
) –X (
Data
) –Y (
Data
) –d (
Data
) –d_contact (
Data
) –y_contact (
Data
) –
- addPDEToTransportProblem((FinleyDomain)arg1, (TransportProblem)arg2, (Data)arg3, (Data)tp, (Data)source, (Data)M, (Data)A, (Data)B, (Data)C, (Data)D, (Data)X, (Data)Y, (Data)d, (Data)y, (Data)d_contact, (Data)y_contact) None : ¶
- Parameters
tp (
AbstractTransportProblem
) –source (
Data
) –M (
Data
) –A (
Data
) –B (
Data
) –C (
Data
) –D (
Data
) –X (
Data
) –Y (
Data
) –d (
Data
) –y (
Data
) –d_contact (
Data
) –y_contact (
Data
) –
- dump((FinleyDomain)arg1, (str)fileName) None : ¶
dumps the mesh to a file with the given name.
- getConnectivityInfo((FinleyDomain)arg1) numpy.ndarray : ¶
- Returns
returns point and connectivity information
- Return type
Data
- getDataShape((FinleyDomain)arg1, (object)functionSpaceCode) object : ¶
- Returns
a pair (dps, ns) where dps=the number of data points per sample, and ns=the number of samples
- Return type
tuple
- getDescription((FinleyDomain)arg1) str : ¶
- Returns
a description for this domain
- Return type
string
- getDim((FinleyDomain)arg1) int : ¶
- Return type
int
- getMPIRank((FinleyDomain)arg1) int : ¶
- Returns
the rank of this process
- Return type
int
- getMPISize((FinleyDomain)arg1) int : ¶
- Returns
the number of processes used for this
Domain
- Return type
int
- getNormal((FinleyDomain)arg1) Data : ¶
- Returns
boundary normals at the quadrature point on the face elements
- Return type
Data
- getNumDataPointsGlobal((FinleyDomain)arg1) int : ¶
- Returns
the number of data points summed across all MPI processes
- Return type
int
- getSize((FinleyDomain)arg1) Data : ¶
- Returns
the element size
- Return type
Data
- getSystemMatrixTypeId((FinleyDomain)arg1, (object)options) int : ¶
- Returns
the identifier of the matrix type to be used for the global stiffness matrix when particular solver options are used.
- Return type
int
- Parameters
options (
SolverBuddy
) –
- getTag((FinleyDomain)arg1, (str)name) int : ¶
- Returns
tag id for
name
- Return type
string
- getTransportTypeId((FinleyDomain)arg1, (object)solver, (object)preconditioner, (object)package, (object)symmetry) int : ¶
- Returns
the identifier of the transport problem type to be used when a particular solver, preconditioner, package and symmetric matrix is used.
- Return type
int
- Parameters
solver (
int
) –preconditioner (
int
) –package (
int
) –symmetry (
int
) –
- getVTKElementType((FinleyDomain)arg1) int : ¶
- Returns
returns the VTK code for this element type
- Return type
Data
- getX((FinleyDomain)arg1) Data : ¶
- Returns
locations in the FEM nodes
- Return type
Data
- isValidTagName((FinleyDomain)arg1, (str)name) bool : ¶
- Returns
True is
name
corresponds to a tag- Return type
bool
- newOperator((FinleyDomain)arg1, (object)row_blocksize, (FunctionSpace)row_functionspace, (object)column_blocksize, (FunctionSpace)column_functionspace, (object)type) Operator : ¶
creates a stiffness matrix and initializes it with zeros
- Parameters
row_blocksize (
int
) –row_functionspace (
FunctionSpace
) –column_blocksize (
int
) –column_functionspace (
FunctionSpace
) –type (
int
) –
- newTransportProblem((FinleyDomain)theta, (object)blocksize, (FunctionSpace)functionspace, (object)type) TransportProblem : ¶
creates a TransportProblem
- Parameters
theta (
float
) –blocksize (
int
) –functionspace (
FunctionSpace
) –type (
int
) –
- onMasterProcessor((FinleyDomain)arg1) bool : ¶
- Returns
True if this code is executing on the master process
- Return type
bool
- print_mesh_info((FinleyDomain)arg1[, (object)full=False]) None : ¶
- Parameters
full (
bool
) –
- setTagMap((FinleyDomain)arg1, (str)name, (object)tag) None : ¶
Give a tag number a name.
- Parameters
name (
string
) – Name for the tagtag (
int
) – numeric id
- Note
Tag names must be unique within a domain
- setX((FinleyDomain)arg1, (Data)arg) None : ¶
assigns new location to the domain
- Parameters
arg (
Data
) –
- showTagNames((FinleyDomain)arg1) str : ¶
- Returns
A space separated list of tag names
- Return type
string
- write((FinleyDomain)arg1, (str)filename) None : ¶
Write the current mesh to a file with the given name.
- esys.finley.GMSHDesign¶
alias of
esys.pycad.gmsh.Design
Functions¶
- esys.finley.Brick(n0=1, n1=1, n2=1, order=1, l0=1.0, l1=1.0, l2=1.0, periodic0=0, periodic1=0, periodic2=0, integrationOrder=- 1, reducedIntegrationOrder=- 1, useElementsOnFace=1, useFullElementOrder=0, optimize=0, **kwargs)¶
- __Brick_driver( (list)params) -> Domain :
Creates a rectangular mesh with n0 x n1 x n2 elements over the brick [0,l0] x [0,l1] x [0,l2].
- param n0
number of elements in direction 0
- type n0
int
- param n1
number of elements in direction 1
- type n1
int
- param n2
number of elements in direction 2
- type n2
int
- param order
=1, =-1 or =2 gives the order of shape function. If -1 macro elements of order 1 are used.
- param l0
length of side 0
- type l0
float
- param l1
length of side 1
- type l1
float
- param l2
length of side 2
- type l2
float
- param periodic0
whether or not boundary conditions are periodic in direction 0
- type periodic0
bool
- param periodic1
whether or not boundary conditions are periodic in direction 1
- type periodic1
bool
- param periodic2
whether or not boundary conditions are periodic in direction 2
- type periodic2
bool
- param integrationOrder
order of the quadrature scheme. If integrationOrder<0 the integration order is selected independently.
- param reducedIntegrationOrder
order of the quadrature scheme. If reducedIntegrationOrder<0 the integration order is selected independently.
- param useElementsOnFace
whether or not to use elements on face
- type useElementsOnFace
int
- param useFullElementOrder
Whether or not to use Hex27 elements
- type useFullElementOrder
bool
- param optimize
Enable optimisation of node labels
- type optimize
bool
- esys.finley.GetMeshFromFile(filename, **kwargs)¶
Reads a mesh from a file, determines the reader to use based on the file extension. All cases require a filename and gmsh files require a number of dimensions (it doesn’t hurt to pass this in all the time). Other keyword args come from the underlying reader functions.
- esys.finley.GlueFaces((list)meshList[, (object)safetyFactor=0.2[, (object)tolerance=1e-08[, (object)optimize=True]]]) Domain : ¶
Detects matching faces in the mesh, removes them from the mesh and joins the elements touched by the face elements.
- esys.finley.JoinFaces((list)meshList[, (object)safetyFactor=0.2[, (object)tolerance=1e-08[, (object)optimize=True]]]) Domain : ¶
Detects matching faces in the mesh and replaces them by joint elements.
- esys.finley.LoadMesh([(str)fileName='file.nc']) Domain : ¶
- Return type
- esys.finley.MakeDomain(design, integrationOrder=- 1, reducedIntegrationOrder=- 1, optimizeLabeling=True, useMacroElements=False)¶
Creates a Finley
Domain
from aesys.pycad.design.Design
object. Currently only gmsh is supported.- Parameters
design (
esys.pycad.design.Design
) – the geometryintegrationOrder (
int
) – integration order. If -1 the default is used.reducedIntegrationOrder (
int
) – reduced integration order. If -1 the default is used.optimizeLabeling (
bool
) – if set the labeling of the mesh nodes is optimizeduseMacroElements (
bool
) – uses macro elements.
- Returns
the Finley domain defined by the design
- Return type
Domain
- esys.finley.Merge((list)meshList) Domain : ¶
Merges a list of meshes into one mesh.
- Return type
Domain
- esys.finley.ReadGmsh(fileName, numDim, integrationOrder=- 1, reducedIntegrationOrder=- 1, optimize=True, useMacroElements=False, **kwargs)¶
- __ReadGmsh_driver( (list)params) -> Domain :
Read a gmsh mesh file
- rtype
Domain
- param fileName
- type fileName
string
- param integrationOrder
order of the quadrature scheme. If integrationOrder<0 the integration order is selected independently.
- type integrationOrder
int
- param reducedIntegrationOrder
order of the quadrature scheme. If reducedIntegrationOrder<0 the integration order is selected independently.
- param optimize
Enable optimisation of node labels
- type optimize
bool
- param useMacroElements
Enable the usage of macro elements instead of second order elements.
- type useMacroElements
bool
- esys.finley.ReadMesh(filename, integrationOrder=- 1, reducedIntegrationOrder=- 1, optimize=True, **kwargs)¶
- __ReadMesh_driver( (list)params) -> Domain :
Read a mesh from a file. For MPI parallel runs fan out the mesh to multiple processes.
- rtype
Domain
- param fileName
- type fileName
string
- param integrationOrder
order of the quadrature scheme. If integrationOrder<0 the integration order is selected independently.
- type integrationOrder
int
- param reducedIntegrationOrder
order of the quadrature scheme. If reducedIntegrationOrder<0 the integration order is selected independently.
- param optimize
Enable optimisation of node labels
- type optimize
bool
- esys.finley.Rectangle(n0=1, n1=1, order=1, l0=1.0, l1=1.0, periodic0=False, periodic1=False, integrationOrder=- 1, reducedIntegrationOrder=- 1, useElementsOnFace=None, useFullElementOrder=0, optimize=0, **kwargs)¶
- __Rectangle_driver( (list)args) -> Domain :
Creates a rectangular mesh with n0 x n1 elements over the brick [0,l0] x [0,l1].
- param n0
- type n0
- param n1
- type n1
- param order
=1, =-1 or =2 gives the order of shape function. If -1 macro elements of order 1 are used.
- param l0
length of side 0
- param l1
- param integrationOrder
order of the quadrature scheme. If integrationOrder<0 the integration order is selected independently.
- param reducedIntegrationOrder
order of the quadrature scheme. If reducedIntegrationOrder<0 the integration order is selected independently.
- param useElementsOnFace
whether or not to use elements on face
- type useElementsOnFace
int
:param periodic0: whether or not boundary conditions are periodic- param periodic1
- param useFullElementOrder
Whether or not to use Rec9 elements
- type useFullElementOrder
bool
- param optimize
Enable optimisation of node labels
- type optimize
bool