SecondPhase.cpp

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_SecondPhase_cpp_
#define model_SecondPhase_cpp_
 
#include <memory>
//#include <LatticePlaneBase.h>
//#include <LatticeVector.h>
//#include <RationalLatticeDirection.h>
#include <SecondPhase.h>
#include <TerminalColors.h>
 
namespace model
{
    
    template<int dim>
    SecondPhase<dim>::SecondPhase(const std::string& _name,
                                  const std::map<const GlidePlaneBase*,std::shared_ptr<GammaSurface>>& _gsMap,
                                  const PlaneNormalContainerType& planes_in) :
    /* init */ name(_name)
    /* init */,gsMap(_gsMap)
    /* init */,planes(planes_in)
    {
        
        std::cout<<greenBoldColor<<"Creating SecondPhase "<<name<<", phaseID= "<<this->sID<<defaultColor<<std::endl;
        
    }
 
    template<int dim>
    double SecondPhase<dim>::misfitEnergy(const VectorDim& b,const GlidePlaneBase* const gpb) const
    {
        const auto gammaIter(gsMap.find(gpb));
        if(gammaIter!=gsMap.end())
        {
            const VectorDim bL(gpb->G2L*b);
            if(std::fabs(bL(2))>FLT_EPSILON)
            {
                throw std::runtime_error("SLIP VECTOR NOT ON GlidePlaneBase");
            }
            return gammaIter->second->operator()(bL.template segment<2>(0));
        }
        else
        {
            return 0.0;
        }
 
    }
 
template<int dim>
double SecondPhase<dim>::misfitEnergy(const VectorDim& b,const size_t& planeID) const
{
    return misfitEnergy(b,planes.at(planeID).get());
}
 
template<int dim>
std::shared_ptr<GammaSurface> SecondPhase<dim>::gammaSurface(const size_t& planeID) const
{
    return gsMap.at(planes.at(planeID).get());
}
    
template struct SecondPhase<3>;
}
#endif