round.h

Go to the documentation of this file.

#ifndef _RHEOLEF_ROUND_H
#define _RHEOLEF_ROUND_H
#include "rheolef/rounder.h"
#include "rheolef/field_expr.h"
 
namespace rheolef {
 
template<class Expr, class T2>
typename
std::enable_if<
  details::is_field_expr_affine_homogeneous<Expr>::value
 ,field_basic<
    typename Expr::scalar_type
   ,typename Expr::memory_type
  >
>::type
round (
 const Expr& expr,
 const T2&   prec)
{
  typedef typename Expr::scalar_type T;
  typedef typename Expr::memory_type M;
  space_basic<T,M> Xh;
  check_macro (expr.have_homogeneous_space (Xh),
    "round(field_expr, prec): field_expr should have homogeneous space. HINT: use round(interpolate(Xh, field_expr), prec)");
  field_basic<T,M> uh = expr; // create a tmp because expr has no expr.end_dof()
  field_basic<T,M> vh (Xh);
  std::transform (uh.begin_dof(), uh.end_dof(), vh.begin_dof(), rounder(prec));
  return vh;
}
 
}//namespace rheolef
#endif // _RHEOLEF_ROUND_H