resulting memory should be allocated by alloc_fusion
thus, the algorithm should indicate where the result is
stored by returning a pointer to the respective
double array