How complete is PETSc.jl for parallel linear solves in FEM/CFD applications?

I’m also looking at these packages for scaling topology optimization as well.