def METHOD_NAME():
    np.random.seed(1111)  # For reproducibility.
    from numpy.random import RandomState
    rnd = RandomState(0)
    # Creating a list of solutions.
    Solutions = [sol0, sol1, sol2, sol3]
    manifold_projection = SVDProjection(Solutions, p="max")
    kernel = ProjectionKernel()
    kernel.calculate_kernel_matrix(manifold_projection.u, manifold_projection.u)
    assert np.round(kernel.kernel_matrix[0, 1], 8) == 6.0