44 auto var = sigma * sigma;
51 auto besselI_0 = [](
double t) ->
double {
52 auto tSquaredOver4 = t * t / 4.0;
58 while(factor > 1.0/1000000.0) {
59 factor *= tSquaredOver4 / (k * k);
66 auto besselI_1 = [](
double t) ->
double {
67 auto tSquaredOver4 = t * t / 4.0;
71 while (factor > 1.0/1000000.0) {
72 factor *= tSquaredOver4 / (k * (k + 1));
82 auto d = std::exp(var);
95 b[n+1] = -(2*n/var) *
b[n] +
b[n-1];
96 gauss[n+1] =
b[n+1] /
d;
VULKAN_HPP_DEFAULT_DISPATCH_LOADER_DYNAMIC_STORAGE auto & d