fel pa num_factors.py