diff --git a/src/bvpy/vforms/elasticity.py b/src/bvpy/vforms/elasticity.py
index 3feae0b9f94a20fcef18dbf8cb55c3d86a7d1463..d611b43b5c527bccdaa626346f1e167cc45150c7 100644
--- a/src/bvpy/vforms/elasticity.py
+++ b/src/bvpy/vforms/elasticity.py
@@ -310,7 +310,7 @@ class HyperElasticForm(ElasticForm):
         elif self.model == 'neo_hookean':
             d = u.geometric_dimension()
             Id = fe.Identity(d)
-            F = Id + fe.nabla_grad(sol)
+            F = Id + fe.grad(sol)
             C = F.T * F
 
             Ic = fe.tr(C)