Internal functions used by VALIB. More...

Average corotation and contrarotation helper functions
VA_DEVICE_FUN void	valib_resvortstrain2D (VA_REAL SO, VA_REAL residual_vorticity, VA_REAL residual_strain, VA_REAL principal_axis)

Numerical integrands on a unit sphere
VA_DEVICE_FUN void	valib_integrand_resvort (VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL avgcorot)
VA_DEVICE_FUN void	valib_integrand_resstrain (VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL S_RAVG)
VA_DEVICE_FUN void	valib_integrand_surface (VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result)

Small helper functions
VA_DEVICE_FUN VA_REAL	valib_sign (VA_REAL a, VA_REAL b)
VA_DEVICE_FUN VA_COMPLEX	valib_makec (VA_REAL a, VA_REAL b)
VA_DEVICE_FUN VA_REAL	valib_creal (VA_COMPLEX a)
VA_DEVICE_FUN VA_REAL	valib_cimag (VA_COMPLEX a)
VA_DEVICE_FUN VA_REAL	valib_cabs (VA_COMPLEX a)
VA_DEVICE_FUN VA_COMPLEX	valib_cadd (VA_COMPLEX a, VA_COMPLEX b)
VA_DEVICE_FUN VA_COMPLEX	valib_cmultiply (VA_COMPLEX a, VA_COMPLEX b)
VA_DEVICE_FUN VA_REAL	valib_cubic_root (VA_REAL x)
VA_DEVICE_FUN VA_REAL	valib_max (VA_REAL a, VA_REAL b)
VA_DEVICE_FUN VA_REAL	valib_min (VA_REAL a, VA_REAL b)
VA_DEVICE_FUN VA_REAL	valib_square (VA_REAL a)
VA_DEVICE_FUN void	valib_swap_values (VA_REAL a, VA_REAL b)
VA_DEVICE_FUN void	valib_solve_cubic_equation (VA_REAL a, VA_REAL b, VA_REAL c, VA_REAL d, VA_REAL x1, VA_COMPLEX x2, VA_COMPLEX x3, int info)

Operations on vectors of length 3
VA_DEVICE_FUN void	valib_vec_zero3 (VA_REAL *v)
VA_DEVICE_FUN VA_REAL	valib_scalar_prod3 (VA_REAL v1, VA_REAL v2)
VA_DEVICE_FUN void	valib_cross_prod3 (VA_REAL v1, VA_REAL v2, VA_REAL *v3)

Operations on 3x3 matrices
VA_DEVICE_FUN void	valib_mat_zero3 (VA_REAL *A)
VA_DEVICE_FUN void	valib_mat_copy3 (VA_REAL A, VA_REAL B)
VA_DEVICE_FUN void	valib_norm_frobenius3 (VA_REAL A, VA_REAL norm)
VA_DEVICE_FUN void	valib_trace3 (VA_REAL A, VA_REAL tr3)
VA_DEVICE_FUN void	valib_second_invariant3 (VA_REAL A, VA_REAL second_invariant)
VA_DEVICE_FUN void	valib_determinant3 (VA_REAL A, VA_REAL det)
VA_DEVICE_FUN void	valib_transpose3 (VA_REAL *A)
VA_DEVICE_FUN void	valib_sym_antisym3 (VA_REAL *A)
VA_DEVICE_FUN void	valib_deviatorise3 (VA_REAL *A)
VA_DEVICE_FUN void	valib_matvec3_prod (int trans_a, VA_REAL A, VA_REAL x, VA_REAL *y)
VA_DEVICE_FUN void	valib_rank1_update3 (VA_REAL alpha, VA_REAL x, VA_REAL y, VA_REAL *A)
VA_DEVICE_FUN void	valib_tatt3 (VA_REAL A, VA_REAL T, VA_REAL *TATT)
VA_DEVICE_FUN void	valib_matmat3_prod (int trans_a, int trans_b, VA_REAL A, VA_REAL B, VA_REAL *C)
VA_DEVICE_FUN void	valib_ttat3 (VA_REAL A, VA_REAL T, VA_REAL *TTAT)
VA_DEVICE_FUN void	valib_gram_schmidt3 (VA_REAL *Q)
VA_DEVICE_FUN void	valib_eigenvalues_sym3 (VA_REAL A, VA_REAL eigval)
VA_DEVICE_FUN int	valib_find_pivot3 (VA_REAL *A, int jcol, int start)
VA_DEVICE_FUN void	valib_find_nullspace3 (VA_REAL A, VA_REAL x, int *info)
VA_DEVICE_FUN void	valib_constructQ3 (VA_REAL sina, VA_REAL cosa, VA_REAL sinb, VA_REAL cosb, VA_REAL sing, VA_REAL cosg, VA_REAL *Q)
VA_DEVICE_FUN void	valib_constructQfromN3 (VA_REAL n, VA_REAL Q)

