sp-mz 3.31; z_solve_LU_sp-mz.D.16