Избыточное копирование объекта при реализации оператора умножения и оператора присваивания

#include <cstdio>
 
class Matrix
{
  private:
    // внутренние данные матрицы
 
  public:
    Matrix ();
    Matrix (const Matrix&);
    ~Matrix ();
 
    Matrix& operator=(const Matrix& a);
    Matrix operator*(const Matrix& a) const;
};
 
Matrix::Matrix ()
{
  printf ("Matrix ()\n");
};
 
Matrix::Matrix (const Matrix&)
{
  printf ("Matrix (const Matrix&)\n");
};
 
Matrix::~Matrix ()
{
  printf ("~Matrix ()\n");
};
 
Matrix& Matrix::operator=(const Matrix& a)
{
  printf ("operator=(const Matrix&)\n");
 
  // тут как бы код, копирущий данные из "a" в "this"
 
  return *this;
}
 
Matrix Matrix::operator*(const Matrix& a) const
{
  printf ("operator*(const Matrix&)\n");
 
  Matrix result;
 
  // тут как бы код, умножающий "a" и "this" и записывающий
  // данные в "result"
 
  return result;
}
 
Matrix b, c;
 
int main (void)
{
  printf ("------------------\n");
  Matrix a = b * c;
  printf ("------------------\n");
}

    a = b;
    a *= c;

Matrix Matrix::operator*(const Matrix& a) const
{
  printf ("operator*(const Matrix&)\n");
 
  Matrix result; //<--- создается объект
 
  // тут как бы код, умножающий "a" и "this" и записывающий
  // данные в "result"
 
  return result;
}

some obj = foo(); //по факту объект построился не "внутри функции", 
   //а сразу же снаружи,
   // что позволяет избежать копирования
// таким образом получаем эффект "перемещения"

#include "math/Math.h"
#include <iostream>
 
ru::math::Matrix4 mat1;
ru::math::Matrix4 mat2;
 
ru::math::Matrix4 getMat()
{
    return mat1 * mat2;
}
 
int main()
{
    ru::math::Matrix4 mat = getMat();
    ru::math::Vector3 srcVec { 10, 20, 30 };
    ru::math::Vector3 resVec = mat * srcVec;
    std::cout << resVec.length() << std::endl;
}

ru::math::Matrix4 mat;
ru::math::Vector3 resVec;
 
int main(void)
{
    mat = getMat();
    ru::math::Vector3 srcVec { 10, 20, 30 };
    resVec = mat * srcVec;
    std::cout << resVec.length() << std::endl;
}

class Matrix;
struct pair_mul
{
    pair_mul(const Matrix &a_, const Matrix &b_)
        : a(a_)
        , b(b_)
    {}
    const Matrix &a;
    const Matrix &b;
};
 
//...
    Matrix& Matrix::operator=(const pair_mul& p)
    {
      printf ("operator=(const Pair_mul&)\n");
     
      //умножаем p.a на p.b, результат пихаем в this
      return *this;
    }
//...
    Matrix a, b, c;
    a = pair_mul(b, c);

a = b * c;

(a = b) *= c;