Approximate mode
Chaining radix-4 FFT stages optically: fewer conversions and more throughput, at a precision of a few effective bits.
What this stage does
Exact arithmetic converts at every stage (chapter 12). Floating-point-style work, an FFT for signal processing or the linear algebra of AI, tolerates small errors, so the values can stay analogue: an operand is one multi-level symbol rather than a stack of bit planes, the twiddles inside a block are fixed phase shifts, and only every S-th stage is converted. This is the mode in which the light carries the arithmetic rather than bit planes of it.
How it’s modelled
A radix-4 decimation-in-frequency FFT is built stage by stage in floating point (checked against NumPy’s FFT to 1.933 × 10⁻¹³). Each block of S stages runs: DAC → S unitary stages, each with 1 dB illustrative of loss and 0.01 rad illustrative rms phase error per waveguide → receiver noise at the design SNR less the block’s loss → ADC → digital rescale. The design SNR is what the receiver sees when one laser feeds as many units as the 6 dB margin of chapter 15 allows. Energy uses the same crossings, laser and phase shifters as exact mode, with the boundary twiddles in logic. The browser does not rerun this Monte Carlo; the animation and charts read its results.
What we found
Fewer conversions help precision, up to a point. With 1 dB per stage the output keeps 2.74 bits effective bits at S = 1 and 3.198 bits at S = 6: each conversion adds receiver noise, so chaining helps. At 3 dB per stage the loss wins and precision falls, from 2.445 bits to 1.268 bits. Better converters barely help: beyond 8 bits the receiver noise sets the floor.
Energy and throughput improve. Energy per butterfly-stage falls from 34.11 pJ to 6.929 pJ (4.923×), throughput rises 1.934× once the power cap stops binding, and the rate is 59× the exact-mode rate on the same units, because an operand is one symbol.
The honest comparison. A digital butterfly that only needs a few bits is tiny. At the precision the optics deliver, logic does the same work for a small fraction of the energy, even with every stage chained.
For one wavelength. These effective bits assume every operand rides the same colour of light. When W wavelengths share a mesh, each channel sees a slightly different transform and the error is the same at every chained stage: the worst of four channels keeps 3.19 bits at S = 6, the worst of 32 only 2.558 bits. See chapter 18.
Data table
| chained stages S | 1 dB/stage | 3 dB/stage |
|---|---|---|
| 1 | 2.767 | 2.445 |
| 2 | 3.091 | 2.45 |
| 3 | 3.212 | 2.253 |
| 4 | 3.199 | 2.094 |
| 5 | 3.146 | 1.721 |
| 6 | 3.22 | 1.268 |
Data table
| converter ENOB (bits) | S = 1 | S = 2 | S = 3 | S = 4 | S = 5 | S = 6 |
|---|---|---|---|---|---|---|
| 4 | 0.8688 | 1.362 | 1.588 | 1.582 | 1.518 | 1.858 |
| 6 | 2.406 | 2.762 | 2.918 | 2.897 | 2.881 | 2.999 |
| 8 | 2.74 | 3.061 | 3.187 | 3.167 | 3.125 | 3.198 |
| 10 | 2.777 | 3.092 | 3.208 | 3.199 | 3.138 | 3.219 |
| 12 | 2.774 | 3.098 | 3.204 | 3.185 | 3.149 | 3.21 |
| 16 | 2.78 | 3.099 | 3.215 | 3.18 | 3.135 | 3.22 |
Data table
| chained stages S | photonic | digital, matched bits |
|---|---|---|
| 1 | 34.11 | 0.1333 |
| 2 | 17.8 | 0.2141 |
| 3 | 12.37 | 0.2141 |
| 4 | 9.648 | 0.2141 |
| 5 | 8.017 | 0.2141 |
| 6 | 6.929 | 0.2141 |
What would need measured data
Stage loss, phase error per stage and the design SNR are illustrative; so are the converters’ energy and the digital comparison’s node scaling. A measured chained mesh’s loss and phase stability, and real ADC/DAC energies at the required rate, would settle where the curves sit; the shape (receiver noise sets a floor of a few bits) follows from the link budget.
What would have to change
Part II’s break-even conditions, together, including one from the wavelengths chapter. None of them is a parameter the simulator could turn up; each is a different chip or a different workload.
- Optics that do the twiddles or the reductionIn exact mode the light saves only a butterfly’s additions; its multiplies, j→ω and reduction mod p stay in logic and cap the gain at the exact-mode ceiling (1.1×).
- Workloads with O(N²) work per conversion that tolerate low precisionApproximate mode carries an operand per symbol and converts once per block; but at the SNR the laser allows the output keeps about 3.198 bits.
- Crossings and recombination far cheaper than the digital operations they replacePer exact transform the chip spends 485.9 pJ against the ASIC’s 32.35 pJ; break-even would need E_conv = −0.1406 pJ/bit.
- Photonic units smaller than the logic they displaceAgainst the equal-silicon ASIC, break-even needs at most 0.06563 mm² per unit; the model assumes 0.4 mm².
- Low-power phase shiftersThermo-optic shifters hold 0.5 W per unit and reach the power cap at U = 1,000; they are most of the chip’s extra energy per bootstrap.
- Enough memory bandwidth that compute bindsWith Part I’s HBM the chip and the ASIC take the same time (1×): no NTT speed-up can show.
- Phases made achromatic, or one mesh per wavelength, if wavelengths are to share the computeA shared mesh’s transform tilts across the band: per stage the wavelength error passes the random phase error beyond ±5.525 nm, and a mesh per wavelength costs 64 mm² against 25.6 mm² at four wavelengths.