diff --git a/README.md b/README.md
index a48be52..b4610e9 100644
--- a/README.md
+++ b/README.md
@@ -139,7 +139,11 @@ Use:
 stress = cutde.fullspace.strain_to_stress(strain, sm, nu)
 ```
 
-to convert from stress to strain assuming isotropic linear elasticity. `sm` is the shear modulus and `nu` is the Poisson ratio.
+to convert from stress to strain assuming isotropic linear elasticity. `sm` is the shear modulus and `nu` is the Poisson ratio. The `strain` array here is expected to have a shape `(N, 6)`. In case you have a matrix or other shaped array, this snippet might help:
+
+```
+strain_to_stress(strain_matrix.reshape((-1, 6)), mu, nu).reshape(strain_matrix.shape)
+```
 
 ## All pairs interactions matrix