PeriodicBoundaryLagrangianForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "MatrixLagrangianForm.hpp"
 
#include <polyfem/basis/ElementBases.hpp>
#include <polyfem/mesh/LocalBoundary.hpp>
#include <polyfem/mesh/Mesh.hpp>
 
namespace polyfem::solver
{
    class PeriodicBoundaryLagrangianForm : public MatrixLagrangianForm
    {
    public:
        struct Mapping
        {
            StiffnessMatrix A;
            Eigen::MatrixXd b;
            RowVectorNd translation;
            std::vector<int> boundary_dofs;
        };
 
        PeriodicBoundaryLagrangianForm(
            int ndof,
            int value_dim,
            const mesh::Mesh &mesh,
            const std::vector<basis::ElementBases> &bases,
            const std::vector<mesh::LocalBoundary> &local_boundary,
            const std::array<int, 2> &boundary_ids,
            double relative_tolerance);
 
        static Mapping build_mapping(
            int ndof,
            int value_dim,
            const mesh::Mesh &mesh,
            const std::vector<basis::ElementBases> &bases,
            const std::vector<mesh::LocalBoundary> &local_boundary,
            const std::array<int, 2> &boundary_ids,
            double relative_tolerance);
 
        std::string name() const override { return "periodic-boundary-lagrangian"; }
 
    private:
        explicit PeriodicBoundaryLagrangianForm(Mapping mapping);
    };
} // namespace polyfem::solver