The unvectorized code reads:
import numpy as np import numpy.ma as ma np.random.seed(42) H = np.random.uniform(0.1, 1.0, s