%--------------------------
 % @Author: Jingqiao Hu
 % @Date: 2021-01-05 20:04:38
 % @LastEditTime: 2021-01-05 20:09:51

 % if contain fixed dofs = 1, else 0
%--------------------------
function fixed = fixed_ele(fixeddofs, edofMat)
    eleNum = size(edofMat, 1);
    fixed  = zeros(eleNum, 1);
    
    for i = 1:eleNum
        edof = edofMat(i, :);
        if sum(intersect(edof, fixeddofs)) > 0
            fixed(i) = 1;
        end
    end
end