[Diff since v0.1.41](v0.1.41...v0.1.42)

**Merged pull requests:**
- bump fill arrays (#128) (@oxinabox)
- modernize the CI (#129) (@oxinabox)
- =drop support for 1.0 and bound LinearAlgebra compat (#130) (@oxinabox)

[Diff since v0.1.40](v0.1.40...v0.1.41)

**Closed issues:**
- `+(::StridedMatrix, ::BlockDiagonal)` fails for non-square blocks (#124)

**Merged pull requests:**
- Correct `+(::BlockDiagonal, ::StridedMatrix)` for non-square blocks (#125) (@mjp98)

[Diff since v0.1.39](v0.1.39...v0.1.40)

**Closed issues:**
- Method ambiguity between `adjoint` and `BlockDiagonal` with different `eltype` (#91)

**Merged pull requests:**
- Fix method ambiguity for `*` between `Adjoint` or `Transpose` and `BlockDiagonal` with different eltype (#122) (@mjp98)

[Diff since v0.1.38](v0.1.38...v0.1.39)

**Closed issues:**
- Sampling from `MvNormal` with `BlockDiagonal` covariance is slow (#116)

**Merged pull requests:**
- Specialise `LinearAlgebra.lmul!` for `LowerTriangular` blockdiagonal matrices (#119) (@mjp98)

[Diff since v0.1.37](v0.1.37...v0.1.38)

**Closed issues:**
- Adding a `BlockDiagonal` and `Diagonal` matrix fails for non-square matrices (#103)

**Merged pull requests:**
- Add tests for `getblock` and `setblock!` (#118) (@mjp98)
- Add fallback for `+(::BlockDiagonal,::Diagonal)` when blocks are not square (#120) (@mjp98)

[Diff since v0.1.36](v0.1.36...v0.1.37)

**Merged pull requests:**
- Correctly handle `size` when `blocks` is empty (#117) (@mjp98)

[Diff since v0.1.35](v0.1.35...v0.1.36)

**Closed issues:**
- Matrix difference extremely slow (#112)

**Merged pull requests:**
- implement faster subtraction (#113) (@mzgubic)

[Diff since v0.1.34](v0.1.34...v0.1.35)

**Merged pull requests:**
- Fix _densification_pullback bug (#111) (@mzgubic)

[Diff since v0.1.33](v0.1.33...v0.1.34)

**Closed issues:**
- `cholesky(BlockDiagonal(Diagonal))` is broken (#109)

**Merged pull requests:**
- cholesky of `BlockDiagonal(Diagonal)` (#110) (@mzgubic)

[Diff since v0.1.32](v0.1.32...v0.1.33)

**Closed issues:**
- Cannot input BlockDiagonal into MvNormal in Distributions.jl (#102)
- BlockDiagonal constructor fails with Vector of BandedMatrices (#108)

**Merged pull requests:**
- Matrixify cholesky components (#106) (@mzgubic)