Made a variable accuracy version of bessela1().  Graphed speeds and accuracies.  It looks like the amos ratio version needs to be at 1 ppm accuracy to reach the performance of the julia or GSL versions.  The speed of GSL is usually faster than my ratio.  My ratio is faster at 1/1000 accuracy than julia and, sometimes, is faster than gsl.  I need more optimization on it.  Looking into that now.