Create Python bindings and build

Add a pybind.cpp file and use it to wrap the basic plasma source
functionality.

Update CMakeLists to build the Python bindings.

Ignore untracked Python artefacts in git.

Author: DanShort12

.gitignore

@@ -33,3 +33,7 @@
 
 # Build path
 build
+
+# Python
+__pycache__
+*.pyc

CMakeLists.txt

@@ -35,3 +35,13 @@ set_target_properties(source_sampling PROPERTIES POSITION_INDEPENDENT_CODE ON)
 target_include_directories(source_sampling PUBLIC ${OPENMC_INC_DIR})
 target_include_directories(source_sampling PUBLIC ${OPENMC_DIR}/vendor/pugixml)
 target_link_libraries(source_sampling ${OPENMC_LIB} gfortran)
+
+# Build plasma_source Python bindings
+list(APPEND plasma_source_pybind_SOURCES
+  ${SRC_DIR}/plasma_source.cpp
+  ${SRC_DIR}/plasma_source_pybind.cpp
+)
+
+add_subdirectory(pybind11)
+
+pybind11_add_module(plasma_source ${plasma_source_pybind_SOURCES})

parametric_plasma_source/plasma_source_pybind.cpp

@@ -0,0 +1,47 @@
+#include <pybind11/pybind11.h>
+#include "plasma_source.hpp"
+
+namespace py = pybind11;
+namespace ps = plasma_source;
+
+PYBIND11_MODULE(plasma_source, m) {
+    m.doc() = "A parametric plasma source";
+
+    py::class_<ps::PlasmaSource>(m, "PlasmaSource")
+        .def(py::init<const double &, const double &, const double &, const double &,
+                      const double &, const double &, const double &, const double &,
+                      const double &, const double &, const double &, const double &,
+                      const double &, const double &, const std::string &, const int &,
+	                  const int &, const double &, const double &>(),
+                          py::arg("ion_density_pedistal")=1.09e20,
+                          py::arg("ion_density_seperatrix")=3e19,
+                          py::arg("ion_density_origin")=1.09e20,
+                          py::arg("ion_temperature_pedistal")=6.09,
+                          py::arg("ion_temperature_seperatrix")=0.1,
+                          py::arg("ion_temperature_origin")=45.9,
+                          py::arg("pedistal_radius")=0.8,
+                          py::arg("ion_density_peaking_factor")=1.0,
+                          py::arg("ion_temperature_peaking_factor")=8.06,
+                          py::arg("minor_radius")=1.5,
+                          py::arg("major_radius")=4.5,
+                          py::arg("elongation")=2.0,
+                          py::arg("triangularity")=0.55,
+                          py::arg("shafranov_shift")=0.0,
+                          py::arg("plasma_type")="plasma",
+                          py::arg("plasma_id")=1,
+                          py::arg("number_of_bins")=100,
+                          py::arg("min_toroidal_angle") = 0.0,
+                          py::arg("max_toridal_angle") = 360.0)
+        .def("ion_density",
+             &ps::PlasmaSource::ion_density,
+             "Calculate the ion density at a specific minor radius",
+             py::arg("minor_radius"))
+        .def("ion_temperature",
+             &ps::PlasmaSource::ion_temperature,
+             "calculate the ion temperature at a specific minor radius",
+             py::arg("minor_radius"))
+        .def("dt_xs",
+             &ps::PlasmaSource::dt_xs,
+             "determine the value of the dt xs cross sections at a specific ion temperature",
+             py::arg("ion_temperature"));
+}