Mustach!

thingiverse

This doesn't appear to be a mathematical problem, but rather an input for a geometric design or computer-aided drafting (CAD) program. The output of this script would depend on the specifics of the program you're using to render it. Here is the revised code with improved comments: ```Python # 3D building dimensions in inches base_width = [4.87671232876713,4.21917808219178,3.23287671232877] heights=[6.02054794520548,4.9054109589041095,4.13027397260274] # Walls have a thickness of 1.6 inches wall_thickness = 1.6 # Floor heights and corresponding dimensions in square meters (not used) square_footage= [ [-34.32719756002873, 21.163265306118226], [-29.065136986301365, 19.130273972602744], [-23.805076412574005, 18.10548162878642 ], [-20.510547945205482, 16.287671232876713 ]] # Walls and columns have no height, while floor levels start from -6m to -3m in increments of 1 elevations = [ # In feet for convenience of user (though input is still metres) but the values are stored in meters internally [-20,-19,-18,-17,-16,-15], # [-10], []] walls= [[1,0,4],[2,3,7,8]] import matplotlib.pyplot as plt import numpy as np # create 2D projections of the 3D shapes def projXZ(xy_list,wall_thickness): points_X = [pt[0] for pt in xy_list] points_Z = [pt[1] for pt in xy_list] wall_list=[] min_height=-height #project along the X axis fig = plt.figure() ax=fig.add_subplot(111) min_val=np.inf max_val =-min_val min_pt=max_val max_pt=min_val for p in points_Z: if (p > -18): pt=np.round(np.array([points_X[min_height-p+wall_thickness],min_height+p-points_X[min_height]])) wall_list.append(list(pt)) #plt.figure(1) ax.add_line(plt.Line2D(xy=zip(*wall_list))) # ax.plot(wall_list,linestyle='-') # plt.legend([0]) for a,b in enumerate(walls): line=ax.add_patch(plt.Rectangle(xy = ((min_val,min_height+(np.array(wall_list[b]))[0][0])+5, max_val +10 ), fill=False,width=np.inf)) if np.mod(a,2) == 1: print(a,0) for h in walls[::2]: if isinstance(h,list): for ind,xpt in enumerate(zip(*[[xy_list[p][0] - xy_list[j][0]for p,j in zip(wall_list[x][3],wall_list[x][2])]] for x in range(0,len(walls)-1,4) )): # ax.set_xlim([np.min(points_X)+.1, np.max(points_X)]) for ind1 in enumerate(h): ax.text(-0.1+h[ind1-1]-xpt,wall_list[h][ind][2]+0,x=" ",verticalalignment='top') #text for support beam and posts h[7-ind]+=0 #make x, y points to be correctly calculated for line if isinstance(h,list) == False: ax.text(wall_list[ind][1],min_height + max(val-list()[3]), x= str(i), ha="left") break wall_list.pop(-5) min_height -=0.5 # min_height +=np.array(points_Z).min()-3 print(max(height-wall_thickness-20),height,wall_list[b]) plt.draw() # plt.show(block=False) return plt # return None for n,ind in enumerate(elevations): if(n!=6): x,y=[item[ind] for item in square_footage[::3]] z=[[1.9*(heights[p]-base_width[h]/4 -x / 8)**2]+1*(h %5 == 0) [ 10* ((n**7+1-1)/4 - 18+( heights[i]+n*5))]*12 for i,h in zip(ind,range(len(base_width)))] xy = [[ind+x*6,n,-10+n]] fig=plt.figure(n) fig = projXZ(xy,wall_thickness=wall_thickness) fig.show() ``` To plot the 3D object and show it: