00001 //============================================================================= 00002 /*! +dgbmatrix operator */ 00003 inline const dgbmatrix& operator+(const dgbmatrix& mat) 00004 {VERBOSE_REPORT; 00005 return mat; 00006 } 00007 00008 //============================================================================= 00009 /*! -dgbmatrix operator */ 00010 inline _dgbmatrix operator-(const dgbmatrix& mat) 00011 {VERBOSE_REPORT; 00012 dgbmatrix newmat(mat.m,mat.n,mat.kl,mat.ku); 00013 for(long i=0; i<(newmat.kl+newmat.ku+1)*newmat.n; i++){ 00014 newmat.array[i]=-mat.array[i]; 00015 } 00016 00017 return _(newmat); 00018 }