docs: def chi_squared_map_from_residual_map_noise_map_and_mask(residual_map, noise_map, mask):
      """Computes the chi-squared map between a masked residual-map and noise-map, where:
  
      Chi_Squared = ((Residuals) / (Noise)) ** 2.0 = ((Data - Model)**2.0)/(Variances)
  
      Although noise-maps should not contain zero values, it is possible that masking leads to zeros which when \
      divided by create NaNs. Thus, nan_to_num is used to replace these entries with zeros.
  
      Parameters
      -----------
      residual_map : np.ndarray
          The residual-map of the model-data fit to the observed data.
      noise_map : np.ndarray
          The noise-map of the observed data.
      mask : np.ndarray
          The mask applied to the residual-map, where *False* entries are included in the calculation.
      """
      return np.square(np.divide(residual_map, noise_map, out=np.zeros_like(residual_map),
                                 where=np.asarray(mask) == 0))