updated least_interpolant to accept a predefined index set. This is producing the write coefficient values but the ordering is messed up