11 if ((arg1->dimensionality()!=1) || arg2->dimensionality()!=1) {
13 <<
"Warning: dimension mismatch in function convolution"
20 _arg1(right._arg1->clone()),
21 _arg2(right._arg2->clone()),
36 const double NDIVISIONS=200.0;
37 double dx = (_x1-_x0)/NDIVISIONS;
39 for (
double x=_x0; x<_x1; x+=dx) {
40 result += (*_arg1)(argument-x)*(*_arg2)(x);