### 17 releases (11 breaking)

0.13.0 | Jul 27, 2023 |
---|---|

0.11.2 | Apr 27, 2023 |

0.11.0 | Mar 16, 2023 |

0.10.0 | Oct 30, 2022 |

0.8.0 | Jul 13, 2022 |

#**8** in Machine learning

**3,627** downloads per month

Used in rurel

**MIT/Apache**

1.5MB

**30K**
SLoC

# dfdx: shape checked deep learning in rust

Ergonomics & safety focused deep learning in Rust.

**Still in pre-alpha state. The next few releases are planned to be breaking releases.**

Features at a glance:

- 🔥 GPU accelerated tensor library with shapes up to 6d!
- Shapes with both compile and runtime sized dimensions. (e.g.

and`Tensor``<``(``usize`, Const`<`10`>``)``>`

)`Tensor``<`Rank2`<`5, 10`>``>` - A large library of tensor operations (including

,`matmul`

, and much more).`conv2d`- All tensor operations shape and type checked at compile time!!

- Ergonomic neural network building blocks (like

,`Linear`

, and`Conv2D`

).`Transformer` - Standard deep learning optimizers such as

,`Sgd`

,`Adam`

,`AdamW`

, and more.`RMSprop`

is on crates.io! Use by adding this to your `dfdx`

:`Cargo .toml`

`dfdx = "0.13.0"
`

See the documentation at docs.rs/dfdx.

[1] https://en.wikipedia.org/wiki/Automatic_differentiation#Reverse_accumulation

## Design Goals

- Ergonomics the whole way down (both frontend interface & internals).
- Check as much at compile time as possible (i.e. don't compile if something is not correct).
- Maximize performance.
- Minimize unsafe code[1]
- Minimize Rc<RefCell> used in internal code[2]

[1] Currently the only unsafe calls are for matrix multiplication.

[2] The only things that use

are tensors to store their data. `Arc`

is used instead of `Arc`

to reduce
allocations when tensors are cloned.`Box`

## GPU acceleration with CUDA

Enable the

feature to start using the `cuda`

device! Requires the installation of nvidia's cuda toolkit. See feature flags docs for more info.`Cuda`

## API Preview

Check examples/ for more details.

- 👌 Simple Neural Networks API, completely shape checked at compile time.

`type` `Mlp` `=` `(`
`(``Linear``<`10, 32`>``,` ReLU`)``,`
`(``Linear``<`32, 32`>``,` ReLU`)``,`
`(``Linear``<`32, 2`>``,` Tanh`)``,`
`)``;`
`fn` `main``(``)`` ``{`
`let` dev`:` Cuda `=` `Default``::`default`(``)``;` `//` or `Cpu`
`let` mlp `=` dev`.``build_module``::``<`Mlp, `f32``>``(``)``;`
`let` x`:` `Tensor``<`Rank1`<`10`>`, `f32`, Cpu`>` `=` dev`.``zeros``(``)``;`
`let` y`:` `Tensor``<`Rank1`<`2`>`, `f32`, Cpu`>` `=` mlp`.``forward``(`x`)``;`
mlp`.``save``(``"`checkpoint.npz`"``)``?``;`
`}`

- 📈 Ergonomic Optimizer API

`type` `Model` `=` `...`
`let` `mut` model `=` dev`.``build_module``::``<`Model, `f32``>``(``)``;`
`let` `mut` grads `=` model`.``alloc_grads``(``)``;`
`let` `mut` sgd `=` `Sgd``::`new`(``&`model`,` SgdConfig `{`
lr`:` 1e`-``2``,`
momentum`:` `Some``(``Momentum``::`Nesterov`(``0.``9``)``)`
`}``)``;`
`let` loss `=` `...`
grads `=` loss`.``backward``(``)``;`
sgd`.``update``(``&``mut` model`,` `&`grads`)``;`

- 💡 Const tensors can be converted to and from normal rust arrays

`let` t0`:` `Tensor``<`Rank0, `f32`, `_``>` `=` dev`.``tensor``(``0.``0``)``;`
`assert_eq!``(`t0`.``array``(``)``,` `&``0.``0``)``;`
`let` t1 `/*`: Tensor<Rank1<3>, f32, _>`*/` `=` dev`.``tensor``(``[``1.``0``,` `2.``0``,` `3.``0``]``)``;`
`assert_eq!``(`t1`.``array``(``)``,` `[``1.``0``,` `2.``0``,` `3.``0``]``)``;`
`let` t2`:` `Tensor``<`Rank2`<`2, 3`>`, `f32`, `_``>` `=` dev`.``sample_normal``(``)``;`
`assert_ne!``(`t2`.``array``(``)``,` `[``[``0.``0``;` `3``]``;` `2``]``)``;`

## Fun/notable implementation details

### Module

`pub` `trait` `Module`<Input> `{`
`type` `Output``;`
`fn` `forward``(``&``self`, `input``:` Input`)`` ``->` `Self``::`Output`;`
`}`

From this flexible trait we get:

- Single & batched inputs (just have multiple impls!)
- Multiple inputs/outputs (multi-headed modules, or rnns)
- Behavior different when tape is present or not (
**not**the .train()/.eval() behavior present in other libraries!).

### Tuples represent feedforward (a.k.a sequential) modules

Since we can implement traits for tuples, which is *not possible in other languages* AFAIK, they provide a very nice frontend
for sequentially executing modules.

`//` no idea why you would do this, but you could!
`type` `Model` `=` `(`ReLU`,` Sigmoid`,` Tanh`)``;`
`let` model `=` dev`.``build_module``::``<`Model, `f32``>``(``)``;`

`type` `Model` `=` `(``Linear``<`10, 5`>``,` Tanh`)`
`let` model `=` dev`.``build_module``::``<`Model, `f32``>``(``)``;`

How implementing Module for a 2-tuple looks:

`impl``<`Input, A, B`>`` ``Module``<`Input`>` `for`` (``A`, `B`)
`where`
Input`:` Tensor,
A`:` `Module``<`Input`>`, `//` A is a module that takes Input
B`:` `Module``<``A``::`Output`>`, `//` B is a module that takes A's Output
`{`
`type` `Output` `=` `B``::`Output`;` `//` the output of this is B's Output
`fn` `forward``(``&``self`, `x``:` Input`)`` ``->` `Self``::`Output `{`
`let` x `=` `self``.``0.``forward``(`x`)``;`
`let` x `=` `self``.``1.``forward``(`x`)``;`
x
`}`
`}`

