Fast solution calculation (vectorized and optimized)

calculate_solution_fast(
  psi,
  rho,
  s_b,
  s_r,
  s_length,
  total_length,
  gamma,
  min.ploidy,
  max.ploidy,
  min.rho,
  max.rho,
  min.goodness,
  distance_value,
  TheoretMaxdist,
  minimise,
  allow100percent,
  skip_zero_check = FALSE
)