domain_decomposition.hpp

Go to the documentation of this file.

#ifndef DOMAIN_DECOMPOSITION_HPP
#define DOMAIN_DECOMPOSITION_HPP
#include "space_settings.hpp"
#include "NamelistReader.hpp"
#ifdef USE_MPI
#  include "my_mpi.hpp"
#endif
#include "field_decomposition.hpp"

// Decomposition class
template<class Device>
class DomainDecomposition {
    public:

    // Constructor
    DomainDecomposition(NLReader::NamelistReader& nlr, int nplanes, int nnodes);

#ifdef USE_MPI
    // MPI communicators etc
    const MyMPI mpi;
#endif

    // Updater
    void update_pol_decomp();

    // Constants
    bool pol_decomp;     
    int pe_per_plane;          
    int plane_totalpe;         
    int n_pids;                
    int my_pid;                
    int my_plane_pid;          

    bool decompose_fields; 
    FieldDecomposition<Device> field_decomp;

    // Variables
    int first_node;  
    int nnodes;      

    // Views
    Kokkos::View<int*,Kokkos::LayoutRight,HostType> gvid0_pid_h; 
    Kokkos::View<int*,Kokkos::LayoutRight,Device> gvid0_pid; 
    Kokkos::View<int*,Kokkos::LayoutRight,HostType> gvid0_pid_old_h; 
};

#include "domain_decomposition.tpp"
#endif