import sys
from PADS.SVG import SVG,colors

size = 500
margin = 10

triangles = [(1j,1,2+1j)]

def subdivide(triangles,fraction):
    x = 0.5*(1-fraction)
    for (p,q,r) in triangles:
        yield (p,(1-x)*p+x*r,q)
        yield (q,x*p+(1-x)*r,r)

for i in range(13):
    fraction = 1./(7*1.2**i)
    triangles = list(subdivide(triangles,fraction))

svg = SVG(size+(size+2*margin)*(1+1j),sys.stdout)

svg.group(style={"stroke":colors.black, "fill":colors.blue})
shift = margin*(1+1j)
for T in triangles:
    points = [p*size+shift for p in T]
    svg.polygon(points)
svg.ungroup()

svg.close()

# The SVG library 
# from http://www.ics.uci.edu/~eppstein/PADS/SVG.py
# by David Eppstein 
# http://www.ics.uci.edu/~eppstein/
"""SVG.py

Simple generation of SVG-format vector graphics files

D. Eppstein, November 2011.
"""

def _coord(x):
    """String representation for coordinate"""
    return ("%.4f" % x).rstrip("0").rstrip(".")

class SVG:
    def __init__(self, bbox, stream,
               standalone=True, prefix=None, indentation=0):
        """Create a new SVG object, to be written to the given stream.
        If standalone is True or omitted, the SVG object becomes a whole
        XML file; otherwise, it becomes an XML object within a larger XML
        file. If the prefix is nonempty, it is used to distinguish svg tags
        from other tags; a reasonable choice for the prefix value would be
        "s" or "svg". If the indentation is nonzero, it gives a number of
        spaces by which every line of the file is indented.
        
        The bbox argument should be a complex number, the farthest visible
        point from the origin in the positive quadrant. The bounding box
        will become the rectangle between the origin and that point.
        All other methods that specify points should do so using
        complex number coordinates."""
        self.stream = stream
        if prefix:
            self.prefix = prefix + ":"
        else:
            self.prefix = ""
        if standalone:
            self.stream.write('''<?xml version="1.0"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN" 
  "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
''')
        self.indentation = indentation
        self.nesting = 0
        br = _coord(bbox.real)
        bi = _coord(bbox.imag)
        self.element('''svg width="%s" height="%s" viewBox="0 0 %s %s"
     xmlns="http://www.w3.org/2000/svg" version="1.1"''' % (br,bi,br,bi),+1)

    def close(self):
        """Output the end of an SVG file."""
        self.element("svg", -1)
        if self.nesting:
            raise Exception("SVG: Unclosed tags")

    def element(self, e, delta=0, unspaced=False, style={}, **morestyle):
        """Output an SVG element.
        The delta argument distinguishes between XML tags that
        open a nested section of the XML file (delta=+1), XML tags
        that close the same section (delta=-1), and XML tags that
        stand alone without anything nested inside them (delta=0).
        Every call with delta=+1 must be matched by a call with delta=-1.
        If the style argument is nonempty, it should be a dictionary
        of style parameters, included within the object; these
        may also be passed as keyword arguments to element.
        If the same keyword is present both in style and as a keyword
        argument, the keyword argument takes priority."""
        if delta < 0:
            self.nesting += delta
        if delta >= 0 or not unspaced:
            output = [" " * (self.indentation + 2*self.nesting), "<"]
        else:
            output = ["<"]
        if delta < 0:
            output.append("/")
        output += [self.prefix, e]
        style = dict(style)
        style.update(morestyle)
        if style:
            output.append(' style="')
            second = False
            for keyword in style:
                if second:
                    output.append("; ")
                second = True
                output += [keyword, ":", style[keyword]]
            output.append('"')
        if delta > 0:
            self.nesting += delta
        elif delta == 0:
            output.append("/")
        output.append(">")
        if delta <= 0 or not unspaced:
            output.append("\n")
        self.stream.write("".join(output))

    def group(self,style={},**morestyle):
        """Start a group of objects, all with the same style"""
        self.element("g", +1, style=style, **morestyle)

    def ungroup(self):
        """End a group of objects"""
        self.element("g", -1)

    def circle(self, center, radius, style={}, **morestyle):
        """Circle with given center and radius"""
        self.element('circle cx="%s" cy="%s" r="%s"' %
                (_coord(center.real), _coord(center.imag), _coord(radius)),
            style=style, **morestyle)

    def rectangle(self, p, q, style={}, **morestyle):
        """Rectangle with corners at points p and q"""
        x = min(p.real,q.real)
        y = min(p.imag,q.imag)
        width = abs((p-q).real)
        height = abs((p-q).imag)
        self.element('rect x="%s" y="%s" width="%s" height="%s"' %
                (_coord(x), _coord(x), _coord(width), _coord(height)),
            style=style, **morestyle)

    def polygon(self, points, style={}, **morestyle):
        """Polygon with corners at the given set of points"""
        pointlist = " ".join(_coord(p.real)+","+_coord(p.imag) for p in points)
        self.element('polygon points="%s"' % pointlist,
                     style=style, **morestyle)

    def segment(self, p, q, style={}, **morestyle):
        """Line segment from p to q"""
        self.element('line x1="%s" y1="%s" x2="%s" y2="%s"' % 
                     (_coord(p.real), _coord(p.imag),
                      _coord(q.real), _coord(q.imag)), style=style, **morestyle)

    def arc(self, p, q, r, large=False, style={}, **morestyle):
        """Circular arc from p to q with radius r.
        If the large flag is set true, the arc will cover more than
        half of a circle. The SVG "sweep-flag" is not provided;
        instead, to achieve the same effect, swap p and q."""
        if large:
            large = "1"
        else:
            large = "0"
        r = _coord(abs(r))
        self.element('path d="M %s,%s A %s,%s 0 %s 0 %s,%s"' %
                     (_coord(p.real),_coord(p.imag),r,r,large,
                      _coord(q.real),_coord(q.imag)), style=style, **morestyle)

    def text(self, label, location, style={}, **morestyle):
        """Text label at the given location.
        Caller is responsible for making the label xml-safe."""
        self.element('text x="%s" y="%s"' %
            (_coord(location.real),_coord(location.imag)),
            delta=1, unspaced=True, style=style, **morestyle)
        self.stream.write(label)
        self.element('text', delta=-1, unspaced=True)

# A small color palette chosen to have high contrast
# even when viewed by color-blind readers
class colors:
    none = "none"
    white = "#FFFFFF"
    black = "#000000"
    red = "#BC1E46"
    blue = "#0081CD"
    green = "#009246"
    yellow = "#FEC200"
    magenta = "#CC33CC"