#
# Copyright 2005,2006,2007,2009,2011 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.
#

# See gnuradio-examples/python/digital for examples

"""
Generic modulation and demodulation.
"""

from gnuradio import gr
from gnuradio.digital.modulation_utils import extract_kwargs_from_options_for_class
from gnuradio.digital.utils import mod_codes
from gnuradio.digital import digital_swig,digital
import math

# default values (used in __init__ and add_options)
_def_samples_per_symbol = 2
_def_excess_bw = 0.35
_def_verbose = False
_def_log = False

# Frequency correction
_def_freq_bw = 2*math.pi/100.0
# Symbol timing recovery
_def_timing_bw = 2*math.pi/100.0
_def_timing_max_dev = 1.5
# Fine frequency / Phase correction
_def_phase_bw = 2*math.pi/100.0
# Number of points in constellation
_def_constellation_points = 16
# Whether differential coding is used.
_def_differential = False

def add_common_options(parser):
    """
    Sets options common to both modulator and demodulator.
    """
    parser.add_option("-p", "--constellation-points", type="int", default=_def_constellation_points,
                      help="set the number of constellation points (must be a power of 2 for psk, power of 4 for QAM) [default=%default]")
    parser.add_option("", "--non-differential", action="store_false",
                      dest="differential",
                      help="do not use differential encoding [default=False]")
    parser.add_option("", "--differential", action="store_true",
                      dest="differential", default=True,
                      help="use differential encoding [default=%default]")
    parser.add_option("", "--mod-code", type="choice", choices=mod_codes.codes,
                      default=mod_codes.NO_CODE,
                      help="Select modulation code from: %s [default=%%default]"
                            % (', '.join(mod_codes.codes),))
    parser.add_option("", "--excess-bw", type="float", default=_def_excess_bw,
                      help="set RRC excess bandwith factor [default=%default]")


# /////////////////////////////////////////////////////////////////////////////
#                             Generic modulator
# /////////////////////////////////////////////////////////////////////////////

class generic_mod(gr.hier_block2):

    def __init__(self, constellation,
                 samples_per_symbol=_def_samples_per_symbol,
                 differential=_def_differential,
                 excess_bw=_def_excess_bw,
                 gray_coded=True,
                 verbose=_def_verbose,
                 log=_def_log):
        """
    Hierarchical block for RRC-filtered differential generic modulation.

    The input is a byte stream (unsigned char) and the
    output is the complex modulated signal at baseband.

    @param constellation: determines the modulation type
    @type constellation: gnuradio.digital.gr_constellation
    @param samples_per_symbol: samples per baud >= 2
    @type samples_per_symbol: float
    @param excess_bw: Root-raised cosine filter excess bandwidth
    @type excess_bw: float
        @param gray_coded: turn gray coding on/off
        @type gray_coded: bool
        @param verbose: Print information about modulator?
        @type verbose: bool
        @param log: Log modulation data to files?
        @type log: bool
    """

    gr.hier_block2.__init__(self, "generic_mod",
                gr.io_signature(1, 1, gr.sizeof_char),       # Input signature
                gr.io_signature(1, 1, gr.sizeof_gr_complex)) # Output signature

        self._constellation = constellation.base()
        self._samples_per_symbol = samples_per_symbol
        self._excess_bw = excess_bw
        self._differential = differential

        #if self._samples_per_symbol < 2:
            #raise TypeError, ("sbp must be >= 2, is %f" % self._samples_per_symbol)

        arity = pow(2,self.bits_per_symbol())

        # turn bytes into k-bit vectors
        self.bytes2chunks = \
          gr.packed_to_unpacked_bb(self.bits_per_symbol(), gr.GR_MSB_FIRST)

        if gray_coded == True:
            self.symbol_mapper = gr.map_bb(self._constellation.pre_diff_code())

        #if differential:
            #self.diffenc = gr.diff_encoder_bb(arity)

        self.chunks2symbols = gr.chunks_to_symbols_bc(self._constellation.points())
    print self._constellation.points()
    #self.constell = digital.digital_constellation(self._constellation.points())

        ## pulse shaping filter
        #nfilts = 32
        #ntaps = nfilts * 11 * int(self._samples_per_symbol)    # make nfilts filters of ntaps each
        #self.rrc_taps = gr.firdes.root_raised_cosine(
            #nfilts,          # gain
            #nfilts,          # sampling rate based on 32 filters in resampler
            #1.0,             # symbol rate
            #self._excess_bw, # excess bandwidth (roll-off factor)
            #ntaps)
        #self.rrc_filter = gr.pfb_arb_resampler_ccf(self._samples_per_symbol,
                                                   #self.rrc_taps)

    # Connect
        blocks = [self, self.bytes2chunks]
        if gray_coded == True:
            blocks.append(self.symbol_mapper)
        #if differential:
            #blocks.append(self.diffenc)
        #blocks += [self.chunks2symbols, self.rrc_filter, self]
    blocks += [self.chunks2symbols, self]
        self.connect(*blocks)

        if verbose:
            self._print_verbage()

        if log:
            self._setup_logging()


    def samples_per_symbol(self):
        return self._samples_per_symbol

    def bits_per_symbol(self):   # static method that's also callable on an instance
        return self._constellation.bits_per_symbol()

    def add_options(parser):
        """
        Adds generic modulation options to the standard parser
        """
        add_common_options(parser)
    add_options=staticmethod(add_options)

    def extract_kwargs_from_options(cls, options):
        """
        Given command line options, create dictionary suitable for passing to __init__
        """
        return extract_kwargs_from_options_for_class(cls, options)
    extract_kwargs_from_options=classmethod(extract_kwargs_from_options)


    def _print_verbage(self):
        print "\nModulator:"
        print "bits per symbol:     %d" % self.bits_per_symbol()
        print "RRC roll-off factor: %.2f" % self._excess_bw

    def _setup_logging(self):
        print "Modulation logging turned on."
        self.connect(self.bytes2chunks,
                     gr.file_sink(gr.sizeof_char, "tx_bytes2chunks.8b"))
        if self._constellation.apply_pre_diff_code():
            self.connect(self.symbol_mapper,
                         gr.file_sink(gr.sizeof_char, "tx_symbol_mapper.8b"))
        if self._differential:
            self.connect(self.diffenc,
                         gr.file_sink(gr.sizeof_char, "tx_diffenc.8b"))
        self.connect(self.chunks2symbols,
                     gr.file_sink(gr.sizeof_gr_complex, "tx_chunks2symbols.32fc"))
        self.connect(self.rrc_filter,
                     gr.file_sink(gr.sizeof_gr_complex, "tx_rrc_filter.32fc"))


