XYCrossFieldSolver.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "MooseTypes.h"
 
#include "libmesh/equation_systems.h"
#include "libmesh/fe_type.h"
#include "libmesh/parallel.h"
#include "libmesh/point_locator_base.h"
#include "libmesh/replicated_mesh.h"
 
#include <complex>
#include <map>
#include <memory>
#include <string>
#include <utility>
#include <vector>
 
namespace libMesh
{
class LinearImplicitSystem;
}
 
class XYCrossFieldSolver
{
public:
  using NodalCrossField = std::map<dof_id_type, std::complex<Real>>;
 
  XYCrossFieldSolver(const libMesh::MeshBase & background_mesh,
                     const std::map<dof_id_type, Real> & boundary_tangent_angles);
 
  void solve();
 
  Real theta(const Point & point) const;
 
  std::pair<Point, Point> crossFrame(const Point & point) const;
 
  const NodalCrossField & nodalCrossField() const { return _nodal_cross_field; }
 
  const std::vector<dof_id_type> & singularNodes() const { return _singular_nodes; }
 
private:
  static void assembleSystem(libMesh::EquationSystems & es, const std::string & system_name);
 
  void assembleLaplace();
 
  void extractNodalCrossField();
 
  std::complex<Real> interpolatedCrossField(const Point & point) const;
 
  libMesh::Parallel::Communicator _communicator;
 
  libMesh::ReplicatedMesh _mesh;
 
  const NodalCrossField _boundary_cross_field;
 
  const libMesh::FEType _fe_type;
 
  std::unique_ptr<libMesh::EquationSystems> _equation_systems;
 
  libMesh::LinearImplicitSystem * _system;
 
  unsigned int _real_variable;
 
  unsigned int _imaginary_variable;
 
  std::unique_ptr<libMesh::PointLocatorBase> _point_locator;
 
  bool _solved;
 
  NodalCrossField _nodal_cross_field;
 
  std::vector<dof_id_type> _singular_nodes;
};