数年来未目睹之while怪现状

onethousand 2006-11-16 09:10:10

VC 6.0 做一个涉及到矩阵运算的问题。
在函数里有 zmax,colmax,minrow等几个变量，还有一个自定义对象 EMatrix.都有值的。
同一函数体内，接下来一个while循环,iter已被赋值true.
while(iter) {
....
}

由于程序运行错误，在debug模式下看，进入循环体后，上面的那几个变量的值被取到-858993460了，当然，我原来的值不是这样的。

第一次碰到这种奇怪现象。ft!

似乎是被重新定义了，因为在Locals下面还能看到另外几个同名变量，是原值。
按理，我没有重新声明，是不会有这种情况的。

...全文

182 7 打赏收藏转发到动态举报

写回复

用AI写文章

7 条回复

切换为时间正序

请发表友善的回复…

发表回复

onethousand 2006-11-16

打赏
举报

关键是while里面出了问题。
有几个变量是类的私有变量。
编译，运行都没问题

onethousand 2006-11-16

打赏
举报

SimMethod::SMsolver()
{
int *basmat=new int[rowmx]; //存放基矩阵的列数，即约束矩阵中哪几列构成基矩阵
long double *xb=new long double[rowmx]; //存放单纯形表右端项
// int xB[rowmx]={0};
// int xN[colmx-rowmx]={0};
long double *w=new long double[rowmx]; //存放单纯形乘子
long double *z=new long double[colmx]; //存放判别数
long double *y=new long double[rowmx]; //存放主列元素
int iter=0; //是否需要继续迭代标志

for(int i=0;i<rowmx;i++) {
basmat[i]=0;
xb[i]=0;
w[i]=0;
y[i]=0;
}
for(i=0;i<colmx;i++) z[i]=0;

//第一次迭代过程
//指定基矩阵构成列
f=0;
cout<<endl;
cout<<"\nPlease appoint the columns of iniational elementary matrix:";
cout<<"\nAttension:the iniational elementary matrix should be a identity matrix."<<endl;
for(i=0;i<rowmx;i++)
cin>>basmat[i];

//计算右端项和单纯形乘子
for(i=0;i<rowmx;i++) {
xb[i]=b[i];
w[i]=c[basmat[i]-1];
}

for(i=0;i<rowmx;i++) {
f+=c[basmat[i]]*xb[i];
}

//计算判别数，并对迭代标志进行置位或复位
long double zmax=-1000000000000;
int colmax=0;
for(i=0;i<colmx;i++) {
int sum=0;
for(int j=0;j<rowmx;j++) sum+=w[j]*ctmatrp[j][i];
z[i]=sum-c[i];
if(z[i]>zmax) { zmax=z[i]; colmax=i;}
if(z[i]>0) iter=1;
}

//计算主列，并选择离基变量
for(i=0;i<rowmx;i++) y[i]=ctmatrp[i][colmax];
long double ymin=1000000000;
int minrow=0;
for(i=0;i<rowmx;i++) {
if(y[i]>0&&xb[i]/y[i]<ymin) {ymin=xb[i]/y[i];minrow=i;}
}
basmat[minrow]=colmax;//修正基矩阵

//构造初始矩阵
EMatrix E(rowmx,minrow,y);

//迭代过程
while(iter) {
f=0;
iter=0;

E*xb;
for(int i=0;i<rowmx;i++) w[i]=c[basmat[i]];
w*E;

for(i=0;i<rowmx;i++) {
f+=c[basmat[i]]*xb[i];
}

long double zmax=-1000000000000;
int colmax=0;
for(i=0;i<colmx;i++) {
int sum=0;
for(int j=0;j<rowmx;j++) sum+=w[j]*ctmatrp[j][i];
z[i]=sum-c[i];
if(z[i]>zmax) { zmax=z[i]; colmax=i;}
if(z[i]>0) iter=1;
}
if(!iter) break;

for(i=0;i<rowmx;i++) y[i]=ctmatrp[i][colmax];
long double ymin=0;
int minrow=0;
for(i=0;i<rowmx;i++) {
if(y[i]>0&&xb[i]/y[i]>ymin) {ymin=xb[i]/y[i];minrow=i;}
}
basmat[minrow]=colmax;

EMatrix E(rowmx,minrow,y);
}

cout<<endl;
cout<<"\nThe solver is:\n";
for(i=0;i<rowmx;i++)
cout<<"X"<<basmat[i]<<"="<<xb[i]<<" ";
cout<<"\nThe optimize value is "<<f;

delete[] basmat;
delete[] xb;
delete[] w;
delete[] z;
delete[] y;
}