trilinos_partitioned_fsi_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Źorrilla
//                   Jordi Cotela
//
 
#if !defined(KRATOS_TRILINOS_PARTITIONED_FSI_UTILITIES_H )
#define  KRATOS_TRILINOS_PARTITIONED_FSI_UTILITIES_H
 
 
 
// System includes
#include <string>
#include <iostream>
 
 
// External includes
#include "Epetra_MpiComm.h"
#include "Epetra_Map.h"
#include "Epetra_Vector.h"
 
// Project includes
#include "includes/define.h"
 
#include "../FSIApplication/custom_utilities/partitioned_fsi_utilities.hpp"
 
namespace Kratos
{
 
 
 
 
 
 
 
  template< class TSpace, class TValueType, unsigned int TDim >
  class TrilinosPartitionedFSIUtilities : public PartitionedFSIUtilities<TSpace,TValueType,TDim>
  {
  public:
 
      KRATOS_CLASS_POINTER_DEFINITION(TrilinosPartitionedFSIUtilities);
 
      typedef typename TSpace::VectorType                     VectorType;
      typedef typename TSpace::MatrixType                     MatrixType;
 
      typedef typename TSpace::VectorPointerType              VectorPointerType;
      typedef typename TSpace::MatrixPointerType              MatrixPointerType;
 
 
      TrilinosPartitionedFSIUtilities(const Epetra_MpiComm& EpetraCommunicator):
        mrEpetraComm(EpetraCommunicator)
      {}
 
      virtual ~TrilinosPartitionedFSIUtilities() override {}
 
 
 
 
 
      VectorPointerType SetUpInterfaceVector(ModelPart& rInterfaceModelPart) override
      {
          // Initialize Epetra_Map for the vector
          int NumLocalInterfaceDofs = rInterfaceModelPart.GetCommunicator().LocalMesh().NumberOfNodes() * TDim;
          int NumGlobalInterfaceDofs = rInterfaceModelPart.GetCommunicator().GetDataCommunicator().SumAll(NumLocalInterfaceDofs);
          int IndexBase = 0; // 0 for C-style vectors, 1 for Fortran numbering
          Epetra_Map InterfaceMap(NumGlobalInterfaceDofs,NumLocalInterfaceDofs,IndexBase,mrEpetraComm);
 
          // Create new vector using given map
          VectorPointerType p_int_vect = Kratos::make_shared<Epetra_FEVector>(InterfaceMap);
 
          // Set interface vector to zero
          p_int_vect->PutScalar(0.0);
 
          return p_int_vect;
      }
 
 
 
 
 
 
 
 
 
 
    protected:
 
 
 
 
 
 
 
      virtual void SetLocalValue(VectorType& rVector, int LocalRow, double Value) const override
      {
          // Note: We don't go through TSpace because TrilinosSpace forces a GlobalAssemble after each SetValue
          rVector.ReplaceMyValue(LocalRow,0,Value);
      }
 
      virtual double GetLocalValue(const VectorType& rVector, int LocalRow) const override
      {
          return rVector[0][LocalRow];
      }
 
 
 
 
 
 
 
 
 
    private:
 
 
 
      const Epetra_MpiComm& mrEpetraComm;
 
 
 
 
 
 
 
 
 
 
 
      TrilinosPartitionedFSIUtilities& operator=(TrilinosPartitionedFSIUtilities const& rOther){}
 
      TrilinosPartitionedFSIUtilities(TrilinosPartitionedFSIUtilities const& rOther){}
 
 
 
    }; // Class TrilinosPartitionedFSIUtilities
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_TRILINOS_PARTITIONED_FSI_UTILITIES_H  defined