Added test function for generating harmonic test data
from eigenvectors and eigenvalues obtained from
Phonopy.

Improved c derivative function in order to make it
completely n dimensional compatible