How do you speed up the linear sparse solver in Zygote?

To compute \nabla_x \Vert A(x)^{-1} b \Vert with current AD tools where A(x) is sparse, you still have the problem of trying to differentiate through the sparse-matrix constructor, no?