[Gretl-devel] Re: Speed of matrix "block" operations

Thanks, Jack. So, given the options I posited, our machines agree on a best chunk size of (25 * 5000 * 8) bytes (=~ 1 MB) for use with memcpy. (5000 being the number of rows in the matrix to be copied, and 8 the number of bytes to represent a double-precision floating point value.) Now, to optimize libgretl's copying of contiguous data, we just have to figure out how that relates to the size of L1 or L2 cache, or whatever is truly the relevant hardware parameter here!