dragonfly-custom.C 125 KB
Newer Older
Nikhil's avatar
Nikhil committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
/*
 * Copyright (C) 2013 University of Chicago.
 * See COPYRIGHT notice in top-level directory.
 *
 */

#include <ross.h>

#include "codes/jenkins-hash.h"
#include "codes/codes_mapping.h"
#include "codes/codes.h"
#include "codes/model-net.h"
#include "codes/model-net-method.h"
#include "codes/model-net-lp.h"
15
#include "codes/net/dragonfly-custom.h"
Nikhil's avatar
Nikhil committed
16
17
18
#include "sys/file.h"
#include "codes/quickhash.h"
#include "codes/rc-stack.h"
19
20
#include <vector>
#include <map>
21
#include <set>
Nikhil's avatar
Nikhil committed
22

23
24
25
26
#ifdef ENABLE_CORTEX
#include <cortex/cortex.h>
#include <cortex/topology.h>
#endif
Nikhil's avatar
Nikhil committed
27

28
#define DUMP_CONNECTIONS 0
29
#define CREDIT_SIZE 8
Misbah Mubarak's avatar
Misbah Mubarak committed
30
#define DFLY_HASH_TABLE_SIZE 4999
31
// debugging parameters
32
#define DEBUG_LP 892
33
#define T_ID 10
34
#define TRACK -1
35
#define TRACK_PKT 0
36
37
38
#define TRACK_MSG -1
#define DEBUG 0
#define MAX_STATS 65536
39
#define SHOW_ADAP_STATS 1
40
41
42
43
44

#define LP_CONFIG_NM_TERM (model_net_lp_config_names[DRAGONFLY_CUSTOM])
#define LP_METHOD_NM_TERM (model_net_method_names[DRAGONFLY_CUSTOM])
#define LP_CONFIG_NM_ROUT (model_net_lp_config_names[DRAGONFLY_CUSTOM_ROUTER])
#define LP_METHOD_NM_ROUT (model_net_method_names[DRAGONFLY_CUSTOM_ROUTER])
Nikhil's avatar
Nikhil committed
45

46
static int BIAS_MIN = 1;
47
static int DF_DALLY = 0;
48
static int adaptive_threshold = 1024;
49

50
51
52
static long num_local_packets_sr = 0;
static long num_local_packets_sg = 0;
static long num_remote_packets = 0;
Nikhil's avatar
Nikhil committed
53
54
55
56
57
58
59
using namespace std;
struct Link {
  int offset, type;
};
struct bLink {
  int offset, dest;
};
60
61
62
63
64
/* Each entry in the vector is for a router id
 * against each router id, there is a map of links (key of the map is the dest
 * router id)
 * link has information on type (green or black) and offset (number of links
 * between that particular source and dest router ID)*/
Nikhil's avatar
Nikhil committed
65
vector< map< int, vector<Link> > > intraGroupLinks;
66
67
/* contains mapping between source router and destination group via link (link
 * has dest ID)*/
Nikhil's avatar
Nikhil committed
68
vector< map< int, vector<bLink> > > interGroupLinks;
69
/*MM: Maintains a list of routers connecting the source and destination groups */
Nikhil's avatar
Nikhil committed
70
71
72
73
74
75
76
77
78
79
vector< vector< vector<int> > > connectionList;

struct IntraGroupLink {
  int src, dest, type;
};

struct InterGroupLink {
  int src, dest;
};

80
81
82
83
84
85
#ifdef ENABLE_CORTEX
/* This structure is defined at the end of the file */
extern "C" {
extern cortex_topology dragonfly_custom_cortex_topology;
}
#endif
Nikhil's avatar
Nikhil committed
86

87
88
89
static int debug_slot_count = 0;
static long term_ecount, router_ecount, term_rev_ecount, router_rev_ecount;
static long packet_gen = 0, packet_fin = 0;
Nikhil's avatar
Nikhil committed
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108

static double maxd(double a, double b) { return a < b ? b : a; }

/* minimal and non-minimal packet counts for adaptive routing*/
static int minimal_count=0, nonmin_count=0;
static int num_routers_per_mgrp = 0;

typedef struct dragonfly_param dragonfly_param;
/* annotation-specific parameters (unannotated entry occurs at the 
 * last index) */
static uint64_t                  num_params = 0;
static dragonfly_param         * all_params = NULL;
static const config_anno_map_t * anno_map   = NULL;

/* global variables for codes mapping */
static char lp_group_name[MAX_NAME_LENGTH];
static int mapping_grp_id, mapping_type_id, mapping_rep_id, mapping_offset;

/* router magic number */
109
static int router_magic_num = 0;
Nikhil's avatar
Nikhil committed
110
111

/* terminal magic number */
112
static int terminal_magic_num = 0;
Nikhil's avatar
Nikhil committed
113

114
115
116
117
/* Hops within a group */
static int num_intra_nonmin_hops = 4;
static int num_intra_min_hops = 2;

118
static FILE * dragonfly_log = NULL;
Nikhil's avatar
Nikhil committed
119

120
121
static int sample_bytes_written = 0;
static int sample_rtr_bytes_written = 0;
Nikhil's avatar
Nikhil committed
122

123
124
static char cn_sample_file[MAX_NAME_LENGTH];
static char router_sample_file[MAX_NAME_LENGTH];
Nikhil's avatar
Nikhil committed
125

126
127
//don't do overhead here - job of MPI layer
static tw_stime mpi_soft_overhead = 0;
128

