From 1cf6cc217db13af0d10b1e84546084c6a59dbd46 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Fri, 10 May 2019 15:08:37 +0200
Subject: [PATCH] one implicit none is enough

---
 src/Lambert.f90 |  4 ---
 src/list.f90    | 14 ----------
 src/math.f90    | 74 -------------------------------------------------
 3 files changed, 92 deletions(-)

diff --git a/src/Lambert.f90 b/src/Lambert.f90
index 40bd04a01..9b7d2ef86 100644
--- a/src/Lambert.f90
+++ b/src/Lambert.f90
@@ -44,7 +44,6 @@ module Lambert
     PI
   use future
 
-  implicit none
   private
   real(pReal), parameter, private :: &
     SPI  = sqrt(PI), &
@@ -77,7 +76,6 @@ function LambertCubeToBall(cube) result(ball)
   use prec, only: &
     dEq0
  
-  implicit none
   real(pReal), intent(in), dimension(3) :: cube
   real(pReal),             dimension(3) :: ball, LamXYZ, XYZ
   real(pReal),             dimension(2) :: T
@@ -141,7 +139,6 @@ pure function LambertBallToCube(xyz) result(cube)
    use math, only: &
     math_clip
  
-  implicit none
   real(pReal), intent(in), dimension(3) :: xyz
   real(pReal),             dimension(3) :: cube, xyz1, xyz3
   real(pReal),             dimension(2) :: Tinv, xyz2
@@ -196,7 +193,6 @@ end function LambertBallToCube
 !--------------------------------------------------------------------------
 pure function GetPyramidOrder(xyz)
 
-  implicit none
   real(pReal),intent(in),dimension(3) :: xyz
   integer,               dimension(3) :: GetPyramidOrder
  
diff --git a/src/list.f90 b/src/list.f90
index 1c4c4a8a0..be80b151d 100644
--- a/src/list.f90
+++ b/src/list.f90
@@ -70,7 +70,6 @@ subroutine add(this,string)
     IO_lc, &
     IO_stringPos
 
-  implicit none
   class(tPartitionedStringList),  target, intent(in) :: this
   character(len=*),                       intent(in) :: string
   type(tPartitionedStringList),   pointer            :: new, temp
@@ -95,7 +94,6 @@ end subroutine add
 !--------------------------------------------------------------------------------------------------
 subroutine show(this)
 
-  implicit none
   class(tPartitionedStringList), target, intent(in) :: this
   type(tPartitionedStringList),  pointer            :: item
 
@@ -114,7 +112,6 @@ end subroutine show
 !--------------------------------------------------------------------------------------------------
 subroutine free(this)
 
-  implicit none
   class(tPartitionedStringList),  intent(inout) :: this
 
   if(associated(this%next)) deallocate(this%next)
@@ -128,7 +125,6 @@ end subroutine free
 !--------------------------------------------------------------------------------------------------
 recursive subroutine finalize(this)
 
-  implicit none
   type(tPartitionedStringList),  intent(inout) :: this
 
   if(associated(this%next)) deallocate(this%next)
@@ -142,7 +138,6 @@ end subroutine finalize
 !--------------------------------------------------------------------------------------------------
 subroutine finalizeArray(this)
 
-  implicit none
   integer :: i
   type(tPartitionedStringList),  intent(inout), dimension(:) :: this
   type(tPartitionedStringList),  pointer :: temp ! bug in Gfortran?
@@ -165,7 +160,6 @@ logical function keyExists(this,key)
   use IO, only: &
     IO_stringValue
 
-  implicit none
   class(tPartitionedStringList), target, intent(in) :: this
   character(len=*),                      intent(in) :: key
   type(tPartitionedStringList),  pointer            :: item
@@ -189,8 +183,6 @@ integer function countKeys(this,key)
   use IO, only: &
     IO_stringValue
 
-  implicit none
-
   class(tPartitionedStringList), target, intent(in) :: this
   character(len=*),                      intent(in) :: key
   type(tPartitionedStringList),  pointer            :: item
@@ -218,7 +210,6 @@ real(pReal) function getFloat(this,key,defaultVal)
     IO_stringValue, &
     IO_FloatValue
 
-  implicit none
   class(tPartitionedStringList), target, intent(in)           :: this
   character(len=*),                      intent(in)           :: key
   real(pReal),                           intent(in), optional :: defaultVal
