TriToQuadConverter.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 "MeshGenerator.h"
#include "MooseEnum.h"
 
#include <array>
#include <vector>
 
class TriToQuadConverter : public MeshGenerator
{
public:
  static InputParameters validParams();
 
  TriToQuadConverter(const InputParameters & parameters);
 
  std::unique_ptr<MeshBase> generate() override;
 
  using QuadCorners = std::array<unsigned int, 4>;
 
  struct RecombineCandidate
  {
    Real eta;
    dof_id_type first_elem_id;
    dof_id_type second_elem_id;
    std::array<dof_id_type, 4> quad_node_ids;
  };
 
  static Real quadQuality(const std::array<Point, 4> & quad_points);
 
  static std::vector<RecombineCandidate>
  greedyMatching(std::vector<RecombineCandidate> & candidates, const Real eta_min);
 
  static std::vector<QuadCorners> triSubdivisionTemplate();
 
  static std::vector<QuadCorners> quadSubdivisionTemplate();
 
protected:
  std::unique_ptr<MeshBase> & _input;
  const MooseEnum _algorithm;
  const Real _eta_min;
  const SubdomainName _tri_subdomain_name_suffix;
  const bool _all_quad;
 
private:
  static RecombineCandidate
  buildCandidate(const Elem & elem, const unsigned int side, const Elem & neighbor);
 
  void subdivide(ReplicatedMesh & mesh) const;
 
  void recombine(ReplicatedMesh & mesh) const;
 
  void moveSurvivingTriangles(ReplicatedMesh & mesh,
                              const subdomain_id_type scratch_subdomain_id) const;
};