open Jig;;
let p = rp 4 in
pzl "permutation" (Array.fold_right (fun j ls -> j::ls) p []);
let (top, embd, [inf; outf]) = mix (
   [[0; 1; 2; 3]], [|[|1.; 0.; 0.|]; [|0.; 1.; 0.|]; [|0.; 0.; 1.|]; [|0.17; -0.21; 0.11|]|],
   [[0; 1; 3]; [0; 2; 3]]) p in
let bl = Glue.morphgen (nz (embd, pyt)) top [] in
(p, bl, (jig embd bl inf outf).lcurl ((if p.(2) < p.(1) then p.(1)-1 else p.(1)),
  (Array.init 3 (fun i -> Array.init 3 (fun j -> if i=j then 1. else 0.)))))
(* Should yield (_, _, [|[|1.; 0.; 0.|]; [|0.; 1.; 0.|]; [|0.; 0.; 1.|]|]) *);;