@@ -255,7 +246,6 @@ integer function getInt(this,key,defaultVal)
     IO_stringValue, &
     IO_IntValue
 
-  implicit none
   class(tPartitionedStringList), target, intent(in)           :: this
   character(len=*),                      intent(in)           :: key
   integer,                               intent(in), optional :: defaultVal
@@ -292,7 +282,6 @@ character(len=65536) function getString(this,key,defaultVal,raw)
     IO_error, &
     IO_stringValue
  
-  implicit none
   class(tPartitionedStringList), target, intent(in)           :: this
   character(len=*),                      intent(in)           :: key
   character(len=*),                      intent(in), optional :: defaultVal
@@ -343,7 +332,6 @@ function getFloats(this,key,defaultVal,requiredSize)
     IO_stringValue, &
     IO_FloatValue
  
-  implicit none
   real(pReal),     dimension(:), allocatable          :: getFloats
   class(tPartitionedStringList), target, intent(in)   :: this
   character(len=*),              intent(in)           :: key
@@ -393,7 +381,6 @@ function getInts(this,key,defaultVal,requiredSize)
     IO_stringValue, &
     IO_IntValue
  
-  implicit none
   integer, dimension(:), allocatable                          :: getInts
   class(tPartitionedStringList), target, intent(in)           :: this
   character(len=*),                      intent(in)           :: key
@@ -443,7 +430,6 @@ function getStrings(this,key,defaultVal,raw)
     IO_error, &
     IO_StringValue
  
-  implicit none
   character(len=65536),dimension(:), allocatable             :: getStrings
   class(tPartitionedStringList),target, intent(in)           :: this
   character(len=*),                     intent(in)           :: key
diff --git a/src/math.f90 b/src/math.f90
index c3aab9936..dcd1aba8f 100644
--- a/src/math.f90
+++ b/src/math.f90
@@ -158,7 +158,6 @@ subroutine math_init
  use numerics, only: &
    randomSeed
 
- implicit none
  integer :: i
  real(pReal), dimension(4) :: randTest
  integer :: randSize
@@ -200,7 +199,6 @@ subroutine math_check
  use prec,     only: tol_math_check
  use IO,       only: IO_error
 
- implicit none
  character(len=64) :: error_msg
 
  ! +++ check vector expansion +++
@@ -230,7 +228,6 @@ end subroutine math_check
 !--------------------------------------------------------------------------------------------------
 recursive subroutine math_sort(a, istart, iend, sortDim)
 
- implicit none
  integer, dimension(:,:), intent(inout) :: a
  integer, intent(in),optional :: istart,iend, sortDim
  integer :: ipivot,s,e,d
@@ -267,7 +264,6 @@ recursive subroutine math_sort(a, istart, iend, sortDim)
  !-------------------------------------------------------------------------------------------------
  integer function qsort_partition(a, istart, iend, sort)
  
-   implicit none
    integer, dimension(:,:), intent(inout) :: a
    integer,                 intent(in)    :: istart,iend,sort
    integer, dimension(size(a,1))          :: tmp
@@ -307,7 +303,6 @@ end subroutine math_sort
 !--------------------------------------------------------------------------------------------------
 pure function math_expand(what,how)
 
- implicit none
  real(pReal),   dimension(:), intent(in) :: what
  integer, dimension(:), intent(in) :: how
  real(pReal), dimension(sum(how)) ::  math_expand
@@ -328,7 +323,6 @@ end function math_expand
 !--------------------------------------------------------------------------------------------------
 pure function math_range(N)
 
- implicit none
  integer, intent(in) :: N                                                                           !< length of range
  integer :: i
  integer, dimension(N) :: math_range
@@ -343,7 +337,6 @@ end function math_range
 !--------------------------------------------------------------------------------------------------
 pure function math_identity2nd(dimen)
 
- implicit none
  integer, intent(in) :: dimen                                                                       !< tensor dimension
  integer :: i
  real(pReal), dimension(dimen,dimen) :: math_identity2nd
@@ -359,7 +352,6 @@ end function math_identity2nd
 !--------------------------------------------------------------------------------------------------
 pure function math_identity4th(dimen)
 
