00001
00003 inline _zgbmatrix operator*(const _zgbmatrix& mat, const double& d)
00004 {
00005 #ifdef CPPL_VERBOSE
00006 std::cerr << "# [MARK] operator*(const _zgbmatrix&, const double&)"
00007 << std::endl;
00008 #endif//CPPL_VERBOSE
00009
00010 zdscal_((mat.KL+mat.KU+1)*mat.N, d, mat.Array, 1);
00011 return mat;
00012 }
00013
00014
00016 inline _zgbmatrix operator/(const _zgbmatrix& mat, const double& d)
00017 {
00018 #ifdef CPPL_VERBOSE
00019 std::cerr << "# [MARK] operator/(const _zgbmatrix&, const double&)"
00020 << std::endl;
00021 #endif//CPPL_VERBOSE
00022
00023 zdscal_((mat.KL+mat.KU+1)*mat.N, 1./d, mat.Array, 1);
00024 return mat;
00025 }