CellDivider.h

Go to the documentation of this file.

#ifndef CELL_DIVIDER_H_INCLUDED
#define CELL_DIVIDER_H_INCLUDED
/*
 * Copyright 2011-2016 Universiteit Antwerpen
 *
 * Licensed under the EUPL, Version 1.1 or  as soon they will be approved by
 * the European Commission - subsequent versions of the EUPL (the "Licence");
 * You may not use this work except in compliance with the Licence.
 * You may obtain a copy of the Licence at: http://ec.europa.eu/idabc/eupl5
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the Licence is distributed on an "AS IS" basis,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the Licence for the specific language governing
 * permissions and limitations under the Licence.
 */
#include "model/ComponentInterfaces.h"
#include "sim/CoreData.h"

#include <array>
#include <functional>
#include <list>
#include <memory>
#include <vector>

namespace SimPT_Sim {

class Cell;
class Edge;
class Mesh;
class Node;
class Sim;
class Wall;
using boost::property_tree::ptree;

class CellDivider
{
public:
        CellDivider(const CoreData& cd);

        CellDivider(Mesh* mesh);

        unsigned int DivideCells();

        void DivideOverAxis(Cell* cell, const std::array<double, 3>& axis);

        std::list<Cell*> GeometricSplitCell(Cell* cell, Node* node1, Node* node2);

        void Initialize(Mesh* mesh);

        void Initialize(const CoreData& cd);

private:
        void AddSharedWallNodes(std::array<std::array<double, 3>, 2>& new_node_vector,
                std::array<Node*, 2>& new_node_index,
                std::vector<Node*>& parent_nodes, std::vector<Node*>& daughter_nodes);

        void AddSplitWalls(Cell* this_cell, Cell* neighbor_cell, int i,
                std::array<Node*, 2>& new_node_index,
                std::array<Wall*, 4>& div_wall, std::array<double, 2>& orig_length,
                std::array<double, 2>& orig_rest_length,
        std::array<double, 2>& orig_rest_length_init);

        void ComputeIntersection(Cell* cell, std::vector<Node*> intersect_indexes,
                const std::array<double, 3>& from, const std::array<double, 3>& to,
                std::array<std::array<double, 3>, 2>& new_node_vector,
                std::array<Node*, 2>& new_node_index, std::array<int, 2>& new_node_flag,
                std::array<Edge, 2>& divided_edges);

        Node* InsertNewNodeAtIntersect(Cell* this_cell, Cell* neighbor_cell,
                std::array<double, 3> node_vec, const Edge& intrsct_edge);

        Cell* SplitCell(Cell* this_cell, std::vector<Node*> intersect_indexes,
                const std::array<double, 3>& from, const std::array<double, 3>& to);

private:
        CoreData                      m_cd;                        
        CellDaughtersComponent        m_cell_daughters;            
        CellSplitComponent            m_cell_split;                
        double                        m_collapse_node_threshold;   
        bool                          m_copy_wall;                 
        Mesh*                         m_mesh;                      
        std::string                   m_model_name;                
        double                        m_target_node_distance;      
        std::function<double()>       m_uniform_generator;         
};

} // namespace

#endif // end_of_include_guard