switch the order of leftMMFT and rightMMFT in callbackfunc in gradientVec. Now the MatrixMatrixFunc is the correct answer