b = zeros((m, p))