The method is based on the gradient descent of Ma et al (2020), with initialization of the universal singular value thresholding as discussed there. The parameter identifiability constraint is the same as in the paper.
References
Z. Ma, Z. Ma, and H. Yuan. Universal latent space model fitting for large networks with edge
covariates. Journal of Machine Learning Research, 21(4):1-67, 2020.