00001 //============================================================================= 00002 /*! +zgbmatrix operator */ 00003 inline const zgbmatrix& operator+(const zgbmatrix& mat) 00004 {VERBOSE_REPORT; 00005 return mat; 00006 } 00007 00008 //============================================================================= 00009 /*! -zgbmatrix operator */ 00010 inline _zgbmatrix operator-(const zgbmatrix& mat) 00011 {VERBOSE_REPORT; 00012 zgbmatrix 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 }