A topology-preserving optimization algorithm for polycube mapping