AxialGeoMultiscaleMapping.hpp

Go to the documentation of this file.

#pragma once
 
#include <vector>
#include "logging/Logger.hpp"
#include "mapping/Mapping.hpp"
 
namespace precice::mapping {

class AxialGeoMultiscaleMapping : public Mapping {
public:
  enum struct MultiscaleType {
    SPREAD,
    COLLECT
  };
  enum struct MultiscaleAxis {
    X = 0,
    Y = 1,
    Z = 2
  };

  enum struct MultiscaleProfile {
    UNIFORM,
    PARABOLIC
  };
  enum struct MultiscaleDimension {
    D1D3,
    D1D2,
    D2D3
  };
  enum struct MultiscaleCrossSection {
    CIRCLE,
    SQUARE
  };

  AxialGeoMultiscaleMapping(Constraint constraint, int dimensions, MultiscaleDimension dimension, MultiscaleType type, MultiscaleAxis axis, double radius, MultiscaleProfile profile = MultiscaleProfile::UNIFORM, MultiscaleCrossSection crossSection = MultiscaleCrossSection::CIRCLE);

  void computeMapping() override;

  void clear() override;
 
  void tagMeshFirstRound() override;
  void tagMeshSecondRound() override;

  std::string getName() const final override;
 
protected:
  void mapConservative(const time::Sample &inData, Eigen::VectorXd &outData) override;

  void mapConsistent(const time::Sample &inData, Eigen::VectorXd &outData) override;
 
private:
  mutable logging::Logger _log{"mapping::AxialGeoMultiscaleMapping"};
 
  // dimensionality of mapping, namely 1D-3D, 1D-2D or 2D-3D
  MultiscaleDimension _dimension;

  MultiscaleType _type;

  MultiscaleAxis _axis;

  double _radius;

  MultiscaleProfile _profile;
 
  // cross-section of the pipe
  MultiscaleCrossSection _crossSection;

  std::vector<double> _vertexDistances;

  std::vector<Eigen::Vector2d> _vertexTransverseCoords;
 
  // Axis of 2D interface line
  int _lineCoord = -1;
 
  // nearest input vertex to output vertices
  std::vector<int> _nearestVertex;
 
  // band of nearest output vertices to input vertex
  std::vector<std::vector<int>> _collectBands;
 
  // Maximum distance of output vertices to each input vertex
  std::vector<double> _maxDistancePerInput;
};
 
} // namespace precice::mapping