129
130
131
typedef struct terminal_custom_message_list terminal_custom_message_list;
struct terminal_custom_message_list {
    terminal_custom_message msg;
Nikhil's avatar
Nikhil committed
132
    char* event_data;
133
134
    terminal_custom_message_list *next;
    terminal_custom_message_list *prev;
Nikhil's avatar
Nikhil committed
135
136
};

137
static void init_terminal_custom_message_list(terminal_custom_message_list *thisO, 
138
    terminal_custom_message *inmsg) {
139
140
141
142
    thisO->msg = *inmsg;
    thisO->event_data = NULL;
    thisO->next = NULL;
    thisO->prev = NULL;
Nikhil's avatar
Nikhil committed
143
144
}

145
146
147
148
static void delete_terminal_custom_message_list(void *thisO) {
    terminal_custom_message_list* toDel = (terminal_custom_message_list*)thisO;
    if(toDel->event_data != NULL) free(toDel->event_data);
    free(toDel);
Nikhil's avatar
Nikhil committed
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
}

struct dragonfly_param
{
    // configuration parameters
    int num_routers; /*Number of routers in a group*/
    double local_bandwidth;/* bandwidth of the router-router channels within a group */
    double global_bandwidth;/* bandwidth of the inter-group router connections */
    double cn_bandwidth;/* bandwidth of the compute node channels connected to routers */
    int num_vcs; /* number of virtual channels */
    int local_vc_size; /* buffer size of the router-router channels */
    int global_vc_size; /* buffer size of the global channels */
    int cn_vc_size; /* buffer size of the compute node channels */
    int chunk_size; /* full-sized packets are broken into smaller chunks.*/
    // derived parameters
    int num_cn;
165
166
    int intra_grp_radix;
    int num_col_chans;
167
    int num_row_chans;
168
169
    int num_router_rows;
    int num_router_cols;
Nikhil's avatar
Nikhil committed
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
    int num_groups;
    int radix;
    int total_routers;
    int total_terminals;
    int num_global_channels;
    double cn_delay;
    double local_delay;
    double global_delay;
    double credit_delay;
    double router_delay;
};

struct dfly_hash_key
{
    uint64_t message_id;
    tw_lpid sender_id;
};

struct dfly_router_sample
{
    tw_lpid router_id;
    tw_stime* busy_time;
    int64_t* link_traffic_sample;
    tw_stime end_time;
    long fwd_events;
    long rev_events;
};

struct dfly_cn_sample
{
   tw_lpid terminal_id;
   long fin_chunks_sample;
   long data_size_sample;
   double fin_hops_sample;
   tw_stime fin_chunks_time;
   tw_stime busy_time_sample;
   tw_stime end_time;
   long fwd_events;
   long rev_events;
};

struct dfly_qhash_entry
{
   struct dfly_hash_key key;
   char * remote_event_data;
   int num_chunks;
   int remote_event_size;
   struct qhash_head hash_link;
};

/* handles terminal and router events like packet generate/send/receive/buffer */
typedef struct terminal_state terminal_state;
typedef struct router_state router_state;

/* dragonfly compute node data structure */
struct terminal_state
{
   uint64_t packet_counter;

   int packet_gen;
   int packet_fin;

   // Dragonfly specific parameters
   unsigned int router_id;
   unsigned int terminal_id;

   // Each terminal will have an input and output channel with the router
   int* vc_occupancy; // NUM_VC
   int num_vcs;
   tw_stime terminal_available_time;
240
241
   terminal_custom_message_list **terminal_msgs;
   terminal_custom_message_list **terminal_msgs_tail;
Nikhil's avatar
Nikhil committed
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
   int in_send_loop;
   struct mn_stats dragonfly_stats_array[CATEGORY_MAX];

   struct rc_stack * st;
   int issueIdle;
   int terminal_length;

   const char * anno;
   const dragonfly_param *params;

   struct qhash_table *rank_tbl;
   uint64_t rank_tbl_pop;

   tw_stime   total_time;
   uint64_t total_msg_size;
   double total_hops;
   long finished_msgs;
   long finished_chunks;
   long finished_packets;

262
   tw_stime * last_buf_full;
Nikhil's avatar
Nikhil committed
263
   tw_stime busy_time;
264
265
266
267
   
   tw_stime max_latency;
   tw_stime min_latency;

Nikhil's avatar
Nikhil committed
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
   char output_buf[4096];
   /* For LP suspend functionality */
   int error_ct;

   /* For sampling */
   long fin_chunks_sample;
   long data_size_sample;
   double fin_hops_sample;
   tw_stime fin_chunks_time;
   tw_stime busy_time_sample;

   char sample_buf[4096];
   struct dfly_cn_sample * sample_stat;
   int op_arr_size;
   int max_arr_size;
283
  
Nikhil's avatar
Nikhil committed
284
285
286
287
288
289
   /* for logging forward and reverse events */
   long fwd_events;
   long rev_events;
};

/* terminal event type (1-4) */
290
typedef enum event_t
Nikhil's avatar
Nikhil committed
291
292
293
294
295
296
297
298
{
  T_GENERATE=1,
  T_ARRIVE,
  T_SEND,
  T_BUFFER,
  R_SEND,
  R_ARRIVE,
  R_BUFFER,
299
} event_t;
Nikhil's avatar
Nikhil committed
300
301
302
303

/* whether the last hop of a packet was global, local or a terminal */
enum last_hop
{
304
   GLOBAL=1,
Nikhil's avatar
Nikhil committed
305
   LOCAL,
306
307
   TERMINAL,
   ROOT
Nikhil's avatar
Nikhil committed
308
309
310
311
312
313
314
};

