Closed kynan closed 12 years ago
The nice thing about this is that once we do Jacobian insertion in the UFL, the rest of the code basically just sees J, inv(J) and det(J) as tensor and scalar fields respectively so kernels taking them as parameters should "just work".
This has been achieved in the generate-transforms branch, merged in e739ee4.
transform_to_physical
will change to only compute the Jacobian determinant and (if necessary) inverse. That meansdetwei
and (if necessary)dShape
need to be computed by the kernel tabulating the local tensor.