Macro lu_solve 

macro_rules! lu_solve {
    ($manager:expr, $matrix_a:expr, $target_ids:expr) => { ... };
}

Macro for LU solve operations that works with any DataManager implementation This macro implements forward and backward substitution for solving Ax = b in-place where A is an LU-decomposed matrix and b are the target vectors.