/* three forms of routing algorithms available, adaptive routing is not
 * accurate and fully functional in the current version as the formulas
 * for detecting load on global channels are not very accurate */
enum ROUTING_ALGO
{
315
    MINIMAL = 1,
Nikhil's avatar
Nikhil committed
316
317
318
319
320
    NON_MINIMAL,
    ADAPTIVE,
    PROG_ADAPTIVE
};

321
322
323
324
325
enum LINK_TYPE
{
    GREEN,
    BLACK,
};
Nikhil's avatar
Nikhil committed
326
327
328
329
330
331
332
333
334
335
336
struct router_state
{
   unsigned int router_id;
   int group_id;
   int op_arr_size;
   int max_arr_size;

   int* global_channel; 
   
   tw_stime* next_output_available_time;
   tw_stime* cur_hist_start_time;
337
   tw_stime** last_buf_full;
Nikhil's avatar
Nikhil committed
338
339
340
341

   tw_stime* busy_time;
   tw_stime* busy_time_sample;

342
343
344
345
   terminal_custom_message_list ***pending_msgs;
   terminal_custom_message_list ***pending_msgs_tail;
   terminal_custom_message_list ***queued_msgs;
   terminal_custom_message_list ***queued_msgs_tail;
Nikhil's avatar
Nikhil committed
346
347
348
   int *in_send_loop;
   int *queued_count;
   struct rc_stack * st;
349
350

   int* last_sent_chan;
Nikhil's avatar
Nikhil committed
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
   int** vc_occupancy;
   int64_t* link_traffic;
   int64_t * link_traffic_sample;

   const char * anno;
   const dragonfly_param *params;

   int* prev_hist_num;
   int* cur_hist_num;
   
   char output_buf[4096];
   char output_buf2[4096];

   struct dfly_router_sample * rsamples;
   
   long fwd_events;
   long rev_events;
};

static short routing = MINIMAL;

static tw_stime         dragonfly_total_time = 0;
static tw_stime         dragonfly_max_latency = 0;


static long long       total_hops = 0;
static long long       N_finished_packets = 0;
static long long       total_msg_sz = 0;
static long long       N_finished_msgs = 0;
static long long       N_finished_chunks = 0;

static int dragonfly_rank_hash_compare(
        void *key, struct qhash_head *link)
{
    struct dfly_hash_key *message_key = (struct dfly_hash_key *)key;
    struct dfly_qhash_entry *tmp = NULL;

    tmp = qhash_entry(link, struct dfly_qhash_entry, hash_link);
    
    if (tmp->key.message_id == message_key->message_id
            && tmp->key.sender_id == message_key->sender_id)
        return 1;

    return 0;
}
static int dragonfly_hash_func(void *k, int table_size)
{
    struct dfly_hash_key *tmp = (struct dfly_hash_key *)k;
399
400
401
402
    uint32_t pc = 0, pb = 0;	
    bj_hashlittle2(tmp, sizeof(*tmp), &pc, &pb);
    return (int)(pc % (table_size - 1));
    /*uint64_t key = (~tmp->message_id) + (tmp->message_id << 18);
Nikhil's avatar
Nikhil committed
403
404
405
    key = key * 21;
    key = ~key ^ (tmp->sender_id >> 4);
    key = key * tmp->sender_id; 
406
    return (int)(key & (table_size - 1));*/
Nikhil's avatar
Nikhil committed
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
}

/* convert GiB/s and bytes to ns */
static tw_stime bytes_to_ns(uint64_t bytes, double GB_p_s)
{
    tw_stime time;

    /* bytes to GB */
    time = ((double)bytes)/(1024.0*1024.0*1024.0);
    /* GiB to s */
    time = time / GB_p_s;
    /* s to ns */
    time = time * 1000.0 * 1000.0 * 1000.0;

    return(time);
}

/* returns the dragonfly message size */
425
int dragonfly_custom_get_msg_sz(void)
Nikhil's avatar
Nikhil committed
426
{
427
	   return sizeof(terminal_custom_message);
Nikhil's avatar
Nikhil committed
428
429
430
431
}

static void free_tmp(void * ptr)
{
432
    struct dfly_qhash_entry * dfly = (dfly_qhash_entry *)ptr; 
433
434
435
436
437
    if(dfly->remote_event_data)
        free(dfly->remote_event_data);
   
    if(dfly)
        free(dfly);
Nikhil's avatar
Nikhil committed
438
}
439

440
441
442
static void append_to_terminal_custom_message_list(  
        terminal_custom_message_list ** thisq,
        terminal_custom_message_list ** thistail,
Nikhil's avatar
Nikhil committed
443
        int index, 
444
        terminal_custom_message_list *msg) {
Nikhil's avatar
Nikhil committed
445
446
447
448
449
450
451
452
453
    if(thisq[index] == NULL) {
        thisq[index] = msg;
    } else {
        thistail[index]->next = msg;
        msg->prev = thistail[index];
    } 
    thistail[index] = msg;
}

454
455
456
static void prepend_to_terminal_custom_message_list(  
        terminal_custom_message_list ** thisq,
        terminal_custom_message_list ** thistail,
Nikhil's avatar
Nikhil committed
457
        int index, 
458
        terminal_custom_message_list *msg) {
Nikhil's avatar
Nikhil committed
459
460
461
462
463
464
465
466
467
    if(thisq[index] == NULL) {
        thistail[index] = msg;
    } else {
        thisq[index]->prev = msg;
        msg->next = thisq[index];
    } 
    thisq[index] = msg;
}

