diff --git a/doc/src/compute_fabric.rst b/doc/src/compute_fabric.rst
index f1a5d3d7f9..b38ffafa48 100644
--- a/doc/src/compute_fabric.rst
+++ b/doc/src/compute_fabric.rst
@@ -146,13 +146,13 @@ m to :math:`M` (inclusive).  A middle asterisk means all types from m to n
 Output info
 """""""""""
 
-This compute calculates a local vector of doubles and a scalar. The vector
-stores the unique components of the first requested tensor in the order
-:math:`xx`, :math:`yy`, :math:`zz`, :math:`xy`, :math:`xz`, :math:`yz`
-followed by the same components for all subsequent tensors.
+This compute calculates a global vector of doubles and a global scalar. The
+vector stores the unique components of the first requested tensor in the
+order :math:`xx`, :math:`yy`, :math:`zz`, :math:`xy`, :math:`xz`,
+:math:`yz` followed by the same components for all subsequent tensors.
 The length of the vector is therefore six times the number of requested
-tensors. The scalar output is the number of pairwise interactions included in
-the calculation of the fabric tensor.
+tensors. The scalar output is the number of pairwise interactions included
+in the calculation of the fabric tensor.
 
 Restrictions
 """"""""""""
diff --git a/doc/src/pair_granular.rst b/doc/src/pair_granular.rst
index 0911a3486a..8a4d44fbf7 100644
--- a/doc/src/pair_granular.rst
+++ b/doc/src/pair_granular.rst
@@ -736,7 +736,7 @@ or
 
 .. math::
 
-   E_{eff,ij} = \frac{E_{ij}}{2(1-\nu_{ij})}
+   E_{eff,ij} = \frac{E_{ij}}{2(1-\nu_{ij}^2)}
 
 These pair styles write their information to :doc:`binary restart files <restart>`, so a pair_style command does not need to be
 specified in an input script that reads a restart file.
diff --git a/src/GRANULAR/gran_sub_mod_normal.cpp b/src/GRANULAR/gran_sub_mod_normal.cpp
index 05639feb9d..f08892d6f2 100644
--- a/src/GRANULAR/gran_sub_mod_normal.cpp
+++ b/src/GRANULAR/gran_sub_mod_normal.cpp
@@ -167,7 +167,7 @@ void GranSubModNormalHertzMaterial::coeffs_to_local()
 
 void GranSubModNormalHertzMaterial::mix_coeffs(double* icoeffs, double* jcoeffs)
 {
-  coeffs[0] = mix_stiffnessE(icoeffs[0], jcoeffs[0],icoeffs[2], jcoeffs[2]);
+  coeffs[0] = mix_stiffnessE(icoeffs[0], jcoeffs[0], icoeffs[2], jcoeffs[2]);
   coeffs[1] = mix_geom(icoeffs[1], jcoeffs[1]);
   coeffs[2] = mix_geom(icoeffs[2], jcoeffs[2]);
   coeffs_to_local();