eval_data.hpp

// FEAT3: Finite Element Analysis Toolbox, Version 3
// Copyright (C) 2010 by Stefan Turek & the FEAT group
// FEAT3 is released under the GNU General Public License version 3,
// see the file 'copyright.txt' in the top level directory for details.
 
#pragma once
 
// includes, FEAT
#include <kernel/trafo/base.hpp>
 
namespace FEAT
{
  namespace Trafo
  {
    template<
      typename EvalTraits_,
      TrafoTags cfg_tags_>
    class EvalData
    {
    public:
      typedef EvalTraits_ EvalTraits;
 
      // Note:
      // The following members are ordered by non-ascending size
      // to avoid unnecessary data alignment padding.
 
      typename EvalTraits::HessianInverseType hess_inv;
      typename EvalTraits::HessianTensorType hess_ten;
      typename EvalTraits::JacobianInverseType jac_inv;
      typename EvalTraits::JacobianMatrixType jac_mat;
      typename EvalTraits::NormalVectorType normal;
      typename EvalTraits::ImagePointType img_point;
      typename EvalTraits::DomainPointType dom_point;
      typename EvalTraits::JacobianDeterminantType jac_det;
 
      static constexpr TrafoTags config = cfg_tags_;
 
#if defined(DEBUG) && !defined(__CUDA_ARCH__)
      EvalData()
      {
        // format all values to NaN
        const auto qnan = Math::nan<typename EvalTraits_::DataType>();
        hess_inv = qnan;
        hess_ten = qnan;
        jac_inv = qnan;
        jac_mat = qnan;
        normal = qnan;
        img_point = qnan;
        dom_point = qnan;
        jac_det = qnan;
      }
#endif // DEBUG
 
      void format()
      {
        typedef typename EvalTraits_::DataType DataType;
        hess_inv = DataType(0);
        hess_ten = DataType(0);
        jac_inv = DataType(0);
        jac_mat = DataType(0);
        normal = DataType(0);
        img_point = DataType(0);
        dom_point = DataType(0);
        jac_det = DataType(0);
      }
    }; // class EvalData<...>
  } // namespace Trafo
} // namespace FEAT