Method operator+=

linalg/multivector.cpp:34–44

  }
  
  MultiVector MultiVector::operator+= (const MultiVectorExpr & expr)
  {
    if (Size() != expr.Size())
      throw Exception("MultiVector assignment-add sizes mismatch, my size = "
                      + ToString(Size()) + " other size = " + ToString(expr.Size()));
    
    Vector<double> ones(Size());
    ones = 1;      
    expr.AddTo(ones, *this);
    return *this;
  }
  
  MultiVector & MultiVector::operator-= (const MultiVectorExpr & expr)
  {