#pragma once

#include <type_traits>

#include <amdis/assembler/FirstOrderTestvecGradTrial.hpp>

namespace AMDiS
{
  /**
   * \addtogroup operators
   * @{
   **/

  namespace tag
  {
    struct gradtest_trialvec {};
  }


  /// first-order operator \f$ \langle\nabla\psi, c\,\Phi\rangle \f$
  template <class GridFct, class QuadCreator>
  class GridFunctionOperator<tag::gradtest_trialvec, GridFct, QuadCreator>
      : public GridFunctionOperator<tag::testvec_gradtrial, GridFct, QuadCreator>
  {
    using Transposed = GridFunctionOperator<tag::testvec_gradtrial, GridFct, QuadCreator>;

  public:
    GridFunctionOperator(tag::gradtest_trialvec, GridFct const& expr, QuadCreator const& quadCreator)
      : Transposed(tag::testvec_gradtrial{}, expr, quadCreator)
    {}

    template <class Context, class QuadratureRule,
              class ElementMatrix, class RowNode, class ColNode, bool sameFE, bool sameNode>
    void calculateElementMatrix(Context const& context,
                                QuadratureRule const& quad,
                                ElementMatrix& elementMatrix,
                                RowNode const& rowNode,
                                ColNode const& colNode,
                                std::integral_constant<bool, sameFE> flagSameFE,
                                std::integral_constant<bool, sameNode> flagSameNode)
    {
      auto elementMatrixTransposed = trans(elementMatrix);
      Transposed::calculateElementMatrix(
        context, quad, elementMatrixTransposed, colNode, rowNode, flagSameFE, flagSameNode);
    }
  };

  /** @} **/

} // end namespace AMDiS