#!/usr/bin/env python

from mapping import *

##################################################################################
#   file   : transpose.py  
#   date   : may 2014
#   author : Alain Greiner
##################################################################################
#  This file describes the mapping of the multi-threaded "transpose" 
#  application on a multi-clusters, multi-processors architecture.
#  There is one thread per processor.
#  The mapping of virtual segments is the following:
#    - There is one shared data vseg in cluster[0][0]
#    - The code vsegs are replicated on all clusters containing processors.
#    - There is one heap vseg per cluster containing processors.
#    - The stacks vsegs are distibuted on all clusters containing processors.
#  This mapping uses 5 platform parameters, (obtained from the "mapping" argument)
#    - x_size    : number of clusters in a row
#    - y_size    : number of clusters in a column
#    - x_width   : number of bits coding x coordinate
#    - y_width   : number of bits coding y coordinate
#    - nprocs    : number of processors per cluster
##################################################################################

######################
def extend( mapping ):

    x_size    = mapping.x_size
    y_size    = mapping.y_size
    nprocs    = mapping.nprocs
    x_width   = mapping.x_width
    y_width   = mapping.y_width

    # define vsegs base & size
    code_base  = 0x10000000
    code_size  = 0x00010000     # 64 Kbytes  (256 Mbytes max)
    
    data_base  = 0x20000000
    data_size  = 0x00010000     # 64 Kbytes  (256 Mbytes max)

    mmap_base  = 0x30000000 
    mmap_size  = 0x10000000     # 256 Mbytes (non mapped) 

    stack_base = 0x40000000 
    stack_size = 0x00010000     # 64 Kbytes per thread  (64 Mbytes max)

    heap_base  = 0x60000000 
    heap_size  = 0x00200000     # 2 Mbytes per cluster  (512 Mbytes max)

    # create vspace
    vspace = mapping.addVspace( name = 'transpose', startname = 'trsp_data', active = False )
    
    # data vseg : shared (only in cluster[0,0])
    mapping.addVseg( vspace, 'trsp_data', data_base , data_size, 
                     'C_WU', vtype = 'ELF', x = 0, y = 0, pseg = 'RAM', 
                     binpath = 'bin/transpose/appli.elf',
                     local = False )

    # mmap vseg : non mapped in physical memory
    mapping.addVseg( vspace, 'trsp_mmap', mmap_base , mmap_size, 
                     'C_WU', vtype = 'MMAP', local = False )

    # code vsegs : local (one copy in each cluster)
    for x in xrange (x_size):
        for y in xrange (y_size):
            cluster_id = (x * y_size) + y
            if ( mapping.clusters[cluster_id].procs ):

                mapping.addVseg( vspace, 'trsp_code_%d_%d' %(x,y), 
                                 code_base , code_size,
                                 'CXWU', vtype = 'ELF', x = x, y = y, pseg = 'RAM', 
                                 binpath = 'bin/transpose/appli.elf',
                                 local = True )

    # stacks vsegs: local (one stack per processor => nprocs stacks per cluster)
    for x in xrange (x_size):
        for y in xrange (y_size):
            cluster_id = (x * y_size) + y
            if ( mapping.clusters[cluster_id].procs ):
                for p in xrange( nprocs ):
                    proc_id = (((x * y_size) + y) * nprocs) + p
                    base    = stack_base + (proc_id * stack_size)
                    size    = stack_size - 4096
                    mapping.addVseg( vspace, 'trsp_stack_%d_%d_%d' % (x,y,p), 
                                     base , size , 'C_WU', vtype = 'BUFFER', 
                                     x = x , y = y , pseg = 'RAM',
                                     local = True )

    # heap vsegs: distributed non local (all heap vsegs can be accessed by all tasks)
    for x in xrange (x_size):
        for y in xrange (y_size):
            cluster_id = (x * y_size) + y
            if ( mapping.clusters[cluster_id].procs ):
                size  = heap_size
                base  = heap_base + (cluster_id * size)

                mapping.addVseg( vspace, 'trsp_heap_%d_%d' % (x,y), base, size, 
                                 'C_WU', vtype = 'HEAP', x = x, y = y, pseg = 'RAM',
                                 local = False, big = True )

    # distribute one thread per processor / main on P[0,0,0]
    for x in xrange (x_size):
        for y in xrange (y_size):
            cluster_id = (x * y_size) + y
            if ( mapping.clusters[cluster_id].procs ):
                for p in xrange( nprocs ):
                    if (x == 0) and (y == 0) and (p == 0) :   # main thread
                        startid = 1
                        is_main = True
                    else :                                    # other threads
                        startid = 0
                        is_main = False

                    mapping.addThread( vspace,
                                       'trsp_%d_%d_%d' % (x,y,p),
                                       is_main,
                                       x, y, p,
                                       'trsp_stack_%d_%d_%d' % (x,y,p),
                                       'trsp_heap_%d_%d' % (x,y),
                                       startid )

    # extend mapping name
    mapping.name += '_transpose'

    return vspace  # useful for test
            
################################ test ##################################################

if __name__ == '__main__':

    vspace = extend( Mapping( 'test', 2, 2, 4 ) )
    print vspace.xml()


# Local Variables:
# tab-width: 4;
# c-basic-offset: 4;
# c-file-offsets:((innamespace . 0)(inline-open . 0));
# indent-tabs-mode: nil;
# End:
#
# vim: filetype=python:expandtab:shiftwidth=4:tabstop=4:softtabstop=4