Reset "inited" in the unsync mode so as to avoid expensive atomic operation by
[charm.git] / NodeHelper.h
1 #ifndef _NODEHELPER_H
2 #define _NODEHELPER_H
3 #include <assert.h>
4
5 #include "charm++.h"
6 #include "NodeHelperAPI.h"
7
8 #define USE_TREE_BROADCAST_THRESHOLD 8
9 #define TREE_BCAST_BRANCH (4)
10 #define CACHE_LINE_SIZE 64
11
12 class FuncSingleHelper;
13
14 class CurLoopInfo{
15     friend class FuncSingleHelper;
16     
17 private:
18     volatile int curChunkIdx;
19     int numChunks;
20     HelperFn fnPtr;
21     int lowerIndex;
22     int upperIndex;
23     int paramNum;
24     void *param;
25     //limitation: only allow single variable reduction of size numChunks!!!
26     void **redBufs;
27         char *bufSpace;
28
29     volatile int finishFlag;
30     
31     //a tag to indicate whether the task for this new loop has been inited
32     //this tag is needed to prevent other helpers to run the old task
33     int inited;
34     
35 public:    
36     CurLoopInfo(int maxChunks):numChunks(0),fnPtr(NULL), lowerIndex(-1), upperIndex(0), 
37     paramNum(0), param(NULL), curChunkIdx(-1), finishFlag(0), redBufs(NULL), bufSpace(NULL), inited(0) 
38         {
39                 redBufs = new void *[maxChunks];
40                 bufSpace = new char[maxChunks * CACHE_LINE_SIZE];
41         for(int i=0; i<maxChunks; i++) redBufs[i] = (void *)(bufSpace+i*CACHE_LINE_SIZE);
42         }
43     
44     ~CurLoopInfo() { 
45                 delete [] redBufs; 
46                 delete [] bufSpace;
47         }
48     
49     void set(int nc, HelperFn f, int lIdx, int uIdx, int numParams, void *p){        /*
50       * WARNING: there's a rare data-racing case here. The current loop is
51       * about to finish (just before setting inited to 0; A helper (say B) 
52       * just enters the stealWork and passes the inited check. The helper 
53       * (say A) is very fast, and starts the next loop, and happens enter
54       * into the middle of this function. Then helper B will face corrupted
55       * task info as it is trying to execute the old loop task!
56       * In reality for user cases, this case happens very rarely!! -Chao Mei
57       */
58         numChunks = nc;
59         fnPtr = f;
60         lowerIndex = lIdx;
61         upperIndex = uIdx;
62         paramNum = numParams;
63         param = p;
64         curChunkIdx = -1;
65         finishFlag = 0;
66         //needs to be set last
67         inited = 1;
68     }
69       
70     void waitLoopDone(int sync){
71         //while(!__sync_bool_compare_and_swap(&finishFlag, numChunks, 0));
72         if(sync) while(finishFlag!=numChunks);
73         //finishFlag = 0;
74         inited = 0;
75     }
76     int getNextChunkIdx(){
77         return __sync_add_and_fetch(&curChunkIdx, 1);
78     }
79     void reportFinished(int counter){
80         if(counter==0) return;
81         __sync_add_and_fetch(&finishFlag, counter);
82     }
83     
84         void **getRedBufs() { return redBufs; }
85         
86     void stealWork();
87 };
88
89 /* FuncNodeHelper is a nodegroup object */
90
91 typedef struct converseNotifyMsg{
92     char core[CmiMsgHeaderSizeBytes];
93     int srcRank;
94     void *ptr;
95 }ConverseNotifyMsg;
96
97 class FuncNodeHelper : public CBase_FuncNodeHelper {
98     friend class FuncSingleHelper;
99         
100 public:
101     static int MAX_CHUNKS;
102 private:    
103     int numHelpers;    
104     FuncSingleHelper **helperPtr; /* ptrs to the FuncSingleHelpers it manages */
105         int useTreeBcast;
106     
107 public:
108         FuncNodeHelper();
109     ~FuncNodeHelper() {
110         delete [] helperPtr;
111     }
112     
113     void parallelizeFunc(HelperFn func, /* the function that finishes a partial work on another thread */
114                         int paramNum, void * param, /* the input parameters for the above func */
115                         int numChunks, /* number of chunks to be partitioned */
116                         int lowerRange, int upperRange, /* the loop-like parallelization happens in [lowerRange, upperRange] */                        
117                         int sync=1, /* whether the flow will continue until all chunks have finished */
118                         void *redResult=NULL, REDUCTION_TYPE type=NODEHELPER_NONE /* the reduction result, ONLY SUPPORT SINGLE VAR of TYPE int/float/double */
119                         );
120     void reduce(void **redBufs, void *redBuf, REDUCTION_TYPE type, int numChunks);
121 };
122
123 void SingleHelperStealWork(ConverseNotifyMsg *msg);
124
125 /* FuncSingleHelper is a chare located on every core of a node */
126 class FuncSingleHelper: public CBase_FuncSingleHelper {
127         friend class FuncNodeHelper;
128 private: 
129     FuncNodeHelper *thisNodeHelper;
130     ConverseNotifyMsg *notifyMsg;
131     CurLoopInfo *curLoop; /* Points to the current loop that is being processed */
132     
133 public:
134     FuncSingleHelper(size_t ndhPtr) {
135         thisNodeHelper = (FuncNodeHelper *)ndhPtr;
136         CmiAssert(thisNodeHelper!=NULL);        
137         int stealWorkHandler = CmiRegisterHandler((CmiHandler)SingleHelperStealWork);
138         curLoop = new CurLoopInfo(FuncNodeHelper::MAX_CHUNKS);
139         
140         notifyMsg = (ConverseNotifyMsg *)malloc(sizeof(ConverseNotifyMsg));
141         if(thisNodeHelper->useTreeBcast){
142             notifyMsg->srcRank = CmiMyRank();
143         }else{
144             notifyMsg->srcRank = -1;
145         }
146         notifyMsg->ptr = (void *)curLoop;
147         CmiSetHandler(notifyMsg, stealWorkHandler);
148         thisNodeHelper->helperPtr[CkMyRank()] = this;
149     }
150
151     ~FuncSingleHelper() {
152         delete curLoop;
153         delete notifyMsg;
154     }
155     
156     FuncSingleHelper(CkMigrateMessage *m) {}            
157 };
158
159 #endif