c8b054bd85ec213d8555997cb218bfafb10118a9
[charm.git] / src / conv-com / pipebroadcastconverse.C
1 #include <math.h>
2 #include "pipebroadcastconverse.h"
3
4 inline int log_of_2 (int i) {
5   int m;
6   for (m=0; i>(1<<m); ++m);
7   return m;
8 }
9
10 //PipeBcastHashKey CODE
11 int PipeBcastHashKey::staticCompare(const void *k1,const void *k2,size_t ){
12     return ((const PipeBcastHashKey *)k1)->
13                 compare(*(const PipeBcastHashKey *)k2);
14 }
15
16 CkHashCode PipeBcastHashKey::staticHash(const void *v,size_t){
17     return ((const PipeBcastHashKey *)v)->hash();
18 }
19
20 void PipeBroadcastConverse::commonInit(){
21   //log_of_2_inv = 1/log((double)2);
22   seqNumber = 0;
23 }
24
25 extern void propagate_handler(void *);
26
27 void propagate_handler_frag(void *message) {
28   int instid = CmiGetXHandler(message);
29   PipeBroadcastConverse *myStrategy = (PipeBroadcastConverse*)ConvComlibGetStrategy(instid);
30   ComlibPrintf("[%d] propagate_handler_frag: calling on instid %d %x\n",CkMyPe(),instid,myStrategy);
31   //CProxy_ComlibManager(CkpvAccess(cmgrID)).ckLocalBranch()->getStrategy(instid);
32   PipeBcastInfo *info = (PipeBcastInfo*)(((char*)message)+CmiReservedHeaderSize);
33   myStrategy->propagate((char*)message, true, info->srcPe, info->chunkSize+CmiReservedHeaderSize+sizeof(PipeBcastInfo), NULL);
34 }
35
36 void PipeBroadcastConverse::propagate(char *env, int isFragmented, int srcPeNumber, int totalSendingSize, setFunction setPeNumber){
37   // find destination processors and send
38   int destination, tmp, k, sizeToSend;
39   int num_pes, *dest_pes;
40   PipeBcastInfo *info = (PipeBcastInfo*)(env+CmiReservedHeaderSize);
41   //int srcPeNumber = isFragmented ? info->srcPe : env->getSrcPe();
42   //int totalSendingSize = isFragmented ? info->chunkSize+CmiReservedHeaderSize+sizeof(PipeBcastInfo) : env->getTotalsize();
43
44   switch (topology) {
45   case USE_LINEAR:
46     if (srcPeNumber == (CkMyPe()+1)%CkNumPes()) break;
47     destination = (CkMyPe()+1) % CkNumPes();
48     ComlibPrintf("[%d] Pipebroadcast sending to %d\n",CkMyPe(), destination);
49     CmiSyncSend(destination, totalSendingSize, env);
50     break;
51   case USE_HYPERCUBE:
52     tmp = srcPeNumber ^ CkMyPe();
53     k = log_of_2(CkNumPes()) + 2;
54     if (tmp) {
55       do {--k;} while (!(tmp>>k));
56     }
57     ComlibPrintf("[%d] tmp=%d, k=%d\n",CkMyPe(),tmp,k);
58     // now 'k' is the last dimension in the hypercube used for exchange
59     if (isFragmented) info->srcPe = CkMyPe();
60     else setPeNumber(env,CkMyPe());  // where the message is coming from
61     dest_pes = (int *)malloc(k*sizeof(int));
62     --k;  // next dimension in the cube to be used
63     num_pes = HypercubeGetBcastDestinations(k, dest_pes);
64
65     /*
66     for ( ; k>=0; --k) {
67       // add the processor destination at level k if it exist
68       dest_pes[num_pes] = CkMyPe() ^ (1<<k);
69       if (dest_pes[num_pes] >= CkNumPes()) {
70         dest_pes[num_pes] &= (-1)<<k;
71         if (CkNumPes()>dest_pes[num_pes]) dest_pes[num_pes] += (CkMyPe() - (CkMyPe() & ((-1)<<k))) % (CkNumPes() - dest_pes[num_pes]);
72       }
73       if (dest_pes[num_pes] < CkNumPes()) {
74         ComlibPrintf("[%d] PipeBroadcast sending to %d\n",CkMyPe(), dest_pes[num_pes]);
75         ++num_pes;
76       }
77     }
78     */
79
80     //CmiSyncListSend(num_pes, dest_pes, env->getTotalsize(), (char *)env);
81     for (k=0; k<num_pes; ++k) {
82       ComlibPrintf("[%d] PipeBroadcast sending to %d\n",CkMyPe(), dest_pes[k]);
83       CmiSyncSend(dest_pes[k], totalSendingSize, env);
84     }
85     //sizeToSend = pipeSize<totalSendingSize ? pipeSize : totalSendingSize;
86     //for (k=0; k<num_pes; ++k) CmiSyncSend(dest_pes[k], sizeToSend, env);
87     free(dest_pes);
88     break;
89
90     // for other strategies
91
92   default:
93     // should NEVER reach here!
94     char *error_msg;
95     sprintf(error_msg, "Error, topology %d not known\n",topology);
96     CmiAbort(error_msg);
97   }
98
99   // deliver messages to local objects (i.e. send it to ComlibManager)
100   storing(env, isFragmented);
101   //CmiSetHandler(env, CmiGetXHandler(env));
102   //CmiSyncSendAndFree(CkMyPe(), env->getTotalsize(), (char *)env);
103
104 }
105
106 void PipeBroadcastConverse::storing(char* fragment, int isFragmented) {
107   char *complete;
108   int isFinished=0;
109   int totalDimension;
110   //ComlibPrintf("isArray = %d\n", (getType() == ARRAY_STRATEGY));
111
112   // check if the message is fragmented
113   if (isFragmented) {
114     // store the fragment in the hash table until completed
115     ComlibPrintf("[%d] deliverer: received fragmented message, storing\n",CkMyPe());
116     PipeBcastInfo *info = (PipeBcastInfo*)(fragment+CmiReservedHeaderSize);
117
118     PipeBcastHashKey key (info->bcastPe, info->seqNumber);
119     PipeBcastHashObj *position = fragments.get(key);
120
121     char *incomingMsg;
122     if (position) {
123       // the message already exist, add to it
124       ComlibPrintf("[%d] adding to an existing message for id %d/%d (%d remaining)\n",CkMyPe(),info->bcastPe,info->seqNumber,position->remaining-1);
125       incomingMsg = position->message;
126       memcpy (incomingMsg+CmiReservedHeaderSize+((pipeSize-CmiReservedHeaderSize-sizeof(PipeBcastInfo))*info->chunkNumber), fragment+CmiReservedHeaderSize+sizeof(PipeBcastInfo), info->chunkSize);
127
128       if (--position->remaining == 0) {  // message completely received
129         isFinished = 1;
130         complete = incomingMsg;
131         totalDimension = position->dimension;
132         // delete from the hash table
133         fragments.remove(key);
134       }
135
136     } else {
137       // the message doesn't exist, create it
138       ComlibPrintf("[%d] creating new message of size %d for id %d/%d; chunk=%d chunkSize=%d\n",CkMyPe(),info->messageSize,info->bcastPe,info->seqNumber,info->chunkNumber,info->chunkSize);
139       incomingMsg = (char*)CmiAlloc(info->messageSize);
140       memcpy (incomingMsg, fragment, CmiReservedHeaderSize);
141       memcpy (incomingMsg+CmiReservedHeaderSize+((pipeSize-CmiReservedHeaderSize-sizeof(PipeBcastInfo))*info->chunkNumber), fragment+CmiReservedHeaderSize+sizeof(PipeBcastInfo), info->chunkSize);
142       int remaining = (int)ceil((double)info->messageSize/(pipeSize-CmiReservedHeaderSize-sizeof(PipeBcastInfo)))-1;
143       if (remaining) {  // more than one chunk (it was not forced to be splitted)
144         PipeBcastHashObj *object = new PipeBcastHashObj(info->messageSize, remaining, incomingMsg);
145         fragments.put(key) = object;
146       } else {  // only one chunk, it was forces to be splitted
147         isFinished = 1;
148         complete = incomingMsg;
149         // nothing to delete from fragments since nothing has been added
150       }
151     }
152     CmiFree(fragment);
153
154   } else {  // message not fragmented
155     ComlibPrintf("[%d] deliverer: received message in single chunk\n",CkMyPe());
156     isFinished = 1;
157     complete = fragment;
158   }
159
160   if (isFinished) {
161     higherLevel->deliverer(complete, totalDimension);
162   }
163 }
164
165 void PipeBroadcastConverse::deliverer(char *msg, int dimension) {
166   ComlibPrintf("{%d} dest = %d, %d, %x\n",CkMyPe(),destinationHandler, dimension,CmiHandlerToInfo(destinationHandler).hdlr);
167   if (destinationHandler) {
168     CmiSetHandler(msg, destinationHandler);
169     CmiSyncSendAndFree(CkMyPe(), dimension, msg);
170   } else {
171     CmiPrintf("[%d] Pipelined Broadcast: message not delivered since destination not set!");
172   }
173 }
174
175 PipeBroadcastConverse::PipeBroadcastConverse(int _topology, int _pipeSize, Strategy *parent) : Strategy(), topology(_topology), pipeSize(_pipeSize) {
176   if (parent) higherLevel = parent;
177   else higherLevel = this;
178   seqNumber = 0;
179   messageBuf = new CkQ<MessageHolder *>;
180   //if (!parent) propagateHandle_frag = CmiRegisterHandler((CmiHandler)propagate_handler_frag);
181   ComlibPrintf("init: %d %d (%x)\n",topology, pipeSize,this);
182   //if (!parent) ComlibPrintf("[%d] registered handler fragmented to %d\n",CkMyPe(),propagateHandle_frag);
183 }
184
185 void PipeBroadcastConverse::insertMessage(MessageHolder *cmsg){
186   ComlibPrintf("[%d] Pipelined Broadcast with converse strategy %d\n",CkMyPe(),topology);
187   messageBuf->enq(cmsg);
188   doneInserting();
189 }
190
191 void PipeBroadcastConverse::doneInserting(){
192   ComlibPrintf("[%d] DoneInserting\n",CkMyPe());
193   while (!messageBuf->isEmpty()) {
194     MessageHolder *cmsg = messageBuf->deq();
195     // modify the Handler to deliver the message to the propagator
196     char *env = cmsg->getMessage();
197     CmiSetHandler(env, CsvAccess(pipeBcastPropagateHandle_frag));
198     conversePipeBcast(env, cmsg->getSize());
199     delete cmsg;
200     //conversePipeBcast(env, env->getTotalsize(), false);
201   }
202 }
203
204 // routine for interfacing with converse.
205 // Require only the converse reserved header if forceSplit is true
206 void PipeBroadcastConverse::conversePipeBcast(char *env, int totalSize) {
207   // set the instance ID to be used by the receiver using the XHandler variable
208   CmiSetXHandler(env, myInstanceID);
209
210   ++seqNumber;
211   // message doesn't fit into the pipe: split it into chunks and propagate them individually
212   ComlibPrintf("[%d] Propagating message in multiple chunks (totalsize=%d)\n",CkMyPe(),totalSize);
213
214   char *sendingMsg;
215   char *nextChunk = env+CmiReservedHeaderSize;
216   int remaining = totalSize-CmiReservedHeaderSize;
217   int reducedPipe = pipeSize-CmiReservedHeaderSize-sizeof(PipeBcastInfo);
218   int sendingMsgSize;
219   ComlibPrintf("reducedPipe = %d, CmiReservedHeaderSize = %d, sizeof(PipeBcastInfo) = %d\n",reducedPipe,CmiReservedHeaderSize,sizeof(PipeBcastInfo));
220   ComlibPrintf("sending %d chunks of size %d, total=%d to handle %d\n",(int)ceil(((double)totalSize-CmiReservedHeaderSize)/reducedPipe),reducedPipe,remaining,CsvAccess(pipeBcastPropagateHandle_frag));
221   CmiSetHandler(env, CsvAccess(pipeBcastPropagateHandle_frag));
222   ComlibPrintf("setting env handler to %d\n",CsvAccess(pipeBcastPropagateHandle_frag));
223   for (int i=0; i<(int)ceil(((double)totalSize-CmiReservedHeaderSize)/reducedPipe); ++i) {
224     sendingMsgSize = reducedPipe<remaining? pipeSize : remaining+CmiReservedHeaderSize+sizeof(PipeBcastInfo);
225     sendingMsg = (char*)CmiAlloc(sendingMsgSize);
226     memcpy (sendingMsg, env, CmiReservedHeaderSize);
227     PipeBcastInfo *info = (PipeBcastInfo*)(sendingMsg+CmiReservedHeaderSize);
228     info->srcPe = CkMyPe();
229     info->bcastPe = CkMyPe();
230     info->seqNumber = seqNumber;
231     info->chunkNumber = i;
232     info->chunkSize = reducedPipe<remaining ? reducedPipe : remaining;
233     info->messageSize = totalSize;
234     memcpy (sendingMsg+CmiReservedHeaderSize+sizeof(PipeBcastInfo), nextChunk, info->chunkSize);
235
236     remaining -= info->chunkSize;
237     nextChunk += info->chunkSize;
238
239     propagate(sendingMsg, true, CkMyPe(), sendingMsgSize, NULL);
240   }
241   CmiFree(env);
242 }
243
244 void PipeBroadcastConverse::pup(PUP::er &p){
245   Strategy::pup(p);
246   ComlibPrintf("[%d] initial of Pipeconverse pup %s\n",CkMyPe(),(p.isPacking()==0)?(p.isUnpacking()?"UnPacking":"sizer"):("Packing"));
247
248   p | pipeSize;
249   p | topology;
250   p | seqNumber;
251
252   ComlibPrintf("[%d] PipeBroadcast converse pupping %s, size=%d, topology=%d\n",CkMyPe(), (p.isPacking()==0)?(p.isUnpacking()?"UnPacking":"sizer"):("Packing"),pipeSize,topology);
253
254   if (p.isUnpacking()) {
255     //log_of_2_inv = 1/log((double)2);
256     messageBuf = new CkQ<MessageHolder *>;
257     //propagateHandle_frag = CmiRegisterHandler((CmiHandler)propagate_handler_frag);
258     ComlibPrintf("[%d] registered handler fragmented to %d\n",CkMyPe(),CsvAccess(pipeBcastPropagateHandle_frag));
259   }
260   if (p.isPacking()) {
261     delete messageBuf;
262   }
263   //p|(*messageBuf);
264   //p|fragments;
265
266 }
267
268 PUPable_def(PipeBroadcastConverse);