2017-02-09 23 views
0

私は、並列コードの作成/実行にかなり慣れています。現在、私はプロセスのための感触を得るために並列コードを書く際の基本的なチュートリアルを試しています。私のコンピュータはMpichとubuntuを使っています。MPI_Send、無効なタグで致命的なエラーが発生しました

私はこのページの「ベクトルを合計するための完全な並列プログラム」と題したコードを実行しようとしています:http://condor.cc.ku.edu/~grobe/docs/intro-MPI.shtml

および/番号を入力するためのプロンプトが表示された後、実行時に次のエラーが発生した午前:

Fatal error in MPI_Send: Invalid tag, error stack: 
MPI_Send(174): MPI_Send(buf=0x7ffeab0f2d3c, count=1, MPI_INT, dest=1, tag=1157242880, MPI_COMM_WORLD) failed 
MPI_Send(101): Invalid tag, value is 1157242880 

コンパイル中に私はまた、警告が表示されます。

sumvecp.f90:41:23: 

    call mpi_send(vector(start_row),num_rows_to_send, mpi_real, an_id, send_data_tag, mpi_comm_world,ierr) 
         1 
Warning: Legacy Extension: REAL array index at (1) 

これは私のコード

です
program sumvecp 

include '/usr/include/mpi/mpif.h' 

parameter (max_rows = 10000000) 
parameter (send_data_tag = 2001, return_data_tag = 2002) 

integer my_id, root_proces, ierr, status(mpi_status_size) 
integer num_procs, an_id, num_rows_to_receive 
integer avg_rows_per_process, num_rows,num_rows_to_send 

real vector(max_rows), vector2(max_rows), partial_sum, sum 


root_process = 0 

call mpi_init(ierr) 

call mpi_comm_rank(mpi_comm_world,my_id,ierr) 
call mpi_comm_size(mpi_comm_world,num_procs,ierr) 

if (my_id .eq. root_process) then 
    print *, "please enter the number of numbers to sum: " 
    read *, num_rows 
    if (num_rows .gt. max_rows) stop "Too many numbers." 

    avg_rows_per_process = num_rows/num_procs 

    do ii = 1,num_rows 
     vector(ii) = float(ii) 
    end do 

    do an_id = 1, num_procs -1 
     start_row = (an_id*avg_rows_per_process) +1 
     end_row = start_row + avg_rows_per_process -1 
     if (an_id .eq. (num_procs - 1)) end_row = num_rows 
     num_rows_to_send = end_row - start_row + 1 

     call mpi_send(num_rows_to_send, 1, mpi_int, an_id, send_data_tag, mpi_comm_world,ierr) 

     call mpi_send(vector(start_row),num_rows_to_send, mpi_real, an_id, send_data_tag, mpi_comm_world,ierr) 
    end do 

    summ = 0.0 
    do ii = 1, avg_rows_per_process 
     summ = summ + vector(ii) 
    end do 

    print *,"sum", summ, "calculated by the root process." 

    do an_id =1, num_procs -1 
     call mpi_recv(partial_sum, 1, mpi_real, mpi_any_source, mpi_any_tag, mpi_comm_world, status, ierr) 

     sender = status(mpi_source) 
     print *, "partial sum", partial_sum, "returned from process", sender 
     summ = summ + partial_sum 
    end do 

    print *, "The grand total is: ", sum 

else 
    call mpi_recv(num_rows_to_receive, 1, mpi_int, root_process, mpi_any_tag, mpi_comm_world,status,ierr) 

    call mpi_recv(vector2,num_rows_to_received, mpi_real,root_process,mpi_any_tag,mpi_comm_world,status,ierr) 

    num_rows_received = num_rows_to_receive 

    partial_sum = 0.0 
    do ii=1,num_rows_received 
     partial_sum = partial_sum + vector2(ii) 
    end do 

    call mpi_send(partial_sum,1,mpi_real,root_process,return_data_tag,mpi_comm_world,ierr) 
endif 

call mpi_finalize(ierr) 
stop 
end 

答えて

1

IMPLICIT NONEが不足しており、宣言されていない変数が多数あります。

send_data_tag = 2001, return_data_tag = 2002 

が暗黙的real変数やないinteger秒であるため、

報告されたエラーです。しかし、おそらくもっと多くの問題があります。

最初にIMPLICIT NONEと宣言または変数を追加します。また、include '/usr/include/mpi/mpif.h'の代わりにuse mpiを入力することを強くお勧めします。これは、より多くの問題を見つけるのに役立つでしょう。


コードは一部のWebサイトからコピーされています。このコードは明らかに間違っているので、私はこのウェブサイトを信頼しません。

+1

今週も同じ基本的な問題で多くの質問があります。 'IMPLICIT NONE'は**必須**です!!! –

+0

IMPLICIT NONEと宣言されていない変数を宣言した後、正常に動作しました。 PythonからFortranに来た。 (それが明らかでない場合)ありがとう! –

関連する問題