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