Sparse array glm fit in pymc3 very long with NUTS

thank you junpenglao. I am incorporating this today will report back.