display_material_mesh_gui.m

this a 3d fdtd code using matlab and it can be used as a reference for other structures. Dear sir,

            [handles] = plot_sigma_m_xy(handles);
        end
        if get(cToggleYZplane,'value')==1 
            [handles] = plot_sigma_m_yz(handles);
        end
        if get(cToggleZXplane,'value')==1 
            [handles] = plot_sigma_m_zx(handles);
        end
    end

    if (get(cCheckEdge,'value') == 0)
        po = findobj(gca,'type','patch');
        set(po,'EdgeColor','none');
    end
    
    if (get(cCheckColorbar,'value') == 1)
        colorbar;
    else
        colorbar off;
    end    
    
    if (get(cCheckAxis,'value') == 1)
        axis on;
    else
        axis off;
    end    
end

function [handles] = plot_eps_r_xy(handles)

    val = get_index(get(cEditK,'string'));
    colormap(jet);
    caxis([handles.eps_r_minval*0.9 handles.eps_r_maxval*1.1]);

    % plot x components;
    pdata = handles.eps_r_x(:,:,val);
    zz = handles.min_z+(val-1)*handles.dz;
    x = handles.vertices_xy_x_x;
    y = handles.vertices_xy_x_y;
    z = y*0+zz;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_xy_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot y components;
    pdata = handles.eps_r_y(:,:,val);
    zz = handles.min_z+(val-1)*handles.dz;
    x = handles.vertices_xy_y_x;
    y = handles.vertices_xy_y_y;
    z = x*0+zz;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_xy_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);

end

function [handles] = plot_eps_r_yz(handles)
    
    val = get_index(get(cEditI,'string'));
    colormap(jet);
    caxis([handles.eps_r_minval*0.9 handles.eps_r_maxval*1.1]);

    % plot y components;
    pdata = squeeze(handles.eps_r_y(val,:,:));
    xx = handles.min_x+(val-1)*handles.dx;
    y = handles.vertices_yz_y_y;
    z = handles.vertices_yz_y_z;
    x = z*0+xx;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_yz_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot z components;
    pdata = squeeze(handles.eps_r_z(val,:,:));
    xx = handles.min_x+(val-1)*handles.dx;
    y = handles.vertices_yz_z_y;
    z = handles.vertices_yz_z_z;
    x = y*0+xx;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_yz_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

end

function [handles] = plot_eps_r_zx(handles)

    val = get_index(get(cEditJ,'string'));
    colormap(jet);
    caxis([handles.eps_r_minval*0.9 handles.eps_r_maxval*1.1]);

    % plot z components;
    pdata = squeeze(handles.eps_r_z(:,val,:));
    yy = handles.min_y+(val-1)*handles.dy;
    z = handles.vertices_zx_z_z;
    x = handles.vertices_zx_z_x;
    y = x*0+yy;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_zx_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot x components;
    pdata = squeeze(handles.eps_r_x(:,val,:));
    yy = handles.min_y+(val-1)*handles.dy;
    z = handles.vertices_zx_x_z;
    x = handles.vertices_zx_x_x;
    y = z*0+yy;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_zx_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_sigma_e_xy(handles)

    val = get_index(get(cEditK,'string'));
    colormap(jet);
    caxis([handles.sigma_e_minval*0.9 handles.sigma_e_maxval*1.1]);

    % plot x components;
    pdata = handles.sigma_e_x(:,:,val);
    zz = handles.min_z+(val-1)*handles.dz;
    x = handles.vertices_xy_x_x;
    y = handles.vertices_xy_x_y;
    z = y*0+zz;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_xy_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot y components;
    pdata = handles.sigma_e_y(:,:,val);
    zz = handles.min_z+(val-1)*handles.dz;
    x = handles.vertices_xy_y_x;
    y = handles.vertices_xy_y_y;
    z = x*0+zz;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_xy_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_sigma_e_yz(handles)
    
    val = get_index(get(cEditI,'string'));
    colormap(jet);
    caxis([handles.sigma_e_minval*0.9 handles.sigma_e_maxval*1.1]);

    % plot y components;
    pdata = squeeze(handles.sigma_e_y(val,:,:));
    xx = handles.min_x+(val-1)*handles.dx;
    y = handles.vertices_yz_y_y;
    z = handles.vertices_yz_y_z;
    x = z*0+xx;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_yz_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot z components;
    pdata = squeeze(handles.sigma_e_z(val,:,:));
    xx = handles.min_x+(val-1)*handles.dx;
    y = handles.vertices_yz_z_y;
    z = handles.vertices_yz_z_z;
    x = y*0+xx;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_yz_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