- implicit none
  integer, intent(in) :: dimen                                                                       !< tensor dimension
  integer :: i,j,k,l
  real(pReal), dimension(dimen,dimen,dimen,dimen) ::  math_identity4th
@@ -380,7 +372,6 @@ end function math_identity4th
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_civita(i,j,k)
 
- implicit none
  integer, intent(in) :: i,j,k
 
  math_civita = 0.0_pReal
@@ -402,7 +393,6 @@ end function math_civita
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_delta(i,j)
 
- implicit none
  integer, intent (in) :: i,j
 
  math_delta = merge(0.0_pReal, 1.0_pReal, i /= j)
@@ -415,7 +405,6 @@ end function math_delta
 !--------------------------------------------------------------------------------------------------
 pure function math_cross(A,B)
 
- implicit none
  real(pReal), dimension(3), intent(in) ::  A,B
  real(pReal), dimension(3) :: math_cross
 
@@ -431,7 +420,6 @@ end function math_cross
 !--------------------------------------------------------------------------------------------------
 pure function math_outer(A,B)
 
- implicit none
  real(pReal), dimension(:), intent(in) ::  A,B
  real(pReal), dimension(size(A,1),size(B,1)) ::  math_outer
  integer :: i,j
@@ -446,7 +434,6 @@ end function math_outer
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_inner(A,B)
 
- implicit none
  real(pReal), dimension(:),         intent(in) :: A
  real(pReal), dimension(size(A,1)), intent(in) :: B
 
@@ -460,7 +447,6 @@ end function math_inner
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_mul33xx33(A,B)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) ::  A,B
  integer :: i,j
  real(pReal), dimension(3,3) ::              C
@@ -476,7 +462,6 @@ end function math_mul33xx33
 !--------------------------------------------------------------------------------------------------
 pure function math_mul3333xx33(A,B)
 
- implicit none
  real(pReal), dimension(3,3) :: math_mul3333xx33
  real(pReal), dimension(3,3,3,3), intent(in) ::  A
  real(pReal), dimension(3,3), intent(in) ::  B
@@ -492,7 +477,6 @@ end function math_mul3333xx33
 !--------------------------------------------------------------------------------------------------
 pure function math_mul3333xx3333(A,B)
 
- implicit none
  integer :: i,j,k,l
  real(pReal), dimension(3,3,3,3), intent(in) ::  A
  real(pReal), dimension(3,3,3,3), intent(in) ::  B
@@ -509,7 +493,6 @@ end function math_mul3333xx3333
 !--------------------------------------------------------------------------------------------------
 pure function math_mul33x33(A,B)
 
- implicit none
  real(pReal), dimension(3,3) ::  math_mul33x33
  real(pReal), dimension(3,3), intent(in) ::  A,B
  integer :: i,j
@@ -524,7 +507,6 @@ end function math_mul33x33
 !--------------------------------------------------------------------------------------------------
 pure function math_mul66x66(A,B)
 
- implicit none
  real(pReal), dimension(6,6) ::  math_mul66x66
  real(pReal), dimension(6,6), intent(in) ::  A,B
  integer :: i,j
@@ -541,7 +523,6 @@ end function math_mul66x66
 !--------------------------------------------------------------------------------------------------
 pure function math_mul99x99(A,B)
 
- implicit none
  real(pReal), dimension(9,9) ::  math_mul99x99
  real(pReal), dimension(9,9), intent(in) ::  A,B
  integer  i,j
@@ -559,7 +540,6 @@ end function math_mul99x99
 !--------------------------------------------------------------------------------------------------
 pure function math_mul33x3(A,B)
 
- implicit none
  real(pReal), dimension(3) ::  math_mul33x3
  real(pReal), dimension(3,3), intent(in) ::  A
  real(pReal), dimension(3),   intent(in) ::  B
@@ -575,7 +555,6 @@ end function math_mul33x3
 !--------------------------------------------------------------------------------------------------
 pure function math_mul33x3_complex(A,B)
 
- implicit none
  complex(pReal), dimension(3) ::  math_mul33x3_complex
  complex(pReal), dimension(3,3), intent(in) ::  A
  real(pReal),    dimension(3),   intent(in) ::  B
@@ -591,7 +570,6 @@ end function math_mul33x3_complex
 !--------------------------------------------------------------------------------------------------
 pure function math_mul66x6(A,B)
 
