In [17], amplifier feedback to the strings was simulated as follows: The sum of all vibrating strings was passed through the cubic nonlinearity, multiplied by a feedback gain, delayed, and summed into the strings. This is easily implemented in Faust. (See freeax.dsp distributed with this module.)