Helper functions for the Q-criterion and its variants.
VA_DEVICE_FUN void	valib_qcriterion_blending (VA_REAL SO, VA_REAL blending_ratio, VA_REAL qcriterion_blended)

Numerical integrators on a sphere
VA_DEVICE_FUN void	valib_integrator_fibonacci (VA_DEVICE_ADDR void integrand(VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result), VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result)
VA_DEVICE_FUN void	valib_integrator_bazant (VA_DEVICE_ADDR void integrand(VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result), VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result)
VA_DEVICE_FUN void	valib_integrator_rectangle (VA_DEVICE_ADDR void integrand(VA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result), VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result, int numIntervals)

Triple decomposition helper functions.
VA_DEVICE_FUN void	valib_get_brf (VA_DEVICE_ADDR VA_REAL A, int num_intervals, VA_REAL alpha_brf, VA_REAL beta_brf, VA_REAL gamma_brf, VA_REAL goal_brf)
VA_DEVICE_FUN void	valib_triple_decomposition_base (VA_DEVICE_ADDR VA_REAL A, int num_intervals, VA_DEVICE_ADDR VA_REAL residual_vorticity_strain, VA_DEVICE_ADDR VA_REAL shear_vorticity_strain, VA_DEVICE_ADDR VA_REAL alpha, VA_DEVICE_ADDR VA_REAL beta, VA_DEVICE_ADDR VA_REAL *gamma)
	Triple Decomposition Method (TDM, a.k.a. residual vorticity and residual strain rate) [14], [10].

Detailed Description

Internal functions used by VALIB.

Function Documentation

◆ valib_cabs()

VA_DEVICE_FUN VA_REAL valib_cabs ( VA_COMPLEX a )

#include "linalg3.h"

Returns absolute value of a complex number a.

Definition at line 98 of file linalg3.h.

◆ valib_cadd()

VA_DEVICE_FUN VA_COMPLEX valib_cadd	(	VA_COMPLEX	a,
		VA_COMPLEX	b )

#include "linalg3.h"

Returns c = a + b for complex numbers.

Definition at line 109 of file linalg3.h.

◆ valib_cimag()

VA_DEVICE_FUN VA_REAL valib_cimag ( VA_COMPLEX a )

#include "linalg3.h"

Returns imaginary part of a complex number a.

Definition at line 75 of file linalg3.h.

◆ valib_cmultiply()

VA_DEVICE_FUN VA_COMPLEX valib_cmultiply	(	VA_COMPLEX	a,
		VA_COMPLEX	b )

#include "linalg3.h"

Returns c = a * b for complex numbers.

Definition at line 118 of file linalg3.h.

◆ valib_constructQ3()

VA_DEVICE_FUN void valib_constructQ3	(	VA_REAL	sina,
		VA_REAL	cosa,
		VA_REAL	sinb,
		VA_REAL	cosb,
		VA_REAL	sing,
		VA_REAL	cosg,
		VA_REAL *	Q )

#include "linalg3.h"

Get orthogonal matrix Q column-wise based on sines and cosines of 3 angles of rotations of coordinate system along coordinate axes by z, y', and z'', see [14].

Definition at line 898 of file linalg3.h.