- implicit none
  real(pReal), dimension(6) ::  math_mul66x6
  real(pReal), dimension(6,6), intent(in) ::  A
  real(pReal), dimension(6),   intent(in) ::  B
@@ -608,7 +586,6 @@ end function math_mul66x6
 !--------------------------------------------------------------------------------------------------
 pure function math_exp33(A,n)
 
- implicit none
  integer :: i
  integer, intent(in), optional :: n
  real(pReal), dimension(3,3), intent(in) :: A
@@ -637,7 +614,6 @@ pure function math_inv33(A)
  use prec, only: &
    dNeq0
 
- implicit none
  real(pReal),dimension(3,3),intent(in)  :: A
  real(pReal) :: DetA
  real(pReal),dimension(3,3) :: math_inv33
@@ -675,7 +651,6 @@ pure subroutine math_invert33(A, InvA, DetA, error)
  use prec, only: &
    dEq0
 
- implicit none
  logical, intent(out) :: error
  real(pReal),dimension(3,3),intent(in)  :: A
  real(pReal),dimension(3,3),intent(out) :: InvA
@@ -713,7 +688,6 @@ function math_invSym3333(A)
  use IO, only: &
    IO_error
 
- implicit none
  real(pReal),dimension(3,3,3,3)            :: math_invSym3333
 
  real(pReal),dimension(3,3,3,3),intent(in) :: A
@@ -744,7 +718,6 @@ end function math_invSym3333
 !--------------------------------------------------------------------------------------------------
 subroutine math_invert2(InvA, error, A)
 
- implicit none
  real(pReal), dimension(:,:), intent(in)  :: A
  
  real(pReal), dimension(size(A,1),size(A,1)), intent(out) :: invA
@@ -762,7 +735,6 @@ end subroutine math_invert2
 !--------------------------------------------------------------------------------------------------
 subroutine math_invert(myDim,A, InvA, error)
 
- implicit none
  integer, intent(in)  :: myDim
  real(pReal), dimension(myDim,myDim), intent(in)  :: A
 
@@ -790,7 +762,6 @@ end subroutine math_invert
 !--------------------------------------------------------------------------------------------------
 pure function math_symmetric33(m)
 
- implicit none
  real(pReal), dimension(3,3) :: math_symmetric33
  real(pReal), dimension(3,3), intent(in) :: m
  
@@ -804,7 +775,6 @@ end function math_symmetric33
 !--------------------------------------------------------------------------------------------------
 pure function math_symmetric66(m)
 
- implicit none
  real(pReal), dimension(6,6) :: math_symmetric66
  real(pReal), dimension(6,6), intent(in) :: m
 
@@ -818,7 +788,6 @@ end function math_symmetric66
 !--------------------------------------------------------------------------------------------------
 pure function math_skew33(m)
 
- implicit none
  real(pReal), dimension(3,3) :: math_skew33
  real(pReal), dimension(3,3), intent(in) :: m
 
@@ -831,7 +800,6 @@ end function math_skew33
 !--------------------------------------------------------------------------------------------------
 pure function math_spherical33(m)
 
- implicit none
  real(pReal), dimension(3,3) :: math_spherical33
  real(pReal), dimension(3,3), intent(in) :: m
 
@@ -845,7 +813,6 @@ end function math_spherical33
 !--------------------------------------------------------------------------------------------------
 pure function math_deviatoric33(m)
 
- implicit none
  real(pReal), dimension(3,3) :: math_deviatoric33
  real(pReal), dimension(3,3), intent(in) :: m
 
@@ -859,7 +826,6 @@ end function math_deviatoric33
 !--------------------------------------------------------------------------------------------------
 pure function math_equivStrain33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
  real(pReal), dimension(3) :: e,s
  real(pReal) :: math_equivStrain33
@@ -880,7 +846,6 @@ end function math_equivStrain33
 !--------------------------------------------------------------------------------------------------
 pure function math_equivStress33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
  real(pReal) :: math_equivStress33
 
@@ -902,7 +867,6 @@ end function math_equivStress33
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_trace33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
 
  math_trace33 = m(1,1) + m(2,2) + m(3,3)
@@ -915,7 +879,6 @@ end function math_trace33
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_det33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
  
  math_det33 = m(1,1)* (m(2,2)*m(3,3)-m(2,3)*m(3,2)) &
