Something like... double * doMathOperation(unsigned N, const double * in, double aux1, double aux2, double * out = 0); (obviously, in C, you would leave

Something like... double * doMathOperation(unsigned N, const double * in, double aux1, double aux2, double * out = 0); (obviously, in C, you would leave the = 0 out and document it instead).