adding vol2diam to aero_data with tests

Author: zdaq12

src/aero_data.F90

@@ -161,4 +161,16 @@ subroutine f_aero_data_diam2vol(ptr_c, diam, vol) bind(C)
 
 end  subroutine
 
+subroutine f_aero_data_vol2diam(ptr_c, vol, diam) bind(C)
+    type(aero_data_t), pointer :: ptr_f => null()
+    type(c_ptr), intent(in) :: ptr_c
+    real(c_double), intent(in) :: vol
+    real(c_double), intent(out) :: diam
+
+    call c_f_pointer(ptr_c, ptr_f)
+
+    diam = aero_data_vol2diam(ptr_f, vol)
+
+end subroutine
+
 end module

src/aero_data.hpp

@@ -23,6 +23,7 @@ extern "C" void f_aero_data_get_prime_radius(const void *ptr, double*) noexcept;
 extern "C" void f_aero_data_rad2vol(const void *ptr, const double*, double*) noexcept;
 extern "C" void f_aero_data_vol2rad(const void *ptr, const double*, double*) noexcept;
 extern "C" void f_aero_data_diam2vol(const void *ptr, const double*, double*) noexcept;
+extern "C" void f_aero_data_vol2diam(const void *ptr, const double*, double*) noexcept;
 
 struct AeroData {
     PMCResource ptr;
@@ -97,5 +98,11 @@ struct AeroData {
         return vol;
     }
 
+    static double vol2diam(const AeroData &self, const double vol) {
+        double diam;
+        f_aero_data_vol2diam(&self.ptr, &vol, &diam);
+        return diam;
+    }
+
 };
 

src/pypartmc.cpp

@@ -77,6 +77,8 @@ PYBIND11_MODULE(_PyPartMC, m) {
             "Convert mass-equivalent volume (m^3) to geometric radius (m)")
         .def("diam2vol", AeroData::diam2vol,
             "Convert geometric diameter (m) to mass-equivalent volume (m^3).")
+        .def("vol2diam", AeroData::vol2diam,
+            "Convert mass-equivalent volume (m^3) to geometric diameter (m).")
     ;
 
     py::class_<AeroParticle>(m, "AeroParticle",

tests/test_aero_data.py

@@ -197,3 +197,38 @@ def test_diam2vol_fractal(aero_data_params:dict):
         # assert
         np.testing.assert_almost_equal(value, (4/3)*np.pi*(sut.prime_radius)**3 *
             (1e-6/sut.prime_radius)**sut.frac_dim / sut.vol_fill_factor)
+
+    @staticmethod
+    def test_vol2diam_sphere():
+        # arrange
+        sut = ppmc.AeroData(AERO_DATA_CTOR_ARG_MINIMAL)
+        vol = 4.19e-18
+        sut.frac_dim = 3.0
+        sut.vol_fill_factor = 1.0
+        sut.prime_radius = 1e-8
+
+        # act
+        value = sut.vol2diam(vol)
+
+        # arrange
+        np.testing.assert_almost_equal(value, 2e-6)
+
+    @staticmethod
+    @pytest.mark.parametrize("aero_data_params", (
+        {"frac_dim": 2.4, "vol_fill_factor": 1.2, 'prime_radius': 1e-7},
+        {"frac_dim": 2.5, "vol_fill_factor": 1.1, "prime_radius": 1e-8},
+        {"frac_dim": 2.2, "vol_fill_factor": 1.3, 'prime_radius': 1e-6}))
+    def test_vol2diam_fractal(aero_data_params:dict):
+        # arrange
+        sut = ppmc.AeroData(AERO_DATA_CTOR_ARG_MINIMAL)
+        vol = 4.19e-18
+        for key,value in aero_data_params.items():
+            setattr(sut, key, value)
+
+        # act
+        value = sut.vol2diam(vol)
+
+        # assert
+        np.testing.assert_almost_equal(value, 2*(sut.prime_radius *
+            (((3*vol/4/np.pi)**(1/3)/sut.prime_radius)**3 * sut.vol_fill_factor)**
+            (1/sut.frac_dim)))