a little cleanup and conformation of usage of CmiChunkHeader. It will prevent in
[charm.git] / src / conv-com / petable.C
1 /*****************************************************************************
2  * $Source$
3  * $Author$
4  * $Date$
5  * $Revision$
6  *****************************************************************************/
7
8 /*********************************************
9  * File : petable.C
10  *
11  * Author: Krishnan V
12  *
13  * The message buffer
14  *********************************************/
15 #include <string.h>
16 #include <stdlib.h>
17 #include <converse.h>
18 #include "comlib.h"
19 #include "petable.h"
20 #include "converse.h"
21
22 #define BIGBUFFERSIZE 65536
23 #define PTPREALLOC    100
24
25 struct AllToAllHdr{
26     char dummy[CmiReservedHeaderSize];
27     int refno;
28     comID id;
29     int ufield;
30     int nmsgs;
31 };
32
33
34 /**************************************************************
35  * Preallocated memory=P*MSGQLEN ptr + 2P ints + 1000 ptrs
36  **************************************************************/
37 PeTable :: PeTable(int n)
38 {
39   NumPes=n;
40   magic=0;
41   PeList = (PTinfo ***)CmiAlloc(sizeof(PTinfo *)*NumPes);
42   //  ComlibPrintf("Pelist[%d][%d]\n", NumPes, MSGQLEN);
43   msgnum=new int[NumPes];
44   MaxSize=new int[NumPes];
45   for (int i=0;i<NumPes;i++) {
46         msgnum[i]=0;
47         MaxSize[i]=MSGQLEN;
48         PeList[i]=(PTinfo **)CmiAlloc(sizeof(PTinfo *)*MSGQLEN);
49         for (int j=0;j<MSGQLEN;j++) PeList[i][j]=0;
50   }
51
52   //ptrlist=(PTinfo **)CmiAlloc(1000*sizeof(PTinfo *));
53   //  FreeList= new GList;
54   //CombBuffer=(char *)CmiAlloc(BIGBUFFERSIZE);
55
56   PTFreeList=NULL;
57   PTFreeChunks=NULL;
58 }
59
60 PeTable :: ~PeTable()
61 {
62   int i;
63   for (i=0;i<NumPes;i++) CmiFree(PeList[i]);
64   CmiFree(PeList);
65   delete msgnum;
66   delete MaxSize;
67   GarbageCollect();
68   //CmiFree(ptrlist);
69   PTinfo *tmp;
70   while (PTFreeChunks) {
71         tmp=PTFreeChunks;
72         PTFreeChunks=PTNEXTCHUNK(tmp);
73         CmiFree(tmp);
74   }
75  // delete FreeList;
76
77 }
78
79 void PeTable:: Purge()
80 {
81   for (int i=0; i<NumPes;i++) {
82         if (msgnum[i]) {
83             // ComlibPrintf("%d Warning: %d Undelivered Messages for %d\n", CkMyPe(), msgnum[i], i);
84           //msgnum[i]=0;
85         }
86   }
87   GarbageCollect();
88   //  ComlibPrintf("combcount = %d\n", combcount);
89   //combcount = 0;
90 }
91
92 void PeTable :: ExtractAndDeliverLocalMsgs(int index)
93 {
94   int j;
95   msgstruct m;
96
97   ComlibPrintf("%d:Delivering %d local messages\n", CkMyPe(), msgnum[index]);
98   for (j=msgnum[index]-1;(j>=0);j--) {
99
100         m.msgsize=PeList[index][j]->msgsize;
101         m.msg=PeList[index][j]->msg;
102
103         if (--(PeList[index][j]->refCount) <=0) {
104             CmiSyncSendAndFree(CkMyPe()/*index*/, m.msgsize, (char*)m.msg);
105             PTFREE(PeList[index][j]);
106         }
107         else {
108             CmiSyncSend(CkMyPe()/*index*/, m.msgsize, (char*)m.msg);
109         }
110         PeList[index][j]=NULL;
111   }
112   msgnum[index]=j+1;
113
114   return;
115 }
116
117
118 #undef PACK
119 #undef PACKMSG
120 //#define PACKINT(data) {((int*)t)[0] = data; t+=sizeof(int);}
121 #define PACK(type,data) {junk=(char *)&(data); memcpy(t, junk, sizeof(type)); t+=sizeof(type);}
122 #define PACKMSG(data, size) { memcpy(p+msg_offset, (data), size); msg_offset += size; }
123
124 /*Used for all to all multicast operations.  Assumes that each message
125   is destined to all the processors, to speed up all to all
126   substantially --Sameer 09/03/03 
127   
128   Protocol:
129   |ref|comid|ufield|nmsgs|size|ref|msg1|size2|ref2|msg2|....
130 */
131
132 char * PeTable ::ExtractAndPackAll(comID id, int ufield, int *length)
133 {
134     int nmsgs = 0, i, j;
135     int index = 0;
136
137     ComlibPrintf("[%d] In Extract And Pack All\n", CkMyPe());
138
139     //Increment magic to detect duplicate messages
140     magic++;
141
142     register int total_msg_size = 0;
143
144     //first compute size
145     for (i=0;i<NumPes;i++) {
146         index = i;
147         for (j=msgnum[index]-1; (j>=0); j--) {
148             if (PeList[index][j]->magic != magic) {                
149                 total_msg_size += ALIGN8(PeList[index][j]->msgsize);
150                 total_msg_size += 2 * sizeof(int);
151                 PeList[index][j]->magic=magic;
152
153                 nmsgs ++;
154             }            
155         }
156     }
157     
158     total_msg_size += ALIGN8(sizeof(AllToAllHdr));
159
160     ComlibPrintf("[%d] Message Size %d, nmsgs %d **%d**\n", CkMyPe(), total_msg_size, nmsgs, sizeof(AllToAllHdr));
161     
162     //poiter to the combined message, UGLY NAME
163     char *p = (char *) CmiAlloc(total_msg_size * sizeof(char));    
164
165     ComlibPrintf("After cmialloc\n");
166
167     //buffer to copy stuff into
168     char *t = p; 
169     char *junk = NULL;
170     
171     int dummy = 0;
172     
173     int refno = 0;
174
175     AllToAllHdr ahdr;
176     ahdr.refno = refno;
177     ahdr.id = id;
178     ahdr.ufield = ufield;
179     ahdr.nmsgs = nmsgs;
180
181     /*
182       PACKINT(refno);    
183       PACK(comID, id);
184       
185       PACKINT(ufield);
186       PACKINT(nmsgs);
187       //    PACKINT(dummy); //Aligning to 8 bytes
188     */
189
190     PACK(AllToAllHdr, ahdr);   
191
192     int msg_offset = ALIGN8(sizeof(AllToAllHdr));
193     
194     //Increment magic again for creating the message
195     magic++;
196     for (i=0;i<NumPes;i++) {
197         index=i;
198         int ref = 0;
199         int size;
200
201         for (j=msgnum[index]-1; (j>=0); j--) {
202             //Check if it is a duplicate
203             if (PeList[index][j]->magic != magic) {                
204                 size = PeList[index][j]->msgsize;
205                 PACKMSG(&size, sizeof(int));
206                 PACKMSG(&ref, sizeof(int));
207                 PeList[index][j]->magic=magic;
208                 PACKMSG(PeList[index][j]->msg, size);
209
210                 msg_offset = ALIGN8(msg_offset);
211             }
212
213             //Free it when all the processors have gotten rid of it
214             if (--(PeList[index][j]->refCount) <=0) {
215                 ComlibPrintf("before cmifree \n");
216                 CmiFree(PeList[index][j]->msg);   
217                 ComlibPrintf("after cmifree \n");
218
219                 PTFREE(PeList[index][j]);
220             }
221             //Assign the current processors message pointer to NULL
222             PeList[index][j] = NULL;
223         }
224         msgnum[index] = 0;
225     }
226     
227     *length = total_msg_size;
228     return p;
229 }
230
231 char * PeTable ::ExtractAndPack(comID id, int ufield, int npe, 
232                                 int *pelist, int *length)
233 {
234     char *junk;
235     int nummsgs, offset, num_distinctmsgs;
236         
237     int tot_msgsize=TotalMsgSize(npe, pelist, &nummsgs, &num_distinctmsgs);
238
239     ComlibPrintf("%d In ExtractAndPack %d, %d\n", CmiMyPe(), npe, nummsgs); 
240
241     if (tot_msgsize ==0) {
242         *length=0;
243         
244         ComlibPrintf("Returning NULL\n");
245         return(NULL);
246     }
247     
248     int msg_offset = CmiReservedHeaderSize + sizeof(comID) 
249         + (npe + 4 + nummsgs) * sizeof(int);  
250
251     msg_offset = ALIGN8(msg_offset);
252     
253     *length = tot_msgsize;
254     *length += msg_offset;
255     char *p;
256     p=(char *)CmiAlloc(*length);
257
258     char *t = p + CmiReservedHeaderSize;
259     int i, j;
260     if (!p) CmiAbort("Big time problem\n");
261     magic++;
262
263     int refno = id.refno;    
264
265     PACK(int, refno);
266     PACK(comID, id);
267     PACK(int, ufield);
268     PACK(int, npe);
269     
270     int lesspe=0;
271     int npacked = 0;
272     for (i=0;i<npe;i++) {
273         int index=pelist[i];
274
275         if (msgnum[index]<=0) {
276             lesspe++;
277             
278             ComlibPrintf("[%d] msgnum[index]<=0 !!!!!\n", CkMyPe());
279             continue;
280         }
281         
282         ComlibPrintf("%d Packing pelist[%d]\n", CkMyPe(), index);
283         register int newval=-1*pelist[i];
284         PACK(int, newval); 
285
286         for (j=0;j<msgnum[index];j++) {
287             if (PeList[index][j]->magic == magic) {
288                 offset=(PeList[index][j]->offset);
289             }
290             else {
291                 npacked ++;
292                 
293                 offset=msg_offset;
294                 PeList[index][j]->magic=magic;
295                 PeList[index][j]->offset=msg_offset;
296                 PTinfo *tempmsg=PeList[index][j];
297                 
298                 CmiChunkHeader hdr;
299                 hdr.size = tempmsg->msgsize;
300                 hdr.ref = -1;
301                 PACKMSG(&hdr, sizeof(CmiChunkHeader));
302                 PACKMSG(tempmsg->msg, tempmsg->msgsize);
303
304                 msg_offset = ALIGN8(msg_offset);
305             }
306             
307             //ComlibPrintf("%d Packing msg_offset=%d\n", CkMyPe(), offset);
308             PACK(int, offset); 
309
310             if (--(PeList[index][j]->refCount) <=0) {
311                 CmiFree(PeList[index][j]->msg);
312                 
313                 PTFREE(PeList[index][j]);
314             }
315             PeList[index][j]=NULL;
316         }
317         msgnum[index]=0;
318     }
319     offset=-1;
320     PACK(int, offset);
321     
322     if (lesspe) {
323         t=p+CmiReservedHeaderSize+2*sizeof(int) + sizeof(comID);
324         npe=npe-lesspe;
325         PACK(int, npe);
326     }
327
328     return(p);
329
330
331 #undef UNPACK
332 #define UNPACK(type,data) {junk=(char *)&(data); memcpy(junk, t, sizeof(type));t+=sizeof(type);}
333 #undef UNPACKMSG
334 #define UNPACKMSG(dest,src, size) { memcpy(dest, src, size); offset += size;}
335
336 int PeTable :: UnpackAndInsert(void *in)
337 {
338   char *junk;
339   char *t =(char *)in + CmiReservedHeaderSize;
340   int i, ufield, npe, pe, tot_msgsize, ptrlistindex=0;
341   comID id;
342   int refno = 0;
343
344   UNPACK(int, refno);
345   
346   //ComlibPrintf("%d UnPacking id\n", CkMyPe());
347   UNPACK(comID, id);
348   UNPACK(int, ufield);
349   UNPACK(int, npe);
350   
351   register int offset;
352   for (i=0;i<npe;i++) {
353         UNPACK(int, pe);
354         pe *= -1;
355
356         UNPACK(int, offset);
357         while (offset > 0) {
358             int tempmsgsize;
359             UNPACKMSG(&(tempmsgsize), (char *)in+offset, sizeof(int));
360             int ptr;
361             UNPACKMSG(&ptr, (char *)in+offset, sizeof(int));
362
363             if (ptr >=0 )  {
364                 if (msgnum[pe] >= MaxSize[pe]) {
365                     REALLOC(PeList[pe], MaxSize[pe]);
366                     MaxSize[pe] *= 2;
367                 }
368                 PeList[pe][msgnum[pe]]=ptrvec[ptr];
369                 (ptrvec[ptr])->refCount++;
370                 msgnum[pe]++;
371
372                 UNPACK(int, offset);
373                 continue;
374             }
375             
376             PTinfo *temp;
377             PTALLOC(temp);
378             temp->msgsize=tempmsgsize;
379             temp->refCount=1;
380             temp->magic=0;
381             temp->offset=0;
382
383             ptrvec.insert(ptrlistindex, temp);
384             memcpy((char *)in+offset-sizeof(int), &ptrlistindex, sizeof(int));
385
386             ptrlistindex++;
387             temp->msg=(void *)((char *)in+offset);
388             if (msgnum[pe] >= MaxSize[pe]) {
389
390                 REALLOC(PeList[pe], MaxSize[pe]);
391                 MaxSize[pe] *= 2;
392             }
393             PeList[pe][msgnum[pe]]=temp;
394             msgnum[pe]++;
395             UNPACK(int, offset);
396         }
397         t -=sizeof(int);
398   }
399   *(int *)((char *)in -sizeof(int))=ptrlistindex; 
400   
401   if (ptrlistindex==0)
402       CmiFree(in);
403   
404   for (i=0;i<ptrlistindex;i++) {
405       char * actualmsg=(char *)(ptrvec[i]->msg);
406       int *rc=(int *)(actualmsg-sizeof(int));
407       *rc=(int)((char *)in-actualmsg);
408       //ComlibPrintf("I am inserting %d\n", *rc);
409   }
410
411   ptrvec.removeAll();
412   
413   return(ufield);
414 }
415
416 /* Unpack and insert an all to all message, the router provides the
417    list of processors to insert into.
418    Same protocol as mentioned earlier.
419 */
420
421 int PeTable :: UnpackAndInsertAll(void *in, int npes, int *pelist){
422   char *junk;
423   char *t =(char *)in /*+CmiReservedHeaderSize*/;
424   int i,  
425       ufield,   // user field or ths stage of the iteration 
426       nmsgs,    // number of messages in combo message
427       refno,    // reference number
428       dummy;    // alignment dummy
429   
430   comID id;
431
432   /*
433     UNPACK(int, refno);      
434     UNPACK(comID, id);
435     
436     UNPACK(int, ufield);
437     UNPACK(int, nmsgs);
438     //UNPACK(int, dummy);
439     int header_size = sizeof(comID) + CmiReservedHeaderSize + 3 *sizeof(int);
440     if(header_size % 8 != 0)
441     t+= 8 - header_size % 8;
442   */
443
444   AllToAllHdr ahdr;
445   UNPACK(AllToAllHdr, ahdr);
446
447   if(sizeof(AllToAllHdr) & 7 != 0)
448       t += 8 - sizeof(AllToAllHdr) & 7;
449
450   refno = ahdr.refno;
451   id = ahdr.id;
452   nmsgs = ahdr.nmsgs;
453   ufield = ahdr.ufield;
454
455   ComlibPrintf("[%d] unpack and insert all %d, %d\n", CkMyPe(), ufield, nmsgs);
456   
457   //Inserting a memory foot print may, change later
458   CmiChunkHeader *chdr= (CmiChunkHeader*)((char*)in - sizeof(CmiChunkHeader));
459
460   int *ref;
461   int size;
462   char *msg;
463   for(int count = 0; count < nmsgs; count++){
464
465     t += sizeof(CmiChunkHeader);
466
467     msg = t;
468     t += ALIGN8(size);
469
470     // Get the size of the message, and set the ref field correctly for CmiFree
471     size = SIZEFIELD(msg);
472     REFFIELD(msg) = (int)((char *)&REFFIELD(in) - (char *)REFFIELD(msg));
473
474     // Do CmiReference(msg), this is done bypassing converse!
475     chdr->ref++;
476
477     /*
478     UNPACK(int, size);
479     ref = (int *)t;
480     t += sizeof(int);
481     
482     *ref = (int)((char *)(&chdr->ref) - (char *)ref);
483     chdr->ref ++;
484
485     ComlibPrintf("ref = %d, global_ref = %d\n", *ref, chdr->ref);
486     
487     msg = t;
488     t += ALIGN8(size);
489     */
490     InsertMsgs(npes, pelist, size, msg);
491   }  
492
493   CmiFree(in);
494   return ufield;
495 }
496
497 PTvectorlist PeTable :: ExtractAndVectorize(comID id, int ufield, int npe, int *pelist) {
498
499 }
500
501 PTvectorlist PeTable :: ExtractAndVectorizeAll(comID id, int ufield) {
502
503 }
504
505 void PeTable :: GarbageCollect()
506 {
507 }
508