Higher-order derivatives for determinant for low rank matrices

[dpfau]

PR #2809 fixes the issue with determinants of low-rank matrices returning NaN. However, higher order derivatives still fail. Given the close relationship between the gradient of the determinant and the matrix inverse, it may be possible to use some of the implicit differentiation tricks from custom_linear_solve for second derivatives of the determinant.