CurrentDensity fast option is not working for non-orthogonal cells

The same problem as before, the derivatives are not converted in the proper coordinate space before been used.