void OperationMultipleEvalPeriodic::mult(DataVector& alpha,
    DataVector& result) {
  AlgorithmDGEMV<SLinearPeriodicBasis> op;
  LinearPeriodicBasis<unsigned int, unsigned int> base;

  op.mult(storage, base, alpha, this->dataset, result);
}
void OperationBLinearBoundary::mult(DataVector& alpha, DataMatrix& data, DataVector& result)
{
	AlgorithmDGEMV<SLinearBoundaryBase> op;
	linearboundaryBase<unsigned int, unsigned int> base;

	op.mult(storage, base, alpha, data, result);
}
void OperationMultipleEvalPolyBoundary::mult(DataVector& alpha,
    DataVector& result) {
  AlgorithmDGEMV<SPolyBoundaryBase> op;

  op.mult(storage, base, alpha, this->dataset, result);
}