The dual numbers can do efficient autodiff!

The codual numbers are a simple method of doing automatic differentiation in reverse mode. They contrast with the dual numbers which provide an easy way of doing automatic differentiation in forward mode. The difference between the two modes is that sometimes one is faster than the other.

The folklore appears to be that forward mode autodiff is easy to implement because it can be done using the beautiful algebra of dual numbers, while the same is assumed to not be the case for reverse mode. This repository presents a counterargument that a variant of the dual numbers called the codual numbers can be used to represent an implementation of reverse mode autodiff that is just as elegant and terse as can be done for forward mode. This idea was first suggested by Sandro Magi (pseudonym: Naasking).

This implementation of the codual numbers differs from Sandro Magis by using simple memoisation to eliminate the exponential worst-case behaviour he encountered. In Magis original implementation, this idea seems obscured, largely because the code was more effectful and therefore the opportunity for memoisation was less apparent. The memoisation is achieved using only one additional line of code!

This implementation should be simpler and more transparent than Magis, I hope. It also suggests that Magis reasoning behind the term codual numbers is perhaps misleading.

Definition of dual number and codual number

The codual numbers are the set

while the codual numbers are a subset of

where the second component is always a linear map.

A notation thats used to write a dual number is , which stands for . Formally, while .

The codual numbers may be represented using exactly the same notation as the dual numbers. They are no different than the dual numbers, except in how theyre represented on a computer! Using lambda calculus notation (which I assume you are familiar with) any dual number can be turned into the codual number , and conversely every codual number can be turned into the dual number . The difference is merely one of data structure; we need a closure to represent the codual numbers.

The definition of an operation on the codual numbers can be inferred from its definition on the dual numbers. We demonstrate this using multiplication. For dual numbers, we may define multiplication by:

For the codual numbers, we may use the correspondence to get:

where by , we mean multiplication of real numbers. Using the fact that and are linear maps, we can rearrange this to:

This is precisely how we define multiplication of codual numbers in the code.

Relationship with other autodiff strategies

It appears that there are three ways of doing reverse-mode autodiff, which correspond directly to the three stages of solving a problem using dynamic programming. See the table below:

This suggests that the tape-based approach can be derived from the coduals.

Exhibit A:

python def fib(n): if n == 0 or n == 1: return n else: return fib(n-1) + fib(n-2) Exhibit B:

```python from functools import cache

@cache def fib(n): if n == 0 or n == 1: return n else: return fib(n-1) + fib(n-2) ``` Exhibit C:

python def fib(n): a, b = 0, 1 for _ in range(n): a, b = b, a + b return a