468
469
470
static terminal_custom_message_list* return_head(
        terminal_custom_message_list ** thisq,
        terminal_custom_message_list ** thistail,
Nikhil's avatar
Nikhil committed
471
        int index) {
472
    terminal_custom_message_list *head = thisq[index];
Nikhil's avatar
Nikhil committed
473
474
475
476
477
478
479
480
481
482
483
484
    if(head != NULL) {
        thisq[index] = head->next;
        if(head->next != NULL) {
            head->next->prev = NULL;
            head->next = NULL;
        } else {
            thistail[index] = NULL;
        }
    }
    return head;
}

485
486
487
static terminal_custom_message_list* return_tail(
        terminal_custom_message_list ** thisq,
        terminal_custom_message_list ** thistail,
Nikhil's avatar
Nikhil committed
488
        int index) {
489
    terminal_custom_message_list *tail = thistail[index];
Nikhil's avatar
Nikhil committed
490
491
492
493
494
495
496
497
498
499
500
501
502
    assert(tail);
    if(tail->prev != NULL) {
        tail->prev->next = NULL;
        thistail[index] = tail->prev;
        tail->prev = NULL;
    } else {
        thistail[index] = NULL;
        thisq[index] = NULL;
    }
    return tail;
}

static void dragonfly_read_config(const char * anno, dragonfly_param *params){
503
504
505
506
507
508
509
510
    /*Adding init for router magic number*/
    uint32_t h1 = 0, h2 = 0; 
    bj_hashlittle2(LP_METHOD_NM_ROUT, strlen(LP_METHOD_NM_ROUT), &h1, &h2);
    router_magic_num = h1 + h2;
    
    bj_hashlittle2(LP_METHOD_NM_TERM, strlen(LP_METHOD_NM_TERM), &h1, &h2);
    terminal_magic_num = h1 + h2;
    
Nikhil's avatar
Nikhil committed
511
512
    // shorthand
    dragonfly_param *p = params;
Nikhil's avatar
Nikhil committed
513
    int myRank;
514
    MPI_Comm_rank(MPI_COMM_CODES, &myRank);
Nikhil's avatar
Nikhil committed
515

516
    int rc = configuration_get_value_int(&config, "PARAMS", "local_vc_size", anno, &p->local_vc_size);
Nikhil's avatar
Nikhil committed
517
518
519
520
521
    if(rc) {
        p->local_vc_size = 1024;
        fprintf(stderr, "Buffer size of local channels not specified, setting to %d\n", p->local_vc_size);
    }

522
523
524
    rc = configuration_get_value_int(&config, "PARAMS", "adaptive_threshold", anno, &adaptive_threshold);
    if(rc) {
    	adaptive_threshold = p->local_vc_size / 8;
525
        printf("\n Setting adaptive threshold to %d ", adaptive_threshold);
526
	}
527
528
529
530
    else
    {
        printf("\n Setting adaptive threshold to %d ", adaptive_threshold);
    }
531

Nikhil's avatar
Nikhil committed
532
533
534
535
536
537
    rc = configuration_get_value_int(&config, "PARAMS", "global_vc_size", anno, &p->global_vc_size);
    if(rc) {
        p->global_vc_size = 2048;
        fprintf(stderr, "Buffer size of global channels not specified, setting to %d\n", p->global_vc_size);
    }

538
539
540
541
    rc = configuration_get_value_int(&config, "PARAMS", "df-dally-vc", anno, &DF_DALLY);
    if(rc) {
        DF_DALLY = 0;
    }
542
543
544
545
546
    
    rc = configuration_get_value_int(&config, "PARAMS", "minimal-bias", anno, &BIAS_MIN);
    if(rc) {
        BIAS_MIN = 0;
    }
547
548
549
    else
	printf("\n Setting minimal bias");

Nikhil's avatar
Nikhil committed
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
    rc = configuration_get_value_int(&config, "PARAMS", "cn_vc_size", anno, &p->cn_vc_size);
    if(rc) {
        p->cn_vc_size = 1024;
        fprintf(stderr, "Buffer size of compute node channels not specified, setting to %d\n", p->cn_vc_size);
    }

    rc = configuration_get_value_int(&config, "PARAMS", "chunk_size", anno, &p->chunk_size);
    if(rc) {
        p->chunk_size = 512;
        fprintf(stderr, "Chunk size for packets is specified, setting to %d\n", p->chunk_size);
    }

    rc = configuration_get_value_double(&config, "PARAMS", "local_bandwidth", anno, &p->local_bandwidth);
    if(rc) {
        p->local_bandwidth = 5.25;
        fprintf(stderr, "Bandwidth of local channels not specified, setting to %lf\n", p->local_bandwidth);
    }

    rc = configuration_get_value_double(&config, "PARAMS", "global_bandwidth", anno, &p->global_bandwidth);
    if(rc) {
        p->global_bandwidth = 4.7;
        fprintf(stderr, "Bandwidth of global channels not specified, setting to %lf\n", p->global_bandwidth);
    }

    rc = configuration_get_value_double(&config, "PARAMS", "cn_bandwidth", anno, &p->cn_bandwidth);
    if(rc) {
        p->cn_bandwidth = 5.25;
        fprintf(stderr, "Bandwidth of compute node channels not specified, setting to %lf\n", p->cn_bandwidth);
    }

580
    rc = configuration_get_value_double(&config, "PARAMS", "router_delay", anno,
Nikhil's avatar
Nikhil committed
581
            &p->router_delay);
582
    if(rc) {
583
584
      p->router_delay = 100;
    }
Nikhil's avatar
Nikhil committed
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601

    configuration_get_value(&config, "PARAMS", "cn_sample_file", anno, cn_sample_file,
            MAX_NAME_LENGTH);
    configuration_get_value(&config, "PARAMS", "rt_sample_file", anno, router_sample_file,
            MAX_NAME_LENGTH);
    
    char routing_str[MAX_NAME_LENGTH];
    configuration_get_value(&config, "PARAMS", "routing", anno, routing_str,
            MAX_NAME_LENGTH);
    if(strcmp(routing_str, "minimal") == 0)
        routing = MINIMAL;
    else if(strcmp(routing_str, "nonminimal")==0 || 
            strcmp(routing_str,"non-minimal")==0)
        routing = NON_MINIMAL;
    else if (strcmp(routing_str, "adaptive") == 0)
        routing = ADAPTIVE;
    else if (strcmp(routing_str, "prog-adaptive") == 0)
Nikhil's avatar
Nikhil committed
602
	      routing = PROG_ADAPTIVE;
Nikhil's avatar
Nikhil committed
603
604
605
606
607
608
609
    else
    {
        fprintf(stderr, 
                "No routing protocol specified, setting to minimal routing\n");
        routing = -1;
    }

610
611
612
613
614
615
616
617
618
619
    // rc = configuration_get_value_int(&config, "PARAMS", "num_vcs_override", anno, &p->num_vcs);
    // if(rc) {
    //     if(routing == PROG_ADAPTIVE)
    //         p->num_vcs = 10;
    //     else
    //         p->num_vcs = 8;
    // }
    // else {
    //     printf("Overriding num_vcs: p->num_vcs=%d\n",p->num_vcs);
    // }
620
621
622
   
if(DF_DALLY == 0) 
{
623
624
625
626
    if(routing == PROG_ADAPTIVE)
        p->num_vcs = 10;
    else
        p->num_vcs = 8;
627
628
629
}
else
{
630
    if(routing == PROG_ADAPTIVE)
631
        p->num_vcs = 5;
632
    else
633
        p->num_vcs = 4;
634
}
Nikhil's avatar
Nikhil committed
635
636
637
    rc = configuration_get_value_int(&config, "PARAMS", "num_groups", anno, &p->num_groups);
    if(rc) {
      printf("Number of groups not specified. Aborting");
638
      MPI_Abort(MPI_COMM_CODES, 1);
Nikhil's avatar
Nikhil committed
639
    }
640
641
642
643
644
    rc = configuration_get_value_int(&config, "PARAMS", "num_col_chans", anno, &p->num_col_chans);
    if(rc) {
//        printf("\n Number of links connecting chassis not specified, setting to default value 3 ");
        p->num_col_chans = 3;
    }
645
646
647
648
649
    rc = configuration_get_value_int(&config, "PARAMS", "num_row_chans", anno, &p->num_row_chans);
    if(rc) {
//        printf("\n Number of links connecting chassis not specified, setting to default value 3 ");
        p->num_row_chans = 1;
    }
650
651
652
653
654
655
656
657
658
659
    rc = configuration_get_value_int(&config, "PARAMS", "num_router_rows", anno, &p->num_router_rows);
    if(rc) {
        printf("\n Number of router rows not specified, setting to 6 ");
        p->num_router_rows = 6;
    }
    rc = configuration_get_value_int(&config, "PARAMS", "num_router_cols", anno, &p->num_router_cols);
    if(rc) {
        printf("\n Number of router columns not specified, setting to 16 ");
        p->num_router_cols = 16;
    }
660
661
662
663
    p->intra_grp_radix = (p->num_router_cols * p->num_row_chans);
    if(p->num_router_rows > 1)
        p->intra_grp_radix += (p->num_router_rows * p->num_col_chans);

664
665
    p->num_routers = p->num_router_rows * p->num_router_cols;
    
666
    rc = configuration_get_value_int(&config, "PARAMS", "num_cns_per_router", anno, &p->num_cn);
Nikhil's avatar
Nikhil committed
667
    if(rc) {
668
669
        printf("\n Number of cns per router not specified, setting to %d ", p->num_routers/2);
        p->num_cn = p->num_routers/2;
Nikhil's avatar
Nikhil committed
670
    }
671

672
673
    rc = configuration_get_value_int(&config, "PARAMS", "num_global_channels", anno, &p->num_global_channels);
    if(rc) {
674
675
        printf("\n Number of global channels per router not specified, setting to 10 ");
        p->num_global_channels = 10;
676
    }
677
    p->radix = p->intra_grp_radix + p->num_global_channels + p->num_cn;
Nikhil's avatar
Nikhil committed
678
679
    p->total_routers = p->num_groups * p->num_routers;
    p->total_terminals = p->total_routers * p->num_cn;
Nikhil's avatar
Nikhil committed
680
681
682
683
    
    // read intra group connections, store from a router's perspective
    // all links to the same router form a vector
    char intraFile[MAX_NAME_LENGTH];
684
    configuration_get_value(&config, "PARAMS", "intra-group-connections", 
Nikhil's avatar
Nikhil committed
685
        anno, intraFile, MAX_NAME_LENGTH);
686
687
    if(strlen(intraFile) <= 0) {
      tw_error(TW_LOC, "Intra group connections file not specified. Aborting");
Nikhil's avatar
Nikhil committed
688
689
    }
    FILE *groupFile = fopen(intraFile, "rb");
690
691
692
    if(!groupFile)
        tw_error(TW_LOC, "intra-group file not found ");

Nikhil's avatar
Nikhil committed
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
    if(!myRank)
      printf("Reading intra-group connectivity file: %s\n", intraFile);

    {
      vector< int > offsets;
      offsets.resize(p->num_routers, 0);
      intraGroupLinks.resize(p->num_routers);
      IntraGroupLink newLink;

      while(fread(&newLink, sizeof(IntraGroupLink), 1, groupFile) != 0) {
        Link tmpLink;
        tmpLink.type = newLink.type;
        tmpLink.offset = offsets[newLink.src]++;
        intraGroupLinks[newLink.src][newLink.dest].push_back(tmpLink);
      }
    }

    fclose(groupFile);

    // read inter group connections, store from a router's perspective
    // also create a group level table that tells all the connecting routers
    char interFile[MAX_NAME_LENGTH];
715
    configuration_get_value(&config, "PARAMS", "inter-group-connections", 
Nikhil's avatar
Nikhil committed
716
        anno, interFile, MAX_NAME_LENGTH);
717
718
    if(strlen(interFile) <= 0) {
      tw_error(TW_LOC, "Inter group connections file not specified. Aborting");
Nikhil's avatar
Nikhil committed
719
720
721
    }
    FILE *systemFile = fopen(interFile, "rb");
    if(!myRank)
722
    {
Nikhil's avatar
Nikhil committed
723
      printf("Reading inter-group connectivity file: %s\n", interFile);
724
725
      printf("\n Total routers %d total groups %d ", p->total_routers, p->num_groups);
    }
Nikhil's avatar
Nikhil committed
726
727
728
729
730
731
732
733
734

    {
      vector< int > offsets;
      offsets.resize(p->total_routers, 0);
      interGroupLinks.resize(p->total_routers);
      connectionList.resize(p->num_groups);
      for(int g = 0; g < connectionList.size(); g++) {
        connectionList[g].resize(p->num_groups);
      }
735
      
Nikhil's avatar
Nikhil committed
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
      InterGroupLink newLink;

      while(fread(&newLink, sizeof(InterGroupLink), 1, systemFile) != 0) {
        bLink tmpLink;
        tmpLink.dest = newLink.dest;
        int srcG = newLink.src / p->num_routers;
        int destG = newLink.dest / p->num_routers;
        tmpLink.offset = offsets[newLink.src]++;
        interGroupLinks[newLink.src][destG].push_back(tmpLink);
        int r;
        for(r = 0; r < connectionList[srcG][destG].size(); r++) {
          if(connectionList[srcG][destG][r] == newLink.src) break;
        }
        if(r == connectionList[srcG][destG].size()) {
          connectionList[srcG][destG].push_back(newLink.src);
        }
      }
    }

    fclose(systemFile);

757
#if DUMP_CONNECTIONS == 1
Nikhil's avatar
Nikhil committed
758
759
760
761
762
763
764
765
    printf("Dumping intra-group connections\n");
    for(int a = 0; a < intraGroupLinks.size(); a++) {
      printf("Connections for router %d\n", a);
      map< int, vector<Link> >  &curMap = intraGroupLinks[a];
      map< int, vector<Link> >::iterator it = curMap.begin();
      for(; it != curMap.end(); it++) {
        printf(" ( %d - ", it->first);
        for(int l = 0; l < it->second.size(); l++) {
766
          // offset is number of local connections
767
          // type is black or green according to Cray architecture 
Nikhil's avatar
Nikhil committed
768
769
770
771
772
773
774
          printf("%d,%d ", it->second[l].offset, it->second[l].type);
        }
        printf(")");
      }
      printf("\n");
    }
#endif
775
#if DUMP_CONNECTIONS == 1
Nikhil's avatar
Nikhil committed
776
777
778
    printf("Dumping inter-group connections\n");
    for(int a = 0; a < interGroupLinks.size(); a++) {
      printf("Connections for router %d\n", a);
779
780
      map< int, vector<bLink> >  &curMap = interGroupLinks[a];
      map< int, vector<bLink> >::iterator it = curMap.begin();
Nikhil's avatar
Nikhil committed
781
      for(; it != curMap.end(); it++) {
782
        // dest group ID 
Nikhil's avatar
Nikhil committed
783
784
        printf(" ( %d - ", it->first);
        for(int l = 0; l < it->second.size(); l++) {
785
786
            // dest is dest router ID
            // offset is number of global connections
Nikhil's avatar
Nikhil committed
787
788
789
790
791
792
793
794
          printf("%d,%d ", it->second[l].offset, it->second[l].dest);
        }
        printf(")");
      }
      printf("\n");
    }
#endif

795
#if DUMP_CONNECTIONS == 1
Nikhil's avatar
Nikhil committed
796
797
798
799
800
801
802
803
804
805
806
807
808
    printf("Dumping source aries for global connections\n");
    for(int g = 0; g < p->num_groups; g++) {
      for(int g1 = 0; g1 < p->num_groups; g1++) {
        printf(" ( ");
        for(int l = 0; l < connectionList[g][g1].size(); l++) {
          printf("%d ", connectionList[g][g1][l]);
        }
        printf(")");
      }
      printf("\n");
    }
#endif
    if(!myRank) {
809
        printf("\n Total nodes %d routers %d groups %d routers per group %d radix %d\n",
Nikhil's avatar
Nikhil committed
810
                p->num_cn * p->total_routers, p->total_routers, p->num_groups,
811
                p->num_routers, p->radix);
Nikhil's avatar
Nikhil committed
812
    }
Nikhil's avatar
Nikhil committed
813

Nikhil's avatar
Nikhil committed
814
815
816
    p->cn_delay = bytes_to_ns(p->chunk_size, p->cn_bandwidth);
    p->local_delay = bytes_to_ns(p->chunk_size, p->local_bandwidth);
    p->global_delay = bytes_to_ns(p->chunk_size, p->global_bandwidth);
817
    p->credit_delay = bytes_to_ns(CREDIT_SIZE, p->local_bandwidth); //assume 8 bytes packet
Nikhil's avatar
Nikhil committed
818
819
}

