Compute the complete Pearson correlation matrix across dataset features with stable handling of numeric inputs.
