I'm working in a solver that should use coordinates generalization (field-aligned or body fitted coordinates0 with OpenFOAM. But I have a question: Does OF already uses a internal generalization procedure to calculate the Jacobian, and all I need is only work on the Mesh (eg: converting that one from cartesian to my generalizes system)? Or I'll need to work on both ways: In the equation code AND in the mesh?