pbar_jacobian.m

matlab实例

function [dx, x_pt] = pBar_jacobian(ksi, p, ElemData, Nodes)
% Jacobian (dx) and global coordinates of this point (x_pt)
% Always linear mapping

XL = Nodes( ElemData.connect(1) ).coord(1);
XR = Nodes( ElemData.connect(p+1) ).coord(1);

dx = (XR - XL)/2.0;
if dx < 1.0e-14
    error('pBar_jacobian: found element with invalid jacobian');
end

% Linear mapping always!
x_pt = XL*(1.0-ksi)*0.5 + XR*(1.0+ksi)*0.5;

% ****************************************************
% ****************************************************