Simple test of the precision of complex_bessel

test_bessel.cpp

#include <complex_bessel.h>
#include <iostream>
#include <iomanip>

int main(int argc, char* argv[])
{
  std::cout << std::setprecision(20) << "complex_bessel: besselI(1,5) = " << sp_bessel::besselI(1,5) << std::endl;
  std::cout << std::setprecision(20) << "besselI(1,5)[complex_bessel]- besselI(1,5)[matlab] = " << sp_bessel::besselI(1,5)-24.335642142450530 << std::endl;
  std::cout << std::setprecision(20) << "besselI(1,5)[complex_bessel]-besselI(1,5)[WolframAlpha] = " << sp_bessel::besselI(1,5)-24.33564214245052719914305045176000846056487436829889815840 << std::endl;
  return 0;
}