FieldPoint.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MODEL is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef _model_FieldPoint_h
#define _model_FieldPoint_h
 
#include <SpatialCellObserver.h>
#include <FieldPointBase.h>
 
 
namespace model {
    
    template<typename Derived,
    /*    */ short unsigned int _dim,
    /*    */ typename ...FieldTypes>
    struct FieldPoint;
    
    /**************************************************************************/
    /**************************************************************************/
    // template specialization: one or more fields
    template<typename Derived,
    /*    */ short unsigned int _dim,
    /*    */ typename FieldType,
    /*    */ typename ...MoreFieldTypes>
    struct FieldPoint<Derived,_dim,FieldType,MoreFieldTypes...> :
    /* inheritance  */ public FieldPoint<Derived,_dim,MoreFieldTypes...>,
    /* inheritance  */ public FieldPointBase<Derived,FieldType>
    {
        
        typedef typename FieldPoint<Derived,_dim,MoreFieldTypes...>::VectorDimD VectorDimD;
        
        /**********************************************************************/
        template<typename...T>
        FieldPoint(const bool& enb, const T&...moreEnb) :
        /* base init */ FieldPoint<Derived,_dim,MoreFieldTypes...>(moreEnb...),
        /* base init */ FieldPointBase<Derived,FieldType>(enb)
        {
        }
        
//        /**********************************************************************/
//        template <typename OtherFieldType>
//        FieldPointBase<Derived,OtherFieldType>& field()
//        {
//            return *static_cast<FieldPointBase<Derived,OtherFieldType>* const>(this);
//        }
        
//        /**********************************************************************/
//        template <typename OtherFieldType>
//        const FieldPointBase<Derived,OtherFieldType>& field() const
//        {
//            return *static_cast<const FieldPointBase<Derived,OtherFieldType>* const>(this);
//        }
        
        /**********************************************************************/
        template <typename OtherFieldType>
        typename OtherFieldType::MatrixType field() const
        {
            return OtherFieldType::get(*static_cast<const FieldPointBase<Derived,OtherFieldType>* const>(this));
        }
        
        /**********************************************************************/
        template <typename OtherFieldType>
        FieldPointBase<Derived,OtherFieldType>& fieldPointBase()
        {
            return *static_cast<FieldPointBase<Derived,OtherFieldType>* const>(this);
        }
        
        /**********************************************************************/
        template <typename OtherFieldType>
        const FieldPointBase<Derived,OtherFieldType>& fieldPointBase() const
        {
            return *static_cast<const FieldPointBase<Derived,OtherFieldType>* const>(this);
        }
        
        
#ifdef _MODEL_MPI_
        
        typedef typename FieldPointBase<Derived,FieldType>::EigenMapType EigenMapType;
        
        void set_mpiID(const size_t& k)
        {
            FieldPoint<Derived,_dim,MoreFieldTypes...>::set_mpiID(k);
            FieldPointBase<Derived,FieldType>::set_mpiID(k);
        }
#endif
        
    };
    
    /**************************************************************************/
    /**************************************************************************/
    // template specialization: zero fields
    template<typename Derived,
    /*    */ short unsigned int _dim>
    struct FieldPoint<Derived,_dim>
    {
        
        typedef SpatialCellObserver<Derived,_dim> SpatialCellObserverType;
        typedef typename SpatialCellObserverType::VectorDimD VectorDimD;
 
#ifdef _MODEL_MPI_
    private:
        size_t _mpiID;
#endif
 
    public:
        
#ifdef _MODEL_MPI_
        /**********************************************************************/
        const size_t& mpiID() const
        {
            return _mpiID;
        }
        
        /**********************************************************************/
        void set_mpiID(const size_t& k)
        {
            _mpiID=k;
        }
#endif
        
        /**********************************************************************/
        template <typename SourcePointType>
        typename SpatialCellObserver<SourcePointType,_dim>::CellMapType neighborCells() const
        {
            return SpatialCellObserver<SourcePointType,_dim>::neighborCells(static_cast<const Derived*>(this)->P);
        }
        
        /**********************************************************************/
        template <typename SourcePointType>
        typename SpatialCellObserver<SourcePointType,_dim>::CellMapType farCells() const
        {
            return SpatialCellObserver<SourcePointType,_dim>::farCells(static_cast<const Derived*>(this)->P);
        }
        
        
    };
        
} // end namespace
#endif