Gauss-Sidel iteration SlaybaughLab#10

wzheng21 · Sep 10, 2017 · 2cfc091 · 2cfc091
1 parent d2d164f
commit 2cfc091
Show file tree

Hide file tree

Showing 6 changed files with 42 additions and 0 deletions.
diff --git a/src/transport/equation_base.cc → src/equation/equation_base.cc b/src/transport/equation_base.cc → src/equation/equation_base.cc
diff --git a/src/transport/equation_base.h → src/equation/equation_base.h b/src/transport/equation_base.h → src/equation/equation_base.h
diff --git a/src/transport/even_parity.cc → src/equation/even_parity.cc b/src/transport/even_parity.cc → src/equation/even_parity.cc
diff --git a/src/transport/even_parity.h → src/equation/even_parity.h b/src/transport/even_parity.h → src/equation/even_parity.h
diff --git a/src/iteration/gauss_sidel.cc b/src/iteration/gauss_sidel.cc
@@ -0,0 +1,24 @@
+#include "mg_base.h"
+
+template <int dim>
+GaussSidel<dim>::GaussSidel (ParameterHandler &prm)
+:
+MGBase<dim> (prm)
+{
+}
+
+template <int dim>
+GaussSidel<dim>::~GaussSidel ()
+{
+}
+
+template <int dim>
+void GaussSidel<dim>::mg_iterations
+(std::vector<Vector<double> > &sflxes_proc,
+ std::vector<std_cxx11::shared_ptr<EquationBase<dim> > > &equ_ptrs)
+{
+  // TODO: fill this up
+}
+
+template class GaussSidel<2>;
+template class GaussSidel<3>;
diff --git a/src/iteration/gauss_sidel.h b/src/iteration/gauss_sidel.h
@@ -0,0 +1,18 @@
+#ifndef __gauss_sidel_h__
+#define __gauss_sidel_h__
+
+#include "mg_base.h"
+
+template <int dim>
+class GaussSidel : public MGBase<dim>
+{
+public:
+  GaussSidel (ParameterHandler);
+  ~GaussSidel ();
+
+  void mg_iterations
+  (std::vector<Vector<double> > &sflxes_proc,
+   std::vector<std_cxx11::shared_ptr<EquationBase<dim> > > equ_ptrs);
+}
+
+#endif //__gauss_sidel_h__