merge_model_parts_process.h

Go to the documentation of this file.

//
//   Project Name:        Kratos
//   Last Modified by:    $Author: anonymous $
//   Date:                $Date: 2007-11-06 12:34:26 $
//   Revision:            $Revision: 1.4 $
//
//  this process save structural elements in a separate list
 
#if !defined(KRATOS_MERGE_MODEL_PARTS_PROCESS_INCLUDED )
#define  KRATOS_MERGE_MODEL_PARTS_PROCESS_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream>
#include <algorithm>
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "processes/process.h"
#include "includes/node.h"
#include "includes/element.h"
#include "includes/model_part.h"
#include "custom_elements/updated_lagrangian_fluid.h"
#include "custom_elements/updated_lagrangian_fluid3D.h"
#include "custom_elements/updated_lagrangian_fluid_inc.h"
#include "custom_elements/updated_lagrangian_fluid3D_inc.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
 
class MergeModelPartsProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(MergeModelPartsProcess);
 
 
    MergeModelPartsProcess()
    //ModelPart& fluid_model_part, ModelPart& structure_model_part, ModelPart& combined_model_part)
    //: mr_fluid_model_part(fluid_model_part), mr_structure_model_part(structure_model_part), mr_combined_model_part(combined_model_part)
    {
    }
 
    ~MergeModelPartsProcess() override
    {
    }
 
 
 
    //  void operator()()
    //  {
    //      MergeParts();
    //  }
 
 
 
    void MergeParts(ModelPart& fluid_model_part, ModelPart& structure_model_part, ModelPart& combined_model_part)
    {
        KRATOS_TRY
        combined_model_part.Elements().clear();
        combined_model_part.Nodes();
        combined_model_part.Nodes().clear();
 
        combined_model_part.Nodes()=fluid_model_part.Nodes();
        combined_model_part.Elements()=fluid_model_part.Elements();
        combined_model_part.Conditions()=fluid_model_part.Conditions();
 
        //sorting and renumbering the fluid elements
        unsigned int id=1;
        for(ModelPart::ElementsContainerType::iterator im = combined_model_part.ElementsBegin() ;
                im != combined_model_part.ElementsEnd() ; ++im)
        {
            im->SetId(id);
//              im->Id() = id;
            id++;
        }
        /*      KRATOS_WATCH("Last FLUID ID is");
                KRATOS_WATCH(id);
 
                KRATOS_WATCH("The size of fluid elements container is");
                KRATOS_WATCH(fluid_model_part.Elements().size());
        */
        //adding the structure elements to the combined part
        //combined_model_part.Elements().push_back(*( structure_model_part.ElementsBegin().base()));
        id += 1000;
        for(ModelPart::ElementsContainerType::iterator im = structure_model_part.ElementsBegin() ;
                im != structure_model_part.ElementsEnd() ; ++im)
        {
            im->SetId(id);
//              im->Id() = id;
            combined_model_part.Elements().push_back(*(im.base()));
            id++;
        }
        //combined_model_part.Elements().Sort();
 
 
        //renumbering
        id=1;
        for(ModelPart::ConditionsContainerType::iterator im = combined_model_part.ConditionsBegin() ;
                im != combined_model_part.ConditionsEnd() ; ++im)
        {
            im->SetId(id);
//              im->Id() = id;
            id++;
        }
        //adding the structure conditions to the combined part
        id += 1000;
        for(ModelPart::ConditionsContainerType::iterator im = structure_model_part.ConditionsBegin() ;
                im != structure_model_part.ConditionsEnd() ; ++im)
        {
            im->SetId(id);
//              im->Id() = id;
            combined_model_part.Conditions().push_back(*(im.base()));
            id++;
        }
        //combined_model_part.Conditions().Sort();
 
        //renumbering
        /*
        unsigned int new_id=1;
        for(ModelPart::ConditionsContainerType::iterator im = structure_model_part.ConditionsBegin() ;
            im != structure_model_part.ConditionsEnd() ; ++im)
        {
            im->Id() =new_id;
            new_id++;
        }
        */
        //WE HAVE TO COPY THE ProcessInfo pointer to the new part, otherwise it is empty
        combined_model_part.SetProcessInfo(fluid_model_part.pGetProcessInfo());
    combined_model_part.SetProperties(fluid_model_part.pProperties());
 
        KRATOS_CATCH("")
    }
 
 
 
 
 
 
 
    std::string Info() const override
    {
        return "MergeModelPartsProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "MergeModelPartsProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    //ModelPart& mr_fluid_model_part;
    //ModelPart& mr_structure_model_part;
    //ModelPart& mr_combined_model_part;
 
 
 
 
 
 
 
 
 
 
//      MergeModelPartsProcess& operator=(MergeModelPartsProcess const& rOther);
 
//      MergeModelPartsProcess(MergeModelPartsProcess const& rOther);
 
 
 
}; // Class MergeModelPartsProcess
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  MergeModelPartsProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const MergeModelPartsProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_MERGE_MODEL_PARTS_PROCESS_INCLUDED  defined