node_configure.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Carlos A. Roig
//
 
#pragma once
 
// System includes
#include <string>
#include <iostream>
#include <limits>
#include <cmath>
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "geometries/point.h"
#include "containers/pointer_vector_set.h"
#include "includes/indexed_object.h"
#include "spatial_containers/tree.h"
#include "spatial_containers/cell.h"
 
namespace Kratos {
 
 
 
 
 
 
class NodeConfigure {
public:
    KRATOS_CLASS_POINTER_DEFINITION(NodeConfigure);
 
    static constexpr auto Epsilon   = std::numeric_limits<double>::epsilon();
    static constexpr auto Dimension = 3;
 
    typedef Point                   PointType;
    typedef Node                 ObjectType;
    typedef ObjectType::Pointer     PointerType;
 
    typedef PointerVectorSet<ObjectType,IndexedObject>  ObjectContainerType;
 
    typedef ObjectContainerType::ContainerType          ContainerType;
    typedef ObjectContainerType::ContainerType          ResultContainerType;
    typedef ContactPair<PointerType>                    ContactPairType;
 
    typedef ContainerType::iterator                     IteratorType;
    typedef ResultContainerType::iterator               ResultIteratorType;
    typedef std::vector<double>::iterator               DistanceIteratorType;
 
    typedef double                                      CoordinateType;
    typedef Tvector<CoordinateType,Dimension>           CoordinateArray;
 
 
    NodeConfigure(){};
 
    virtual ~NodeConfigure(){}
 
 
 
    static inline void CalculateBoundingBox(
        const PointerType& rObject, 
        PointType& rLowPoint, 
        PointType& rHighPoint
        ) 
    {
        rHighPoint = rLowPoint = *rObject;
    }
 
    static inline void CalculateBoundingBox(
        const PointerType& rObject, 
        PointType& rLowPoint, 
        PointType& rHighPoint, 
        const double Radius
        ) 
    {
        auto radiusExtension = PointType(Radius, Radius, Radius);
 
        rLowPoint  = PointType{*rObject - radiusExtension};
        rHighPoint = PointType{*rObject + radiusExtension};
    }
 
    static inline void CalculateCenter(
        const PointerType& rObject, 
        PointType& rCentralPoint
        ) 
    {
        rCentralPoint = *rObject;
    }
 
    static inline bool Intersection(
        const PointerType& rObj_1, 
        const PointerType& rObj_2
        ) 
    {
        for(std::size_t i = 0; i < Dimension; i++) {
            if(std::fabs((*rObj_1)[i] - (*rObj_2)[i]) > Epsilon) {
                return false;
            }
        }
 
        return true;
    }
 
    static inline bool Intersection(
        const PointerType& rObj_1, 
        const PointerType& rObj_2, 
        double Radius
    ) 
    {
        for(std::size_t i = 0; i < Dimension; i++) {
            if(std::fabs((*rObj_1)[i] - (*rObj_2)[i]) > Epsilon + Radius) {
                return false;
            }
        }
 
        return true;
    }
 
    static inline bool IntersectionBox(
        const PointerType& rObject, 
        const PointType& rLowPoint, 
        const PointType& rHighPoint
        ) 
    {
        for(std::size_t i = 0; i < Dimension; i++) {
            if( (*rObject)[i] < rLowPoint[i] - Epsilon || (*rObject)[i] > rHighPoint[i] + Epsilon) {
                return false;
            }
        }
 
        return true;
    }
 
    static inline bool IntersectionBox(
        const PointerType& rObject, 
        const PointType& rLowPoint, 
        const PointType& rHighPoint, 
        const double Radius
        ) 
    {
        for(std::size_t i = 0; i < Dimension; i++) {
            if( ((*rObject)[i] + Radius) < rLowPoint[i] - Epsilon || ((*rObject)[i] - Radius) > rHighPoint[i] + Epsilon) {
                return false;
            }
        }
 
        return true;
    }
 
    static inline void Distance(
        const PointerType& rObj_1, 
        const PointerType& rObj_2, 
        double& distance
        ) 
    {
        double pwdDistance = 0.0f;
 
        for(std::size_t i = 0; i < Dimension; i++) {
            pwdDistance += std::pow((*rObj_1)[i] - (*rObj_2)[i], 2);
        }
 
        distance = std::sqrt(pwdDistance);
    }
 
    static inline double GetObjectRadius(
        const PointerType& rObject, 
        const double Radius
        ) 
    {
        return 0.0f;
    }
 
 
 
 
    virtual std::string Info() const {
        return "Spatial Containers Configure for 'Points'";
    }
 
    virtual std::string Data() const {
        return "Dimension: " + std::to_string(Dimension);
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const {
        rOStream << Info() << std::endl;
    }
 
    virtual void PrintData(std::ostream& rOStream) const {
        rOStream << Data() << Dimension << std::endl;
    }
 
 
 
protected:
 
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
    NodeConfigure& operator=(NodeConfigure const& rOther);
 
    NodeConfigure(NodeConfigure const& rOther);
 
 
}; // Class NodeConfigure
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream, NodeConfigure& rThis){
    return rIStream;
}
 
inline std::ostream& operator << (std::ostream& rOStream, const NodeConfigure& rThis){
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
  return rOStream;
}
 
 
} // namespace Kratos.