# /////////////////////////////////////////////////////////////////////////////
#                             Generic demodulator
#
#      Differentially coherent detection of differentially encoded generically
#      modulated signal.
# /////////////////////////////////////////////////////////////////////////////

class generic_demod(gr.hier_block2):

    def __init__(self, constellation,
                 samples_per_symbol=_def_samples_per_symbol,
                 differential=_def_differential,
                 excess_bw=_def_excess_bw,
                 gray_coded=True,
                 freq_bw=_def_freq_bw,
                 timing_bw=_def_timing_bw,
                 phase_bw=_def_phase_bw,
                 verbose=_def_verbose,
                 log=_def_log):
        """
    Hierarchical block for RRC-filtered differential generic demodulation.

    The input is the complex modulated signal at baseband.
    The output is a stream of bits packed 1 bit per byte (LSB)

    @param constellation: determines the modulation type
    @type constellation: gnuradio.digital.gr_constellation
    @param samples_per_symbol: samples per symbol >= 2
    @type samples_per_symbol: float
    @param excess_bw: Root-raised cosine filter excess bandwidth
    @type excess_bw: float
        @param gray_coded: turn gray coding on/off
        @type gray_coded: bool
        @param freq_bw: loop filter lock-in bandwidth
        @type freq_bw: float
        @param timing_bw: timing recovery loop lock-in bandwidth
        @type timing_bw: float
        @param phase_bw: phase recovery loop bandwidth
        @type phase_bw: float
        @param verbose: Print information about modulator?
        @type verbose: bool
        @param debug: Print modualtion data to files?
        @type debug: bool
    """

    gr.hier_block2.__init__(self, "generic_demod",
                gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
                gr.io_signature(1, 1, gr.sizeof_char))       # Output signature

        self._constellation = constellation.base()
        self._samples_per_symbol = samples_per_symbol
        self._excess_bw = excess_bw
        self._phase_bw = phase_bw
        self._freq_bw = freq_bw
        self._timing_bw = timing_bw
        self._timing_max_dev= _def_timing_max_dev
        self._differential = differential

        #if self._samples_per_symbol < 2:
            #raise TypeError, ("sbp must be >= 2, is %d" % self._samples_per_symbol)

        arity = pow(2,self.bits_per_symbol())

        #nfilts = 32
        #ntaps = 11 * int(self._samples_per_symbol*nfilts)

        # Automatic gain control
        #self.agc = gr.agc2_cc(0.6e-1, 1e-3, 1, 1, 100)

        ## Frequency correction
        #fll_ntaps = 55
        #self.freq_recov = digital_swig.fll_band_edge_cc(self._samples_per_symbol, self._excess_bw,
                                                        #fll_ntaps, self._freq_bw)

        # symbol timing recovery with RRC data filter
        #taps = gr.firdes.root_raised_cosine(nfilts, nfilts*self._samples_per_symbol,
                                            #1.0, self._excess_bw, ntaps)
        #self.time_recov = gr.pfb_clock_sync_ccf(self._samples_per_symbol,
                                                #self._timing_bw, taps,
                                                #nfilts, nfilts//2, self._timing_max_dev)

        fmin = -0.25
        fmax = 0.25

        self.receiver = digital_swig.constellation_receiver_cb(
            self._constellation, self._phase_bw,
            fmin, fmax)

        self.receiver2 = digital_swig.constellation_decoder_cb(self._constellation)
        print self._constellation
        # Do differential decoding based on phase change of symbols
        #if differential:
            #self.diffdec = gr.diff_decoder_bb(arity)

        if gray_coded:
            self.symbol_mapper = gr.map_bb(
                mod_codes.invert_code(self._constellation.pre_diff_code()))

        # unpack the k bit vector into a stream of bits
        self.unpack = gr.unpack_k_bits_bb(self.bits_per_symbol()) #ber mode
    self.unpack2 = gr.unpacked_to_packed_bb(self.bits_per_symbol(), gr.GR_MSB_FIRST)

        #if verbose:
            #self._print_verbage()

        #if log:
            #self._setup_logging()

        # Connect and Initialize base class
        #blocks = [self, self.agc, self.freq_recov,
                  #self.time_recov, self.receiver]
        blocks = [self, self.receiver2,self.symbol_mapper,self]
        #if differential:
            #blocks.append(self.diffdec)
        #if self._constellation.apply_pre_diff_code():
            #blocks.append(self.symbol_mapper)
        #blocks += [self]
        self.connect(*blocks)

    def samples_per_symbol(self):
        return self._samples_per_symbol

    def bits_per_symbol(self):   # staticmethod that's also callable on an instance
        return self._constellation.bits_per_symbol()

    #def _print_verbage(self):
        #print "\nDemodulator:"
        #print "bits per symbol:     %d"   % self.bits_per_symbol()
        #print "RRC roll-off factor: %.2f" % self._excess_bw
        #print "FLL bandwidth:       %.2e" % self._freq_bw
        #print "Timing bandwidth:    %.2e" % self._timing_bw
        #print "Phase bandwidth:     %.2e" % self._phase_bw

    #def _setup_logging(self):
        #print "Modulation logging turned on."
        #self.connect(self.agc,
                     #gr.file_sink(gr.sizeof_gr_complex, "rx_agc.32fc"))
        #self.connect((self.freq_recov, 0),
                     #gr.file_sink(gr.sizeof_gr_complex, "rx_freq_recov.32fc"))
        #self.connect((self.freq_recov, 1),
                     #gr.file_sink(gr.sizeof_float, "rx_freq_recov_freq.32f"))
        #self.connect((self.freq_recov, 2),
                     #gr.file_sink(gr.sizeof_float, "rx_freq_recov_phase.32f"))
        #self.connect((self.freq_recov, 3),
                     #gr.file_sink(gr.sizeof_float, "rx_freq_recov_error.32f"))
        #self.connect((self.time_recov, 0),
                     #gr.file_sink(gr.sizeof_gr_complex, "rx_time_recov.32fc"))
        #self.connect((self.time_recov, 1),
                     #gr.file_sink(gr.sizeof_float, "rx_time_recov_error.32f"))
        #self.connect((self.time_recov, 2),
                     #gr.file_sink(gr.sizeof_float, "rx_time_recov_rate.32f"))
        #self.connect((self.time_recov, 3),
                     #gr.file_sink(gr.sizeof_float, "rx_time_recov_phase.32f"))
        #self.connect((self.receiver, 0),
                     #gr.file_sink(gr.sizeof_char, "rx_receiver.8b"))
        #self.connect((self.receiver, 1),
                     #gr.file_sink(gr.sizeof_float, "rx_receiver_error.32f"))
        #self.connect((self.receiver, 2),
                     #gr.file_sink(gr.sizeof_float, "rx_receiver_phase.32f"))
        #self.connect((self.receiver, 3),
                     #gr.file_sink(gr.sizeof_float, "rx_receiver_freq.32f"))
        #if self._differential:
            #self.connect(self.diffdec,
                         #gr.file_sink(gr.sizeof_char, "rx_diffdec.8b"))
        #if self._constellation.apply_pre_diff_code():
            #self.connect(self.symbol_mapper,
                         #gr.file_sink(gr.sizeof_char, "rx_symbol_mapper.8b"))
        #self.connect(self.unpack,
                     #gr.file_sink(gr.sizeof_char, "rx_unpack.8b"))

    #def add_options(parser):
        #"""
        #Adds generic demodulation options to the standard parser
        #"""
        ## Add options shared with modulator.
        #add_common_options(parser)
        ## Add options specific to demodulator.
        #parser.add_option("", "--freq-bw", type="float", default=_def_freq_bw,
                          #help="set frequency lock loop lock-in bandwidth [default=%default]")
        #parser.add_option("", "--phase-bw", type="float", default=_def_phase_bw,
                          #help="set phase tracking loop lock-in bandwidth [default=%default]")
        #parser.add_option("", "--timing-bw", type="float", default=_def_timing_bw,
                          #help="set timing symbol sync loop gain lock-in bandwidth [default=%default]")
    #add_options=staticmethod(add_options)

    #def extract_kwargs_from_options(cls, options):
        #"""
        #Given command line options, create dictionary suitable for passing to __init__
        #"""
        #return extract_kwargs_from_options_for_class(cls, options)
    #extract_kwargs_from_options=classmethod(extract_kwargs_from_options)