ITHACAgeometry.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
     ██╗████████╗██╗  ██╗ █████╗  ██████╗ █████╗       ███████╗██╗   ██╗
     ██║╚══██╔══╝██║  ██║██╔══██╗██╔════╝██╔══██╗      ██╔════╝██║   ██║
     ██║   ██║   ███████║███████║██║     ███████║█████╗█████╗  ██║   ██║
     ██║   ██║   ██╔══██║██╔══██║██║     ██╔══██║╚════╝██╔══╝  ╚██╗ ██╔╝
     ██║   ██║   ██║  ██║██║  ██║╚██████╗██║  ██║      ██║      ╚████╔╝
     ╚═╝   ╚═╝   ╚═╝  ╚═╝╚═╝  ╚═╝ ╚═════╝╚═╝  ╚═╝      ╚═╝       ╚═══╝
 
 * In real Time Highly Advanced Computational Applications for Finite Volumes
 * Copyright (C) 2017 by the ITHACA-FV authors
-------------------------------------------------------------------------------
License
    This file is part of ITHACA-FV
    ITHACA-FV is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
    ITHACA-FV is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
    GNU Lesser General Public License for more details.
    You should have received a copy of the GNU Lesser General Public License
    along with ITHACA-FV. If not, see <http://www.gnu.org/licenses/>.
File
    ITHACAgeometry
Description
    set of methods for geometrical parametrization
SourceFiles
    ITHACAgeometry.C
\*---------------------------------------------------------------------------*/
 
 
#ifndef ITHACAgeometry_H
#define ITHACAgeometry_H
 
#include "fvCFD.H"
#include "IOmanip.H"
#include "freestreamFvPatchField.H"
#include <sys/stat.h>
#include <unistd.h>
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wold-style-cast"
#include <Eigen/Eigen>
#pragma GCC diagnostic pop
#include <functional>
#include "./colormod.H"
#include "polyMeshTools.H"
#include <chrono>
#include "mixedFvPatchFields.H"
#include "fvMeshSubset.H"
using namespace std::placeholders;
#include "Foam2Eigen.H"
 
namespace ITHACAutilities
{
labelList getIndicesFromBox(const fvMesh& mesh, List<label> indices,
                            Eigen::MatrixXd Box, List<vector>& points2Move);
 
//--------------------------------------------------------------------------
List<label> getIndices(const fvMesh& mesh, int index, int layers);
 
//--------------------------------------------------------------------------
List<label> getIndices(const fvMesh& mesh, int index_row, int index_col,
                       int layers);
 
//--------------------------------------------------------------------------
void getPointsFromPatch(fvMesh& mesh, label ind, List<vector>& points,
                        labelList& indices);
 
//--------------------------------------------------------------------------
List<vector> displacedSegment(List<vector> x0, double mux1, double mux2,
                              double muy1, double muy2);
//----------------------------------------------------------------------
vector displacePolabel(vector x0, vector x_low, vector x_up,
                       double mux_low, double mux_up,
                       double muy_low, double muy_up, double muz_low = 0.0, double muz_up = 0.0);
 
//--------------------------------------------------------------------------
Field<vector> rotateMesh(fvMesh& mesh, double r1, double r2, vector axis,
                         double alpha, labelList movingPointsIDs,
                         List<double> radii, word angleVariationMethod = "Linear", double v = 3);
 
//--------------------------------------------------------------------------
Eigen::MatrixXd rotationMatrix(vector AxisOfRotation,
                               double AngleOfRotation);
//--------------------------------------------------------------------------
Eigen::VectorXd boudaryFaceToCellDistance(fvMesh& mesh, label BC_ind);
 
//--------------------------------------------------------------------------
List<label> getIndicesFromDisc(const fvMesh& mesh, double radius,
                               vector origin, vector axis, List<double>& radii);
 
//--------------------------------------------------------------------------
template<typename type_f>
List<label> getIndicesFromBox(
    GeometricField<type_f, fvPatchField, volMesh>& field, Eigen::MatrixXd Box);
 
//--------------------------------------------------------------------------
template<typename type_f>
fvMeshSubset* getSubMeshFromBox(
    GeometricField<type_f, fvPatchField, volMesh>& field, Eigen::MatrixXd Box);
 
//--------------------------------------------------------------------------
volScalarField meshNonOrtho(fvMesh& mesh, volScalarField& NonOrtho);
 
 
//--------------------------------------------------------------------------
List<vector> rotatePoints(const List<vector>& originalPoints,
                          vector AxisOfRotation, double AngleOfRotation);
 
}
 
#endif