ＦＣＣ（面心立方）の位置に粒子を配置する方法

質問

main関数の中で以下のように書きました．よりよくするにはどうすればいいでしょうか．（計算速度，見た目）

i=0;
for(j=0;j<M;++j){
for(k=0;k<M;++k){
for(l=0;l<M;++l){
if(i<M*M*M){
rx[i]=0.01+0.5*L+a*(j-M/2);
ry[i]=0.01+0.5*L+a*(k-M/2);
rz[i]=0.01+0.5*L+a*(l-M/2);
}
++i;
}
}
}
for(j=0;j<M;++j){
for(k=0;k<M;++k){
for(l=0;l<M;++l){
if(i<2*M*M*M){
rx[i]=0.01+0.5*L+a*(j-M/2)+a/2;
ry[i]=0.01+0.5*L+a*(k-M/2)+a/2;
rz[i]=0.01+0.5*L+a*(l-M/2);
}
++i;
}
}
}
for(j=0;j<M;++j){
for(k=0;k<M;++k){
for(l=0;l<M;++l){
if(i<3*M*M*M){
rx[i]=0.01+0.5*L+a*(j-M/2)+a/2;
ry[i]=0.01+0.5*L+a*(k-M/2);
rz[i]=0.01+0.5*L+a*(l-M/2)+a/2;
}
++i;
}
}
}
for(j=0;j<M;++j){
for(k=0;k<M;++k){
for(l=0;l<M;++l){
if(i<4*M*M*M){
rx[i]=0.01+0.5*L+a*(j-M/2);
ry[i]=0.01+0.5*L+a*(k-M/2)+a/2;
rz[i]=0.01+0.5*L+a*(l-M/2)+a/2;
}
++i;
}
}
}

回答No.2

Tacosan · ベストアンサー

最初に 1回だけ計算して配列に入れておく, ってこと. 同じ値を何回も計算するのはバカらしいので.

お礼

なるほど．最もですね．ありがとうございます．

回答No.1

Tacosan

if の存在意義がまったくわからない.

rx とかの式の右辺を先に計算して表にすればいいのに.

お礼

回答ありがとうございます．
プログラム初心者なもので，よくわかりません．
表にすればいいというのはどういうことでしょうか．