diff --git a/src/math.f90 b/src/math.f90
index 315da2642..e592ae1a9 100644
--- a/src/math.f90
+++ b/src/math.f90
@@ -452,7 +452,6 @@ pure function math_tensorproduct(A,B)
  implicit none
  real(pReal), dimension(:), intent(in) ::  A,B
  real(pReal), dimension(size(A,1),size(B,1)) ::  math_tensorproduct
-
  integer(pInt) :: i,j
 
  forall (i=1_pInt:size(A,1),j=1_pInt:size(B,1)) math_tensorproduct(i,j) = A(i)*B(j)