Make every header self-contained

include/integratorxx/batch/spherical_micro_batcher.hpp

@@ -7,6 +7,14 @@
 #include <iomanip>
 #include <vector>
 #include <array>
+#include <algorithm>
+#include <cmath>
+#include <cstddef>
+#include <iterator>
+#include <memory>
+#include <stdexcept>
+#include <tuple>
+#include <utility>
 
 
 namespace IntegratorXX {

include/integratorxx/composite_quadratures/product_quadrature.hpp

@@ -1,6 +1,8 @@
 #pragma once
 
 #include <integratorxx/quadrature.hpp>
+#include <memory>
+#include <vector>
 
 namespace IntegratorXX {
 

include/integratorxx/composite_quadratures/pruned_spherical_quadrature.hpp

@@ -6,6 +6,9 @@
 #include <integratorxx/type_traits.hpp>
 
 #include <utility>
+#include <cassert>
+#include <memory>
+#include <vector>
 
 namespace IntegratorXX {
 

include/integratorxx/composite_quadratures/spherical_quadrature.hpp

@@ -3,6 +3,7 @@
 #include <integratorxx/composite_quadratures/product_quadrature.hpp>
 #include <integratorxx/types.hpp>
 #include <integratorxx/type_traits.hpp>
+#include <memory>
 
 namespace IntegratorXX {
 

include/integratorxx/composite_quadratures/sub_quadrature.hpp

@@ -2,6 +2,8 @@
 
 #include <integratorxx/quadrature.hpp>
 #include <cassert>
+#include <tuple>
+#include <utility>
 
 namespace IntegratorXX {
 

include/integratorxx/quadrature.hpp

@@ -4,6 +4,7 @@
 #include <cstdint>
 #include <cstddef>
 #include <cassert>
+#include <utility>
 #include <integratorxx/type_traits.hpp>
 
 namespace IntegratorXX {

include/integratorxx/quadratures/primitive/gausschebyshev1.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadrature.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/primitive/gausschebyshev2.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadrature.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/primitive/gausschebyshev2modified.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadrature.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/primitive/gausschebyshev3.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadrature.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/primitive/gausslegendre.hpp

@@ -3,6 +3,10 @@
 #include <integratorxx/quadrature.hpp>
 #include <integratorxx/util/legendre.hpp>
 #include <vector>
+#include <cmath>
+#include <limits>
+#include <stdexcept>
+#include <tuple>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/primitive/gausslobatto.hpp

@@ -3,6 +3,10 @@
 #include <integratorxx/quadrature.hpp>
 #include <integratorxx/util/legendre.hpp>
 #include <vector>
+#include <cmath>
+#include <limits>
+#include <stdexcept>
+#include <tuple>
 
 namespace IntegratorXX {
 
@@ -47,7 +51,7 @@ struct quadrature_traits<GaussLobatto<PointType, WeightType>> {
     for(size_t idx = 1; idx < mid; ++idx) {
       // Initial guess
       const point_type i = static_cast<point_type>(npts - 1 - idx);
-      point_type z = cos (i * M_PI / ( npts - 1.0));
+      point_type z = std::cos(i * M_PI / (npts - 1.0));
 
       // Old value of root
       point_type z_old = -2.0;

include/integratorxx/quadratures/radial/becke.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadratures/primitive/gausschebyshev2.hpp>
 #include <integratorxx/quadratures/radial/radial_transform.hpp>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/radial/mhl.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadratures/primitive/uniform.hpp>
 #include <integratorxx/quadratures/radial/radial_transform.hpp>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/radial/muraknowles.hpp

@@ -3,6 +3,7 @@
 #include <integratorxx/quadratures/primitive/uniform.hpp>
 #include <integratorxx/quadratures/radial/radial_transform.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/radial/radial_transform.hpp

@@ -1,6 +1,8 @@
 #pragma once
 
 #include <memory>
+#include <cstddef>
+#include <tuple>
 
 #include <integratorxx/quadrature.hpp>
 

include/integratorxx/quadratures/radial/treutlerahlrichs.hpp

@@ -2,6 +2,7 @@
 
 #include <integratorxx/quadratures/primitive/gausschebyshev2.hpp>
 #include <integratorxx/quadratures/radial/radial_transform.hpp>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/s2/ahrens_beylkin.hpp

@@ -4,6 +4,7 @@
 #include <integratorxx/quadratures/s2/ahrens_beylkin/ahrens_beylkin_grids.hpp>
 #include <integratorxx/util/copy_grid.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/s2/delley.hpp

@@ -4,6 +4,7 @@
 #include <integratorxx/util/copy_grid.hpp>
 #include <integratorxx/quadratures/s2/delley/delley_grids.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/s2/lebedev_laikov.hpp

@@ -4,6 +4,7 @@
 #include <integratorxx/quadratures/s2/lebedev_laikov/lebedev_laikov_grids.hpp>
 #include <integratorxx/util/copy_grid.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/quadratures/s2/womersley.hpp

@@ -5,6 +5,7 @@
 #include <integratorxx/util/create_array.hpp>
 #include <integratorxx/quadratures/s2/womersley/womersley_grids.hpp>
 #include <vector>
+#include <cmath>
 
 namespace IntegratorXX {
 

include/integratorxx/util/bound_transform.hpp

@@ -1,4 +1,6 @@
 #pragma once
+#include <cassert>
+#include <tuple>
 
 #include <limits>
 

include/integratorxx/util/copy_grid.hpp

@@ -1,4 +1,5 @@
 #pragma once
+#include <algorithm>
 
 namespace IntegratorXX {
 

include/integratorxx/util/create_array.hpp

@@ -1,4 +1,7 @@
 #pragma once
+#include <array>
+#include <cstddef>
+#include <utility>
 
 namespace IntegratorXX {
 

include/integratorxx/util/legendre.hpp

@@ -1,4 +1,5 @@
 #pragma once
+#include <cstddef>
 #include <tuple>
 
 namespace IntegratorXX {