import math

def pytha(c,b):
    value = math.sqrt(c*c - b*b)
    return value

#location of slot
xpos = 40
ypos = 0

# size in mm
mill = 0.6
cir1 = 5.3
cir2 = 6.3
cir3 = 2.3
ycr3 = 4.0
min1 = 0.05

step = 4 # degrees

# position
xhole1 = 0.0
yhole1 = 6.0
xhole2 = -4.0 # top mill
yhole2 = 3.0

# size solder holes
xhole = 1.5
yhole = 0.8

#pcb height
zmax = -22

#diode
dxsize = 2.5
dysize = 7.5 - dxsize


def pytha(c,b):
    value = math.sqrt(c*c - b*b)
    return value


def pcbmillcherry(xpos,ypos):

    c1 = ((cir1 - min1)/ 2) - (mill /2)
    c2 = ((cir2 - min1)/ 2) - (mill /2)
    c3 = ((cir3 - min1)/ 2) - (mill /2)
    y3 = ycr3
    x =  round (math.degrees(math.atan(c3/(pytha(c2, c3)))))

    print ('G00 Z5.0')
    print ('G00 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + pytha(c1, c3)))
    for z in range (-2, zmax, -2):
        print ('G01 F50')
        print ('G01 Z{0:.4f}' .format(z/10))

        print ('G01 F50')

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + pytha(c1, c3)))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + y3))

        for i in range (-90,90,step*3):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c3 * math.sin(math.radians(i)) , ypos + y3 + c3 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos + y3))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos + pytha(c1, c3)  ))

        for i in range (15,75,step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c1, c3), ypos + c3))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c2, c3), ypos + c3))

        for i in range (90-x,90+x, step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c2, c3), ypos - c3))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c1, c3), ypos - c3))

        for i in range (15+90,75+90,step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos - pytha(c1, c3)))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos - pytha(c2, c3)))

        for i in range (180-x,180+x, step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos - pytha(c2, c3)))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos - pytha(c1, c3)))

        for i in range (15+180,75+180,step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c1, c3), ypos - c3))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c2, c3), ypos - c3))

        for i in range (270-x,270+x, step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c2, c3), ypos + c3))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c1, c3), ypos + c3))

        for i in range (15+270,75+270,step):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + pytha(c1, c3)  ))

    c1 = (cir1 / 2) - (mill /2)
    c2 = (cir2 / 2) - (mill /2)
    c3 = (cir3 / 2) - (mill /2)
    y3 = ycr3
    x =  round (math.degrees(math.atan(c3/(pytha(c2, c3)))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + pytha(c1, c3)))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + y3))

    for i in range (-90,90,step*3):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c3 * math.sin(math.radians(i)) , ypos + y3 + c3 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos + y3))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos + pytha(c1, c3)  ))

    for i in range (15,75,step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c1, c3), ypos + c3))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c2, c3), ypos + c3))

    for i in range (90-x,90+x, step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c2, c3), ypos - c3))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + pytha(c1, c3), ypos - c3))

    for i in range (15+90,75+90,step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos - pytha(c1, c3)))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + c3, ypos - pytha(c2, c3)))

    for i in range (180-x,180+x, step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos - pytha(c2, c3)))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos - pytha(c1, c3)))

    for i in range (15+180,75+180,step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c1, c3), ypos - c3))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c2, c3), ypos - c3))

    for i in range (270-x,270+x, step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c2 * math.sin(math.radians(i)) , ypos + c2 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c2, c3), ypos + c3))
    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - pytha(c1, c3), ypos + c3))

    for i in range (15+270,75+270,step):
        print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + c1 * math.sin(math.radians(i)) , ypos + c1 * math.cos(math.radians(i))))

    print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - c3, ypos + pytha(c1, c3)  ))
    print ('G00 Z5.0')

    xh = (xhole - mill)/2
    yh = (yhole - mill)/2

    print ('G00 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
    for z in range (-2, zmax, -2):
        print ('G01 F50')
        print ('G01 Z{0:.4f}' .format(z/10))
        print ('G01 F50')
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
    print ('G00 Z5.0')


    for z in range (-2,-6, -2):
        print ('G01 F50')
        print ('G01 Z{0:.4f}' .format(z/10))
        print ('G01 F50')

        xh = (mill * 0.5) + (xhole - mill)/2
        yh = (mill * 0.5) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))

        xh = (mill * 1.0) + (xhole - mill)/2
        yh = (mill * 1.0) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))

        xh = (mill * 1.5) + (xhole - mill)/2
        yh = (mill * 1.5) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 - xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole1 + xh, ypos + yhole1 + yh))
    print ('G00 Z5.0')




    xh = (xhole - mill)/2
    yh = (yhole - mill)/2

    print ('G00 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
    for z in range (-2, zmax, -2):
        print ('G01 F50')
        print ('G01 Z{0:.4f}' .format(z/10))
        print ('G01 F50')
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 - yh))
    print ('G00 Z5.0')

    for z in range (-2,-6, -2):
        print ('G01 F50')
        print ('G01 Z{0:.4f}' .format(z/10))
        print ('G01 F50')

        xh = (mill * 0.5) + (xhole - mill)/2
        yh = (mill * 0.5) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))

        xh = (mill * 1.0) + (xhole - mill)/2
        yh = (mill * 1.0) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))

        xh = (mill * 1.5) + (xhole - mill)/2
        yh = (mill * 1.5) + (yhole - mill)/2
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 + yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 - xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 - yh))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + xhole2 + xh, ypos + yhole2 + yh))
    print ('G00 Z5.0')


