I added the testBessels script in /tests.  It looks like I should switch all my bessels to GSL's scaled versions.  However, on small N, my functions are still faster.  Should check with a housemate on this.