qcriterion_common.h

Go to the documentation of this file.

 
#ifndef VA_QCRITERION_COMMON_H
#define VA_QCRITERION_COMMON_H
 
#include "common_defs.h"
#include "linalg3.h"
 
/***************************************************************************/
/***************************************************************************/
inline VA_DEVICE_FUN void valib_qcriterion_blending(VA_REAL *SO,
                                                    VA_REAL blending_ratio,
                                                    VA_REAL *qcriterion_blended)
{
    // prepare symmetric and antisymmetric parts of SO
    valib_sym_antisym3(SO);
 
    // evaluate Q = 1/2 (||Omega||_F^2 - ||S||_F^2)
    *qcriterion_blended = 0.5
                        * (2.0*( SO[1]*SO[1]         // )
                                +SO[2]*SO[2]         // } Omega^2
                                +SO[5]*SO[5])        // )
                           - blending_ratio *
                             ( SO[0]*SO[0]           // )
                              +SO[4]*SO[4]           // } diagonal of S^2
                              +SO[8]*SO[8]           // )
                              +2.0*( SO[3]*SO[3]     // )
                                    + SO[6]*SO[6]    // } off-diagonal of S^2
                                    + SO[7]*SO[7])));// )
}
 
/***************************************************************************/
#endif // VA_QCRITERION_COMMON_H