◆ valib_constructQfromN3()

VA_DEVICE_FUN void valib_constructQfromN3	(	VA_REAL *	n,
		VA_REAL *	Q )

#include "linalg3.h"

Get orthogonal matrix Q column-wise based on components of a normal vector n which corresponds to the last row of Q. Vector n need not be normalized to \( \|n\|_2 = 1\).

Definition at line 920 of file linalg3.h.

◆ valib_creal()

VA_DEVICE_FUN VA_REAL valib_creal ( VA_COMPLEX a )

#include "linalg3.h"

Returns real part of a complex number a.

Definition at line 52 of file linalg3.h.

◆ valib_cross_prod3()

VA_DEVICE_FUN void valib_cross_prod3	(	VA_REAL *	v1,
		VA_REAL *	v2,
		VA_REAL *	v3 )

#include "linalg3.h"

Function for cross product of two vectors

\[ v_3 = v_1 \times v_2 \]

.

Definition at line 328 of file linalg3.h.

◆ valib_cubic_root()

VA_DEVICE_FUN VA_REAL valib_cubic_root ( VA_REAL x )

#include "linalg3.h"

Function for finding the real cubic root.

Definition at line 134 of file linalg3.h.

◆ valib_determinant3()

VA_DEVICE_FUN void valib_determinant3	(	VA_REAL *	A,
		VA_REAL *	det )

#include "linalg3.h"

Compute the determinant of a 3x3 matrix (3rd tensor invariant).

Definition at line 410 of file linalg3.h.

◆ valib_deviatorise3()

VA_DEVICE_FUN void valib_deviatorise3 ( VA_REAL * A )

#include "linalg3.h"

Find the deviatoric part of a 3x3 matrix, i.e. shift the diagonal such that the matrix has zero trace,

\[ A_D = A - \frac{1}{3}\mbox{trace}(A)\ I. \]

Definition at line 461 of file linalg3.h.

◆ valib_eigenvalues_sym3()

VA_DEVICE_FUN void valib_eigenvalues_sym3	(	VA_REAL *	A,
		VA_REAL *	eigval )

#include "linalg3.h"

Computing eigenvalues. Uses the method given in [17].

Definition at line 626 of file linalg3.h.

◆ valib_find_nullspace3()

VA_DEVICE_FUN void valib_find_nullspace3	(	VA_REAL *	A,
		VA_REAL *	x,
		int *	info )

#include "linalg3.h"

Find a nontrivial solution to the equation Ax = 0. This function is used for searching the eigenvector corresponding to real eigenvalue in case of simple eigenvalues. It expects one-dimensional nullspace, otherwise returns an error.

Definition at line 709 of file linalg3.h.

◆ valib_find_pivot3()

VA_DEVICE_FUN int valib_find_pivot3	(	VA_REAL *	A,
		int	jcol,
		int	start )

#include "linalg3.h"

Find a the index of the pivot in the column of a matrix. It searches the jcol-th column of matrix A from index start to 3, i.e. index of the maximum of absolute value of A(start:,jcol).

Definition at line 682 of file linalg3.h.

◆ valib_get_brf()

VA_DEVICE_FUN void valib_get_brf	(	VA_DEVICE_ADDR VA_REAL *	A,
		int *	num_intervals,
		VA_REAL *	alpha_brf,
		VA_REAL *	beta_brf,
		VA_REAL *	gamma_brf,
		VA_REAL *	goal_brf )

#include "triple_decomposition.h"

Function for determining the Basic Reference Frame (BRF) introduced in [14]. This is the coordinate frame in which the shear is maximized and in which it is later eliminated. The BRF is found as an optimization problem, searching the maxima in formula (10) of [14], i.e.

\[ \max_{\alpha \in [0,\pi], \\ \beta \in [0,\pi], \\ \gamma \in [0,\pi/2]} \left( |S_{12}\Omega_{12}| + |S_{23}\Omega_{23}| + |S_{31}\Omega_{31}| \right), \]

