qdcriterion.h

 
#ifndef VA_QDCRITERION_H
#define VA_QDCRITERION_H
 
#include "qcriterion_common.h"
 
/***************************************************************************/
VA_DEVICE_FUN void valib_qdcriterion(VA_REAL *A, VA_REAL *qd_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_D = 1/2 (||Omega||_F^2 - ||S_D||_F^2)
    valib_qcriterion_blending(SO, 1.0, qd_criterion);
}
 
#endif // VA_QDCRITERION_H