10 lines
319 B
Mathematica
10 lines
319 B
Mathematica
|
function[jacobian]=Jacobian(shapeFunction,elemCoordinates)
|
|||
|
|
% Compute Jacobian and calculate global cordinates
|
|||
|
|
|
|||
|
|
jacobian=struct();
|
|||
|
|
|
|||
|
|
jacobian.matrix=elemCoordinates'*shapeFunction.dfun;
|
|||
|
|
jacobian.invmatrix=inv(jacobian.matrix);
|
|||
|
|
jacobian.globalDerivatives=(jacobian.invmatrix*shapeFunction.dfun')';%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ƫ<EFBFBD><EFBFBD>
|
|||
|
|
|
|||
|
|
return
|