qcriterion.h

 
#ifndef VA_QCRITERION_H
#define VA_QCRITERION_H
 
#include "qcriterion_common.h"
 
/***************************************************************************/
VA_DEVICE_FUN void valib_qcriterion(VA_REAL *A, VA_REAL *qcriterion)
{
    VA_REAL SO[9];    // symmetric (upper triangle) and antisymmetric
                      // (lower triangle) parts of velocity gradient
 
    // copy A to SO
    valib_mat_copy3(A, SO);
 
    // evaluate Q = 1/2 (||Omega||_F^2 - ||S||_F^2)
    valib_qcriterion_blending(SO, 1.0, qcriterion);
}
 
#endif // VA_QCRITERION_H