Go to the documentation of this file.00001
00002
00003 inline zgbmatrix& zgbmatrix::operator*=(const double& d)
00004 {VERBOSE_REPORT;
00005 zdscal_((kl+ku+1)*n, d, array, 1);
00006 return *this;
00007 }
00008
00009
00010
00011 inline zgbmatrix& zgbmatrix::operator/=(const double& d)
00012 {VERBOSE_REPORT;
00013 zdscal_((kl+ku+1)*n, 1./d, array, 1);
00014 return *this;
00015 }
00016
00017
00018
00019
00020
00021
00022
00023 inline _zgbmatrix operator*(const zgbmatrix& mat, const double& d)
00024 {VERBOSE_REPORT;
00025 zgbmatrix newmat(mat.m, mat.n, mat.kl, mat.ku);
00026 for(long i=0; i<(newmat.kl+newmat.ku+1)*newmat.n; i++){
00027 newmat.array[i] =mat.array[i]*d;
00028 }
00029
00030 return _(newmat);
00031 }
00032
00033
00034
00035 inline _zgbmatrix operator/(const zgbmatrix& mat, const double& d)
00036 {VERBOSE_REPORT;
00037 zgbmatrix newmat(mat.m, mat.n, mat.kl, mat.ku);
00038 for(long i=0; i<(newmat.kl+newmat.ku+1)*newmat.n; i++){
00039 newmat.array[i] =mat.array[i]/d;
00040 }
00041
00042 return _(newmat);
00043 }