diff --git a/src/homogenization_mechanical_RGC.f90 b/src/homogenization_mechanical_RGC.f90
index 71d1542c8..f0c2900e2 100644
--- a/src/homogenization_mechanical_RGC.f90
+++ b/src/homogenization_mechanical_RGC.f90
@@ -652,7 +652,7 @@ module function RGC_updateState(P,F,avgF,dt,dPdF,ce) result(doneAndHappy)
     real(pReal), dimension(6,6) :: C
 
 
-    C = phase_homogenizedC(material_phaseID(grainID,ce),material_phaseEntry(grainID,ce))
+    C = phase_homogenizedC66(material_phaseID(grainID,ce),material_phaseEntry(grainID,ce))
     equivalentMu = lattice_equivalent_mu(C,'voigt')
 
   end function equivalentMu
diff --git a/src/phase.f90 b/src/phase.f90
index a36fc26b9..9c4212c71 100644
--- a/src/phase.f90
+++ b/src/phase.f90
@@ -230,10 +230,10 @@ module phase
     end function phase_mechanical_constitutive
 
     !ToDo: Merge all the stiffness functions
-    module function phase_homogenizedC(ph,en) result(C)
+    module function phase_homogenizedC66(ph,en) result(C)
       integer, intent(in) :: ph, en
       real(pReal), dimension(6,6) :: C
-    end function phase_homogenizedC
+    end function phase_homogenizedC66
     module function phase_damage_C(C_homogenized,ph,en) result(C)
       real(pReal), dimension(3,3,3,3), intent(in)  :: C_homogenized
       integer,                         intent(in)  :: ph,en
@@ -299,7 +299,7 @@ module phase
 
   public :: &
     phase_init, &
-    phase_homogenizedC, &
+    phase_homogenizedC66, &
     phase_f_phi, &
     phase_f_T, &
     phase_K_phi, &
diff --git a/src/phase_damage.f90 b/src/phase_damage.f90
index 49b9f9528..8e4b66ae1 100644
--- a/src/phase_damage.f90
+++ b/src/phase_damage.f90
@@ -417,7 +417,7 @@ function phase_damage_deltaState(Fe, ph, en) result(broken)
    sourceType: select case (phase_damage(ph))
 
     case (DAMAGE_ISOBRITTLE_ID) sourceType
-      call isobrittle_deltaState(phase_homogenizedC(ph,en), Fe, ph,en)
+      call isobrittle_deltaState(phase_homogenizedC66(ph,en), Fe, ph,en)
       broken = any(IEEE_is_NaN(damageState(ph)%deltaState(:,en)))
       if (.not. broken) then
         myOffset = damageState(ph)%offsetDeltaState
diff --git a/src/phase_mechanical_elastic.f90 b/src/phase_mechanical_elastic.f90
index 4d7037830..c9e5b6490 100644
--- a/src/phase_mechanical_elastic.f90
+++ b/src/phase_mechanical_elastic.f90
@@ -166,7 +166,7 @@ module subroutine phase_hooke_SandItsTangents(S, dS_dFe, dS_dFi, &
     i, j
 
 
-  C = math_66toSym3333(phase_homogenizedC(ph,en))
+  C = math_66toSym3333(phase_homogenizedC66(ph,en))
   C = phase_damage_C(C,ph,en)
 
   E = 0.5_pReal*(matmul(transpose(Fe),Fe)-math_I3)                                                  !< Green-Lagrange strain in unloaded configuration
@@ -184,7 +184,7 @@ end subroutine phase_hooke_SandItsTangents
 !> @brief returns the homogenized elasticity matrix
 !> ToDo: homogenizedC66 would be more consistent
 !--------------------------------------------------------------------------------------------------
-module function phase_homogenizedC(ph,en) result(C)
+module function phase_homogenizedC66(ph,en) result(C)
 
   real(pReal), dimension(6,6) :: C
   integer,      intent(in)    :: ph, en
@@ -197,7 +197,7 @@ module function phase_homogenizedC(ph,en) result(C)
      C = elastic_C66(ph,en)
   end select plasticType
 
-end function phase_homogenizedC
+end function phase_homogenizedC66
 
 
 end submodule elastic