where \(\alpha\), \(\beta\), and \(\gamma\) are the three angles of rotation of the coordinate frame. The optimal value is determined by an exhaustive search over all possible values of the angles. The resolution is determined as \( h = \pi / n_i \), with \(n_i\) the number of intervals (the num_intervals variable). For example, num_intervals = 180 corresponds to 1 degree resolution. This optimization problem can become costly for large number of intervals.

Definition at line 38 of file triple_decomposition.h.

◆ valib_gram_schmidt3()

VA_DEVICE_FUN void valib_gram_schmidt3 ( VA_REAL * Q )

#include "linalg3.h"

Gram-Schmidt orthogonalization of a 3x3 matrix.

Definition at line 593 of file linalg3.h.

◆ valib_integrand_resstrain()

VA_DEVICE_FUN void valib_integrand_resstrain	(	VA_REAL *	n,
		VA_REAL *	weight,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	S_RAVG )

#include "corotation_common.h"

Integrand to be used with different integration rules for determining averaged contrarotation in a given integration point with a weight, compute the local contribution and add it to the average contrarotation.

Definition at line 163 of file corotation_common.h.

◆ valib_integrand_resvort()

VA_DEVICE_FUN void valib_integrand_resvort	(	VA_REAL *	n,
		VA_REAL *	weight,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	avgcorot )

#include "corotation_common.h"

Integrand to be used with different integration rules for determining averaged corotation in a given integration point with a weight, compute the local contribution and add it to the average corotation vector.

Definition at line 94 of file corotation_common.h.

◆ valib_integrand_surface()

VA_DEVICE_FUN void valib_integrand_surface	(	VA_REAL *	n,
		VA_REAL *	weight,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	result )

#include "corotation_common.h"

Integrand to be used for debugging purposes. Integrator should return value of the surface of a unit sphere.

Definition at line 269 of file corotation_common.h.

◆ valib_integrator_bazant()

VA_DEVICE_FUN void valib_integrator_bazant	(	VA_DEVICE_ADDR void	integrandVA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	result )

#include "spherical_integrators.h"

Integrator on a unit sphere following [1].

Definition at line 74 of file spherical_integrators.h.

◆ valib_integrator_fibonacci()

VA_DEVICE_FUN void valib_integrator_fibonacci	(	VA_DEVICE_ADDR void	integrandVA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	result )

#include "spherical_integrators.h"

Integrator on a unit sphere using Fibonacci numbers, [7]. Programmed following Appendix A from [5].

Definition at line 20 of file spherical_integrators.h.

◆ valib_integrator_rectangle()

VA_DEVICE_FUN void valib_integrator_rectangle	(	VA_DEVICE_ADDR void	integrandVA_REAL n, VA_REAL weight, VA_DEVICE_ADDR VA_REAL A, VA_DEVICE_ADDR VA_REAL result,
		VA_DEVICE_ADDR VA_REAL *	A,
		VA_DEVICE_ADDR VA_REAL *	result,
		int	numIntervals )

#include "spherical_integrators.h"

Integrator on a unit sphere using the simplest rectangular rule. Corresponds to the approach described in [13].

Definition at line 116 of file spherical_integrators.h.

◆ valib_makec()

VA_DEVICE_FUN VA_COMPLEX valib_makec	(	VA_REAL	a,
		VA_REAL	b )

#include "linalg3.h"

Creates a complex number from two real numbers.

Definition at line 29 of file linalg3.h.

◆ valib_mat_copy3()

VA_DEVICE_FUN void valib_mat_copy3	(	VA_REAL *	A,
		VA_REAL *	B )

#include "linalg3.h"

Copy a 3x3 matrix A to a 3x3 matrix B.

Definition at line 360 of file linalg3.h.

◆ valib_mat_zero3()

VA_DEVICE_FUN void valib_mat_zero3 ( VA_REAL * A )

#include "linalg3.h"

Zero a 3x3 matrix.

Definition at line 346 of file linalg3.h.

◆ valib_matmat3_prod()

VA_DEVICE_FUN void valib_matmat3_prod	(	int	trans_a,
		int	trans_b,
		VA_REAL *	A,
		VA_REAL *	B,
		VA_REAL *	C )

