BEGIN_C_DECLS
/*============================================================================
* User function definitions
*============================================================================*/
/*----------------------------------------------------------------------------
* Define periodic faces.
*
* This is done by calling one of the cs_join_perio_add_*() functions for
* each periodicity to add.
*
* The first arguments to cs_join_perio_add_() are the same as for
* mesh joining:
* sel_criteria <-- boundary face selection criteria string
* fraction <-- value of the fraction parameter;
* the initial tolerance radius associated to each vertex
* is equal to the lenght of the shortest incident edge,
* multiplied by this fraction.
* plane <-- value of the plane parameter;
* when subdividing faces, 2 faces are considered
* coplanar and may be joined if angle between their
* normals (in degrees) does not exceed this parameter.
* verbosity <-- level of verbosity required
*
* The last arguments depend on the type of periodicity to define,
* and are described below.
*
* The function returns a number (1 to n) associated with the
* new joining. This number may optionnally be used to assign advanced
* parameters to the joining.
*----------------------------------------------------------------------------*/
void
cs_user_periodicity(void)
{
int join_num;
int verbosity = 1, visualization = 1;
float fraction = 0.10, plane = 0.25;
fraction = 0.10;
plane = 25.0;
verbosity = 1; /* processor-local files if > 1, debug level if >= 3 */
visualization = 1; /* debug level if >= 3 */
/* Example 1: define a periodicity of translation */
/* ---------------------------------------------- */
{
const double translation[3] = {1.0, 0.0, 0.0}; /* Translation vector */
join_num = cs_join_perio_add_translation("98 or 99",
fraction,
plane,
verbosity,
visualization,
translation);
}
/* Example 2: define a periodicity of rotation */
/* ------------------------------------------- */
{
double theta = 20; /* angle in degrees */
double axis[3] = {1.0, 0, 0}; /* axis of rotation */
double invariant[3] = {0, 0, 0}; /* invariant point */
/* change default values */
fraction = 0.2;
verbosity = 2;
join_num = cs_join_perio_add_rotation("3",
fraction,
plane,
verbosity,
visualization,
theta,
axis,
invariant);
/* restore default values */
fraction = 0.1;
verbosity = 1;
}
/* Example 3: define a general periodicity */
/* --------------------------------------- */
/* We define a general transformation using a homogeneous matrix:
*
* We define the first 3 rows of a 4x4 matrix:
* _ _
* | r11 r12 r13 tx | t(x,y,z) : translation vector
* | r21 r22 r23 ty | r(i,j) : rotation matrix
* | r31 r32 r33 tz |
* |_ 0 0 0 1 _|
*
* Transformations may be combined using matrix multiplication,
* so this be used for helecoidal transformations for instance. */
{
double matrix[3][4] = {{1., 0., 0., 0.5},
{0., 1., 0., 0.},
{0., 0., 1., 0.}};
join_num = cs_join_perio_add_mixed("all[]",
fraction,
plane,
verbosity,
visualization,
matrix);
}
/*--------------------------------------------------------------------------*/
/* Example with advanced parameters;
Advanced parameters may be modified to solve errors during the
joining step or to get a better mesh quality. */
{
/* Merge tolerance factor:
* used to locally modify the tolerance associated to each
* vertex BEFORE the merge step.
* = 0 => no vertex merge;
* < 1 => vertex merge is more strict. It may increase the number
* of tolerance reduction and so define smaller subset of
* vertices to merge together but it can drive to loose
* intersections;
* = 1 => no change;
* > 1 => vertex merge is less strict. The subset of vertices able
* to be merged together is greater. */
double mtf = 1.00;
/* Pre-merge factor:
* this parameter is used to define a bound under which two vertices
* are merged before the merge step.
* Tolerance limit for the pre-merge = pmf * fraction. */
double pmf = 0.10;
/* Tolerance computation mode:
*
* 1: (default) tol = min. edge length related to a vertex * fraction
* 2: tolerance is computed like in mode 1 with in addition, the
* multiplication by a coefficient equal to the max sin(e1, e2)
* where e1 and e2 are two edges sharing the same vertex V for which
* we want to compute the tolerance.
* 11: as 1 but taking into account only the selected faces
* 12: as 2 but taking into account only the selected faces */
int tcm = 1;
/* Intersection computation mode:
* 1: (default) Original algorithm.
* Try to snap intersection to extremity.
* 2: New intersection algorithm.
* Avoid snapping intersection to extremity. */
int icm = 1;
/* Maximum number of equivalence breaks which is
* enabled during the merge step. */
int max_break = 500;
/* Maximum number of sub-faces when splitting a selected face. */
int max_sub_face = 100;
/* tml, tmb and tmr are parameters of the searching algorithm for
* face intersections between selected faces (octree structure).
* Useful to adjust speed vs. memory consumption. */
/* Tree Max Level:
* deepest level reachable when building the tree */
int tml = 30;
/* Tree Max Boxes:
* max. number of bounding boxes which can be linked to a leaf of the tree
* (not necessary true for the deepest level) */
int tmb = 25;
/* Tree Max. Ratio:
* stop refining the tree structure when
* number of bounding boxes > tmr * number of faces to locate.
* In parallel, a separate (usually lower) value may be set for
* the initial coarse tree used to determine distribution.
* Reducing this will help reduce memory consumption. */
double tmr = 5.0;
double tmr_distrib = 2.0;
/* Set advanced parameters */
cs_join_set_advanced_param(join_num,
mtf, pmf, tcm, icm,
max_break, max_sub_face,
tml, tmb, tmr, tmr_distrib);
}
}
void
cs_gui_mesh_define_periodicities(void)
{
int perio_id;
double angle, trans[3], axis[3], invariant[3], matrix[3][4];
int n_perio = 0;
int n_modes = 0;
char **modes = NULL;
char *path = NULL;
if (!cs_gui_file_is_loaded())
return;
n_perio
= cs_gui_get_tag_number("/solution_domain/periodicity/face_periodicity",
1);
if (n_perio == 0)
return;
/* Get modes associated with each periodicity */
path = cs_xpath_init_path();
cs_xpath_add_elements(&path, 3,
"solution_domain",
"periodicity", "face_periodicity");
cs_xpath_add_attribute(&path, "mode");
modes = cs_gui_get_attribute_values(path, &n_modes);
if (n_modes != n_perio)
bft_error(__FILE__, __LINE__, 0,
_("Number of periodicities (%d) and modes (%d) do not match."),
n_perio, n_modes);
BFT_FREE(path);
/* loop on periodicities */
for (perio_id = 0; perio_id < n_perio; perio_id++) {
char *selector_s = _get_periodicity("selector", perio_id+1);
char *fraction_s = _get_periodicity("fraction", perio_id+1);
char *plane_s = _get_periodicity("plane", perio_id+1);
char *verbosity_s = _get_periodicity("verbosity", perio_id+1);
char *visu_s = _get_periodicity("visualization", perio_id+1);
double fraction = (fraction_s != NULL) ? atof(fraction_s) : 0.1;
double plane = (plane_s != NULL) ? atof(plane_s) : 25.0;
int verbosity = (verbosity_s != NULL) ? atoi(verbosity_s) : 1;
int visualization = (visu_s != NULL) ? atoi(visu_s) : 1;
if (!strcmp(modes[perio_id], "translation")) {
_get_periodicity_translation(perio_id+1, trans);
cs_join_perio_add_translation(selector_s,
fraction,
plane,
verbosity,
visualization,
trans);
}
else if (!strcmp(modes[perio_id], "rotation")) {
_get_periodicity_rotation(perio_id+1, &angle, axis, invariant);
cs_join_perio_add_rotation(selector_s,
fraction,
plane,
verbosity,
visualization,
angle,
axis,
invariant);
}
else if (!strcmp(modes[perio_id], "mixed")) {
_get_periodicity_mixed(perio_id+1, matrix);
cs_join_perio_add_mixed(selector_s,
fraction,
plane,
verbosity,
visualization,
matrix);
}
else
bft_error(__FILE__, __LINE__, 0,
_("Periodicity mode \"%s\" unknown."), modes[perio_id]);
#if _XML_DEBUG_
bft_printf("==> cs_gui_mesh_define_periodicities\n");
bft_printf("--selector = %s\n", selector_s);
bft_printf("--fraction = %s\n", fraction_s);
bft_printf("--plane = %s\n", plane_s);
bft_printf("--verbosity = %s\n", verbosity_s);
bft_printf("--visualization = %s\n", visu_s);
#endif
BFT_FREE(selector_s);
BFT_FREE(fraction_s);
BFT_FREE(plane_s);
BFT_FREE(verbosity_s);
BFT_FREE(visu_s);
}
for (perio_id = 0; perio_id < n_perio; perio_id++)
BFT_FREE(modes[perio_id]);
BFT_FREE(modes);
}