Switch to polar coordinates in the solver.

1 files changed, 6 insertions, 0 deletions

diff --git a/internal.h b/internal.h
index 501d390..372cc0d 100644
--- a/internal.h
+++ b/internal.h
@@ -38,8 +38,14 @@
  */
 #define EPS 1E-08
 
+enum CoordSystem {
+    COORD_SYSTEM_CYLINDRICAL,
+    COORD_SYSTEM_RADIAL,
+};
 
 typedef struct BDPriv {
+    int coord_system;
+
     BasisSetContext *basis[2];
     QFuncContext    *qfunc;