A new multidimensional slow continued fraction algorithm and stepped surface

We give a new algorithm of slow continued fraction expansion related to any real cubic number field as a 2-dimensional version of the Farey map. Using our algorithm, we can find the generators of dual substitutions (so-called tiling substitutions) for any stepped surface for any cubic direction.