add function in(), thus parallelize feas_di_2D()