Fix sphinx and doxygen, add test debug statement, change tolerance

This commit fixes all sphinx or doxygen warnings, adds a statement to
debug the python reactor test, test_tolerances, and changes the problem
tolerance from 1e-6 to 1e-5. This test should still test what intended.
The change fixes the test for the MoleReactor class on Windows. Thought
about how tolerances are set for mole reactors is needed.

doc/sphinx/cython/zerodim.rst

@@ -42,6 +42,10 @@ Reactor
 ^^^^^^^
 .. autoclass:: Reactor
 
+MoleReactor
+^^^^^^^^^^^
+.. autoclass:: MoleReactor(contents=None, *, name=None, energy='on')
+
 IdealGasReactor
 ^^^^^^^^^^^^^^^
 .. autoclass:: IdealGasReactor(contents=None, *, name=None, energy='on')
@@ -54,6 +58,10 @@ ConstPressureReactor
 ^^^^^^^^^^^^^^^^^^^^
 .. autoclass:: ConstPressureReactor(contents=None, *, name=None, energy='on')
 
+ConstPressureMoleReactor
+^^^^^^^^^^^^^^^^^^^^^^^^
+.. autoclass:: ConstPressureMoleReactor(contents=None, *, name=None, energy='on')
+
 IdealGasConstPressureReactor
 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 .. autoclass:: IdealGasConstPressureReactor(contents=None, *, name=None, energy='on')

include/cantera/zeroD/MoleReactor.h

@@ -38,11 +38,11 @@ class MoleReactor : public Reactor
 
 protected:
     //! Get moles of the system from mass fractions stored by thermo object
-    //! @param y[out] vector for system moles to be put into
+    //! @param y vector for moles to be put into
     virtual void getMoles(double* y);
 
     //! Set internal mass variable based on moles given
-    //! @param y[in] vector of moles of the system
+    //! @param y vector of moles of the system
     virtual void setMassFromMoles(double* y);
 
     virtual void evalSurfaces(double* LHS, double* RHS, double* sdot);

test/python/test_reactor.py

@@ -220,10 +220,9 @@ def integrate(atol, rtol):
 
         n_baseline = integrate(1e-10, 1e-20)
         n_rtol = integrate(5e-7, 1e-20)
-        n_atol = integrate(1e-10, 1e-6)
-
-        self.assertTrue(n_baseline > n_rtol)
-        self.assertTrue(n_baseline > n_atol)
+        n_atol = integrate(1e-10, 1e-5)
+        self.assertTrue(n_baseline > n_rtol, f"{n_baseline} is less than {n_rtol}")
+        self.assertTrue(n_baseline > n_atol, f"{n_baseline} is less than {n_atol}")
 
     def test_advance_limits(self):
         P0 = 10 * ct.one_atm