Starting the implementation of Gauss Legendre quad

src/demeter/quadrature.hpp

@@ -0,0 +1,3 @@
+#pragma once
+
+#include "solve/gauss_legendre.hpp"

src/demeter/quadrature/gauss_legendre.cpp

@@ -0,0 +1,16 @@
+#include "gauss_legendre.hpp"
+
+#include <cmath>
+#include <iostream>
+#include <Eigen/Dense>
+
+namespace Demeter {
+
+void GaussLegendre::compute_angles() {
+  // to do
+} 
+
+void GaussLegendre::compute_weights() {
+  // to do
+} 
+} // namespace Demeter

src/demeter/quadrature/gauss_legendre.hpp

@@ -0,0 +1,24 @@
+#pragma once
+
+
+#include <Eigen/Core>
+#include <Eigen/Dense>
+
+namespace Demeter {
+
+class GaussLegendre {
+ public:
+
+  GaussLegendre(const int& nphi, const int& nteta) 
+     : nphi_(nphi), nteta_(nteta) {
+  }
+
+  void compute_angles();
+  void compute_weights();
+
+ private:
+  int nphi_;  // nber of azimutal angles
+  int nteta_; // nber of polar    angles
+
+};
+} // namespace Demeter

src/demeter/solve/linear_solver.hpp

@@ -12,7 +12,6 @@ namespace Demeter {
 class LinearSolver {
  public: 
 
-  // Constructeur vide
   LinearSolver(const Eigen::Matrix2f& A, const Eigen::Matrix2f& b) 
       : A_(A), b_(b) {
       // Initialisation par défaut