from zencad import *

p_1=[(-6.5,2.045,0),
     (-6,0.278,0),
     (-5.5,-0.68,0),
     (-5,-1.346,0),
     (-4.5,-1.85,0),
     (-4,-2.248,0),
     (-3.5,-2.566,0),
     (-3,-2.821,0),
     (-2.5,-3.021,0),
     (-2,-3.172,0),
     (-1.5,-3.276,0),
     (-1,-3.338,0),
     (-0.5,-3.357,0),
     (0,-3.335,0),
     (0.5,-3.274,0),
     (1,-3.176,0),
     (1.5,-3.041,0),
     (2,-2.872,0),
     (2.5,-2.67,0),
     (3,-2.436,0),
     (3.5,-2.173,0),
     (4,-1.883,0),
     (4.5,-1.567,0),
     (5,-1.226,0),
     (5.5,-0.862,0),
     (6,-0.477,0),
     (6.5,-0.072,0),
     (7,0.351,0),
     (7.5,0.792,0),
     (8,1.249,0),
     (8.5,1.721,0),
     (9,2.207,0),
     (9.5,2.706,0),
     (10,3.217,0),
     (10.5,3.74,0),
     (11,4.246,0),#
     (11.5,4.79,0),
     (11,4.704,0),
     (10.5,4.412,0),
     (10,4.029,0),
     (9.5,3.657,0),
     (9,3.298,0),
     (8.5,2.952,0),
     (8,2.621,0),
     (7.5,2.305,0),
     (7,2.006,0),
     (6.5,1.726,0),
     (6,1.46,0),
     (5.5,1.226,0),
     (5,1.009,0),
     (4.5,0.816,0),
     (4,0.647,0),
     (3.5,0.505,0),
     (3,0.39,0),
     (2.5,0.305,0),
     (2,0.249,0),
     (1.5,0.223,0),
     (1,0.229,0),
     (0.5,0.268,0),
     (0,0.339,0),
     (-0.5,0.444,0),
     (-1,0.583,0),
     (-1.5,0.757,0),
     (-2,0.97,0),
     (-2.5,1.224,0),
     (-3,1.526,0),
     (-3.5,1.884,0),
     (-4,2.316,0),
     (-4.5,2.848,0),
     (-5,3.532,0),
     (-5.5,4.485,0),
     (-6,3.782,0),#
     (-6.5,2.045,0),
     ]

p_2=[(-6.5,2.211,15),
     (-6,0.632,15),
     (-5.5,-0.267,15),
     (-5,-0.905,15),
     (-4.5,-1.395,15),
     (-4,-1.79,15),
     (-3.5,-2.111,15),
     (-3,-2.371,15),
     (-2.5,-2.578,15),
     (-2,-2.736,15),
     (-1.5,-2.856,15),
     (-1,-2.909,15),
     (-0.5,-2.926,15),
     (0,-2.898,15),
     (0.5,-2.826,15),
     (1,-2.711,15),
     (1.5,-2.555,15),
     (2,-2.359,15),
     (2.5,-2.126,15),
     (3,-1.859,15),
     (3.5,-1.558,15),
     (4,-1.226,15),
     (4.5,-0.865,15),
     (5,-0.477,15),
     (5.5,-0.064,15),
     (6,0.373,15),
     (6.5,0.832,15),
     (7,1.311,15),
     (7.5,1.81,15),
     (8,2.326,15),
     (8.5,2.859,15),
     (9,3.39,15),
     (9.5,3.93,15),
     (10,4.483,15),#
     #(10.5,5.02,15),
     (10.5,5.351,15),#
     (10,5.059,15),
     (9.5,4.606,15),
     (9,4.165,15),
     (8.5,3.736,15),
     (8,3.321,15),
     (7.5,2.92,15),
     (7,2.537,15),
     (6.5,2.171,15),
     (6,1.824,15),
     (5.5,1.499,15),
     (5,1.197,15),
     (4.5,0.92,15),
     (4,0.669,15),
     (3.5,0.447,15),
     (3,0.256,15),
     (2.5,0.097,15),
     (2,-0.027,15),
     (1.5,-0.116,15),
     (1,-0.168,15),
     (0.5,-0.182,15),
     (0,-0.157,15),
     (-0.5,-0.094,15),
     (-1,0.008,15),
     (-1.5,0.148,15),
     (-2,0.325,15),
     (-2.5,0.541,15),
     (-3,0.796,15),
     (-3.5,1.097,15),
     (-4,1.452,15),
     (-4.5,1.878,15),
     (-5,2.409,15),
     (-5.5,3.113,15),
     (-6,4.173,15),
     (-6.5,2.211,15),
     ]

p_3=[(-6.5,-2.282,35),
     (-6,-2.704,35),
     (-5.5,-2.986,35),
     (-5,-3.159,35),
     (-4.5,-3.246,35),
     (-4,-3.262,35),
     (-3.5,-3.219,35),
     (-3,-3.123,35),
     (-2.5,-2.982,35),
     (-2,-2.801,35),
     (-1.5,-2.583,35),
     (-1,-2.333,35),
     (-0.5,-2.051,35),
     (0,-1.741,35),
     (0.5,-1.405,35),
     (1,-1.045,35),
     (1.5,-0.661,35),
     (2,-0.256,35),
     (2.5,0.171,35),
     (3,0.617,35),
     (3.5,1.081,35),
     (4,1.564,35),
     (4.5,2.064,35),
     (5,2.58,35),
     (5.5,3.112,35),
     (6,3.659,35),
     (6.5,4.221,35),
     (7,4.785,35),
     (7.5,5.35,35),
     (8,5.825,35),#
     (8.5,6.48,35),
     #(8.5,6.984,35),
     (8,6.353,35),#
     (7.5,6.004,35),
     (7,5.518,35),
     (6.5,5.038,35),
     (6,4.564,35),
     (5.5,4.098,35),
     (5,3.643,35),
     (4.5,3.198,35),
     (4,2.767,35),
     (3.5,2.349,35),
     (3,1.947,35),
     (2.5,1.561,35),
     (2,1.192,35),
     (1.5,0.842,35),
     (1,0.51,35),
     (0.5,0.197,35),
     (0,-0.095,35),
     (-0.5,-0.368,35),
     (-1,-0.621,35),
     (-1.5,-0.854,35),
     (-2,-1.066,35),
     (-2.5,-1.259,35),
     (-3,-1.429,35),
     (-3.5,-1.575,35),
     (-4,-1.693,35),
     (-4.5,-1.777,35),
     (-5,-1.813,35),
     (-5.5,-1.783,35),
     (-6.4,-2.068,35),#
     (-6.5,-2.282,35),
     ]

POINTS=points2([p_1,p_2,p_3])

##m=interpolate2(POINTS,degmin=1,degmax=2)
##m1=fill(interpolate(p_1))
##m2=fill(interpolate(p_2))
m1=interpolate(p_1)
m2=interpolate(p_2)
m3=interpolate(p_3)
m=loft([m1,m2,m3],smooth=True,shell=False,maxdegree=3)
#disp(m)
##disp(m1)
##disp(m2)
nodes,triangles=triangulate(m,0.05)
mm=polyhedron(nodes,triangles)
disp(polyhedron(nodes,triangles))
disp(POINTS,color=color.red)
to_stl(m,'c:/1/02.stl',0.05)

show()