So far these changes are included:

Added

• Add a warning in the docs that C**(1/2) equals to C**0. The square root operator must be used instead.

Changed

• Evaluate the i-th power of the inverse of a tensor for negative exponents in C**-3.
• Initialize all tensor components to zero.

Fixed

• Fix assignments of scalars to tensors.
• Ensure major- and minor-symmetric result in cdya(): C(i,j,k,l) = (A(i,k) B(j,l) + A(i,l) B(j,k) + B(i,k) A(j,l) + B(i,l) A(j,k))/4.

Removed

• Remove ln, exp, dexp functions which were based on isotropic tensor function-approximations. This could be misleading if one assumes these are analytic (exact) functions.
• Remove unused internal variables and comments.
• Remove unused archive_functions.md from the docs.
• Remove incorrect implementation of libinnercrossdyadic.f.

If nothing very important pops up in the next two weeks, this will be released as v2.0.0