@@ -930,7 +893,6 @@ end function math_det33
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_detSym33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
  
   math_detSym33 = -(m(1,1)*m(2,3)**2 + m(2,2)*m(1,3)**2 + m(3,3)*m(1,2)**2) &
@@ -944,7 +906,6 @@ end function  math_detSym33
 !--------------------------------------------------------------------------------------------------
 pure function math_33to9(m33)
 
- implicit none
  real(pReal), dimension(9)               :: math_33to9
  real(pReal), dimension(3,3), intent(in) :: m33
  
@@ -960,7 +921,6 @@ end function math_33to9
 !--------------------------------------------------------------------------------------------------
 pure function math_9to33(v9)
 
- implicit none
  real(pReal), dimension(3,3)           :: math_9to33
  real(pReal), dimension(9), intent(in) :: v9
  
@@ -979,7 +939,6 @@ end function math_9to33
 !--------------------------------------------------------------------------------------------------
 pure function math_sym33to6(m33,weighted)
 
- implicit none
  real(pReal), dimension(6)               :: math_sym33to6
  real(pReal), dimension(3,3), intent(in) :: m33
  logical,     optional,       intent(in) :: weighted
@@ -1006,7 +965,6 @@ end function math_sym33to6
 !--------------------------------------------------------------------------------------------------
 pure function math_6toSym33(v6,weighted)
 
- implicit none
  real(pReal), dimension(3,3)           :: math_6toSym33
  real(pReal), dimension(6), intent(in) :: v6
  logical,     optional,     intent(in) :: weighted
@@ -1033,7 +991,6 @@ end function math_6toSym33
 !--------------------------------------------------------------------------------------------------
 pure function math_3333to99(m3333)
 
- implicit none
  real(pReal), dimension(9,9)                 :: math_3333to99
  real(pReal), dimension(3,3,3,3), intent(in) :: m3333
 
@@ -1050,7 +1007,6 @@ end function math_3333to99
 !--------------------------------------------------------------------------------------------------
 pure function math_99to3333(m99)
 
- implicit none
  real(pReal), dimension(3,3,3,3)         :: math_99to3333
  real(pReal), dimension(9,9), intent(in) :: m99
 
@@ -1070,7 +1026,6 @@ end function math_99to3333
 !--------------------------------------------------------------------------------------------------
 pure function math_sym3333to66(m3333,weighted)
 
- implicit none
  real(pReal), dimension(6,6)                 :: math_sym3333to66
  real(pReal), dimension(3,3,3,3), intent(in) :: m3333
  logical,     optional,           intent(in) :: weighted
@@ -1098,7 +1053,6 @@ end function math_sym3333to66
 !--------------------------------------------------------------------------------------------------
 pure function math_66toSym3333(m66,weighted)
 
- implicit none
  real(pReal), dimension(3,3,3,3)            :: math_66toSym3333
  real(pReal), dimension(6,6),    intent(in) :: m66
  logical,     optional,          intent(in) :: weighted
@@ -1127,7 +1081,6 @@ end function math_66toSym3333
 !--------------------------------------------------------------------------------------------------
 pure function math_Voigt66to3333(m66)
 
- implicit none
  real(pReal), dimension(3,3,3,3) :: math_Voigt66to3333
  real(pReal), dimension(6,6), intent(in) :: m66
  integer :: i,j
@@ -1147,7 +1100,6 @@ end function math_Voigt66to3333
 !--------------------------------------------------------------------------------------------------
 pure function math_qMul(A,B)
 
- implicit none
  real(pReal), dimension(4) ::  math_qMul
  real(pReal), dimension(4), intent(in) ::  A, B
 
@@ -1164,7 +1116,6 @@ end function math_qMul
 !--------------------------------------------------------------------------------------------------
 pure function math_qConj(Q)
 
- implicit none
  real(pReal), dimension(4) ::  math_qConj
  real(pReal), dimension(4), intent(in) ::  Q
 
@@ -1178,7 +1129,6 @@ end function math_qConj
 !--------------------------------------------------------------------------------------------------
 pure function math_qRot(Q,v)
 
- implicit none
  real(pReal), dimension(4), intent(in) :: Q
  real(pReal), dimension(3), intent(in) :: v
  real(pReal), dimension(3) :: math_qRot
