] cp.cnstr_rhs = np.zeros((len(cp.cnstr_lhs),), dtype=float) # @todo - renaming of methods # @todo - projection on the caf - to get the initial vector # @todo - gemetry transformator # @todo - derivatives of caf for the current position. # @todo - rthombus generator with cut-away elements # @todo - time step counting - save the initial step separately from the time history X0 = cp.generate_X0() X_fc = cp.solve(X0 + 1e-6) print "nodes", cp.get_new_nodes(X_fc) # # print 'nodes' # new_nodes = cp.get_new_nodes(X_fc) # cp2 = CreasePattern(nodes = new_nodes, # crease_lines = cp.crease_lines, # facets = cp.facets, # n_steps = 1, # show_iter = True, # z0_ratio = 0.1, # MAX_ITER = 200) # # face_z_t = CnstrTargetFace(F = [r_, s_, 0]) # # cp2.tf_lst = [(face_z_t, n_arr)]
] cp.cnstr_rhs = np.zeros((len(cp.cnstr_lhs),), dtype=float) # @todo - renaming of methods # @todo - projection on the caf - to get the initial vector # @todo - gemetry transformator # @todo - derivatives of caf for the current position. # @todo - rthombus generator with cut-away elements # @todo - time step counting - save the initial step separately from the time history X0 = cp.generate_X0() X_fc = cp.solve(X0 + 1e-6) print 'nodes', cp.get_new_nodes(X_fc) my_model = CreasePatternView(data=cp, ff_resolution=30, show_cnstr=True) my_model.configure_traits() # al = InfocadLink(data = cp, n_split = 3) # al.model_name = 'HP_shell' # al.build_inp() # al = InfocadLink(data=cp, n_split=3)
connectivity = [(vertex, neighbors) for vertex, neighbors in zip(cp.interior_vertices, cp.cycled_neighbors.T)] print "connectivity", connectivity uf = EqConsDevelopability(cp, connectivity=connectivity) cp.eqcons["uf"] = uf # the derivatives are not correct cp.use_G_du = False u_unfoldable = cp.solve(u0 + 1e-6, acc=1e-4) # =========================================================================== # Unfolding # =========================================================================== # new_nodes = cp.get_new_nodes(u_unfoldable) cp2 = CreasePattern( nodes=new_nodes, crease_lines=cp.crease_lines, facets=cp.facets, n_steps=1, show_iter=True, z0_ratio=0.1, time_arr=np.linspace(1, 0, 10), MAX_ITER=200, ) cp2.tf_lst = [(face_z_t, n_arr)] cp2.cnstr_lhs = [ [(n_h[1, 0], 0, 1.0)], # 0