FEAT: enable dipole6

Author: aaronjridley

src/quadtree.cpp

@@ -18,11 +18,14 @@ int64_t iProcQuery = -1;
 Quadtree::Quadtree(std::string shape) {
   if (shape == "cubesphere")
     nRootNodes = 6;
-
+  if (shape == "sphere")
+    nRootNodes = 1;
+  if (shape == "dipole")
+    nRootNodes = 1;
   if (shape == "dipole2")
     nRootNodes = 2;
-  else
-    nRootNodes = 1;
+  if (shape == "dipole6")
+    nRootNodes = 6;
 }
 
 // --------------------------------------------------------------------------
@@ -73,6 +76,13 @@ void Quadtree::build(std::string gridtype) {
     IsSphere = true;
   }
 
+  if (grid_input.shape == "dipole6") {
+    origins = Dipole6::ORIGINS;
+    rights = Dipole6::RIGHTS;
+    ups = Dipole6::UPS;
+    IsSphere = true;
+  }
+  
   arma_vec o(3), r(3), u(3);
 
   // This captures the limits of the sphere, independent of what the