end

function [handles] = plot_sigma_e_zx(handles)

    val = get_index(get(cEditJ,'string'));
    colormap(jet);
    caxis([handles.sigma_e_minval*0.9 handles.sigma_e_maxval*1.1]);

    % plot z components;
    pdata = squeeze(handles.sigma_e_z(:,val,:));
    yy = handles.min_y+(val-1)*handles.dy;
    z = handles.vertices_zx_z_z;
    x = handles.vertices_zx_z_x;
    y = x*0+yy;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.faces_zx_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot x components;
    pdata = squeeze(handles.sigma_e_x(:,val,:));
    yy = handles.min_y+(val-1)*handles.dy;
    z = handles.vertices_zx_x_z;
    x = handles.vertices_zx_x_x;
    y = z*0+yy;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.faces_zx_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);

end

function [handles] = plot_mu_r_xy(handles)

    val = get_index(get(cEditK,'string'));
    colormap(jet);
    caxis([handles.mu_r_minval*0.9 handles.mu_r_maxval*1.1]);

    % plot x components;
    pdata = handles.mu_r_x(:,:,val);
    zz = handles.min_z+(val-0.5)*handles.dz;
    x = handles.mu_vertices_xy_x_x;
    y = handles.mu_vertices_xy_x_y;
    z = y*0+zz;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_xy_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot y components;
    pdata = handles.mu_r_y(:,:,val);
    zz = handles.min_z+(val-0.5)*handles.dz;
    x = handles.mu_vertices_xy_y_x;
    y = handles.mu_vertices_xy_y_y;
    z = x*0+zz;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_xy_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_mu_r_yz(handles)
    val = get_index(get(cEditI,'string'));
    colormap(jet);
    caxis([handles.mu_r_minval*0.9 handles.mu_r_maxval*1.1]);

    % plot y components;
    pdata = squeeze(handles.mu_r_y(val,:,:));
    xx = handles.min_x+(val-0.5)*handles.dx;
    y = handles.mu_vertices_yz_y_y;
    z = handles.mu_vertices_yz_y_z;
    x = z*0+xx;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_yz_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot z components;
    pdata = squeeze(handles.mu_r_z(val,:,:));
    xx = handles.min_x+(val-0.5)*handles.dx;
    y = handles.mu_vertices_yz_z_y;
    z = handles.mu_vertices_yz_z_z;
    x = y*0+xx;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_yz_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_mu_r_zx(handles)

    val = get_index(get(cEditJ,'string'));
    colormap(jet);
    caxis([handles.mu_r_minval*0.9 handles.mu_r_maxval*1.1]);

    % plot z components;
    pdata = squeeze(handles.mu_r_z(:,val,:));
    yy = handles.min_y+(val-0.5)*handles.dy;
    z = handles.mu_vertices_zx_z_z;
    x = handles.mu_vertices_zx_z_x;
    y = x*0+yy;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_zx_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot x components;
    pdata = squeeze(handles.mu_r_x(:,val,:));
    yy = handles.min_y+(val-0.5)*handles.dy;
    z = handles.mu_vertices_zx_x_z;
    x = handles.mu_vertices_zx_x_x;
    y = z*0+yy;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_zx_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_sigma_m_xy(handles)

    val = get_index(get(cEditK,'string'));
    colormap(jet);
    caxis([handles.sigma_m_minval*0.9 handles.sigma_m_maxval*1.1]);

    % plot x components;
    pdata = handles.sigma_m_x(:,:,val);
    zz = handles.min_z+(val-0.5)*handles.dz;
    x = handles.mu_vertices_xy_x_x;
    y = handles.mu_vertices_xy_x_y;
    z = y*0+zz;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_xy_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot y components;
    pdata = handles.sigma_m_y(:,:,val);
    zz = handles.min_z+(val-0.5)*handles.dz;
    x = handles.mu_vertices_xy_y_x;
    y = handles.mu_vertices_xy_y_y;
    z = x*0+zz;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_xy_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_sigma_m_yz(handles)
    
    val = get_index(get(cEditI,'string'));
    colormap(jet);
    caxis([handles.sigma_m_minval*0.9 handles.sigma_m_maxval*1.1]);

    % plot y components;
    pdata = squeeze(handles.sigma_m_y(val,:,:));
    xx = handles.min_x+(val-0.5)*handles.dx;
    y = handles.mu_vertices_yz_y_y;
    z = handles.mu_vertices_yz_y_z;
    x = z*0+xx;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_yz_y,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot z components;
    pdata = squeeze(handles.sigma_m_z(val,:,:));
    xx = handles.min_x+(val-0.5)*handles.dx;
    y = handles.mu_vertices_yz_z_y;
    z = handles.mu_vertices_yz_z_z;
    x = y*0+xx;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_yz_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [handles] = plot_sigma_m_zx(handles)

    val = get_index(get(cEditJ,'string'));
    colormap(jet);
    caxis([handles.sigma_m_minval*0.9 handles.sigma_m_maxval*1.1]);

    % plot z components;
    pdata = squeeze(handles.sigma_m_z(:,val,:));
    yy = handles.min_y+(val-0.5)*handles.dy;
    z = handles.mu_vertices_zx_z_z;
    x = handles.mu_vertices_zx_z_x;
    y = x*0+yy;
    pdata = reshape(pdata,1,[]);
    patch('faces',handles.mu_faces_zx_z,'vertices',[x y z],'facecolor','flat','cdata',pdata);

    % plot x components;
    pdata = squeeze(handles.sigma_m_x(:,val,:));
    yy = handles.min_y+(val-0.5)*handles.dy;
    z = handles.mu_vertices_zx_x_z;
    x = handles.mu_vertices_zx_x_x;
    y = z*0+yy;
    pdata = reshape(pdata.',1,[]);
    patch('faces',handles.mu_faces_zx_x,'vertices',[x y z],'facecolor','flat','cdata',pdata);
end

function [val] = get_index(str)
    [token, remain] = strtok(str,'=');
    remain = strrep(remain, '=', ' ');
    [token, remain] = strtok(remain,'/');
    val = str2num(token);
end
function [] = set_maxval(handles)
    
    val = get_index(get(cEditJ,'string'));
    maxval = handles.ny;
    if get(cRadioPermittivity,'value')==1 || get(cRadioEconductivity,'value')==1
        maxval = maxval + 1;
    end
    if val>maxval
        val=maxval;
    end
    set(cEditJ,'string',['j = ' num2str(val) ' / ' num2str(maxval)]);

    val = get_index(get(cEditI,'string'));
    maxval = handles.nx;
    if get(cRadioPermittivity,'value')==1 || get(cRadioEconductivity,'value')==1
        maxval = maxval + 1;
    end
    if val>maxval
        val=maxval;
    end
    set(cEditI,'string',['i = ' num2str(val) ' / ' num2str(maxval)]);

    val = get_index(get(cEditK,'string'));
    maxval = handles.nz;
    if get(cRadioPermittivity,'value')==1 || get(cRadioEconductivity,'value')==1
        maxval = maxval + 1;
    end
    if val>maxval
        val=maxval;
    end
    set(cEditK,'string',['k = ' num2str(val) ' / ' num2str(maxval)]);
end

end