shared_points_mapper.h

Go to the documentation of this file.

//
//   Project Name:        Kratos
//   Last Modified by:    $Author: pooyan $
//   Date:                $Date: 2006-11-27 16:07:42 $
//   Revision:            $Revision: 1.1.1.1 $
//
//
 
 
#if !defined(KRATOS_SHARED_POINTS_MAPPER_H_INCLUDED )
#define  KRATOS_SHARED_POINTS_MAPPER_H_INCLUDED
 
 
 
// System includes
#include <set>
 
 
// External includes
#include "boost/smart_ptr.hpp"
 
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "containers/pointer_vector.h"
#include "includes/node.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
class SharedPointsMapper
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(SharedPointsMapper);
 
 
    //**************************************************************************************************
    //**************************************************************************************************
    SharedPointsMapper(
        const ModelPart::NodesContainerType& OriginNodes,
        const ModelPart::NodesContainerType& DestinationNodes,
        double tol = 1e-9)
    {
        KRATOS_TRY
 
        if( OriginNodes.size()!=0 && DestinationNodes.size()!=0)
        {
            //if( OriginNodes.size() != DestinationNodes.size() )
            //    KRATOS_THROW_ERROR(std::logic_error,"wrong number of nodes","");
 
            mOriginNodes.reserve( OriginNodes.size() );
            mDestinationNodes.reserve( DestinationNodes.size() );
 
            for(ModelPart::NodesContainerType::const_iterator origin_it = OriginNodes.begin(); origin_it != OriginNodes.end(); origin_it++)
            {
                for(ModelPart::NodesContainerType::const_iterator destination_it = DestinationNodes.begin(); destination_it != DestinationNodes.end(); destination_it++)
                {
                    if  (
                        fabs(origin_it->X() - destination_it->X() ) < tol  &&
                        fabs(origin_it->Y() - destination_it->Y() ) < tol  &&
                        fabs(origin_it->Z() - destination_it->Z() ) < tol
                    )
                    {
                        mOriginNodes.push_back(  *(origin_it.base() ) );
                        mDestinationNodes.push_back(  *(destination_it.base() ) );
                    }
                }
 
            }
        }
 
        KRATOS_CATCH("")
    }
 
 
    virtual ~SharedPointsMapper() {}
 
 
 
 
    //**************************************************************************************************
    //**************************************************************************************************
    void ScalarMap( const Variable<double>& rOriginVariable,
                    const Variable<double>& rDestinationVariable)
    {
        KRATOS_TRY
 
        PointerVector< Node >::iterator it_origin = mOriginNodes.begin();
        PointerVector< Node >::iterator it_destination = mDestinationNodes.begin();
 
        for(unsigned int i = 0 ; i < mOriginNodes.size() ; ++i)
        {
            (it_destination++ )->FastGetSolutionStepValue(rDestinationVariable) =
                (it_origin++ )->FastGetSolutionStepValue(rOriginVariable);
 
        }
 
        KRATOS_CATCH("")
    }
 
    void InverseScalarMap( const Variable<double>& rOriginVariable,
                    const Variable<double>& rDestinationVariable)
    {
        KRATOS_TRY
 
        PointerVector< Node >::iterator it_origin = mOriginNodes.begin();
        PointerVector< Node >::iterator it_destination = mDestinationNodes.begin();
 
        for(unsigned int i = 0 ; i < mOriginNodes.size() ; ++i)
        {
            (it_origin++ )->FastGetSolutionStepValue(rOriginVariable) =
                (it_destination++ )->FastGetSolutionStepValue(rDestinationVariable);
 
        }
 
        KRATOS_CATCH("")
    }
 
    //**************************************************************************************************
    //**************************************************************************************************
    void VectorMap( const Variable<array_1d<double,3> >& rOriginVariable,
                    const Variable<array_1d<double,3> >& rDestinationVariable)
    {
        KRATOS_TRY
 
        PointerVector< Node >::iterator it_origin = mOriginNodes.begin();
        PointerVector< Node >::iterator it_destination = mDestinationNodes.begin();
 
        for(unsigned int i = 0 ; i < mOriginNodes.size() ; ++i)
        {
            (it_destination++ )->FastGetSolutionStepValue(rDestinationVariable) =
                (it_origin++ )->FastGetSolutionStepValue(rOriginVariable);
 
        }
 
        KRATOS_CATCH("")
    }
 
    //**************************************************************************************************
    //**************************************************************************************************
    void InverseVectorMap( const Variable<array_1d<double,3> >& rOriginVariable,
                    const Variable<array_1d<double,3> >& rDestinationVariable)
    {
        KRATOS_TRY
 
        PointerVector< Node >::iterator it_origin = mOriginNodes.begin();
        PointerVector< Node >::iterator it_destination = mDestinationNodes.begin();
 
        for(unsigned int i = 0 ; i < mOriginNodes.size() ; ++i)
        {
            noalias((it_origin++ )->FastGetSolutionStepValue(rOriginVariable)) =
                 (it_destination++ )->FastGetSolutionStepValue(rDestinationVariable);
 
 
        }
 
        KRATOS_CATCH("")
    }
 
 
 
 
 
 
 
 
    virtual void PrintInfo(std::ostream& OStream) const
    {
    }
 
 
    virtual void PrintData(std::ostream& OStream) const
    {
    }
 
 
 
 
 
protected:
 
 
    PointerVector< Node > mOriginNodes;
 
    PointerVector< Node > mDestinationNodes;
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
 
 
 
 
 
 
    SharedPointsMapper& operator=(const SharedPointsMapper& Other);
 
    //SharedPointsMapper(const SharedPointsMapper& Other);
 
 
 
}; // Class SharedPointsMapper
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SHARED_POINTS_MAPPER_H_INCLUDED  defined