[Gretl-devel] Re: MPI data transfer
On 17.10.2019 21:41, Allin Cottrell wrote:
Now in git I've arranged for mpibcast, mpisend and mpirecv to work
directly for bundles -- and as a by-product, enabled MPI transfer for
arrays. I've done a fair amount of testing of "direct" mpibcast of
bundles, and that's currently the default in git.
That is great!!!
The default for mpisend/mpirecv is still the old XML method (since I
just coded the new variant today and it's barely tested) but if you're
building gretl yourself you can try out the new variants thus: look
for the definition of ALT_BUNSEND in lib/src/gretl_mpi.c and change it
from 0 to 1.
I wanted to test new version of BACE, but I get the following error:
--------------
Incompatible types in assignment: series = matrix
*** error in function BACE_initial_model, line 26
series YyY = M.Y_mat
--------------
The strange thing is, that the following code works:
<hansl>
open arma
matrix U = u
series u_new = U
</hansl>
I'll do more testing myself, but if anyone else wants to join
in,
please do. My results so far are very encouraging. Here's some
comparative output from a script that transfers a bundle holding an
array of 20 big matrices:
My results:
old XML version:
gretlmpi 2019d-git
Bieżąca sesja: 2019-10-18 08:03
have 20 matrices
bcast time: 0,0124994s
got 20 matrices
send/recv time: 1,23351s
diff_arrays gave 8,88178e-16
new direct version:
gretlmpi 2019d-git
Bieżąca sesja: 2019-10-18 08:05
have 20 matrices
bcast time: 0,0125704s
got 20 matrices
send/recv time: 0,00509449s
diff_arrays gave 0
Marcin
--
Marcin Błażejowski