[未答复] 循环次数比较大,如何实现并行运算?

[复制链接]
天天想放假123 发表于 5 天前
load mprofile.mat
load ec_gene_names.mat
mprofile_size=size(mprofile);
cell_num=mprofile_size(2);
node_HP=zeros(mprofile_size(2),mprofile_size(1));
gene_names=ec_gene_names;
for c=107:113
    clear adjacent_network;
    fid=fopen(['tmp_network',num2str(c),'_adj_edges.txt']);
    adjacent_network={};
    j=0;
    while ~feof(fid)
        tline=fgetl(fid);
        j=j+1;
        adjacent_network{j}=regexp(tline, '\t', 'split');
    end
    fclose(fid);
    total_node_num=j;

    for na=1:total_node_num
        center=na;
        cell_gene_name{c}{na}=gene_names{center};
        gx=mprofile(center,:);
        e=0;
        rcc=0;
        clear p_rho;
        if (length(adjacent_network{na})==1)&&(str2num(adjacent_network{na}{1})==na)
            node_HP(c,na)=0;
        else
            for n=1:length(adjacent_network{na})
                nei=adjacent_network{na}{n};
                e=e+1;
                gy=mprofile(str2num(nei),:);
                curr_weight= abs(edge_weight(gx,gy));
                node_p_rho(n)=abs(curr_weight(c)*mprofile(str2num(nei),c));
            end            
        end
    end
end

您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表