Modules implemented for Tuples up to 6 elements, but *you can arbitrarily nest them*!

### No `Rc``<`RefCells`<`T`>``>`

used - Gradient tape is not kept behind a cell!

`Rc``<`RefCells`<`T`>``>`Other implementations may store a reference to the gradient tape directly on tensors, which requires mutating tensors or using Rc/Refcells all over the place.

We've figured out an elegant way to avoid this, reducing references and dynamic borrow checks to 0!

Since all operations result in exactly 1 child, we can always move the gradient tape to the child of the last operation. Additionally, no model parameters (all tensors) will ever own the gradient tape because they will never be the result of any operation. This means we know exactly which tensor owns the gradient tape, and the tensors that have it will always be intermediate results that don't need to be maintained across gradient computation.

*All of this together gives users unprecedented control/precision over what tensors are recorded on the gradient tape!*

One advanced use case requires that tensors be re-used multiple times in a computation graph. This can be handled by cloning the tensor, and manually moving the gradient tape around.

### Type checked backward

tl;dr: If you forget to include a call to

or `trace``(``)`

, the program won't compile!`traced``(``)`

`-`let pred = module.forward(x);
`+`let pred = module.forward(x.traced(grads));
let loss = (y - pred).square().mean();
let gradients = loss.backward();

Since we know exactly what tensors own the gradient tape, we can require the tensor passed into

to own the gradient tape!
And further, we can require it be moved into `.``backward``(``)`

, so it can destruct the tape and construct the gradients!`.``backward``(``)`

**All of this can be checked at compile time 🎉**

### 📄 Validated against pytorch

All functions & operations are tested against behavior shown by similar code in pytorch.

# License

Dual-licensed to be compatible with the Rust project.

Licensed under the Apache License, Version 2.0 http://www.apache.org/licenses/LICENSE-2.0 or the MIT license http://opensource.org/licenses/MIT, at your option. This file may not be copied, modified, or distributed except according to those terms.

#### Dependencies

~1–2.6MB

~53K SLoC