Question

我最近在这里询问了一个相关的问题：

https://stackoverflow.com/questions/21171836/storing-matrix-output-in-higher-dimensional-matrix

我现在只想在并行配置中运行代码，但是当我这样做时，我收到以下错误：

???使用==＆gt;时出错parallel_function at 598 ==＆gt;中的错误OIRE在136 未定义的函数或变量“best_index”。

==＆gt;中的错误OIRE_MSE_test在73 parfor t = 1：nsims

==＆gt;中的错误OIRE_MSE_test在95 [b_OIRE OIRE_opt_b（：，：，t）] = OIRE（y，x，iter）;

只要我放弃matlab池打开/关闭和parfor命令，代码就可以正常工作。为什么这不会在并行工具箱中运行？

clc;
n=100;
p=7; 
alpha1=0.999999;
error_vol=0.1;

iter=1000;
nsims=200;

OIRE_opt_b=zeros(7,3,nsims);
OIRE_opt_MSE=zeros(1,3,nsims);
OIRE_opt_index=zeros(1,3,nsims);

GIREI_opt_b=zeros(7,3,nsims);
GIREI_opt_MSE=zeros(1,3,nsims);
GIREI_opt_index=zeros(1,3,nsims);

GIREII_opt_b=zeros(7,3,nsims);
GIREII_opt_MSE=zeros(1,3,nsims);
GIREII_opt_index=zeros(1,3,nsims);

LRRE_opt_b=zeros(7,3,nsims);
LRRE_opt_MSE=zeros(1,3,nsims);
LRRE_opt_index=zeros(1,3,nsims);



 matlabpool open


x=zeros(n,p);

for i=1:n
    z_i4=normrnd(0,1);
    x(i,p)=z_i4;
    for j=1:p-1
        x(i,j)=x(i,j)+alpha1*z_i4;
        x(i,j)= x(i,j)+(1-alpha1^2)^(0.5)*normrnd(0,1);
    end  
end

b_act=[5;1;10;-20;200;30;-2];

parfor t=1:nsims

    residuals=normrnd(0,error_vol,n,1);
    y=x*b_act + residuals;
    y_store(:,t)=y;
    y=y_store(:,t);

    [b_OIRE OIRE_opt_b(:,:,t) OIRE_opt_MSE(:,:,t) OIRE_opt_index(:,:,t)]=OIRE(y,x,iter);

end

调用函数

function [b_OIRE OIRE_opt_b OIRE_opt_MSE OIRE_opt_index]=OIRE(y,x,iter)  
dim=1;
pool=[10,100,1000,10000,1000000,10000000];
count=0;

[n, p]=size(x);

b=x\y; 
b_OIRE = b; % [#1] initialize b_LRRE as b
sigma_sq=((y-x*b)'*(y-x*b))/(n-p); %'

b_act=[1;0;1;1;0;1;1];

econFlag=0;
[U,sigma,V] = svd(x,econFlag);


U1=U(:,1:p);
d=zeros(p,1);
d=diag(d);

alpha=V'*b_OIRE; %'
Delta=sigma.^1;        
Delta=diag(Delta);
f=Delta.*alpha;
F=diag(f);
Theta=sum(f);
c=p^2*sigma_sq+p*Theta^2;
g=Theta*sum(alpha);
I=ones(p,1);
a=sigma_sq*I+Theta*f;
b=F*alpha;
k=zeros(p,1);
A=sigma_sq*eye(p)+F.^2;
varRho=(g-a'*pinv(A)*b)*pinv(c-a'*pinv(A)*a); 
k=pinv(A)*b-varRho*pinv(A)*a;
K=diag(k);
D=varRho*I*I';

b_OIRE= V*(K+D)*U1'*y;


MSE=(k'*A*k)+(2*varRho*a'*k)-(2*b'*k)+(c*varRho^2)-(2*g*varRho)+(alpha'*alpha);


best_OIRE_MSE=MSE;
best_b_OIRE=b_OIRE;


for jj=1:iter           % [## "iter" denotes the iteration number]

    alpha=V'*b_OIRE;        %'
    Delta=sigma.^1;         % [Error! not sigma.^2 but sigma.^1]
    Delta=diag(Delta);
    f=Delta.*alpha;
    F=diag(f);
    Theta=sum(f);
    c=p^2*sigma_sq+p*Theta^2;
    g=Theta*sum(alpha);
    I=ones(p,1);
    a=sigma_sq*I+Theta*f;
    b=F*alpha;
    k=zeros(p,1);
    A=sigma_sq*eye(p)+F.^2;
    varRho=(g-a'*pinv(A)*b)*pinv(c-a'*pinv(A)*a);
    k=pinv(A)*b-varRho*pinv(A)*a;
    K=diag(k);
    D=varRho*I*I';

    b_OIRE= V*(K+D)*U1'*y;

    MSE=(k'*A*k)+(2*varRho*a'*k)-(2*b'*k)+(c*varRho^2)-(2*g*varRho)+(alpha'*alpha);

    if(MSE<best_OIRE_MSE)

        best_b_OIRE=b_OIRE;
        best_OIRE_MSE=MSE;
        best_index=jj+1;
    end

    if( any(jj == pool))
        count=count+1;
        OIRE_opt_b(:,count)=best_b_OIRE;
        OIRE_opt_MSE(count)=best_OIRE_MSE;
        OIRE_opt_index(count)=best_index;
    end


end


end

matlabpool close

Answer 1

您的问题是best_index永远不会在MSE<best_OIRE_MSE if语句之外的任何地方定义。

这是我怀疑正在发生的事情。当您运行代码而没有 parfor时，它将按顺序迭代（即jj==1将始终发生在jj==2之前等）。在没有查看代码的复杂性的情况下，我怀疑这意味着MSE<best_OIRE_MSE在any(jj==pool)之前（或在同一次迭代中）始终为真。这意味着，当您的代码到达any(jj == pool) if语句时，best_index将始终设置。

问题在于parfor无法保证循环的执行顺序。 jj==10可能会在jj==2之前运行。这意味着无法保证best_index在到达any(jj == pool) if语句时定义 - 因此您的错误。

我的建议是将主best_index循环之外的parfor定义为-1，然后忽略OIRE_opt_index为负数的任何情况。

程序在Matlab中正常运行，但在Parallel Toolbox中不起作用

1 个答案: