GalerkinClimbSolver.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_GalerkinClimbSolver_h_
#define model_GalerkinClimbSolver_h_
 
//#include <PlanarDislocationNode.h>
 
 
#ifndef NDEBUG
#define VerboseGalerkinClimbSolver(N,x) if(verboseGalerkinClimbSolver>=N){std::cout<<x;}
#else
#define VerboseGalerkinClimbSolver(N,x)
#endif
 
#include <vector>
 
#include <DislocationDynamicsModule.h>
#include <DislocationClimbSolverFactory.h>
 
namespace model
{
    template<int dim> class DislocationNetwork;
    
    template <int dim>
    class GalerkinClimbSolver : public DislocationClimbSolverBase<dim>
    {
        typedef typename DislocationNetwork<dim>::VectorDim   VectorDim;
        typedef typename DislocationNetwork<dim>::MatrixDim   MatrixDim;
        typedef typename DislocationNetwork<dim>::NetworkNodeType   NetworkNodeType;
        typedef typename DislocationNetwork<dim>::NetworkLinkType NetworkLinkType;
        typedef typename DislocationNetwork<dim>::LoopNodeType   LoopNodeType;
        typedef typename DislocationNetwork<dim>::LoopType   LoopType;
        constexpr static int NdofXnode=NetworkNodeType::NdofXnode;
        static constexpr int mSize=ClusterDynamicsParameters<dim>::mSize;     // Cv, Ci, C2i, C3i
        static constexpr int iSize=ClusterDynamicsParameters<dim>::iSize;  // Nc, Na1, Na2, Na3, rc, ra1, ra2, ra3
        typedef Eigen::Matrix<double,2,mSize> ForceVectorMatrixType;
        typedef Eigen::Matrix<double,2*mSize,2> StiffnessMatrixType;
        typedef Eigen::SparseMatrix<double> SparseMatrixType;
        typedef std::deque<Eigen::Triplet<double> > TripletContainerType;
        
        std::vector<Eigen::VectorXd> Fc;
        std::vector<Eigen::VectorXd> KKc;
 
        ForceVectorMatrixType clusterForceVector(const NetworkLinkType& networkLink) const;
        ForceVectorMatrixType clusterForceKernel(const int& k,const NetworkLinkType& networkLink) const;
        StiffnessMatrixType clusterStiffnessKernel(const int& k,const NetworkLinkType& fieldSegment,const NetworkLinkType& sourceSegment) const;
        StiffnessMatrixType clusterStiffnessMatrix(const NetworkLinkType& fieldSegment,const NetworkLinkType& sourceSegment) const;
        
        
        Eigen::VectorXd getNodeVelocitiesBulk() const;
        Eigen::VectorXd getNodeVelocitiesPipe() const;
 
        void computeClimbScalarVelocitiesBulk(const bool& updateStiffness);
 
        
        public:
        
        GalerkinClimbSolver(const DislocationNetwork<dim>&,const ClusterDynamics<dim>* const );
        Eigen::VectorXd getNodeVelocities() const override;
        void computeClimbScalarVelocities(const bool& updateStiffness) override;
        
    };
    
}
#endif