Performance for circular correlation coefficient

Hi all,

I'm trying to calculate the correlation coefficient between two waves containing angles. From what I gathered, StatsCircularCorrelationTest() provides the functionality I want. However, I run into issues when the waves become very long (Igor locking up, beach ball etc). I want to do the calculation for waves containing >10k points.

Here is a minimal example:

Function test()
	Make/O/N=1000 one, two //barely works for 1000, locks up for large numbers, e.g. 10000
	setscale/I x, 0, 1, one, two
	one = gnoise(1)
	two = gnoise(1)
	
	StatsCircularCorrelationTest/PAA/Q one, two
	Wave W_StatsCircularCorrelationTest
	print W_StatsCircularCorrelationTest[0]
End

This is on Igor 8.04B01 (Mac).

Is there a better way to only get the correlation coefficient?