X_norm = (X_norm / (X_norm.max(axis=0, keepdims=True)
                            + np.spacing(0)))