differential equations – Obtaining more values from `NDEigenvalues`

When I run this code

diskLevels = 
  NDEigenvalues[
                {-Laplacian[u[x, y], {x, y}], 
                DirichletCondition[u[x, y] == 0, True]}, 
                u[x, y], {x, y} [Element] Disk[], 1000
         ];

Mathematica objects that it can only supply me with 961 eigenvalues:

NDEigenvalues::maxeigen: A maximum number of 961 eigenvalues and functions can be computed for this discretized system.

Is there a way to get more values? I tried messing around with some options but any time I managed to prevent the error message the calculation would not terminate.