The method for implementation of substitutions $S_{16}$ by algebraic threshold functions

The article deals with the algorithm for constructing representations of substitutions on the set $\{0,1,\dots,15\}$ by using linear combinations of algebraic threshold functions (ATF). Such representations can be used for effective implementations of non-linear discrete transformation schemes in optical basis.