Emulate directed rounding in default rounding mode

Zygote-differentiable way of flattening and unflattening data structures

Implementation of Takahashi''s inverse subset algorithm for sparse matrices

-