Thank you @mohamed82008 for this package. Can it solve https://github.com/JuliaDiff/ForwardDiff.jl/issues/602? I wonder if we can easily use the factorizations implemented here with Turing.jl models.