Wish: Increase the upper bound for the base p for prime powers in create_galois_field

[ugroempi]

Hi Rob,
your function create_galois_field is very helpful for me (thanks again for exporting it!), and I now have an application for which relatively large values of q are needed. There seems to be a bound on the base p at 50 (lhs::create_galois_field(47^2) works, lhs::create_galois_field(53^2) doesn't). The error message says: GF(2809) = GF(53^2) is not included in this program. To add it, consider modifying gfields.c. This sounds as though it is not much work. If that assumption is correct, it would help me (for a planned package on covering arrays) if the bound could be increased to something somewhat larger, e.g., 120.
Thanks and regards,
Ulrike