43 #ifndef CORE_DATATYPES_DENSECOLMAJMATRIXMULTIPLICATION_H
44 #define CORE_DATATYPES_DENSECOLMAJMATRIXMULTIPLICATION_H
57 if (beg == -1) beg = 0;
58 if (end == -1) end = this->nrows_;
62 for (i=beg; i<end; i++)
65 for (j=0; j<this->ncols_; j++)
67 sum += iget(i, j) * x[j];
74 for (i=beg; i<end; i++) b[i] = 0.0;
75 for (j=0; j<this->ncols_; j++)
79 for (i=beg; i<end; i++)
81 b[i] += iget(i, j) * x[j];
97 if (beg == -1) beg = 0;
98 if (end == -1) end = this->ncols_;
102 for (i=beg; i<end; i++)
105 for (j=0; j<this->nrows_; j++)
107 sum += iget(j, i) * x[j];
114 for (i=beg; i<end; i++) b[i] = 0.0;
115 for (j=0; j<this->nrows_; j++)
119 for (i=beg; i<end; i++)
121 b[i] += iget(j, i) * x[j];
#define ASSERTEQ(c1, c2)
Definition: Assert.h:98
#define ASSERT(condition)
Definition: Assert.h:110
Definition: ColumnMatrix.h:55
long long index_type
Definition: Types.h:39
virtual void mult(const ColumnMatrix &x, ColumnMatrix &b, index_type beg=-1, index_type end=-1, int spVec=0) const
Definition: DenseColMajMatrixMultiplication.h:50
virtual void mult_transpose(const ColumnMatrix &x, ColumnMatrix &b, index_type beg=-1, index_type end=-1, int spVec=0) const
Definition: DenseColMajMatrixMultiplication.h:90
size_type nrows() const
Definition: Matrix.h:132