static void flux_jacobian_eigen_values(const MODEL::Properties& p, const GV& direction, EV& Dv, OP& op ) { Dv = p.v * direction; Dv = Dv.unaryExpr( op ); }