empty_like -> zeros_like

corbass/evaluation.py

@@ -291,7 +291,7 @@ def coeffs(posterior, mu_coeffs, cov_coeffs, r_ref, r_at=REARTH):
 
     err_16, err_84 = np.percentile(ens, (16, 84), axis=1)
     mean = (mu_coeffs.T * gm_weights.flatten()).sum(axis=1)
-    var = np.empty_like(mean)
+    var = np.zeros_like(mean)
 
     for it in range(gm_weights.size):
         var += cov_coeffs[it].diagonal() \
@@ -343,7 +343,7 @@ def spectrum(posterior, mu_coeffs, cov_coeffs, r_ref, r_at=REARTH):
 
     # 1st two moments of mixture
     mean = np.sum(mu_coeffs * gm_weights.flatten()[:, None], axis=0)
-    cov = np.empty_like(cov_coeffs[0])
+    cov = np.zeros_like(cov_coeffs[0])
 
     for it in range(gm_weights.size):
         cov += cov_coeffs[it] \