/* -*- c++ -*- */
/*
 * Copyright 2012 Free Software Foundation, Inc.
 *
 * This file is part of GNU Radio
 *
 * GNU Radio is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3, or (at your option)
 * any later version.
 *
 * GNU Radio is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with GNU Radio; see the file COPYING.  If not, write to
 * the Free Software Foundation, Inc., 51 Franklin Street,
 * Boston, MA 02110-1301, USA.
 */

#ifdef HAVE_CONFIG_H
#include "config.h"
#endif

#include <gr_vector_map.h>
#include <gr_io_signature.h>
#include <string.h>

std::vector<int>
get_in_sizeofs(size_t item_size, std::vector<size_t> in_vlens)
{
  std::vector<int> in_sizeofs;
  for(unsigned int i = 0; i < in_vlens.size(); i++) {
    in_sizeofs.push_back(in_vlens[i]*item_size);
  }
  return in_sizeofs;
}

std::vector<int>
get_out_sizeofs(size_t item_size,
		std::vector< std::vector< std::vector<size_t> > > mapping)
{
  std::vector<int> out_sizeofs;
  for(unsigned int i = 0; i < mapping.size(); i++) {
    out_sizeofs.push_back(mapping[i].size()*item_size);
  }
  return out_sizeofs;
} 

gr_vector_map_sptr
gr_make_vector_map (size_t item_size, std::vector<size_t> in_vlens,
		    std::vector< std::vector< std::vector<size_t> > > mapping)
{
  return gnuradio::get_initial_sptr(new gr_vector_map(item_size,
						      in_vlens,
						      mapping));
}

gr_vector_map::gr_vector_map(size_t item_size, std::vector<size_t> in_vlens,
			     std::vector< std::vector< std::vector<size_t> > > mapping)
  : gr_sync_block("vector_map",
		  gr_make_io_signaturev(in_vlens.size(), in_vlens.size(),
					get_in_sizeofs(item_size, in_vlens)),
		  gr_make_io_signaturev(mapping.size(), mapping.size(),
					get_out_sizeofs(item_size, mapping))),
    d_item_size(item_size), d_in_vlens(in_vlens)
{
  set_mapping(mapping);
}

void
gr_vector_map::set_mapping(std::vector< std::vector< std::vector<size_t> > > mapping) {
  // Make sure the contents of the mapping vectors are possible.
  for(unsigned int i=0; i<mapping.size(); i++) {
    for(unsigned int j=0; j<mapping[i].size(); j++) {
      if(mapping[i][j].size() != 2) {
	throw std::runtime_error("Mapping must be of the form (out_mapping_stream1, out_mapping_stream2, ...), where out_mapping_stream1 is of the form (mapping_element1, mapping_element2, ...), where mapping_element1 is of the form (input_stream, input_element).  This error is raised because a mapping_element vector does not contain exactly 2 items.");
      }
      unsigned int s = mapping[i][j][0];
      unsigned int index = mapping[i][j][1];
      if(s >= d_in_vlens.size()) {
	throw std::runtime_error("Stream numbers in mapping must be less than the number of input streams.");
      }
      if((index < 0) || (index >= d_in_vlens[s])) {
	throw std::runtime_error ("Indices in mapping must be greater than 0 and less than the input vector lengths.");
      }
    }
  }
  gruel::scoped_lock guard(d_mutex);
  d_mapping = mapping;
}

int
gr_vector_map::work(int noutput_items,
		    gr_vector_const_void_star &input_items,
		    gr_vector_void_star &output_items)
{
  const char **inv = (const char **) &input_items[0];
  char **outv = (char **) &output_items[0];

  for(unsigned int n = 0; n < (unsigned int)(noutput_items); n++) {
    for(unsigned int i = 0; i < d_mapping.size(); i++) {
      unsigned int out_vlen = d_mapping[i].size();
      for(unsigned int j = 0; j < out_vlen; j++) {
	unsigned int s = d_mapping[i][j][0];
	unsigned int k = d_mapping[i][j][1];
	memcpy(outv[i] + out_vlen*d_item_size*n +
	       d_item_size*j, inv[s] + d_in_vlens[s]*d_item_size*n +
	       k*d_item_size, d_item_size);
      }
    }
  }

  return noutput_items;
}