@@ -1209,7 +1159,6 @@ end function math_qRot
 !--------------------------------------------------------------------------------------------------
 pure function math_RtoEuler(R)
 
- implicit none
  real(pReal), dimension (3,3), intent(in) :: R
  real(pReal), dimension(3) :: math_RtoEuler
  real(pReal) :: sqhkl, squvw, sqhk
@@ -1241,7 +1190,6 @@ end function math_RtoEuler
 !--------------------------------------------------------------------------------------------------
 pure function math_RtoQ(R)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: R
  real(pReal), dimension(4)   :: absQ, math_RtoQ
  real(pReal) :: max_absQ
@@ -1298,7 +1246,6 @@ end function math_RtoQ
 !--------------------------------------------------------------------------------------------------
 pure function math_EulerToR(Euler)
 
- implicit none
  real(pReal), dimension(3), intent(in) :: Euler
  real(pReal), dimension(3,3) :: math_EulerToR
  real(pReal) :: c1, C, c2, s1, S, s2
@@ -1337,7 +1284,6 @@ end function math_EulerToR
 !--------------------------------------------------------------------------------------------------
 pure function math_axisAngleToR(axis,omega)
 
- implicit none
  real(pReal), dimension(3,3) :: math_axisAngleToR
  real(pReal), dimension(3), intent(in) :: axis
  real(pReal), intent(in) :: omega
@@ -1379,7 +1325,6 @@ pure function math_qToRodrig(Q)
  use prec, only: &
    tol_math_check
 
- implicit none
  real(pReal), dimension(4), intent(in) :: Q
  real(pReal), dimension(3) :: math_qToRodrig
 
@@ -1395,7 +1340,6 @@ real(pReal) function math_sampleGaussVar(meanvalue, stddev, width)
  use prec, only: &
    tol_math_check
 
- implicit none
  real(pReal), intent(in) ::            meanvalue, &                                                 ! meanvalue of gauss distribution
                                        stddev                                                       ! standard deviation of gauss distribution
  real(pReal), intent(in), optional ::  width                                                        ! width of considered values as multiples of standard deviation
@@ -1426,7 +1370,6 @@ end function math_sampleGaussVar
 !--------------------------------------------------------------------------------------------------
 subroutine math_eigenValuesVectorsSym(m,values,vectors,error)
 
- implicit none
  real(pReal), dimension(:,:),                  intent(in)  :: m
  real(pReal), dimension(size(m,1)),            intent(out) :: values
  real(pReal), dimension(size(m,1),size(m,1)),  intent(out) :: vectors
@@ -1453,7 +1396,6 @@ end subroutine math_eigenValuesVectorsSym
 !--------------------------------------------------------------------------------------------------
 subroutine math_eigenValuesVectorsSym33(m,values,vectors)
  
- implicit none
  real(pReal), dimension(3,3),intent(in)  :: m
  real(pReal), dimension(3),  intent(out) :: values
  real(pReal), dimension(3,3),intent(out) :: vectors
@@ -1506,7 +1448,6 @@ end subroutine math_eigenValuesVectorsSym33
 !--------------------------------------------------------------------------------------------------
 function math_eigenvectorBasisSym(m)
 
- implicit none
  real(pReal), dimension(:,:),      intent(in)  :: m
  real(pReal), dimension(size(m,1))             :: values
  real(pReal), dimension(size(m,1),size(m,1))   :: vectors
@@ -1531,7 +1472,6 @@ end function math_eigenvectorBasisSym
 !--------------------------------------------------------------------------------------------------
 pure function math_eigenvectorBasisSym33(m)
 
- implicit none
  real(pReal), dimension(3,3)              :: math_eigenvectorBasisSym33
  real(pReal), dimension(3)                :: invariants, values
  real(pReal), dimension(3,3), intent(in) :: m
@@ -1596,7 +1536,6 @@ end function math_eigenvectorBasisSym33
 !--------------------------------------------------------------------------------------------------
 pure function math_eigenvectorBasisSym33_log(m)
 
- implicit none
  real(pReal), dimension(3,3)              :: math_eigenvectorBasisSym33_log
  real(pReal), dimension(3)                :: invariants, values
  real(pReal), dimension(3,3), intent(in) :: m
