ogprojection.h

// This file is part of Hermes2D.
//
// Hermes2D is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 2 of the License, or
// (at your option) any later version.
//
// Hermes2D is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with Hermes2D.  If not, see <http://www.gnu.org/licenses/>.

#ifndef __H2D_OGPROJECTION_H
#define __H2D_OGPROJECTION_H

#include "../function/solution.h"
#include "../forms.h"
#include "../weakform/weakform.h"

namespace Hermes
{
  namespace Hermes2D
  {
    template<typename Scalar>
    class HERMES_API OGProjection : public Hermes::Mixins::Loggable
    {
    public:
      static void project_global(SpaceSharedPtr<Scalar> space,
        MatrixFormVol<Scalar>* custom_projection_jacobian,
        VectorFormVol<Scalar>* custom_projection_residual,
        Scalar* target_vec);

      static void project_global(SpaceSharedPtr<Scalar> space,
        MatrixFormVol<Scalar>* custom_projection_jacobian,
        VectorFormVol<Scalar>* custom_projection_residual,
        MeshFunctionSharedPtr<Scalar> target_sln);

      static void project_global(const std::vector<SpaceSharedPtr<Scalar> > spaces,
        const std::vector<MatrixFormVol<Scalar>*>& custom_projection_jacobian,
        const std::vector<VectorFormVol<Scalar>*>& custom_projection_residual,
        Scalar* target_vec);

      static void project_global(const std::vector<SpaceSharedPtr<Scalar> > spaces,
        const std::vector<MatrixFormVol<Scalar>*>& custom_projection_jacobian,
        const std::vector<VectorFormVol<Scalar>*>& custom_projection_residual,
        const std::vector<MeshFunctionSharedPtr<Scalar> >& target_slns);

      static void project_global(SpaceSharedPtr<Scalar> space, MeshFunctionSharedPtr<Scalar> source_meshfn,
        Scalar* target_vec, NormType proj_norm = HERMES_UNSET_NORM);

      static void project_global(SpaceSharedPtr<Scalar> space, MeshFunctionSharedPtr<Scalar> source_meshfn,
        Hermes::Algebra::Vector<Scalar>* target_vec, NormType proj_norm = HERMES_UNSET_NORM);

      static void project_global(SpaceSharedPtr<Scalar> space,
        MeshFunctionSharedPtr<Scalar> source_sln, MeshFunctionSharedPtr<Scalar>  target_sln,
        NormType proj_norm = HERMES_UNSET_NORM);

      static void project_global(std::vector<SpaceSharedPtr<Scalar> > spaces, std::vector<MeshFunctionSharedPtr<Scalar> > source_meshfns,
        Scalar* target_vec, std::vector<NormType> proj_norms = std::vector<NormType>());

      static void project_global(std::vector<SpaceSharedPtr<Scalar> > spaces, std::vector<MeshFunctionSharedPtr<Scalar> > source_meshfns,
        Hermes::Algebra::Vector<Scalar>* target_vec, std::vector<NormType> proj_norms = std::vector<NormType>());

      static void project_global(std::vector<SpaceSharedPtr<Scalar> > spaces,
        std::vector<MeshFunctionSharedPtr<Scalar> > source_slns, std::vector<MeshFunctionSharedPtr<Scalar> > target_slns,
        std::vector<NormType> proj_norms = std::vector<NormType>(), bool delete_old_mesh = false);

    protected:
      static void project_internal(SpaceSharedPtr<Scalar> space, WeakFormSharedPtr<Scalar> proj_wf, Scalar* target_vec);
    };
  }
}
#endif