Line 66 and 163 of lib/linalg/NNLS.cpp use hard-coded values of n_proc_max_for_partial_matrix_ = 15 and nb = 3, which seem to be used mainly for distributing the matrix.
There is no sufficient description as to why these values are chosen.
For better scalability of this solver, shouldn't these values be determined by the users?
Line 66 and 163 of
lib/linalg/NNLS.cppuse hard-coded values ofn_proc_max_for_partial_matrix_ = 15andnb = 3, which seem to be used mainly for distributing the matrix.