libsequence  1.9.5
auxillary.cc
1 #include <cstdint>
2 #include <cmath>
3 
4 namespace Sequence
5 {
6  namespace summstats_aux
7  {
8  double
9  a_sub_n(const std::uint32_t nsam)
10  {
11  double rv = 0.0;
12  for (std::uint32_t i = 1; i < nsam; ++i)
13  {
14  rv += 1.0 / static_cast<double>(i);
15  }
16  return rv;
17  }
18 
19  double
20  b_sub_n(const std::uint32_t nsam)
21  {
22  double rv = 0.0;
23  for (std::uint32_t i = 1; i < nsam; ++i)
24  {
25  rv += 1.0 / std::pow(static_cast<double>(i), 2.0);
26  }
27  return rv;
28  }
29 
30  double
31  b_sub_n_plus1(const std::uint32_t nsam)
32  {
33  double rv = 0.0;
34  for (std::uint32_t i = 1; i < nsam + 1; ++i)
35  {
36  rv += 1.0 / std::pow(static_cast<double>(i), 2.0);
37  }
38  return rv;
39  }
40  } // namespace summstats_aux
41 } // namespace Sequence
The namespace in which this library resides.