Drazin inverse/group inverse in Julia?

Is there an implementation of the Drazin inverse (aka group inverse) in Julia somewhere? It’s basically the unique pseudoinverse of a matrix x that commutes with x, if we slightly weaken the notion of pseudoinverse to higher powers of x.