Checked computation of hardening matrix

In case of constitutive_pheno.f90, the computation can be optimized.
Mattex file slightly modified

trunk/constitutive_dislo.f90

@@ -720,7 +720,7 @@ do i=1,material_maxN
    do k=1,constitutive_maxNstatevars
 !* Hardening type *
       do l=1,constitutive_maxNstatevars
-	     if (crystal_SlipIntType(j,k,l)==l) then
+	     if (crystal_SlipIntType(j,k,i)==l) then
 		    K_inter=material_SlipIntCoeff(l,i)
 		 else
 		    K_inter=0.0_pReal

trunk/constitutive_pheno.f90

@@ -691,10 +691,12 @@ do i=1,material_maxN
    do k=1,constitutive_maxNstatevars
 !* Hardening type *
       do l=1,constitutive_maxNstatevars
-	     if (crystal_SlipIntType(j,k,l)==l) then
-		    K_inter=material_SlipIntCoeff(l,i)
-		 else
+	     if (crystal_SlipIntType(j,k,i)==0) then
 		    K_inter=0.0_pReal
+		 elseif (crystal_SlipIntType(j,k,i)==1) then
+		    K_inter=material_SlipIntCoeff(1,i)
+		 else
+		    K_inter=material_SlipIntCoeff(2,i)
 	     endif
       enddo
       constitutive_HardeningMatrix(j,k,i)=K_inter

trunk/mattex.mpie

@@ -5,15 +5,19 @@ Nslip				12
 C11				106.75e3
 C12				60.41e3
 C44				28.34e3
-# Parameters for phenomenological modeling
+
+## Parameters for phenomenological modeling
+# Unit in [MPa]
 s0_slip				31.0
 gdot0_slip			0.001
 n_slip				20
 h0				75
 s_sat				63
 w0				2.25
+# Self and latent hardening coefficients
 hardening_coefficients        	1.0 1.4
-# Parameters for dislocation-based modeling
+
+## Parameters for dislocation-based modeling
 # Initial dislocation density [m]²
 rho0				              1.5e11 
 # Burgers vector [m]