generator.py

Go to the documentation of this file.

from .constants import HERMITE, LAGUERRE, computeSize
 
__all__ = ['IndexGenerator', 'HermiteIndexGenerator', 'LaguerreIndexGenerator']
 
 
class IndexGenerator:
    """Base class for shapelet index generators.
 
    Parameters
    ----------
    order: `int`
       Order of the function indices to generate.
    """
 
    __slots__ = "order", "size"
 
    def __init__(self, order):
        self.orderorder = order
        self.sizesize = computeSize(self.orderorder)
 
    @staticmethod
    def make(self, order, basisType):
        if basisType == HERMITE:
            return HermiteIndexGenerator(order)
        elif basisType == LAGUERRE:
            return LaguerreIndexGenerator(order)
 
    def __len__(self):
        return self.sizesize
 
 
class HermiteIndexGenerator(IndexGenerator):
    """Iterable that generates tuples of indices.
 
    Yields
    ------
    i : `int`
        the overall coefficient index for a 2-d shapelet expansion (just counts from zero)
    nx : `int`
        the order of the x expansion
    ny : `int`
        the order of the y expansion
    """
 
    def __iter__(self):
        i = 0
        for n in range(0, self.orderorder+1):
            for nx in range(0, n+1):
                yield (i, nx, n - nx)
                i += 1
 
 
class LaguerreIndexGenerator(IndexGenerator):
    """Iterable that generates tuples of indices.
 
    Yields
    ------
    i : `int`
        the overall coefficient index for a 2-d shapelet expansion (just counts from zero)
    p, q : `int`
        indices of the polar shapelet expansion (see BasisTypeEnum).
    re : `bool`
        True, if this is the real part of the coefficient.
    """
 
    def __iter__(self):
        i = 0
        for n in range(0, self.orderorder+1):
            p = n
            q = 0
            while p > q:
                yield (i, p, q, True)
                i += 1
                yield (i, p, q, False)
                i += 1
                p -= 1
                q += 1
            if p == q:
                yield (i, p, q, True)
                i += 1