• 締切済み

ほぼ同じプログラムで結果が違う問題

akayoroshiの回答

回答No.1

最適化をしないあるいは最低限にしてコンパイルしたらどうなりますか こういうとき、原因として最適化機能に不具合があることが多いので

この回答がついた質問に戻る
  1. カテゴリ
  2. [技術者向] コンピューター
  3. プログラミング・開発
  4. C・C++・C#

関連するQ&A

  • fortran 配列受け渡し時の次元の一致

    fortran90、コンパイラはifortです。 普通メインプログラムとサブルーチン間での配列の受け渡しは、次元を揃えて渡すと思います。 とあるコード（以後コードA）を読んでいると、2次元配列を渡し、1次元配列で受け取っていました。 例 program test1 integer :: a(3,3) call testsub(a) end program subroutine testsub(b) integer :: b(9) end subroutine これが受け取り側でどう処理されているのかわからず、調べるために適当なテストコードを書きました。 a 123 456 789 ↓ b 123456789 になるとか 結果、コンパイルは通ったのですがサブルーチン側では全て0で置き換えられてしまいました。 コードAはpgiかなんかでコンパイルしていたようなのでコンパイラの違いでしょうか？ よくわらなかったので質問させて頂きました。 質問をまとめますと、 (1)次元の異なる配列の受け渡しができるかどうか (2)その場合中身はどうなるか よろしくお願いします。 ---以下テストコード--- program testa implicit none integer :: a(3,3),i,j do i=1,3 do j=1,3 a(i,j)=j+(i-1)*3 enddo enddo do i=1,3 do j=1,3 write(6,*) a(i,j) enddo enddo call sub1(a) end program subroutine sub1(b) integer :: b(9),i do j=1,9 write(6,*) b(i),'sub' enddo end subroutine

  • 院試の問題(FORTRAN)で。

    FORTRANを昨日勉強し始めたばかりの初心者です。 とある大学院の試験問題でわからない部分があるので質問させていただきます。 問題文とソースは以下のとおりです。 以下のFORTRANプログラムについて、標準出力への出力を答えよ。 program main intger :: m(3,3),i,j,k,n n=3 m(1,1)=2; m(1,2)=1; m(1,3)=2 m(2,1)=2; m(2,2)=3; m(2,3)=12 m(3,1)=8; m(3,2)=-6; m(3,3)=2 do i=1,n-1 do j=i+1,n m(j,i)=m(j,i)/m(i,i) ← do k=i+1,n m(j,k)=m(j,k)-m(j,i)*m(i,k) enddo enddo enddo write(6,*) m(2,1),m(2,3),m(3,2),m(3,3) end これを計算していくとm(2,2)=0(計算ミスだったらすみません)となってしまって 上のソースの矢印の部分でエラーが出てしまうと思い、 それだと出力出来ないのではないかと思ったんですが、 それでは答えにならないので質問させていただきました。 コンパイラを使って実行してみたいのもやまやまなんですが、 なにしろ試験まで時間があまりないもので・・・ どなたか回答できる方いらっしゃいましたらよろしくお願いします！

  • fortran　途中まで考えたのですが。。。。

    エラトステネスのふるい（素数の倍数を除いていって残ったのが素数）で３桁の素数を求めて表示したいです。 私が考えたのは、 １・２～９９までの数を素数かどうか調べて、素数を配列に入れていく ２・１００～９９９まで素数の配列の中の数で割って、割り切れたらおしまい。割り切れなかったら表示していく ということです。 しかし下のプログラムではうまく素数配列ができていないようなのです。 ここまでかなり時間がかかったのでこのプログラムに手をいれて これ以外におかしくなるところもどこを直せばいいのか教えてくださるとうれしいです。 C C q223.f C PROGRAM q223 C IMPLICIT NONE C INTEGER N,i,K,s,l REAL a(99),b(99),c(99),X,Y C real M C a(1)=2 a(2)=3 l=2 C DO N=2,99,1 M=N**(0.5) S=M DO i=2,S,1 K=MOD(N,i) IF(K ==0)THEN exit ELSE IF(K /=0)THEN l=l+1 a(l)=N ENDIF ENDDO ENDDO C do N=100,999 do l=1,99 X=a(l) Y=N/X if(Y==0)then exit else if(Y/=0)then write(*,*)N end if end do end do c end よろしくおねがいします

  • この問題の意味が分かりません。

    配列aの全要素を配列bに対して逆順にコピーするプログラムを作成せよ。 なお、二つの配列の要素数は同一であると仮定して良い。 この問題の意味なんですが、要素数はｎであるとして、 (1)a[n-1]をb[n-1]に、a[n-2]をb[n-2]に・・・a[0]をb[0]にコピーしていくということですか？ 結果、配列aの全要素はbと同じ。 (2)それとも、a[0]をb[n-1]に、a[1]をb[n-2]に、・・・a[n-1]をb[0]にコピーしていうことですか？ 結果、配列bはaの全要素を逆順にした要素が入ってる。 ちなみに、どちらの場合もプログラムは分かります。 やっぱり(2)の意味でしょうか？

    • ベストアンサー
    • Java
  • fortran 行列ベクトル積

    行列ベクトル積を計算するプログラムを下のように書いたのですが、実行した結果の答えが実際計算した答えと異なります。初期の要素の設定がおかしいのでしょうか？教えて下さい。よろしくお願いします。 program list2_14 implicit none integer , parameter :: n = 2 real(8) a(n,n), x(n), y(n) integer i, j, k, l a(1,1:2) = (/1.2d0,3.4d0/) a(2,1:2) = (/5.6d0,7.8d0/) x(:) = (/9.0d0,10.0d0/) do i = 1, n y(i) = 0.0d0 do j = 1, n y(i) = y(i) + a(i,j) * x(j) enddo enddo do k = 1, n write(*,*) (a(k,l), l = 1, n) enddo write(*,*) x(:) write(*,*) y(n) end program list2_14 実行結果 1.2 3.4 5.6 7.8 9. 10. 128.4

  • fortranプログラムについての質問です。

    fortranプログラムについての次の問題について質問させてください。 ｉ＝１、２、・・・、５０に対し（i,a500(i)）(i,a2000(i))を次の演算に従い計算せよという問題です。わかりにくいですが ５００と２０００、さらに次の式のaのあとのnやn+1,n-1などは添字です。 an+1(i)=0.5{an(i+1)+an(i-1)-2.0d0an(i)}+an(i)　　(i=2,3,・・・49),　(n=1,2,・・・,500,・・・,2000) ただしa1(i)=0.0(i=1,2,・・・,49), a1(50)=1.0, an(1)=0.0, an(50)=1.0 (n=1,2,3,・・・,500,・・・,2000)　である。 というものです。自分は添字のｎなどについてはプログラム上では無視し、一次元配列ａとｂを宣言して、「ＤＯループも用いて ｉ＝１～４９を回してまずa2（これがb）を求め、そのa2（すなわち）を元の漸化式にの右辺に結果を入れ、それを５００回と２０００回 別々に回す」という操作を考えたのですがうまくいきません。i-1やi+1をどう扱うかや、そもそもこの考え方は合っているのか、 どうするほうが良いのかなどについてアドバイスをいただきたいです。よろしくお願いします。

  • アドバイスお願いします

    今，以下の様なプログラムを作っています． 正の整数ｎを上限とする数字の集合を考える． そのうち a*i+b*j で表現できない値の個数を出力する． a,bは事前に与えられた正の整数であり，i,jは正の整数の変数である． またa,b,nの取り得る最大値は1000000とする． とりあえずfortranで総当たりの方法で作ってみました． integer a,b,n,ans(1000000),c read(*,*)a,b,n 　do i=1,n 　　ans(i)=i 　enddo icount=0 do i=0,n 　do j=0,n 　　c=a*i+b*j 　　if(c.le.n)then 　　　do k=1,n 　　　　if(ans(k).eq.c)then 　　　　　ans(k)=-1 　　　　endif 　　　enddo 　　endif 　enddo enddo do i=1,n 　if(ans(i).gt.0)then 　　icount=icount+1 　endif enddo write(*,*)icount end ｎの値があまり大きくなければ問題なく動くのですが， 最大値である1000000を入力して走らすと，とんでもなく 時間がかかってしまいます． もっとスマートな方法は無いものでしょうか？ アドバイスお願い致します．

  • 配列のサイズが明示的でない計算について

    Fortranを例にとった場合、 行列A,Bの積Cを計算する場合、例えば以下のようなループで計算していました。(Fortran77) do i=1,10 do j=1,10 sum=0.0 do m=1,10 sum=sum+a(i,m)*b(m,j) enddo c(i,j)=sum enddo enddo それが今では、以下1行です。(Fortran95) c=matmul(a,b) 配列１つ１つの成分について計算する必要がないと言えそうです。Pythonなんかは全部こういう思想ですね。 プログラム言語はすべてこういう方向に進んでいるように思えます。 これがやや不安に思えるところがあります。配列のサイズを指定していないのでどこまで使って計算しているかわからないからです。 古い方（上段）では配列のサイズや上限を指定しているので、指定しているところより先がどうなっていてもとりあえず被害はありません。 しかし、サイズについてのケアがないと思いがけないような計算をしてくるんじゃないかと思うのですが。計算のバグ取のときここ大丈夫かな？となって疑心暗鬼ということにならないかと思うのですが。 配列の最大値なども関数１つと言う感じですが、どこまでの範囲での最大なの？と聞きたくなるわけですが。 いかがでしょうか。

  • Fortranコードの疑問

    FORTRANのコードに関して教えてください。 いま作っているプログラムを走らせるとCommand terminated abnormally.と言われて途中で止まってしまいます。これは小さな計算（2000格子点）くらいだと上手くいくのですが、16000格子点くらいになるとダメです。これくらいだとメモリ的にもまだまだ余裕があるので、色々と他の理由を探してみましたが万策尽きました・・・　どうやら write(*,*)n…繰り返し数 n=2000 write(*,*)n do i=1,n -- code ----- enddo としたとき、出力は普通、 0 2000 となってDo文を実行すると思うのですが、上記の場合のように大きな（？）計算をさせようとすると、 1081265151 2000 などと、定義してもいないのにとんでもない値がｎに代入されてしまうようです。何故このような現象が起こるのかご存知の方がいらっしゃいましたら、どうかご教示頂けると幸いです。

  • 並列化プログラムの作成方法について

    物理計算の並列化についてお尋ねします。計算機環境(GPU, メモリ分散（ハブ）、メモリ共有（CPU内での並列））、コンパイラの選択とかいろいろ条件があると思います。また、作りこみの度合い（本気度）によっても違うと思います。特に並列計算は並列化の度合いによって速度が急激に伸びるなどいろいろ特性があると思います。 現時点では、 通常のパソコン（Core-i7)で、コンパイラはintel fortran, OSはWindows10というごく普通の環境です。Cygwinのgfortran, Visual studioなどもあります。 物理プログラムは、構成として以下のようです。 ２次元空間の時間発展の計算です。 do n=1,nt do i=1,nx ! x方向格子へのループ do j=1,ny ! y要綱格子へのループ enddo ! x方向格子ループの端 enddo ! y方向格子ループの橋 enddo ! 時間方向ループの終端 stop end 時間発展のループの中に、空間方向のループが何回も出てくるという構成です。計算が遅くなる理由は空間方向の格子数が極端に増えるからです。 このような場合、計算の速度を向上させるための並列化の方法についてゼロベースで教えて頂きたいのですが。ゼロベースなので、何から手を付けるべきかということですが。 よろしくお願いします。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する
質問する