) |
) |
self) |
In the current implementation this is always 3, though you can define 2D Formices by given only two coordinates: the third will automatically be set to zero.
) |
This is the product of the nelems() and nplex().
) |
The shape of a Formex is the shape of its data array, i.e. a tuple (nelems, nplex, ndim).
) |
) |
) |
) |
) |
) |
) |
) |
The bounding box is the smallest rectangular volume in global coordinates, such that no points of the Formex are outside the box. It is returned as a [2,3] array: the first row holds the minimal coordinates and the second one the maximal.
) |
) |
) |
) |
) |
The bounding sphere is the smallest sphere with center in the center() of the Formex, and such that no points of the Formex are lying outside the sphere. It is not necessarily the smallest sphere surrounding all points of the Formex.
) |
The centroid of an element is the point whose coordinates are the mean values of all points of the element. The return value is a plex-1 Formex.
p,n) |
p is a point specified by 3 coordinates. n is the normal vector to a plane, specified by 3 components.
The return value is a (nelems(),nplex()) shaped array with the distance of each point to the plane containing the point p and having normal n. Distance values are positive if the point is on the side of the plane indicated by the positive normal.
p,q) |
p and q are two points specified by 3 coordinates.
The return value is a (nelems(),nplex()) shaped array with the distance of each point to the line through p and q. All distance values are positive or zero.
p) |
p is a point specified by 3 coordinates.
The return value is a (nelems(),nplex()) shaped array with the distance of each point to the line through p and q. All distance values are positive or zero.
nodesperbox=1,repeat=True,rtol=1.e-5,atol=1.e-5) |
A tuple of two arrays is returned. The first is a float array with the coordinates of the unique nodes of the Formex. The second is an integer array with the node numbers connected by each element. The elements come in the same order as they are in the Formex, but the order of the nodes is unspecified. By the way, the reverse operation of coords,elems = feModel(F) is accomplished by F = Formex(coords[elems]). There is a (very small) probability that two very close nodes are not equivalenced by this procedure. Use it multiple times with different parameters to check.
rtol and atol are the relative, resp. absolute tolerances used to decide whether any nodal coordinates are considered to be equal.