from dppy.beta_ensembles import CircularEnsemble
circular = CircularEnsemble(
beta=2) # beta must be in {0,1,2,4}, default beta=2
# 1. Plot the eigenvalues, they lie on the unit circle
circular.sample_full_model(size_N=30, haar_mode='QR') # Sample
circular.plot() # Plot of the eigenvalues
# 2. Histogram of the angle of more points, should look uniform on [0,2pi]
circular.flush_samples() # Flush previous sample
circular.sample_full_model(size_N=1000, haar_mode='QR') # Sample
circular.hist()
from dppy.beta_ensembles import CircularEnsemble
circular = CircularEnsemble(beta=2) # beta must be >=0 integer, default beta=2
# See the cristallization of the configuration as beta increases
for b in [0, 1, 5, 10]:
circular.beta = b
circular.sample_banded_model(size_N=30)
circular.plot()
# circular.hist()