local_refine_tetrahedra_mesh_only_on_boundaries.hpp

Go to the documentation of this file.

// KRATOS  __  __ _____ ____  _   _ ___ _   _  ____
//        |  \/  | ____/ ___|| | | |_ _| \ | |/ ___|
//        | |\/| |  _| \___ \| |_| || ||  \| | |  _
//        | |  | | |___ ___) |  _  || || |\  | |_| |
//        |_|  |_|_____|____/|_| |_|___|_| \_|\____| APPLICATION
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Pablo Becker
//
 
#pragma once
 
// NOTE: Before compute the remeshing it is necessary to compute the neighbours
 
// System includes
 
// External includes
 
// Project includes
#include "utilities/parallel_utilities.h"
#include "utilities/variable_utils.h"
 
// Application includes
#include "custom_utilities/local_refine_tetrahedra_mesh_parallel_to_boundaries.hpp"
 
namespace Kratos
{
class LocalRefineTetrahedraMeshOnlyOnBoundaries : public LocalRefineTetrahedraMeshParallelToBoundaries
{
public:
 
 
    explicit LocalRefineTetrahedraMeshOnlyOnBoundaries(ModelPart& rModelPart)
     : LocalRefineTetrahedraMeshParallelToBoundaries(rModelPart)
    {
    }
 
    ~LocalRefineTetrahedraMeshOnlyOnBoundaries() //TODO maybe {}
    = default;
 
    void SearchEdgeToBeRefined(
        ModelPart& rThisModelPart,
        compressed_matrix<int>& rCoord
        ) override
    {
        KRATOS_TRY;
        for (auto& r_elem: rThisModelPart.Elements()) {
            if (r_elem.GetValue(SPLIT_ELEMENT)) {
                Element::GeometryType& r_geom = r_elem.GetGeometry(); // Nodes of the element
                for (unsigned int i = 0; i < r_geom.size(); i++) {
                    int index_i = mMapNodeIdToPos[r_geom[i].Id()];
                    bool is_boundary_i = r_geom[i].Is(BOUNDARY);
                    for (unsigned int j = 0; j < r_geom.size(); j++) {
                        int index_j = mMapNodeIdToPos[r_geom[j].Id()];
                        bool is_boundary_j = r_geom[j].Is(BOUNDARY);
                        if (index_j > index_i && (is_boundary_j&&is_boundary_i)) {
                            rCoord(index_i, index_j) = -2;
                        }
                    }
                }
            }
        }
 
 
        KRATOS_CATCH("");
    }
 
 
protected:
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
};
 
} // namespace Kratos.