820
void dragonfly_custom_configure(){
Nikhil's avatar
Nikhil committed
821
822
823
    anno_map = codes_mapping_get_lp_anno_map(LP_CONFIG_NM_TERM);
    assert(anno_map);
    num_params = anno_map->num_annos + (anno_map->has_unanno_lp > 0);
824
    all_params = (dragonfly_param *)malloc(num_params * sizeof(*all_params));
Nikhil's avatar
Nikhil committed
825
826
827
828
829
830
831
832

    for (int i = 0; i < anno_map->num_annos; i++){
        const char * anno = anno_map->annotations[i].ptr;
        dragonfly_read_config(anno, &all_params[i]);
    }
    if (anno_map->has_unanno_lp > 0){
        dragonfly_read_config(NULL, &all_params[anno_map->num_annos]);
    }
833
834
835
#ifdef ENABLE_CORTEX
	model_net_topology = dragonfly_custom_cortex_topology;
#endif
Nikhil's avatar
Nikhil committed
836
837
838
}

/* report dragonfly statistics like average and maximum packet latency, average number of hops traversed */
839
void dragonfly_custom_report_stats()
Nikhil's avatar
Nikhil committed
840
841
842
843
844
845
{
   long long avg_hops, total_finished_packets, total_finished_chunks;
   long long total_finished_msgs, final_msg_sz;
   tw_stime avg_time, max_time;
   int total_minimal_packets, total_nonmin_packets;
   long total_gen, total_fin;
846
   long total_local_packets_sr, total_local_packets_sg, total_remote_packets;
Nikhil's avatar
Nikhil committed
847

848
849
850
851
852
853
854
   MPI_Reduce( &total_hops, &avg_hops, 1, MPI_LONG_LONG, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &N_finished_packets, &total_finished_packets, 1, MPI_LONG_LONG, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &N_finished_msgs, &total_finished_msgs, 1, MPI_LONG_LONG, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &N_finished_chunks, &total_finished_chunks, 1, MPI_LONG_LONG, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &total_msg_sz, &final_msg_sz, 1, MPI_LONG_LONG, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &dragonfly_total_time, &avg_time, 1,MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_CODES);
   MPI_Reduce( &dragonfly_max_latency, &max_time, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_CODES);
Nikhil's avatar
Nikhil committed
855
   
856
   MPI_Reduce( &packet_gen, &total_gen, 1, MPI_LONG, MPI_SUM, 0, MPI_COMM_CODES);
857
   MPI_Reduce(&packet_fin, &total_fin, 1, MPI_LONG, MPI_SUM, 0, MPI_COMM_CODES);
858
859
    MPI_Reduce( &num_local_packets_sr, &total_local_packets_sr, 1, MPI_LONG, MPI_SUM, 0, MPI_COMM_CODES);
    MPI_Reduce( &num_local_packets_sg, &total_local_packets_sg, 1, MPI_LONG, MPI_SUM, 0, MPI_COMM_CODES);
860
   MPI_Reduce( &num_remote_packets, &total_remote_packets, 1, MPI_LONG, MPI_SUM, 0, MPI_COMM_CODES);
861
   if(routing == ADAPTIVE || routing == PROG_ADAPTIVE || SHOW_ADAP_STATS)
Nikhil's avatar
Nikhil committed
862
    {
863
864
	MPI_Reduce(&minimal_count, &total_minimal_packets, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_CODES);
 	MPI_Reduce(&nonmin_count, &total_nonmin_packets, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_CODES);
Nikhil's avatar
Nikhil committed
865
866
867
868
869
870
871
    }

   /* print statistics */
   if(!g_tw_mynode)
   {	
      printf(" Average number of hops traversed %f average chunk latency %lf us maximum chunk latency %lf us avg message size %lf bytes finished messages %lld finished chunks %lld \n", 
              (float)avg_hops/total_finished_chunks, avg_time/(total_finished_chunks*1000), max_time/1000, (float)final_msg_sz/total_finished_msgs, total_finished_msgs, total_finished_chunks);
872
     if(routing == ADAPTIVE || routing == PROG_ADAPTIVE || SHOW_ADAP_STATS)
Nikhil's avatar
Nikhil committed
873
874
875
              printf("\n ADAPTIVE ROUTING STATS: %d chunks routed minimally %d chunks routed non-minimally completed packets %lld \n", 
                      total_minimal_packets, total_nonmin_packets, total_finished_chunks);
 
876
      printf("\n Total packets generated %ld finished %ld Locally routed- same router %ld different-router %ld Remote (inter-group) %ld \n", total_gen, total_fin, total_local_packets_sr, total_local_packets_sg, total_remote_packets);
Nikhil's avatar
Nikhil committed
877
878
879
880
881
882
883
   }
   return;
}


