[Gretl-devel] Re: MPI data transfer
On Fri, 18 Oct 2019, Marcin Błażejowski wrote:
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!!!
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
--------------
Here's an indication of how to debug this sort of thing:
<hansl>
open data4-10 -q
# matrix m = mnormal(51,1) # works
matrix m = mnormal(50,1) # error
catch series s = m
if $error
printf "nobs=%d, tmax=%d, matrix is %d x %d\n",
$nobs, $tmax, rows(m), cols(m)
endif
</hansl>
The matrix must have a single column, and a number of rows equal to
either $nobs or $tmax.
> 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
That's the new, fast bcast time (current git default). If necessary
you can revert to the XML-based version by setting NEW_BUNCAST to 0 in
lib/src/gretl_mpi.c.
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
Allin