//=============================================================================
/*! double*zhsmatrix operator */
inline _zhsmatrix operator*(const double& d, const zhsmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] operator*(const double&, const zhsmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
zhsmatrix newmat(mat.N, mat.CAP);
for(long c=0; c<mat.VOL; c++){
newmat.fput(mat.Indx[c], mat.Jndx[c], d*mat.Array[c]);
}
return _(newmat);
}