Chain 1 failed due to Mass matrix contains zeros on the diagonal

[chigface]

But I replaced init=‘adapt_diag’ with init=‘map’ and it work but cann’t really differentiate the two.