Inverse of the 4th rank tensor
I am looking for an algorithm to get the inverse of a 4th rank tensor (e.g. the compliance tensor S_(ijkl) from elastic stiffness tensor C_(ijkl)) S_(ijkl)=C_(ijkl)^(-1)
I am programming in FORTRAN, and for this purpose I wasn't able to find neither any algorithm nor any existing subroutine.
If anyone at this forum has any idea about this inversion, kindly guide me.