MeshedDislocationLoop.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the 
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_MeshedDislocationLoop_h_
#define model_MeshedDislocationLoop_h_
 
 
 
 
#ifndef NDEBUG
#define VerboseMeshedDislocationLoop(N,x) if(verboseMeshedDislocationLoop>=N){std::cout<<x;}
#else
#define VerboseMeshedDislocationLoop(N,x)
#endif
 
#include <deque>
#include <vector>
#include <Plane.h>
#include <DefectiveCrystal/DefectiveCrystalBase.h>
#include <DislocationNetwork.h>
#include <DislocationDynamics/DislocationLoopPatches.h>
 
namespace model
{
    struct MeshedDislocationLoop //: public TriangularMesh
    {
        typedef Eigen::Matrix<double,3,1> VectorDim;
        typedef Eigen::Matrix<double,3,3> MatrixDim;
 
        typedef CompareVectorsByComponent<double,2,float> CompareType;
        std::map<Eigen::Matrix<double,2,1>,int,CompareType> uniquePointsIDs;
 
        
        VectorDim burgers;
        std::vector<VectorDim> periodicShifts;
        std::deque<VectorDim> points;
        std::deque<VectorDim> heightPoints;
        std::deque<VectorDim> displacements;
        std::deque<VectorDim> heightDisplacements;
 
        std::deque<MatrixDim> invVertexMatrix;
        std::deque<MatrixDim> defGradients;
 
        std::deque<Eigen::Vector3i> triangles;
        
    //    MeshedDislocationLoop(const VectorDim& burgers_in,const DefectiveCrystalBase<3>& ddBase,const GlidePlane<3>& plane_in,const std::vector<Eigen::Matrix<double,3,1>>& globalBndPts,const double& meshSize,const double& localMeshSize);
//        MeshedDislocationLoop(const std::shared_ptr<PeriodicPlanePatch<3>>& patch,const std::vector<Eigen::Matrix<double,3,1>>& globalBndPts,const DislocationNetwork<3,0>& DN,const double& meshSize,const double& localMeshSize);
        MeshedDislocationLoop(const VectorDim& burgers_in,const GlidePlane<3>& plane,const std::vector<Eigen::Matrix<double,3,1>>& globalBndPts,const DislocationNetwork<3>& DN,const double& meshSize,const double& localMeshSize);
 
        void update();
        
        VectorDim plasticDisplacementKernel(const Eigen::Ref<const VectorDim>& x) const;
        Eigen::Matrix<double,Eigen::Dynamic,3> plasticDisplacement(Eigen::Ref<const Eigen::Matrix<double,Eigen::Dynamic,3>>) const;
        VectorDim triangleAreaVector(const Eigen::Vector3i&) const;
        double solidAngle(const VectorDim& x,const size_t& triID) const;
        
        template <typename T>
        static int sgn(const T& val)
        {
            return (val > T(0)) - (val < T(0));
        }
 
    };
    
}
#endif