def pcbmilldiode(xpos,ypos):

    dx = (dxsize - mill ) / 2.0
    dy = (dysize - mill ) / 2.0
    dc = dx

    print ('G00 Z5.0')
    print ('G00 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos + dy))
    for z in range (-2, zmax, -2):

        print ('G01 Z{0:.4f}' .format(z/10))

        for i in range (-90,90,step*3):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + dc * math.sin(math.radians(i)) , ypos + dy + dc * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + dx, ypos + dy))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + dx, ypos - dy))

        for i in range (90,270,step*3):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + dc * math.sin(math.radians(i)) , ypos - dy + dc * math.cos(math.radians(i))))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos - dy))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos + dy))
    print ('G00 Z5.0')

def pcbmilldioderotate (xpos,ypos):

    dx = (dysize - mill ) / 2.0 # rotated
    dy = (dxsize - mill ) / 2.0
    dc = dy

    print ('G00 Z5.0')
    print ('G00 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos - dy))
    for z in range (-2, zmax, -2):

        print ('G01 Z{0:.4f}' .format(z/10))

        for i in range (180,360,step*3):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos - dx + dc * math.sin(math.radians(i)) , ypos + dc * math.cos(math.radians(i))))

        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos + dy))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + dx, ypos + dy))

        for i in range (0,180,step*3):
            print ('G01 X{0:.4f} Y{1:.4f}' .format(xpos + dx + dc * math.sin(math.radians(i)) , ypos + dc * math.cos(math.radians(i))))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos + dx, ypos - dy))
        print ('G01 X{0:.4f} Y{1:.4f}' .format( xpos - dx, ypos - dy))
    print ('G00 Z5.0')


print ('G00 F50')
print ('G01 F50')

pcbmillcherry(21.43125,88.10625)
pcbmillcherry(40.48125,88.10625)
pcbmillcherry(59.53125,88.10625)
pcbmillcherry(78.58125,88.10625)
pcbmillcherry(97.63125,88.10625)
pcbmillcherry(116.68125,88.10625)
pcbmillcherry(140.49375,88.10625)
pcbmillcherry(183.35625,50.00625)
pcbmillcherry(11.90625,69.05625)
pcbmillcherry(30.95625,69.05625)
pcbmillcherry(50.00625,69.05625)
pcbmillcherry(69.05625,69.05625)
pcbmillcherry(88.10625,69.05625)
pcbmillcherry(107.15625,69.05625)
pcbmillcherry(126.20625,69.05625)
pcbmillcherry(145.25625,69.05625)
pcbmillcherry(16.66875,50.00625)
pcbmillcherry(35.71875,50.00625)
pcbmillcherry(54.76875,50.00625)
pcbmillcherry(73.81875,50.00625)
pcbmillcherry(92.86875,50.00625)
pcbmillcherry(111.91875,50.00625)
pcbmillcherry(138.1125,50.00625)
pcbmillcherry(164.30625,50.00625)
pcbmillcherry(26.19375,30.95625)
pcbmillcherry(45.24375,30.95625)
pcbmillcherry(64.29375,30.95625)
pcbmillcherry(83.34375,30.95625)
pcbmillcherry(102.39375,30.95625)
pcbmillcherry(121.35,30.95625)
pcbmillcherry(145.35,30.95625)
pcbmillcherry(164.30625,30.95625)
pcbmillcherry(183.35625,30.95625)
pcbmillcherry(23.8125,11.90625)
pcbmillcherry(52.3875,11.90625)
pcbmillcherry(76.2,11.90625)
pcbmillcherry(100.0125,11.90625)
pcbmillcherry(123.825,11.90625)
pcbmillcherry(145.25625,11.90625)
pcbmillcherry(164.30625,11.90625)
pcbmillcherry(183.35625,11.90625)





pcbmilldiode (30.95625,84.93125)
pcbmilldiode (50.00625,84.93125)
pcbmilldiode (69.05625,84.93125)
pcbmilldiode (88.10625,84.93125)
pcbmilldiode (107.15625,84.93125)
pcbmilldiode (126.20625,84.93125)
pcbmilldiode (150.01875,84.93125)
pcbmilldiode (21.43125,65.88125)
pcbmilldiode (40.48125,65.88125)
pcbmilldiode (59.53125,65.88125)
pcbmilldiode (78.58125,65.88125)
pcbmilldiode (97.63125,65.88125)
pcbmilldiode (116.68125,65.88125)
pcbmilldiode (135.73125,65.88125)
pcbmilldiode (154.78125,65.88125)
pcbmilldiode (26.19375,46.83125)
pcbmilldiode (45.24375,46.83125)
pcbmilldiode (64.29375,46.83125)
pcbmilldiode (83.34375,46.83125)
pcbmilldiode (102.39375,46.83125)
pcbmilldiode (121.44375,46.83125)
pcbmilldiode (147.6375,46.83125)
pcbmilldiode (35.71875,27.78125)
pcbmilldiode (54.76875,27.78125)
pcbmilldiode (73.81875,27.78125)
pcbmilldiode (92.86875,27.78125)
pcbmilldiode (111.91875,27.78125)
pcbmilldiode (33.3375,8.73125)
pcbmilldiode (61.9125,8.73125)
pcbmilldiode (85.725,8.73125)
pcbmilldiode (109.5375,8.73125)
pcbmilldiode (133.35,8.73125)
pcbmilldiode (154.78125,8.73125)
pcbmilldiode (173.83125,8.73125)
pcbmilldiode (173.83125,46.83125)
pcbmilldiode (135.825,27.78125)
pcbmilldiode (154.875,27.78125)
pcbmilldiode (173.83125,27.78125)
pcbmilldiode (154.78125,46.83125)


pcbmilldioderotate(183.35625,21.43125)



print ('G00 X0 Y0 Z5')