diff --git a/src/coreComponents/constitutive/docs/CompositionalMultiphaseFluid.rst b/src/coreComponents/constitutive/docs/CompositionalMultiphaseFluid.rst
index 3a58539fd0..bd6874e706 100644
--- a/src/coreComponents/constitutive/docs/CompositionalMultiphaseFluid.rst
+++ b/src/coreComponents/constitutive/docs/CompositionalMultiphaseFluid.rst
@@ -179,13 +179,16 @@ The flash calculation process is as follows:
 #. These phase compositions are then used to calculate the component fugacities :math:`\phi_{iL}` and :math:`\phi_{iV}` in the liquid and vapor phases using the equation of state.
 
 #. Convergence is reached when the fugacities are equal for all components. The convergence criterion is defined as:
-.. math::
-    \sum_{i=1}^{N_c} \left( \phi_{iL} - \phi_{iV} \right)^2 < \varepsilon
+
+   .. math::
+       \sum_{i=1}^{N_c} \left( \phi_{iL} - \phi_{iV} \right)^2 < \varepsilon
+   
    where :math:`\varepsilon` is the convergence tolerance.
 
 #. If convergence is not achieved, successive substitution is used to update the set of K-values for the next iteration. The new K-values at iteration  :math:`t+1` are given by:
-.. math::
-    K_i^{(t+1)} = K_i^{(t)} \frac{\phi_{iL}}{\phi_{iV}}
+
+   .. math::
+       K_i^{(t+1)} = K_i^{(t)} \frac{\phi_{iL}}{\phi_{iV}}
 
 Parameters
 =========================
diff --git a/src/coreComponents/constitutive/fluid/multifluid/compositional/functions/StabilityTest.hpp b/src/coreComponents/constitutive/fluid/multifluid/compositional/functions/StabilityTest.hpp
index 9eda854596..22a1485f5f 100644
--- a/src/coreComponents/constitutive/fluid/multifluid/compositional/functions/StabilityTest.hpp
+++ b/src/coreComponents/constitutive/fluid/multifluid/compositional/functions/StabilityTest.hpp
@@ -52,7 +52,7 @@ struct StabilityTest
    * @param[out] kValues the k-values estimated from the stationary points
    * @return a flag indicating that 2 stationary points have been found
    */
-  template< integer USD1 >
+  template< integer USD1, integer USD2 >
   GEOS_HOST_DEVICE
   static bool compute( integer const numComps,
                        real64 const pressure,
@@ -61,7 +61,7 @@ struct StabilityTest
                        ComponentProperties::KernelWrapper const & componentProperties,
                        EquationOfStateType const & equationOfState,
                        real64 & tangentPlaneDistance,
-                       arraySlice1d< real64 > const & kValues )
+                       arraySlice1d< real64, USD2 > const & kValues )
   {
     constexpr integer numTrials = 2;    // Trial compositions
     stackArray2d< real64, 4*maxNumComps > workSpace( 4, numComps );