parallel_unstructured_interface.h

Go to the documentation of this file.

// ------------------------------------------------------
//
// SPDX-License-Identifier: GPL-2.0-or-later
//
// SPDX-FileCopyrightText: Copyright (C) 2011-2026 by the ASPECT authors.
//
// This file is part of ASPECT.
//
// Detailed license information governing the source code and
// contributions can be found in the folder LICENSES and in
// LICENSE and CONTRIBUTING.md at the top level directory.
//
// ------------------------------------------------------


#ifndef _aspect_mesh_deformation_parallel_unstructured_interface_h
#define _aspect_mesh_deformation_parallel_unstructured_interface_h

#include <aspect/mesh_deformation/interface.h>


namespace aspect
{
  namespace MeshDeformation
  {
    template <int dim>
    class ParallelUnstructuredInterface : public Interface<dim>, public SimulatorAccess<dim>
    {
      public:
        virtual
        void
        compute_velocity_constraints_on_boundary(const DoFHandler<dim> &mesh_deformation_dof_handler,
                                                 AffineConstraints<double> &mesh_velocity_constraints,
                                                 const std::set<types::boundary_id> &boundary_ids) const override;

        virtual
        std::vector<Tensor<1,dim>>
        compute_updated_velocities_at_points (const std::vector<std::vector<double>> &current_solution_at_points) const = 0;


      protected:
        // Helper functions to be used in derived classes:

        void
        set_evaluation_points (const std::vector<Point<dim>> &evaluation_points);

        std::vector<std::vector<double>>
        evaluate_aspect_solution_at_points () const;

        LinearAlgebra::Vector
        interpolate_external_vector_field_to_surface_support_points (const std::vector<Tensor<1,dim>> &velocities) const;

        std::vector<Point<dim>> evaluation_points;

        std::unique_ptr<Utilities::MPI::RemotePointEvaluation<dim, dim>> remote_point_evaluator;

        struct DofToEvalPointData
        {
          types::global_dof_index dof_index;
          unsigned int            evaluation_point_rank;
          unsigned int            evaluation_point_index;
          unsigned int            component;
          double                  squared_distance;

          template <class Archive>
          void
          serialize(Archive &ar, const unsigned int /*version*/)
          {
            ar &dof_index;
            ar &evaluation_point_rank;
            ar &evaluation_point_index;
            ar &component;
            ar &squared_distance;
          }
        };

        std::vector<DofToEvalPointData> map_dof_to_eval_point;
    };
  }
}

#endif