/* initialize a dragonfly compute node terminal */
void 
884
terminal_custom_init( terminal_state * s, 
Nikhil's avatar
Nikhil committed
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
	       tw_lp * lp )
{
    s->packet_gen = 0;
    s->packet_fin = 0;

    int i;
    char anno[MAX_NAME_LENGTH];

    // Assign the global router ID
    // TODO: be annotation-aware
    codes_mapping_get_lp_info(lp->gid, lp_group_name, &mapping_grp_id, NULL,
            &mapping_type_id, anno, &mapping_rep_id, &mapping_offset);
    if (anno[0] == '\0'){
        s->anno = NULL;
        s->params = &all_params[num_params-1];
    }
    else{
        s->anno = strdup(anno);
        int id = configuration_get_annotation_index(anno, anno_map);
        s->params = &all_params[id];
    }

   int num_lps = codes_mapping_get_lp_count(lp_group_name, 1, LP_CONFIG_NM_TERM,
           s->anno, 0);

910
   s->terminal_id = codes_mapping_get_lp_relative_id(lp->gid, 0, 0);
911
   s->router_id=(int)s->terminal_id / (s->params->num_cn);
Nikhil's avatar
Nikhil committed
912
913
   s->terminal_available_time = 0.0;
   s->packet_counter = 0;
914
915
916
   s->min_latency = INT_MAX;
   s->max_latency = 0;   

Nikhil's avatar
Nikhil committed
917
918
919
920
921
922
923
924
925
926
927
928
929
930
   s->finished_msgs = 0;
   s->finished_chunks = 0;
   s->finished_packets = 0;
   s->total_time = 0.0;
   s->total_msg_size = 0;

   s->busy_time = 0.0;

   s->fwd_events = 0;
   s->rev_events = 0;

   rc_stack_create(&s->st);
   s->num_vcs = 1;
   s->vc_occupancy = (int*)malloc(s->num_vcs * sizeof(int));
931
   s->last_buf_full = (tw_stime*)malloc(s->num_vcs * sizeof(tw_stime));
Nikhil's avatar
Nikhil committed
932
933
934

   for( i = 0; i < s->num_vcs; i++ )
    {
935
      s->last_buf_full[i] = 0.0;
Nikhil's avatar
Nikhil committed
936
937
938
939
      s->vc_occupancy[i]=0;
    }


940
   s->rank_tbl = NULL;
Nikhil's avatar
Nikhil committed
941
   s->terminal_msgs = 
942
       (terminal_custom_message_list**)malloc(s->num_vcs*sizeof(terminal_custom_message_list*));
Nikhil's avatar
Nikhil committed
943
   s->terminal_msgs_tail = 
944
       (terminal_custom_message_list**)malloc(s->num_vcs*sizeof(terminal_custom_message_list*));
Nikhil's avatar
Nikhil committed
945
946
947
948
949
950
951
952
953
954
955
   s->terminal_msgs[0] = NULL;
   s->terminal_msgs_tail[0] = NULL;
   s->terminal_length = 0;
   s->in_send_loop = 0;
   s->issueIdle = 0;

   return;
}

