%--------------------------
 % @Author: Jingqiao Hu
 % @Date: 2021-04-16 21:35:35
 % @LastEditTime: 2021-04-18 15:30:58
%--------------------------
function plot_deformation_overlap(def_cell, microx, nelx, nely, lx, global_rho)
    
    [~, ~, edofMat_mi, ~] = forAssemble(microx, microx);
%     [~,~, cdofs, ~,~,~,~] = microPBC(microx, microx,lx, lx);
%     idofs = setdiff(1:2*(microx+1)^2, cdofs);

%     dx = lx / microx;
%     jac = dx^2/4;        

%     globalx = microx * nelx;
    globaly = nely * microx;
%     [~, ~, edofMat_full, ~] = forAssemble(globalx, globaly);
    
%     [ref_full, ~, ~, ~] = GenerateMesh(nelx*lx, nely*lx, globalx, globaly);
%     ref0 = reshape(ref_full', [], 1);
    
    figure;

    for i = 1:nelx
        for j = 1:nely
            umicro = reshape(def_cell{j,i}', [], 1);

            def_cx = zeros(4, microx^2);
            def_cy = zeros(4, microx^2);
            info_local = zeros(4, microx^2);

            for ii = 1:microx
                for jj = 1:microx
                    % which micro-ele in this subdomain
                    micro_ele_local = jj + (ii-1)*microx;
                    edof_mi = edofMat_mi(micro_ele_local, :);
                    ue   = umicro(edof_mi,:); % [8,9]

                    % which micro element in whole domain
                    last_col_full = (i-1)*microx + ii-1;
                    this_col_full = (j-1)*microx + jj;
                    micro_ele_full = last_col_full * globaly + this_col_full;
%                     edof_full = edofMat_full(micro_ele_full,:);

                    def_coor = ue;
                    def_cx(:, micro_ele_local) = def_coor(1:2:end);
                    def_cy(:, micro_ele_local) = def_coor(2:2:end);
                    info_local(:, micro_ele_local) = repmat(global_rho(micro_ele_full), 4, 1);
                end
            end
            h = fill(def_cx, def_cy, 1-info_local); axis equal; hold on;
            set(h, {'LineStyle'}, {'none'});
            set(h, 'facealpha', 0.6);
        end
    end
end