home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #30 / NN_1992_30.iso / spool / comp / lang / cplus / 18155 < prev    next >
Encoding:
Text File  |  1992-12-18  |  2.6 KB  |  85 lines
  1. Newsgroups: comp.lang.c++
  2. Path: sparky!uunet!usc!zaphod.mps.ohio-state.edu!caen!uvaarpa!murdoch!virginia.edu!gs4t
  3. From: gs4t@virginia.edu (Gnanasekaran Swaminathan)
  4. Subject: Re: Matrix Addition  Was Re: Complex Addition.
  5. Message-ID: <1992Dec17.204049.20009@murdoch.acc.Virginia.EDU>
  6. Keywords: operator, memory leak
  7. Sender: usenet@murdoch.acc.Virginia.EDU
  8. Organization: University of Virginia
  9. References: <1992Dec11.170849.18525@rs6000.bham.ac.uk> <1992Dec11.232954.7071@borland.com> <36223@sophia.inria.fr>
  10. Date: Thu, 17 Dec 1992 20:40:49 GMT
  11. Lines: 72
  12.  
  13. gchow@hawai.inria.fr (Gloria Chow) writes:
  14. : MATRIX& MATRIX::operator+( const MATRIX& addend )
  15. : {
  16. : register int    i,j;
  17. : MATRIX          sum(nrows, ncols);
  18. : 
  19. :     for( i=0 ; i<nrows ; i++ )
  20. :         for ( j=0 ; j<ncols ; j++ )
  21. :             sum.entry[i][j] = entry[i][j] + addend.entry[i][j];
  22. : 
  23. :     return sum;
  24. : }
  25. : 
  26. : 
  27. : Okay, the problem is when I tried to use the "+" operator as follows:
  28. : 
  29. :     matrix_c = matrix_a + matrix_b;
  30. : 
  31. : the "+" operation works and sum contains the proper matrix, that is until
  32. : "+" finishes and goes out of scope.  The "return" statement makes a temporary
  33. : copy of sum (ie. the values of its nrows, ncols, and the pointer entry), and
  34. : then proceed to destroy sum automatically.  So, now the entry pointer of this
  35. : temporary is no longer pointing to a valid block of memory!!   
  36.  
  37. You are returning a reference to a local variable sum.
  38. That is, no temporary copy is made in your case.
  39. So, when the function operator+ goes out of scope, sum
  40. gets destroyed and hence, the reference refers to junk.
  41.  
  42. One technique advocated several times in this forum and
  43. by several good books is to use reference counting and 
  44. return by value in case of big objects like Matrix.
  45.  
  46. For example,
  47.  
  48. class Matrix {
  49.     struct MatRep {
  50.         int    ref;
  51.         int      rows, cols;
  52.         double*    d;
  53.     
  54.         MatRep (int r, int c)
  55.         : ref(1), rows(r), cols(c), d(new double[rows*cols]) {}
  56.         ~MatRep () { delete [] d; }
  57.     }* rep;
  58. public:
  59.     Matrix (int r, int c): rep(new MatRep(r,c)) {}
  60.     Matrix (const Matrix& m): rep (m.rep) { ++rep->ref; }
  61.     ~Matrix () { if (--rep->ref == 0) delete rep; }
  62.  
  63.     double    operator () (int i, int j) { return rep->d[i*rep->rows+j]; }
  64.     Matrix& operator =  (const Matrix& m) {
  65.         if (this != &m) {
  66.             this->Matrix::~Matrix (); // self destruct
  67.             rep = m.rep; ++rep->ref;
  68.         }
  69.         return *this;
  70.     }
  71.  
  72.     friend Matrix  operator + (const Matrix& a, const Matrix& b);
  73. };
  74.  
  75. // note the return by value here
  76. Matrix operator + (const Matrix& a, const Matrix& b)
  77. {
  78.     Matrix sum (a.rep->rows, a.rep->cols);
  79.     for (int i=0; i < sum.rep->rows * sum.rep->cols;  i++)
  80.         sum.rep->d[i] = a.rep->d[i] + b.rep->d[i];
  81.     return sum;
  82. }
  83.     
  84.     
  85.