qwcriterion.h

 
#ifndef VA_QWCRITERION_H
#define VA_QWCRITERION_H
 
#include "qcriterion_common.h"
 
/***************************************************************************/
VA_DEVICE_FUN void valib_qwcriterion(VA_REAL *A, VA_REAL *qw_criterion)
{
    VA_REAL SO[9];    // symmetric (upper triangle) and antisymmetric
                      // (lower triangle) parts of velocity gradient
 
    // copy A to SO
    valib_mat_copy3(A, SO);
 
    // find the deviatoric part of A
    valib_deviatorise3(SO);
 
    // evaluate Q_M = 1/2 (||Omega||_F^2 - 6/5 ||S_D||_F^2)
    valib_qcriterion_blending(SO, 1.2, qw_criterion);
}
 
#endif // VA_QWCRITERION_H