#include "linalg3.h"

Multiplication of two 3x3 matrices.

\[ C = A B + C \]

.

Definition at line 546 of file linalg3.h.

◆ valib_matvec3_prod()

VA_DEVICE_FUN void valib_matvec3_prod	(	int	trans_a,
		VA_REAL *	A,
		VA_REAL *	x,
		VA_REAL *	y )

#include "linalg3.h"

Multiplication of a 3x3 matrix with a vector.

\[ y = A x \]

.

Definition at line 478 of file linalg3.h.

◆ valib_max()

VA_DEVICE_FUN VA_REAL valib_max	(	VA_REAL	a,
		VA_REAL	b )

#include "linalg3.h"

Fortran-like MAX function. Returns the larger of the two numbers.

Definition at line 144 of file linalg3.h.

◆ valib_min()

VA_DEVICE_FUN VA_REAL valib_min	(	VA_REAL	a,
		VA_REAL	b )

#include "linalg3.h"

Fortran-like MIN function. Returns the smaller of the two numbers.

Definition at line 154 of file linalg3.h.

◆ valib_norm_frobenius3()

VA_DEVICE_FUN void valib_norm_frobenius3	(	VA_REAL *	A,
		VA_REAL *	norm )

#include "linalg3.h"

Compute Frobenius norm of a 3x3 matrix.

Definition at line 374 of file linalg3.h.

◆ valib_qcriterion_blending()

VA_DEVICE_FUN void valib_qcriterion_blending	(	VA_REAL *	SO,
		VA_REAL	blending_ratio,
		VA_REAL *	qcriterion_blended )

inline

#include "qcriterion_common.h"

Function common for \(Q\), \(Q_D\), \(Q_W\), and \(Q_M\) criterions with prescribed blending factor \(\beta\) of the norms \( \| \Omega \|_F^2 \) and \( \| S_X \|_F^2\). Evaluates

\[ Q = \frac{1}{2} ( \| \Omega \|_F^2 - \beta \| S_X \|_F^2 ) > 0, \]

where \( \Omega = \frac{1}{2} ( \nabla u - ( \nabla u )^T ) \) is the antisymmetric part and \( S = \frac{1}{2} ( \nabla u + ( \nabla u )^T ) \) is the symmetric part of the velocity gradient \( \nabla u \).

For the \(Q\)-criterion [8] \( S_X = S \) and \(\beta = 1\).

For the \(Q_D\)-criterion [15] \( S_X = S_D \) and \(\beta = 1\).

For the \(Q_M\)-criterion [11] \( S_X = S_D \) and \(\beta = 1.5\).

For the \(Q_W\)-criterion [12] \( S_X = S_D \) and \(\beta = 1.2\).

Here, \( S_D \) is the deviatoric part of \( S \), which is the same as the symmetric part of the deviatoric part of the velocity gradient,

\[( \nabla u )_D = \nabla u - \frac{1}{3} \mbox{trace}(\nabla u)\ I. \]

Definition at line 50 of file qcriterion_common.h.

◆ valib_rank1_update3()

VA_DEVICE_FUN void valib_rank1_update3	(	VA_REAL	alpha,
		VA_REAL *	x,
		VA_REAL *	y,
		VA_REAL *	A )

#include "linalg3.h"

Rank-1 update of a matrix A of dimensions 3x3, computes

\[ A = A + \alpha x y^T \]

.

Definition at line 500 of file linalg3.h.

◆ valib_resvortstrain2D()

VA_DEVICE_FUN void valib_resvortstrain2D	(	VA_REAL *	SO,
		VA_REAL *	residual_vorticity,
		VA_REAL *	residual_strain,
		VA_REAL *	principal_axis )

#include "corotation_common.h"

Function for determining residual vorticity and residual strain in 2D. The 2D case is seen as the leading principal 2x2 block of the 3x3 matrix. The function implements formulas (16)-(27) from [14].

Definition at line 23 of file corotation_common.h.

◆ valib_scalar_prod3()

