私はC ++にかなり慣れていないので、これを正しく説明しているかどうかわかりません。
プログラム全体を実行しても終了しないコード (関数) があります。
この関数の目的は、最初に 0 から 1 までの乱数を生成することです。次に、その乱数が 0.6 より大きい場合、別の関数に入ります。その関数 virus_death は完全に機能します。乱数が 0.6 未満の場合、ウィルス ポイントは、reached と呼ばれるベクトルに押し戻されます。また、ポイントをファイルに書き込みます。ただし、ステートメントは停止しません。
void replication(vector<Virus>& virus, int a, int b, double& times, int& iv,
ofstream& data_file, int p, vector<Direction>& direction_v,
vector<Virus>&replicate, int step)
{
point R,P;
double replicates = ((double) rand() / (RAND_MAX));
//how long the virus takes to reach the necleus
if( replicates>0.6)
{
cout<<"No Replication "<<endl;
virus_death(times, iv, virus,a,b,data_file,p,direction_v, step);
}
else if( replicates<=0.6)
{
cout<<"Replication"<<endl;
ofstream tom;
tom.open("Timing_nuc.dat");
cout<<"Points"<<endl;
R.x = virus[p].getx();
R.y = virus[p].gety();
R.t = virus[p].gett();
cout<<"Write to File"<<endl;
tom<<p+1<<"\t"<<R.x<<"\t"<<R.y<<"\t"<<R.t<<"\n";
cout<<"New Points"<<endl;
P.x = a;
P.y = b;
P.t = 0;
cout<<P.x<<"\t"<<P.y<<"\t"<<P.t<<"\n";
replicate.push_back(P);
int aspa = replicate.size();
cout<<"Replication Size:"<<aspa<<endl;
}
}
ポイントがベクトルにプッシュされると、関数を離れない理由を誰かが判断するのを手伝ってくれませんか。
この関数を呼び出すコードは次のとおりです。
void location_of_virus(vector<Virus>&virus,double& times, int& iv,
ofstream& data_file, vector<Direction>& direction_v,
vector<Virus>&virus_engulf, circle&cell_c,
circle&nuclei_c, vector<Virus>&replicate, int step)
{
int a;
int b;
//int nc,cc,cmc,nci;
float dis_c, dis_n;
int virus_size= virus.size();
//type = "virus";
for(int i =0; i<virus_size; i++)
{
a = virus[i].getx();
b = virus[i].gety();
//to check where the virus is in location to the environment
dis_c = cell_c.distance(a,b);
dis_n = nuclei_c.distance(a,b);
if(dis_n < 5)
{
location =" Inside Nucleus";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<"\n";
replication(virus,a,b,times,iv,data_file,i,direction_v,replicate,step);
}
else if(dis_n ==5)
{
location = "Nucleus";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<"\n";
nucleus(virus,a,b,times,iv,data_file,i,direction_v,step);
}
else if(dis_c <25 && dis_n>5)
{
location = "Cytoplasm";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<"\n";
cytoplasm(virus,a,b,times,iv,data_file,p,direction_v,step);
}
else if(dis_c ==25)
{
location = "Cell Membrane";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<"\n";
endocytosis(virus,a,b,times,iv, data_file,i,direction_v,virus_engulf,cell_c,nuclei_c);
}
else if (((a ==0 || a == 50 || a>0 || a<50)&& (b ==0 || b == 50 || b>0 || b<50)) && dis_c>15)
{
location ="Extracellular Matrix";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<direction_v[i]<<"\t"<<location<<"\n";
}
else
{
location ="Left the System";
cout<<i+1<<"\t"<<virus[i]<<"\t"<<location<<endl;
data_file<<i+1<<"\t"<<virus[i]<<"\t"<<location<<"\n";
}
}
}