Add parallelization; find r from range
Add parallelization; find r from range
Add penalization phase