If we want to describe a static spherically symmetric star we can use a metric which matches the Schwarzschild solution with correct mass on the outside of the star but differs from Schwartzschild in the inside of the matter distribution.
Basically we solve the Einstein equations with a source $T_{\mu\nu}$, for instance $$T_{\mu\nu}=(\rho+p)u_{\mu}u_{\nu}+p\,g_{\mu\nu}$$ where $u_{\mu}$ has zero spatial components, meaning it is the velocity in a static fluid (this can also be seen as a consequence of Einstein equations).
Can we do something similar for a rotating star using the metric for a Kerr black hole?
I heard that it is a much more difficult problem and I would like to understand how difficult it is (Is it possible?) and what makes it so difficult.