3 ms·
He's proposing to skip the whole computation of the parts[] array and within the kernel, replace for (int i = 0; i < parts[bx].length; i++) { // bx is th
by tromp 3y ago
He's proposing to skip the whole computation of the parts[] array and within the kernel, replace
for (int i = 0; i < parts[bx].length; i++) { // bx is the global thread index
char c = buffer[parts[bx].offset-buffer_offset + i];
by something like
long long split_size = size / num_parts;
long long offset = bx * split_size;
while (buffer[offset++] != '\n')
;
for (int i = 0; buffer[offset+i] != '\n'; i++) {
char c = buffer[offset + i];
(ignoring how to deal with buffer_offset for now).
- tspeterkim 3y agofor (int i = 0; buffer[offset+i] != '\n'; i++) { This would only process the current line, though. Here, each thread processes ~split_size bytes (multiple lines). Even if were to read multiple lines, how would a thread know when to stop? (at which offset?) And when it does, it should communicate where it stopped with the other threads to prevent re-reading the buffer. My brain's hurting now.
- tromp 3y agoIt can compute both the starting offset for itself and next_offset for (bx+1)...
- kevincox 3y agoYeah, this is a bug in their quick example code. Basically think that you simply give each thread a start and stop offset with simple byte division. Then each thread handles the first line that starts after the start offset up to but excluding the first line that starts after its end offset. So there will be a small amount of "overread" at the boundaries but this is probably negligible. let chunk_size = buf.length / threads; let start_byte = chunk_size * bx; let end_byte = start_byte + chunk_size; let i = start_byte; while (buffer[i++] != '\n') {} // Skip until the first new row. while (i < end_byte) { let row_start = i; let row_end = i; while (buffer[row_end++] != '\n') {} process_row(buffer, row_start, row_end); i = row_end + 1; } Basically you first roughly split the data with simple byte division. Then each thread aligns itself to the underling data chunks. This alignment can be done in parallel across all threads rather than being part of a serial step that examines every byte before the parallel work starts. You need to take care that the alignment each thread does doesn't skip or duplicate any rows, but for simple data formats like this I don't think that should be a major difficulty.
- tspeterkim 3y agoI tried this out today. While it works (no longer a pre-split step required), it makes the CUDA kernel run ridiculously slow. I believe it's because of the while loop: while (i < end_byte) { Comparing it to my original solution, 50X divergent branches are introduced! (ncu profiling) The only difference between the two is that the for loop could deterministically iterate, yet this while loop iterates for an unknown amount (at kernel launch time). I admit, I don't perfectly understand the reason. But this is the most likely culprit.