Add numx_stats_correlation (Pearson's r) to the stats module, operating on two equal-length arrays.
Follow CONTRIBUTING.md exactly: full Doxygen header, NULL-checks, numx_status_t return, all four deliverables (header, impl, test, docs page). Use Kahan (compensated) summation for the accumulation, matching the pattern already used in numx_vec_dot (src/linalg.c) after #57.
Add numx_stats_correlation (Pearson's r) to the stats module, operating on two equal-length arrays.
Follow CONTRIBUTING.md exactly: full Doxygen header, NULL-checks, numx_status_t return, all four deliverables (header, impl, test, docs page). Use Kahan (compensated) summation for the accumulation, matching the pattern already used in numx_vec_dot (src/linalg.c) after #57.