Армадилло поэлементное умножение каждого столбца в матрице1, с каждым столбцом в матрице2, хранить в матрице3

Question

Следующий цикл генерирует желаемый результат, но я чувствую, что есть более похожий на броненосец способ получения того же результата.Возможно с .each_col (), но я не могу заставить его работать должным образом.

mat1 = nx 8, mat2 = nx 8, mat3 = nx 64, arma v6.1, c ++ 98

for (unsigned i= 0; i< 8; i++)
{
    for (unsigned j= 0; j < 8; j++)
    {
        mat3.col(i * 8 + j) = mat1.col(i) % mat2.col(j);
    }
}