clear_contact_conditions_mesher_process.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosContactMechanicsApplication $
//   Created by:          $Author:              JMCarbonell $
//   Last modified by:    $Co-Author:                       $
//   Date:                $Date:                  July 2016 $
//   Revision:            $Revision:                    0.0 $
//
//
 
#if !defined(KRATOS_CLEAR_CONTACT_CONDITIONS_MESHER_PROCESS_H_INCLUDED )
#define  KRATOS_CLEAR_CONTACT_CONDITIONS_MESHER_PROCESS_H_INCLUDED
 
 
// System includes
 
// External includes
 
// Project includes
 
#include "includes/model_part.h"
#include "custom_utilities/mesher_utilities.hpp"
#include "custom_processes/mesher_process.hpp"
 
//Data:
//StepData:
//Flags:    (checked)
//          (set)
//          (modified)
//          (reset)
 
 
namespace Kratos
{
 
 
  typedef  ModelPart::NodesContainerType NodesContainerType;
  typedef  ModelPart::ElementsContainerType ElementsContainerType;
  typedef  ModelPart::ConditionsContainerType ConditionsContainerType;
 
 
 
 
 
  class ClearContactConditionsMesherProcess
    : public MesherProcess
  {
  public:
 
    KRATOS_CLASS_POINTER_DEFINITION( ClearContactConditionsMesherProcess );
 
 
    ClearContactConditionsMesherProcess(ModelPart& rModelPart,
                  int EchoLevel = 0)
      : mrModelPart(rModelPart)
    {
      mEchoLevel = EchoLevel;
    }
 
    virtual ~ClearContactConditionsMesherProcess()
    {
    }
 
 
 
    void operator()()
    {
      Execute();
    }
 
 
 
    void Execute() override
    {
      KRATOS_TRY
 
      this->ClearContactConditions();
 
      KRATOS_CATCH(" ")
    };
 
 
 
 
 
 
 
    std::string Info() const override
    {
      return "ClearContactConditionsMesherProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
      rOStream << "ClearContactConditionsMesherProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
  protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
  private:
 
 
    ModelPart& mrModelPart;
 
    int mEchoLevel;
 
 
 
    //**************************************************************************
    //**************************************************************************
 
    void ClearContactConditions()
    {
 
      KRATOS_TRY
 
      //Clear contact conditions from contact domain
    if( mrModelPart.IsNot(CONTACT) )
      std::cout<<" ModelPart Supplied do not corresponds to the Contact Domain: ("<<mrModelPart.Name()<<")"<<std::endl;
 
 
      ClearContactConditions(mrModelPart);
 
 
      KRATOS_CATCH( "" )
 
    }
 
    //**************************************************************************
    //**************************************************************************
 
    void ClearContactConditions(ModelPart& rModelPart)
    {
 
      KRATOS_TRY
 
      //*******************************************************************
      //clearing contact conditions
      //
 
      if( mEchoLevel >= 1 ){
    std::cout<<" ["<<rModelPart.Name()<<" :: CONDITIONS [OLD:"<<rModelPart.NumberOfConditions();
      }
 
      ModelPart::ConditionsContainerType PreservedConditions;
 
      for(ModelPart::ConditionsContainerType::iterator ic = rModelPart.ConditionsBegin(); ic!= rModelPart.ConditionsEnd(); ++ic)
    {
      if(ic->IsNot(CONTACT)){
        PreservedConditions.push_back(*(ic.base()));
      }
    }
 
      rModelPart.Conditions().swap(PreservedConditions);
 
      if( mEchoLevel >= 1 ){
    std::cout<<" / PRE:"<<rModelPart.NumberOfConditions();
      }
 
      //rModelPart.Conditions().Sort();
      //rModelPart.Conditions().Unique();
 
      if( mEchoLevel >= 1 ){
    std::cout<<" / NEW:"<<rModelPart.NumberOfConditions()<<"] "<<std::endl;
      }
 
      KRATOS_CATCH( "" )
 
    }
 
 
 
 
 
 
 
    ClearContactConditionsMesherProcess& operator=(ClearContactConditionsMesherProcess const& rOther);
 
    //ClearContactConditionsMesherProcess(ClearContactConditionsMesherProcess const& rOther);
 
 
 
  }; // Class ClearContactConditionsMesherProcess
 
 
 
 
 
 
  inline std::istream& operator >> (std::istream& rIStream,
                    ClearContactConditionsMesherProcess& rThis);
 
  inline std::ostream& operator << (std::ostream& rOStream,
                    const ClearContactConditionsMesherProcess& rThis)
  {
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
  }
 
 
}  // namespace Kratos.
 
#endif // CONTACT_CLEAR_CONTACT_CONDITIONS_MESHER_PROCESS_H_INCLUDED  defined