Performs the outer product between mat1 (2D tensor) and mat2 (2D tensor), multiply the resulting matrix by the scalar value and add the result to self (4D tensor). In other words,

self_ijkl = self_ijkl + value * mat1_ij * mat2_kl

If mat1 is a n x m matrix, mat2 a p x q matrix, self must be a n x m x p x q tensor.