RuntimeWarning during the processing from get_j/statistics.py: divide by zero encountered in log p = np.polyfit(np.log(dias[ind]), np.log(number_distribution[ind]), 1)

[arsalanmstn]

Describe the bug
A clear and concise description of what the bug is.

To Reproduce
Detailed steps to reproduce the behavior.

Expected behavior
A clear and concise description of what you expected to happen.

Screenshots
If applicable add screenshots to help explain your problem.

Additional context
Add any other context about the problem here.

[arsalanmstn]

Suggestion: update the following line in the get_j function:
p = np.polyfit(np.log(dias[ind]), np.log(number_distribution[ind]), 1)
to:
p = np.polyfit(np.log(dias[ind]), np.log(number_distribution[ind], where= number_distribution[ind]> 0), 1)