/* sets up the router virtual channels, global channels, 
 * local channels, compute node channels */
956
void router_custom_setup(router_state * r, tw_lp * lp)
Nikhil's avatar
Nikhil committed
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
{
    
    char anno[MAX_NAME_LENGTH];
    codes_mapping_get_lp_info(lp->gid, lp_group_name, &mapping_grp_id, NULL,
            &mapping_type_id, anno, &mapping_rep_id, &mapping_offset);

    if (anno[0] == '\0'){
        r->anno = NULL;
        r->params = &all_params[num_params-1];
    } else{
        r->anno = strdup(anno);
        int id = configuration_get_annotation_index(anno, anno_map);
        r->params = &all_params[id];
    }

    // shorthand
    const dragonfly_param *p = r->params;

975
    num_routers_per_mgrp = codes_mapping_get_lp_count (lp_group_name, 1, "modelnet_dragonfly_custom_router",
Nikhil's avatar
Nikhil committed
976
977
978
979
980
981
982
            NULL, 0);
    int num_grp_reps = codes_mapping_get_group_reps(lp_group_name);
    if(p->total_routers != num_grp_reps * num_routers_per_mgrp)
        tw_error(TW_LOC, "\n Config error: num_routers specified %d total routers computed in the network %d "
                "does not match with repetitions * dragonfly_router %d  ",
                p->num_routers, p->total_routers, num_grp_reps * num_routers_per_mgrp);

983
   r->router_id = codes_mapping_get_lp_relative_id(lp->gid, 0, 0);
Nikhil's avatar
Nikhil committed
984
   r->group_id=r->router_id/p->num_routers;
985
986
    
   //printf("\n Local router id %d global id %d ", r->router_id, lp->gid);
Nikhil's avatar
Nikhil committed
987
988
989
990

   r->fwd_events = 0;
   r->rev_events = 0;

991

Nikhil's avatar
Nikhil committed
992
993
994
995
996
997
998
   r->global_channel = (int*)malloc(p->num_global_channels * sizeof(int));
   r->next_output_available_time = (tw_stime*)malloc(p->radix * sizeof(tw_stime));
   r->cur_hist_start_time = (tw_stime*)malloc(p->radix * sizeof(tw_stime));
   r->link_traffic = (int64_t*)malloc(p->radix * sizeof(int64_t));
   r->link_traffic_sample = (int64_t*)malloc(p->radix * sizeof(int64_t));
   r->cur_hist_num = (int*)malloc(p->radix * sizeof(int));
   r->prev_hist_num = (int*)malloc(p->radix * sizeof(int));
999
1000
  
   r->last_sent_chan = (int*) malloc(p->num_router_rows * sizeof(int));
For faster browsing, not all history is shown. View entire blame