public void computeGaussianPerplexity_larger()
{
var points = yinyang.Submatrix(null, 0, 1).ToJagged();
double perplexity = 0.5;
int N = points.Rows();
int D = points.Columns();
var X = points.ToMatrix();
double[,] expected = new double[N, N];
TSNEWrapper.computeGaussianPerplexity(X, N, D, expected, perplexity);
double[][] actual = Jagged.Zeros(N, N);
TSNE.computeGaussianPerplexity(points, N, D, ref actual, perplexity);
Assert.IsTrue(actual.IsEqual(expected, rtol: 1e-5));
}