In this paper we describe a general, computationally feasible strategy to
deduce a family of interpolatory non-stationary subdivision schemes from a
symmetric non-stationary, non-interpolatory one satisfying quite mild
assumptions. To achieve this result we extend our previous work [C.Conti,
L.Gemignani, L.Romani, Linear Algebra Appl. 431 (2009), no. 10, 1971-1987] to
full generality by removing additional assumptions on the input symbols. For
the so obtained interpolatory schemes we prove that they are capable of
reproducing the same exponential polynomial space as the one generated by the
original approximating scheme. Moreover, we specialize the computational
methods for the case of symbols obtained by shifted non-stationary affine
combinations of exponential B-splines, that are at the basis of most
non-stationary subdivision schemes. In this case we find that the associated
family of interpolatory symbols can be determined to satisfy a suitable set of
generalized interpolating conditions at the set of the zeros (with reversed
signs) of the input symbol. Finally, we discuss some computational examples by
showing that the proposed approach can yield novel smooth non-stationary
interpolatory subdivision schemes possessing very interesting reproduction
properties