On: Let the callers of your API control allocations

[cozzyd]

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).