VA_DEVICE_FUN VA_REAL valib_scalar_prod3	(	VA_REAL *	v1,
		VA_REAL *	v2 )

#include "linalg3.h"

Function for scalar product of two vectors

\[ result = v_1 \cdot v_2 \]

.

Definition at line 313 of file linalg3.h.

◆ valib_second_invariant3()

VA_DEVICE_FUN void valib_second_invariant3	(	VA_REAL *	A,
		VA_REAL *	second_invariant )

#include "linalg3.h"

Compute the second invariant of a 3x3 matrix (2nd tensor invariant).

Definition at line 399 of file linalg3.h.

◆ valib_sign()

VA_DEVICE_FUN VA_REAL valib_sign	(	VA_REAL	a,
		VA_REAL	b )

#include "linalg3.h"

Fortran-like SIGN function. Returns absolute value of a with sign of b.

Definition at line 20 of file linalg3.h.

◆ valib_solve_cubic_equation()

VA_DEVICE_FUN void valib_solve_cubic_equation	(	VA_REAL	a,
		VA_REAL	b,
		VA_REAL	c,
		VA_REAL	d,
		VA_REAL *	x1,
		VA_COMPLEX *	x2,
		VA_COMPLEX *	x3,
		int *	info )

#include "linalg3.h"

Function for solving cubic equation using Cardano formulas based on K. Rektorys: Prehled uzite matematiky, p. 39

\[ a*x^3 + b*x^2 + c*x + d = 0 \]

If discriminant is positive, it finds the three roots, info = 0. If discriminant is negative, it cannot solve the problem, info = -1.

Definition at line 188 of file linalg3.h.

◆ valib_square()

VA_DEVICE_FUN VA_REAL valib_square ( VA_REAL a )

#include "linalg3.h"

Square of a number.

Definition at line 163 of file linalg3.h.

◆ valib_swap_values()

VA_DEVICE_FUN void valib_swap_values	(	VA_REAL *	a,
		VA_REAL *	b )

#include "linalg3.h"

Simple function for swapping two values.

Definition at line 172 of file linalg3.h.

◆ valib_sym_antisym3()

VA_DEVICE_FUN void valib_sym_antisym3 ( VA_REAL * A )

#include "linalg3.h"

Find the symmetric and the antisymmetric parts of a 3x3 matrix. Store the symmetric part in the upper triangle, and the antisymmetric part in the lower triangle.

Definition at line 437 of file linalg3.h.

◆ valib_tatt3()

VA_DEVICE_FUN void valib_tatt3	(	VA_REAL *	A,
		VA_REAL *	T,
		VA_REAL *	TATT )

#include "linalg3.h"

Transformation of a square 3x3 matrix A

\[ TATT = T A T^T \]

.

Definition at line 523 of file linalg3.h.

◆ valib_trace3()

VA_DEVICE_FUN void valib_trace3	(	VA_REAL *	A,
		VA_REAL *	tr3 )

#include "linalg3.h"

Compute the trace of a 3x3 matrix (1st tensor invariant).

Definition at line 390 of file linalg3.h.

◆ valib_transpose3()

VA_DEVICE_FUN void valib_transpose3 ( VA_REAL * A )

#include "linalg3.h"

Transpose A inplace.

Definition at line 424 of file linalg3.h.

◆ valib_triple_decomposition_base()

VA_DEVICE_FUN void valib_triple_decomposition_base	(	VA_DEVICE_ADDR VA_REAL *	A,
		int *	num_intervals,
		VA_DEVICE_ADDR VA_REAL *	residual_vorticity_strain,
		VA_DEVICE_ADDR VA_REAL *	shear_vorticity_strain,
		VA_DEVICE_ADDR VA_REAL *	alpha,
		VA_DEVICE_ADDR VA_REAL *	beta,
		VA_DEVICE_ADDR VA_REAL *	gamma )

#include "triple_decomposition.h"

Triple Decomposition Method (TDM, a.k.a. residual vorticity and residual strain rate) [14], [10].

