discrete_problem_linear.h

#ifndef __H2D_DISCRETE_PROBLEM_LINEAR_H
#define __H2D_DISCRETE_PROBLEM_LINEAR_H

#include "discrete_problem.h"

namespace Hermes
{
  namespace Hermes2D
  {
    template<typename Scalar>
    class HERMES_API DiscreteProblemLinear : public DiscreteProblem<Scalar>
    {
    public:
      DiscreteProblemLinear(const WeakForm<Scalar>* wf, Hermes::vector<const Space<Scalar> *> spaces);

      DiscreteProblemLinear(const WeakForm<Scalar>* wf, const Space<Scalar>* space);

      DiscreteProblemLinear();

      virtual ~DiscreteProblemLinear();

      virtual void assemble(SparseMatrix<Scalar>* mat, Vector<Scalar>* rhs = NULL, bool force_diagonal_blocks = false,
        Table* block_weights = NULL);

    protected:
      virtual void assemble_matrix_form(MatrixForm<Scalar>* form, int order, Func<double>** base_fns, Func<double>** test_fns, Func<Scalar>** ext, Func<Scalar>** u_ext,
      AsmList<Scalar>* current_als_i, AsmList<Scalar>* current_als_j, Traverse::State* current_state, int n_quadrature_points, Geom<double>* geometry, double* jacobian_x_weights);

      template<typename T> friend class KellyTypeAdapt;
      template<typename T> friend class NewtonSolver;
      template<typename T> friend class PicardSolver;
      template<typename T> friend class RungeKutta;
      template<typename T> friend class LinearSolver;
    };
  }
}
#endif