While this is a trivial example, I have a few other models where avoiding linear algebra operations are difficult. In case anyone else runs into these issues, I’ve had success with gradiant free samplers like DEMetropolis (or Metropolis) step methods.