public Vector row(int r) { Vector result = new Vector(columns_); for (int i = 0; i < columns_; i++) result[i] = this[r, i]; return result; }
public virtual double value(double t) { Matrix m = param_.diffusion(t, new Vector()); double u = 0; m.row(i_).ForEach((ii, vv) => u += vv * m.row(j_)[ii]); return(u); }