Re: [Gretl-devel] MPI doc
On Mon, 16 Feb 2015, Sven Schreiber wrote:
Advancing line-by-line here... Next error below the sample script:
<hansl>
mpi --np=2
if $mpirank == 0 # parent process
matrix b = ones(1000,1)
matrix a = zeros(1000,1)
endif
# -- up to here no errors --
# divide up the input data
mpiscatter(&b, byrows)
The trouble here is that the matrices a and b exist only in the rank 0
process; they need to exist everywhere.
<hansl>
mpi --np=2
matrix a
matrix b
if $mpirank == 0
b = ones(1000,1)
a = zeros(1000,1)
endif
mpiscatter(&b, byrows)
</hansl>
Allin