  We describe recent progress in the numerical study of the structure of
rapidly rotating superfluid neutron star models in full general relativity. The
superfluid neutron star is described by a model of two interpenetrating and
interacting fluids, one representing the superfluid neutrons and the second
consisting of the remaining charged particles (protons, electrons, muons). We
consider general stationary configurations where the two fluids can have
different rotation rates around a common rotation axis. The previously
discovered existence of configurations with one fluid in a prolate shape is
confirmed.
