Fix memory leak in cyclic parderiv solver

Allocating 2D arrays with cmatrix, but freeing with delete[] (!)
Now changed to free_cmatrix