Template creation 2 store GaussLegendre weights

src/demeter/quadrature/:

@@ -0,0 +1,73 @@
+#include "gauss_legendre.hpp"
+
+#include <cmath>
+#include <iostream>
+#include <Eigen/Dense>
+
+// To do
+// 1 - Definir une structure qui permettent d'aller chercher
+//     les bons tableaux selon le nombre d'angles polaires
+// 2 - Simplifier la structure car plusieurs fois les memes 
+//     valeurs a un signe pres
+// 3 - Automatiser le calcul des racines et des poids pour
+//     aller jusqua 48
+// Tous les poids sont a diviser par deux
+// teta est dans [0,pi] pas besoin des angles teta
+// dans [pi,2pi] on fait teta_k+pi avec 
+// w_{k+pi} = w_k  
+
+namespace Demeter {
+
+
+// N = 2
+template <>
+const auto Eigen::ArrayXd::GaussLegendre<2>::test_ << 2.000000; 
+// nteta = 1
+ //ArrayXd  costeta(1);
+ //ArrayXd  weights(1);
+
+//  costeta << 0.00000;
+//  weights << 2.00000;
+
+//// nteta = 2
+// ArrayXd  costeta(2);
+// ArrayXd  weights(2);
+//
+// costeta << -0.57735, 0.57735;
+// weights <<  1.00000, 1.00000;
+//
+//// nteta = 4
+// ArrayXd  costeta(3);
+// ArrayXd  weights(3);
+//
+// costeta << -0.339981, 0.339981,
+//            -0.861136, 0.861136;
+// weights <<  0.652145, 0.652145,
+//             0.347855, 0.347855;
+
+
+// nteta = 6
+// nteta = 8
+// nteta = 10
+// nteta = 12
+// nteta = 14
+// nteta = 16
+// nteta = 18
+// nteta = 20
+// nteta = 22
+// nteta = 24
+// nteta = 26
+// nteta = 28
+// nteta = 30
+// nteta = 32
+// nteta = 34
+// nteta = 36
+// nteta = 38
+// nteta = 40
+// nteta = 42
+// nteta = 44
+// nteta = 46
+// nteta = 48
+
+} // namespace Demeter
+

src/demeter/quadrature/:!

@@ -0,0 +1,20 @@
+#pragma once
+
+
+#include <Eigen/Core>
+#include <Eigen/Dense>
+
+#include <array>
+#include <span>
+
+namespace Demeter {
+
+template <std::size_t N>
+class GaussLegendre {
+ public:
+  std::span<const double> wgt() const { return {weights_.begin(), weights_.end()}; }
+
+ private:
+  static const std::array<double, N / 2> weights_;
+};
+} // namespace Demeter

src/demeter/quadrature/gauss_legendre.cpp

@@ -5,28 +5,29 @@
 #include <Eigen/Dense>
 
 // To do
