Parallel matrix col/row operations give incorrect results

General Usage
,
2

I don’t think you update different columns in the threads. j2 in add_col! is i1 in add_colrow!, which is the same k in mod_col! in all the threads.