context.h

Go to the documentation of this file.

#ifndef GHOST_CONTEXT_H
#define GHOST_CONTEXT_H

#include "config.h"
#include "types.h"
#include "error.h"
#include "map.h"

typedef struct ghost_context ghost_context;

typedef enum
{
    GHOST_PERMUTATION_ORIG2PERM,
    GHOST_PERMUTATION_PERM2ORIG
}
ghost_permutation_direction;

typedef enum {
    GHOST_CONTEXT_DEFAULT = 0, 
    GHOST_CONTEXT_DIST_NZ = 4,
    GHOST_CONTEXT_DIST_ROWS = 8

} ghost_context_flags_t;

#ifdef __cplusplus
inline ghost_context_flags_t operator|(const ghost_context_flags_t &a,
        const ghost_context_flags_t &b)
{
    return static_cast<ghost_context_flags_t>(
            static_cast<int>(a) | static_cast<int>(b));
}

inline ghost_context_flags_t operator&(const ghost_context_flags_t &a,
        const ghost_context_flags_t &b)
{
    return static_cast<ghost_context_flags_t>(
            static_cast<int>(a) & static_cast<int>(b));
}
#endif


typedef enum{
      GHOST_KACZ_METHOD_MC,
      GHOST_KACZ_METHOD_BMC_RB,
      GHOST_KACZ_METHOD_BMC_one_sweep,
      GHOST_KACZ_METHOD_BMC_two_sweep,
      GHOST_KACZ_METHOD_BMC,
      GHOST_KACZ_METHOD_BMCshift,
      GHOST_KACZ_METHOD_BMCNORMAL
}
ghost_kacz_method;

//TODO zone ptr can be moved here 
typedef struct {
      
      ghost_kacz_method kacz_method;
      ghost_lidx active_threads;
}
ghost_kacz_setting;
    

struct ghost_context
{
    ghost_gidx gnnz;
    ghost_lidx nnz;
    ghost_context_flags_t flags;
    double weight;
    ghost_mpi_comm mpicomm;
    ghost_mpi_comm mpicomm_parent;
    ghost_map *row_map;
    ghost_map *col_map;
    ghost_lidx * wishes; 
    ghost_lidx ** wishlist; 
    ghost_lidx * dues; 
    ghost_lidx ** duelist;
    ghost_lidx ** cu_duelist;
    ghost_lidx* hput_pos; 
    int *duepartners;
    int nduepartners;
    int *wishpartners;
    int nwishpartners;
    ghost_lidx *entsInCol;
    ghost_lidx nChunkAvg;
    ghost_lidx *avg_ptr;
    ghost_lidx nElemAvg;
    ghost_lidx *mapAvg; 
    ghost_lidx *mappedDuelist; //This might be removed in a future version
    int *nrankspresent;

    int nmats;
    ghost_gidx lowerBandwidth;
    ghost_gidx upperBandwidth;
    ghost_gidx bandwidth;
    ghost_gidx maxColRange; 
    ghost_lidx ncolors;
    ghost_lidx *color_ptr;
    ghost_lidx nzones;
    ghost_lidx *zone_ptr;
    ghost_kacz_setting kacz_setting;
    double kaczRatio;
 
};

extern const ghost_context GHOST_CONTEXT_INITIALIZER;


#ifdef __cplusplus
extern "C" {
#endif

    ghost_error ghost_context_create(ghost_context **context, ghost_gidx gnrows, ghost_gidx gncols, ghost_context_flags_t flags, ghost_mpi_comm comm, double weight); 
    
    ghost_error ghost_context_string(char **str, ghost_context *context);
    void ghost_context_destroy(ghost_context *ctx);
    
    char * ghost_context_workdist_string(ghost_context_flags_t flags);
    ghost_error ghost_global_perm_inv(ghost_gidx *toPerm, ghost_gidx *fromPerm, ghost_context *context);
    ghost_map *ghost_context_map(const ghost_context *ctx, const ghost_maptype mt);
    ghost_map *ghost_context_other_map(const ghost_context *ctx, const ghost_maptype mt);
    ghost_map *ghost_context_max_map(const ghost_context *ctx);
    ghost_error ghost_context_set_map(ghost_context *ctx, ghost_maptype which, ghost_map *map);
    ghost_error ghost_context_comm_string(char **str, ghost_context *ctx, int root);

#ifdef __cplusplus
} //extern "C"
#endif

#endif