def init_random(seed=0):
    """
    Initialize and seed the random number generator.
    """
    assert seed >= 0, "Seed must be a non-negative integer."
    EigenMatrix.rnd_state = rnd_struct()
    EigenMatrix.rnd_state_p = ct.pointer(EigenMatrix.rnd_state)
    _eigenmat.init_random(EigenMatrix.rnd_state_p, ct.c_int(seed+1))

    Initialize and seed the random number generator.
    
eigenmat/eigenmat.py