heaviside.m

efg code with matlab

% *************************************************************************
%                 TWO DIMENSIONAL ELEMENT FREE GALERKIN CODE
%                            Nguyen Vinh Phu
%                        LTDS, ENISE, Juillet 2006
% *************************************************************************
function [f,dfdx,dfdy] = heaviside(distance)
% Compute the Heaviside function
% Input : distance is the signed distance to the crack line

f = sign(distance) ;

dfdx = 0 ;
dfdy = 0 ;