qmcriterion.h

 
#ifndef VA_QMCRITERION_H
#define VA_QMCRITERION_H
 
#include "qcriterion_common.h"
 
/***************************************************************************/
VA_DEVICE_FUN void valib_qmcriterion(VA_REAL *A, VA_REAL *qm_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 - 3/2 ||S_D||_F^2)
    valib_qcriterion_blending(SO, 1.5, qm_criterion);
}
 
#endif // VA_QMCRITERION_H