Changes in new version of the dataset (8th April, 2019)

.for files:
	- 3D elements have been added (see example 2b-SingleEdgeNotchSpecimen-Tensile-3D)!
	- an Even-odd iteration split is introduced to completely "freeze" the solution estimates between the staggered iterations - detailed explanation in AlgorithmExplanation.pdf
	- size of the common block significantly reduced (nsvinttot reduced from 9 to 5)
	- nodal phase-field solution rewrite in layer 3 moved before the integration point loop to save CPU time
	- fixed incorrect iteration count after unsuccessful increment
	- indentation for easier navigation through the code
	
.inp file:
	- "extrapolation = no" is used so the first solution estimate in the 1st iteration is not linearly extrapolated from the previous converged increment - very useful for multi-step and restart analyses  
	- constant user-defined tolerance parameter for convergence is substituted with the non-constant tolerance parameters, i.e., residuals are not checked against some constant user-defined tolerance as was the case before, but the largest residual has to be smaller than the time average force (time average moment for the "rotational" DOFs 4 and 5, i.e. phase-field and dummy DOFs) by some parameter. However, one can easily use the constant user-definted tolerance criterion as before - explanation of the convergence parameter used in .inp file

To be added in the next versions:
	- deformation energy decompositions
	- thread-parallel computation for increased efficiency
	- large displacement theory