Gismo related documentation

gismo.GenerateMicrostructures_fun.generate(size_1: int, size_2: int, size_3: int, BOX_LENGTH: float, BOX_HEIGHT: float, EPS: float, INLET_BOUNDARY_ID: int, OUTLET_BOUNDARY_ID: int, knots_y: ~numpy.ndarray, TILING: list, CLOSING_FACE: str, MICROTILE: str) -> (<class 'splinepy.microstructure.microstructure.Microstructure'>, <class 'list'>)

Generate a microstructure

For a given size_1, size_2, and size_3, generate a microstructure with the given parameters.

Parameters:

  • size_1 (int) – Size of microtile on lowest level

  • size_2 (int) – Size of microtile on middle level

  • size_3 (int) – Size of microtile on highest level

  • BOX_LENGTH (float) – Length of the box

  • BOX_HEIGHT (float) – Height of the box

  • EPS (float) – Precision

  • INLET_BOUNDARY_ID (int) – Identifier for inlet boundary

  • OUTLET_BOUNDARY_ID (int) – Identifier for outlet boundary

  • knots_y (np.ndarray) – Knots for y direction

  • TILING (list) – Tiling of the microstructure [nx, ny] with nx, ny integers the number of tiles in x and y direction

  • CLOSING_FACE (str) – Closing face of the microstructure, either “x” or “y”

  • MICROTILE (str) – Type of microtile. See splinepy.microstructure.Microstructure for details

Returns:

  • sp.microstructure.Microstructure – Microstructure

  • list – List of additional blocks for gismo export

Raises:

ValueError – If the parametric function doesnt find the according y-value

This code is copied from: https://github.com/markriegler/splinepy/blob/ft-gismo-bc-export/splinepy/io/gismo.py as the installation from this splinepy version did not work properly.

class gismo.gismo_export.AdditionalBlocks

Bases: object

Class to manage additional XML blocks for gismo export, including boundary conditions, function definitions, and assembly options. Each block is stored as a Python dictionary, which is then later handled by the export function.

add_assembly_options(block_id, dirichlet_strategy=11, dirichlet_values=101, interface_strategy=1, bda=2, bdb=1, bdo=0.333, qua=1, qub=1, qurule=1, comment=None)

Create Python dictionary of g+smo assembly options to be used in gismo export function

Parameters:

  • block_id (int) – ID of the OptionList-block

  • dirichlet_strategy (int) – Method for enforcement of Dirichlet BCs [11..14]

  • dirichlet_values (int) – Method for computation of Dirichlet DoF values [100..103]

  • interface_strategy (int) – Method of treatment of patch interfaces [0..3]

  • bda (int) – Estimated nonzeros per column of the matrix: bdA*deg + bdB

  • bdb (int) – Estimated nonzeros per column of the matrix: bdA*deg + bdB

  • bdo (int) – Overhead of sparse mem. allocation: (1+bdO)(bdA*deg + bdB) [0..1]

  • qua (int) – Number of quadrature points: quA*deg + quB

  • qub (int) – Number of quadrature points: quA*deg + quB

  • qurule (int) – Quadrature rule [1:GaussLegendre,2:GaussLobatto]

  • comment (str) – Comment above OptionList-block

add_boundary_conditions(block_id, dim, function_list, bc_list, cv_list=None, unknown_id=0, multipatch_id=0, comment=None)

Create Python dictionary of boundary function’s block to be used in gismo export function.

Parameters:

  • block_id (int) – ID number of boundaryConditions-block

  • dim (int) – Dimension of physical space

  • function_list (list<str/tuple<str,str,...>>) – List of functions where entries are strings of the function of the BC. If tuple is given, it will be interpreted as component-wise functions (e.g. for vector valued functions)

  • bc_list (list<tuple<str, str, int>>) – List of which boundary attains which type of boundary condition. First entry is the boundary name, second the type of condition (e.g. “Dirichlet”), third is on which variable/unknown it should be applied

  • cv_list (list<tuple<str, str, str, str>>) – List of which corner values of each patch should have a boundary condition. First entry is the variable/unknown, second the patch number, third is on which corner it should be applied and fourth is the function or value

  • unknown_id (int (0)) – ID of the unknown/variable on which this boundaryConditions block should be applied

  • multipatch_id (int (0)) – ID of the multipatch

  • comment (str) – Comment above boundaryFunctions-block

add_function(dim, block_id, function_string, comment=None)

Create Python dictionary of custom function’s block to be used in gismo export function.

Parameters:

  • dim (int) – Physical dimension of the problem

  • id (int) – ID of function block

  • function_string (str or tuple<str,str>) – Function as string. If tuple is given, function will be vector-valued.

  • comment (str) – Comment above function block

to_list()

Returns list of created XML-blocks

Returns:

blocks – List of XML-blocks

Return type:

list<dict>

gismo.gismo_export.export(fname, multipatch=None, indent=True, labeled_boundaries=True, additional_blocks=None, export_fields=False, field_mask=None, as_base64=False)

Export as gismo readable xml geometry file Use gismo-specific xml-keywords to export (list of) splines. All Bezier patches are excluded as their respective non uniform counterpart

Parameters:

  • fname (string) – name of the output file

  • spline_list (Multipatch (preferred)) – (list of) Spline-Types in splinepy format

  • indent (bool) – Appends white spaces using xml.etree.ElementTree.indent, if possible.

  • labeled_boundaries (bool) – Writes boundaries with labels into the MultiPatch part of the XML

  • additional_blocks (list) – List of dictionaries, that specify model related options, like boundary conditions, assembler options, etc.. The dictionaries must have the following keys, ‘tag’->string, ‘text’->string (optional), ‘attributes’->dictionary (optional), ‘children’->list in the same format (optional)

  • export_fields (bool) – Export fields to separate file ending with fields.xml, e.g., filename.xml.fields.xml. Only non-zero splines are exported to save memory

  • field_mask (list) – Selection of active fields that are exported (to save memory and speed up export)

  • as_base64 (bool (False)) – If set, coordinates, weights and knot-vectors are exported in B64 encoding. Available in gismo, once gismo/gismo#634 merges.

Return type:

None