Feedback delay networks were first suggested for artificial reverberation by Gerzon , who proposed an ``orthogonal matrix feedback reverberation unit''. He noted that individual feedback comb filters yielded poor quality, but that several such filters could sound good when cross-coupled. An ``orthogonal matrix feedback'' around a parallel bank of delay lines was suggested as a means of obtaining maximally rich cross-coupling. He was especially concerned with good stereo spreading of the reverberation at a time when most artificial reverberators sought merely to decorrelate the reverberation in each output channel.
Later, and apparently independently, Stautner and Puckette  suggested a specific four-channel FDN reverberator and gave general stability conditions for the FDN. They proposed the feedback matrix
which can be understood as a permutation (with one row sign inversion) of a Hadamard matrix (see §3.7.2).
More recently, Jot [218,217] developed a systematic FDN design methodology allowing largely independent setting of reverberation time in different frequency bands. Using Jot's methodology, FDN reverberators can be polished to a high degree of quality, and they are presently considered to be among the best choices for high-quality artificial reverberation.
Jot's early work was concerned only with single-input, single-output (SISO) reverberators. Later work  with Jullien and others at IRCAM was concerned also with spatializing the reverberation.
An example FDN reverberator using three delay lines is shown in Fig.3.10. It can be seen as an FDN (introduced in §2.7), plus an additional low-order filter applied to the non-direct signal. This filter is called a ``tonal correction'' filter by Jot, and it serves to equalize modal energy irrespective of the reverberation time in each band. In other words, if the decay time is made very short in some band, will have a large gain in that band so that the total energy in the band's impulse-response is unchanged. This is another example of orthogonalization of reverberation parameters: In this case, adjustments in reverberation time, in any frequency band, do not alter total signal energy in the impulse response in that band.