delta.h

Go to the documentation of this file.

 
#ifndef VA_DELTA_H
#define VA_DELTA_H
 
#include "common_defs.h"
#include "linalg3.h"
 
/***************************************************************************/
VA_DEVICE_FUN void valib_delta(VA_REAL *A, VA_REAL *delta)
{
    VA_REAL Q, R;
 
    // find second invariant
    valib_second_invariant3(A, &Q);
 
    // find third invariant (determinant)
    valib_determinant3(A, &R);
 
    // determine delta
    *delta = Q*Q*Q / (VA_REAL) 27.  + R*R / (VA_REAL) 4.;
}
 
#endif // VA_DELTA_H