Example #1
0
  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 );
  }