def initialize_pose_cell_network(self, pc_shape):

        pcnet = PoseCellNetwork(pc_shape)
        xy_pc = pc_shape[0]/2+1 
        th_pc = pc_shape[2]/2+1 
        
        pcnet.posecells[xy_pc, xy_pc, th_pc] = 1

        v_temp = self.view_templates[0]
        pcnet.update(self.visual_odometer.delta, v_temp)
        #pcmax = pcnet.get_pc_max(pcnet.avg_xywrap, pcnet.avg_thwrap)

        return pcnet
    def initialize_pose_cell_network(self, pc_shape):

        pcnet = PoseCellNetwork(pc_shape)
        xy_pc = pc_shape[0]/2+1 
        th_pc = pc_shape[2]/2+1 
        
        pcnet.posecells[xy_pc, xy_pc, th_pc] = 1 # initalize the matrix with shape

        v_temp = self.view_templates[0]
        pcnet.update(self.visual_odometer.delta, v_temp) #determines the effect of visual template on poe cell network
        #pcmax = pcnet.get_pc_max(pcnet.avg_xywrap, pcnet.avg_thwrap)
		# I think so it doesn't plot all the activites because it only reutrns the mzximum value of pose cell index 
		# and we need to get it return all the pose cell activity so that I can understand the pattern.
        return pcnet