-// 1 - Definir une structure qui permettent d'aller chercher
-//     les bons tableaux selon le nombre d'angles polaires
-// 2 - Simplifier la structure car plusieurs fois les memes 
-//     valeurs a un signe pres
-// 3 - Automatiser le calcul des racines et des poids pour
-//     aller jusqua 48
+// 1 - Automatiser le calcul des racines et des poids pour
+//     aller jusqua 48 (python)
 // Tous les poids sont a diviser par deux
 // teta est dans [0,pi] pas besoin des angles teta
 // dans [pi,2pi] on fait teta_k+pi avec 
 // w_{k+pi} = w_k  
 
 namespace Demeter {
 
-void GaussLegendre::compute_angles() {
-  // to do
-} 
+// nteta = 2
+template <>
+Eigen::ArrayXd GaussLegendre<2>::weights_(1); 
+template <>
+Eigen::ArrayXd GaussLegendre<2>::costeta_(1); 
+template <>
+Eigen::ArrayXd GaussLegendre<2>::teta_(1); 
+template <>
+void GaussLegendre<2>::initialization() {
+    weights_ << 2.00000;
+    costeta_ << 0.00000;
+    teta_    << 0.00000;
+}
 
-void GaussLegendre::compute_weights() {
-  // to do
-} 
-
-
 // nteta = 1
  //ArrayXd  costeta(1);
  //ArrayXd  weights(1);
@@ -75,4 +76,3 @@ void GaussLegendre::compute_weights() {
 // nteta = 48
 
 } // namespace Demeter
-

src/demeter/quadrature/gauss_legendre.hpp

@@ -1,24 +1,26 @@
 #pragma once
 
-
-#include <Eigen/Core>
 #include <Eigen/Dense>
 
 namespace Demeter {
 
+template <std::size_t nteta>
 class GaussLegendre {
  public:
+  // Define array initializer
+  static void initialization();
 
-  GaussLegendre(const int& nphi, const int& nteta) 
-     : nphi_(nphi), nteta_(nteta) {
-  }
+  // Define getter methods
+  const Eigen::ArrayXd& getweights() const { return weights_; }
+  const Eigen::ArrayXd& getcosteta() const { return costeta_; }
+  const Eigen::ArrayXd& getteta() const { return teta_; }
 
-  void compute_angles();
-  void compute_weights();
 
  private:
-  size_t nphi_;  // nber of azimutal angles
-  size_t nteta_; // nber of polar    angles
-
+  static Eigen::ArrayXd weights_;
+  static Eigen::ArrayXd costeta_;
+  static Eigen::ArrayXd teta_;
 };
+
+
 } // namespace Demeter

src/demeter/quadrature/w

@@ -0,0 +1,77 @@
+#include "gauss_legendre.hpp"
+
+#include <cmath>
+#include <iostream>
+#include <Eigen/Dense>
+
+// To do
+// 1 - Definir une structure qui permettent d'aller chercher
+//     les bons tableaux selon le nombre d'angles polaires
+// 2 - Simplifier la structure car plusieurs fois les memes 
+//     valeurs a un signe pres
+// 3 - Automatiser le calcul des racines et des poids pour
+//     aller jusqua 48
+// Tous les poids sont a diviser par deux
+// teta est dans [0,pi] pas besoin des angles teta
+// dans [pi,2pi] on fait teta_k+pi avec 
+// w_{k+pi} = w_k  
+
+namespace Demeter {
+
+
+// N = 2
+template <>
+
+const auto GaussLegendre<2>::weights_ = flux_ = ArrayXd::Ones(2);
+const std::array<double, 1> GaussLegendre<2>::wgt_ = {
+    1.00000000000000000000000000000000000e+00};
+
+// nteta = 1
+ //ArrayXd  costeta(1);
+ //ArrayXd  weights(1);
+
+//  costeta << 0.00000;
+//  weights << 2.00000;
+
+//// nteta = 2
+// ArrayXd  costeta(2);
+// ArrayXd  weights(2);
+//
+// costeta << -0.57735, 0.57735;
+// weights <<  1.00000, 1.00000;
+//
+//// nteta = 4
+// ArrayXd  costeta(3);
+// ArrayXd  weights(3);
+//
+// costeta << -0.339981, 0.339981,
+//            -0.861136, 0.861136;
+// weights <<  0.652145, 0.652145,
+//             0.347855, 0.347855;
+
+
+// nteta = 6
+// nteta = 8
+// nteta = 10
+// nteta = 12
+// nteta = 14
+// nteta = 16
+// nteta = 18
+// nteta = 20
+// nteta = 22
+// nteta = 24
+// nteta = 26
+// nteta = 28
+// nteta = 30
+// nteta = 32
+// nteta = 34
+// nteta = 36
+// nteta = 38
+// nteta = 40
+// nteta = 42
+// nteta = 44
+// nteta = 46
+// nteta = 48
+
+} // namespace Demeter
+