This code is written to efficiently compute Kij = Dkl * Iika * Ijla where D is a two dimensional matrix and I is a three dimensional matrix. It uses parallelism with MPI and multithreading to achieve this purpose.