Перепишите двойной цикл более красивым и, возможно, более коротким способом

Мне интересно, можно ли написать следующий код более приятным способом. По сути, я хочу вычислить z = f (x, y) для сетки (x, y) .

a = linspace(0, xr, 100)                                                                  
b = linspace(0, yr, 100)                                                                  

for i in xrange(100):
   for j in xrange(100):
      z[i][j] = f(a[i],b[j])