I want to deal with this problem using manopt. However, the feasible region is not a manifold owing to (4). So how can i transform this problem? I try to move (4) into the objective function as a penalty term, but the obtained solution may not satisfy (4). Can you help me? thank you.