1. A method of modeling an unknown face recognition algorithm, the method comprising: providing an unknown face recognition algorithm; providing a training set (X) comprising a plurality of face images; providing a score matrix (D) containing the scores between pairs of the plurality of face images of the training set produced by the unknown face recognition algorithm, wherein a score is the distance between a pair of the face images; converting the score matrix (D) to an equivalent Euclidean matrix (DE); and using stress minimization with iterative majorization to determine a set of coordinates for the training set from the equivalent Euclidean matrix (DE).