When I tried to move cumsum into fortran (as opposed to using slower method using numpy outside fortran interface), I forgot to initialize the first element of the output with the input. So, the first element was a pseudo random number (whatever was in memory).
When I tried to move cumsum into fortran (as opposed to using slower method using numpy outside fortran interface), I forgot to initialize the first element of the output with the input. So, the first element was a pseudo random number (whatever was in memory).