view_transform_generators 8 1.0 0.0 0.0 0 0.0 1.0 0.0 0 0.0 0.0 1.0 0 0 0 0 1 -1.0 0.0 0.0 0 0.0 1.0 0.0 0 0.0 0.0 1.0 0 0 0 0 1 1.0 0.0 0.0 0 0.0 -1.0 0.0 0 0.0 0.0 1.0 0 0 0 0 1 -1.0 0.0 0.0 0 0.0 -1.0 0.0 0 0.0 0.0 1.0 0 0 0 0 1 0.0 1.0 0.0 0 1.0 0.0 0.0 0 0.0 0.0 -1.0 0 0 0 0 1 0.0 1.0 0.0 0 -1.0 0.0 0.0 0 0.0 0.0 -1.0 0 0 0 0 1 0.0 -1.0 0.0 0 1.0 0.0 0.0 0 0.0 0.0 -1.0 0 0 0 0 1 0.0 -1.0 0.0 0 -1.0 0.0 0.0 0 0.0 0.0 -1.0 0 0 0 0 1 SYMMETRIC_CONTENT SYMMETRY_GROUP "cubocta" PARAMETER angle = 119.99999999999999 #define WALLT (-cos((angle/2.0)*pi/180)) // virtual tension of facet on plane constraint 1 formula: 1.0*x +0.0*y +0.0*z =0.0 energy: // for dihedral angle e1: WALLT * (0.0*x +-0.0*y +0.0*z +-0.0) e2: WALLT * (0.0*x +0.0*y +-0.5*z +-0.0) e3: WALLT * (-0.0*x +0.5*y +0.0*z +-0.0) constraint 2 formula: 0.0*x +1.0*y +-0.0*z =0.0 constraint 3 formula: -0.57735026918962573*x +0.57735026918962573*y +-0.57735026918962573*z =0.0 energy: // for dihedral angle e1: WALLT * (0.0*x +0.57735026918962573*y +0.57735026918962573*z +-0.0) e2: WALLT * (-0.57735026918962573*x +0.0*y +0.57735026918962573*z +-0.0) e3: WALLT * (-0.57735026918962573*x +-0.57735026918962573*y +0.0*z +-0.0) constraint 4 formula: 0.70710678118654757*x +-0.70710678118654757*y +0.0*z =0.0 #define axis axial_point constraints 3 4 vertices 1 -0.35355339059327373 -0.35355339059327373 0.0 axis 2 0.0 -0.7071067811865475 -0.7071067811865475 constraints 1 3 3 0.0 0.0 -0.7071067811865475 constraints 1 2 edges 2 1 3 3 1 2 constraint 3 color green wrap 36 4 2 3 constraint 1 color red 5 3 2 constraint 2 color blue wrap 36 facets 1 -3 2 -4 2 -2 3 -5 bodies 1 1 2 volume 0.09428090415820633563 read raw_cells showall := { transform_expr "abcdefgh"; transforms on; } showall groom:= {t 0.02; w 5e-4; V; u; l 0.02; u; t 0.02; w 5e-4; V; u} gogo := {w 0.1; t 0.1; u; g 10; l 0.1; u; w 1e-2; t 1e-2; V; u; g 100; u; l 0.1; V; u; w 1e-2; t 1e-2; V; u; g 200; w 1e-3; t 1e-3; u; l 0.02; V; u; g 100; g 500; groom; groom; U; g 500; groom; groom; g 500; groom; groom; g 500; groom; groom; hessian_seek; hessian_seek; hessian_seek; hessian_seek; hessian; hessian; hessian; hessian; hessian; hessian; hessian; hessian; hessian; U}