discrete_problem_helpers.h

#ifndef __H2D_DISCRETE_PROBLEM_HELPERS_H
#define __H2D_DISCRETE_PROBLEM_HELPERS_H

#include "hermes_common.h"
#include "weakform/weakform.h"

namespace Hermes
{
  namespace Hermes2D
  {
    namespace Mixins
    {
      template<typename Scalar>
      class HERMES_API DiscreteProblemRungeKutta
      {
      protected:
        DiscreteProblemRungeKutta();

        virtual void set_RK(int original_spaces_count, bool force_diagonal_blocks = nullptr, Table* block_weights = nullptr);

        bool rungeKutta;

        int RK_original_spaces_count;

        bool force_diagonal_blocks;

        Table* block_weights;

        // Return scaling coefficient.
        double block_scaling_coeff(MatrixForm<Scalar>* form) const;
        double block_scaling_coeff(MatrixFormDG<Scalar>* form) const;
      };

      template<typename Scalar>
      class HERMES_API DiscreteProblemWeakForm
      {
      public:
        WeakFormSharedPtr<Scalar> get_weak_formulation() const;

      protected:
        DiscreteProblemWeakForm(WeakFormSharedPtr<Scalar> wf = WeakFormSharedPtr<Scalar>(nullptr));

        virtual void set_weak_formulation(WeakFormSharedPtr<Scalar> wf);

        WeakFormSharedPtr<Scalar> wf;
      };

      template<typename Scalar>
      class HERMES_API DiscreteProblemMatrixVector
      {
      protected:
        DiscreteProblemMatrixVector();

        virtual bool set_matrix(SparseMatrix<Scalar>* mat);
        virtual bool set_rhs(Vector<Scalar>* rhs);

        SparseMatrix<Scalar>* current_mat;
        Vector<Scalar>* current_rhs;
      };
    }

    HERMES_API unsigned char init_geometry_points_allocated(RefMap** reference_mapping, unsigned short reference_mapping_count, int order, GeomVol<double>& geometry, double* jacobian_x_weights);
    HERMES_API unsigned char init_geometry_points_allocated(RefMap* rep_reference_mapping, int order, GeomVol<double>& geometry, double* jacobian_x_weights);
    HERMES_API unsigned char init_surface_geometry_points_allocated(RefMap** reference_mapping, unsigned short reference_mapping_count, int& order, unsigned char isurf, int marker, GeomSurf<double>& geometry, double* jacobian_x_weights);
    HERMES_API unsigned char init_surface_geometry_points_allocated(RefMap* rep_reference_mapping, int& order, unsigned char isurf, int marker, GeomSurf<double>& geometry, double* jacobian_x_weights);
  }
}
#endif