@@ -1665,7 +1604,6 @@ function math_rotationalPart33(m)
  use IO, only: &
    IO_warning
 
- implicit none
  real(pReal), intent(in), dimension(3,3) :: m
  real(pReal), dimension(3,3) :: math_rotationalPart33
  real(pReal), dimension(3,3) :: U , Uinv
@@ -1691,7 +1629,6 @@ function math_eigenvaluesSym(m)
  use, intrinsic :: &
    IEEE_arithmetic
  
- implicit none
  real(pReal), dimension(:,:),                  intent(in)  :: m
  real(pReal), dimension(size(m,1))                         :: math_eigenvaluesSym
  real(pReal), dimension(size(m,1),size(m,1))               :: vectors
@@ -1716,7 +1653,6 @@ end function math_eigenvaluesSym
 !--------------------------------------------------------------------------------------------------
 function math_eigenvaluesSym33(m)
 
- implicit none
  real(pReal), intent(in), dimension(3,3) :: m
  real(pReal), dimension(3) :: math_eigenvaluesSym33,invariants
  real(pReal) :: P, Q, rho, phi
@@ -1747,7 +1683,6 @@ end function math_eigenvaluesSym33
 !--------------------------------------------------------------------------------------------------
 pure function math_invariantsSym33(m)
 
- implicit none
  real(pReal), dimension(3,3), intent(in) :: m
  real(pReal), dimension(3) :: math_invariantsSym33
 
@@ -1764,7 +1699,6 @@ end function math_invariantsSym33
 !--------------------------------------------------------------------------------------------------
 integer pure function math_factorial(n)
 
- implicit none
  integer, intent(in) :: n
  integer :: i
  
@@ -1778,7 +1712,6 @@ end function math_factorial
 !--------------------------------------------------------------------------------------------------
 integer pure function math_binomial(n,k)
 
- implicit none
  integer, intent(in) :: n, k
  integer :: i, j
  
@@ -1793,7 +1726,6 @@ end function math_binomial
 !--------------------------------------------------------------------------------------------------
 integer pure function math_multinomial(alpha)
 
- implicit none
  integer, intent(in), dimension(:) :: alpha
  integer :: i
  
@@ -1810,7 +1742,6 @@ end function math_multinomial
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_volTetrahedron(v1,v2,v3,v4)
 
- implicit none
  real(pReal), dimension (3), intent(in) :: v1,v2,v3,v4
  real(pReal), dimension (3,3) :: m
 
@@ -1828,7 +1759,6 @@ end function math_volTetrahedron
 !--------------------------------------------------------------------------------------------------
 real(pReal) pure function math_areaTriangle(v1,v2,v3)
 
- implicit none
  real(pReal), dimension (3), intent(in) :: v1,v2,v3
 
  math_areaTriangle = 0.5_pReal * norm2(math_crossproduct(v1-v2,v1-v3))
@@ -1841,7 +1771,6 @@ end function math_areaTriangle
 !--------------------------------------------------------------------------------------------------
 pure function math_rotate_forward33(tensor,rot_tensor)
 
- implicit none
 
  real(pReal), dimension(3,3) ::  math_rotate_forward33
  real(pReal), dimension(3,3), intent(in) :: tensor, rot_tensor
@@ -1856,7 +1785,6 @@ end function math_rotate_forward33
 !--------------------------------------------------------------------------------------------------
 pure function math_rotate_backward33(tensor,rot_tensor)
 
- implicit none
  real(pReal), dimension(3,3) ::  math_rotate_backward33
  real(pReal), dimension(3,3), intent(in) :: tensor, rot_tensor
 
@@ -1870,7 +1798,6 @@ end function math_rotate_backward33
 !--------------------------------------------------------------------------------------------------
 pure function math_rotate_forward3333(tensor,rot_tensor)
 
- implicit none
  real(pReal), dimension(3,3,3,3) ::  math_rotate_forward3333
  real(pReal), dimension(3,3), intent(in) :: rot_tensor
  real(pReal), dimension(3,3,3,3), intent(in) :: tensor
@@ -1895,7 +1822,6 @@ real(pReal) pure elemental function math_clip(a, left, right)
  use, intrinsic :: &
    IEEE_arithmetic
 
- implicit none
  real(pReal), intent(in) :: a
  real(pReal), intent(in), optional :: left, right