在C++中如何实现一维变带宽存储的矩阵与向量相乘

矩阵K存储在一维变带宽数组中，知道主对角元位置数组kd[].如何实现K与一位数组y[n]相乘。