Evaluates the quantities appearing in the triple decomposition of the velocity gradient following [14],

\[ \nabla u = S_{RES} + \Omega_{RES} + (\nabla u)_{SH} \]

where the symmetric tensor \( S_{RES} \) corresponds to the irrotational straining motion (residual strain rate), \( \Omega_{RES} \) is the antisymmetric tensor of rigid-body rotation (residual vorticity) stored as a vector, and \( (\nabla u)_{SH} \) contains the effective pure shearing motion. In the Basic Reference Frame (BRF), the last tensor has a strictly asymmetric form, i.e.

\[ u_{i,j} = 0\ \mbox{or}\ u_{j,i} = 0\ \mbox{for all}\ i,j. \]

Parameters

[in]	A	3x3 matrix of velocity gradient stored column-wise
[in]	num_intervals	number of steps, i.e. division of the interval \( [0,\pi] \)
[out]	residual_vorticity_strain	jointly the antisymmetric tensor of rigid rotation \( \omega_{RES} \) (lower triangle) and the symmetric tensor of irrotational straining \( S_{RES} \) (upper triangle)
[out]	shear_vorticity_strain	jointly the antisymmetric tensor of rotation due to shear \( \omega_{SH} \) (lower triangle) and the symmetric tensor of straining due to shear \( S_{SH} \) (upper triangle)
[out]	alpha	rotation angle alpha of the basic reference frame (BRF) in radians
[out]	beta	rotation angle beta of the basic reference frame (BRF) in radians
[out]	gamma	rotation angle gamma of the basic reference frame (BRF) in radians

Definition at line 191 of file triple_decomposition.h.

◆ valib_ttat3()

VA_DEVICE_FUN void valib_ttat3	(	VA_REAL *	A,
		VA_REAL *	T,
		VA_REAL *	TTAT )

#include "linalg3.h"

Transformation of a square 3x3 matrix A

\[ TTAT = T^T A T \]

.

Definition at line 575 of file linalg3.h.

◆ valib_vec_zero3()

VA_DEVICE_FUN void valib_vec_zero3 ( VA_REAL * v )

#include "linalg3.h"

Zero a vector of length 3.

Definition at line 301 of file linalg3.h.

Average corotation and contrarotation helper functions

Numerical integrands on a unit sphere

Small helper functions

Operations on vectors of length 3

Operations on 3x3 matrices

Helper functions for the Q-criterion and its variants.

Numerical integrators on a sphere

Triple decomposition helper functions.

Detailed Description

Function Documentation

◆ valib_cabs()

◆ valib_cadd()

◆ valib_cimag()

◆ valib_cmultiply()

◆ valib_constructQ3()

◆ valib_constructQfromN3()

◆ valib_creal()

◆ valib_cross_prod3()

◆ valib_cubic_root()

◆ valib_determinant3()

◆ valib_deviatorise3()

◆ valib_eigenvalues_sym3()

◆ valib_find_nullspace3()

◆ valib_find_pivot3()

◆ valib_get_brf()

◆ valib_gram_schmidt3()

◆ valib_integrand_resstrain()

◆ valib_integrand_resvort()

◆ valib_integrand_surface()

◆ valib_integrator_bazant()

◆ valib_integrator_fibonacci()

◆ valib_integrator_rectangle()

◆ valib_makec()

◆ valib_mat_copy3()

◆ valib_mat_zero3()

◆ valib_matmat3_prod()

◆ valib_matvec3_prod()

◆ valib_max()

◆ valib_min()

◆ valib_norm_frobenius3()

◆ valib_qcriterion_blending()

◆ valib_rank1_update3()

◆ valib_resvortstrain2D()

◆ valib_scalar_prod3()

◆ valib_second_invariant3()

◆ valib_sign()

◆ valib_solve_cubic_equation()

◆ valib_square()

◆ valib_swap_values()

◆ valib_sym_antisym3()

◆ valib_tatt3()

◆ valib_trace3()

◆ valib_transpose3()

◆ valib_triple_decomposition_base()

◆ valib_ttat3()

◆ valib_vec_zero3()