Turned easy in several steps:
a[1, :] *= 2 a[2, :] *= 2 result = np.sum(a, axis=0) result = const * result