Понимание возврата функции

Я начинающий программист и только кратко рассмотрел анатомию вызова функции (настройка стека и т. Д.). Я могу написать функцию двумя разными способами, и мне интересно, какой (если любой) более эффективен. Это для программы с конечными элементами, поэтому эту функцию можно вызывать несколько тысяч раз. Он использует библиотеку линейной алгебры Aramdillo.

Первый способ:

void Q4::stiffness(mat &stiff)
{
stiff.zeros; // sets all elements of the matrix to zero
// a bunch of linear algebra calculations
// ...
stiff *= h;
}

int main()
{
mat elementStiffness(Q4__DOF, Q4__DOF);
mat globalStiffness(totalDOF, totalDOF);
for (int i = 0; i < reallyHugeNumber; i++)
{
elements[i].stiffness(&elementStiffness, PSTRESS);
assemble(&globalStiffness, &elementStiffness);
}
return 0;
}

Второй способ:

mat Q4::stiffness()
{
mat stiff(Q4__DOF, Q4__DOF); // initializes element stiffness matrix
// a bunch of linear algebra calculations
// ...
return stiff *= h;
}

int main()
{
mat elementStiffness(Q4__DOF, Q4__DOF);
mat globalStiffness(totalDOF, totalDOF);
for (int i = 0; i < reallyHugeNumber; i++)
{
elementStiffness = elements[i].stiffness(PSTRESS);
assemble(&globalStiffness, &elementStiffness);
}
return 0;
}

Я думаю, что я спрашиваю: использование второго способа mat stiff толкнул в стек, а затем скопировал в elementStiffness? Потому что я представляю, что матрица, помещаемая в стек, а затем копируемая, намного дороже, чем передача эталонной матрицы и установка ее элементов в ноль.