Use that r = m/p.
Small update to diagfrob.c
Small changes to diagfrob/test/t-points.
Improve diagfrob test code, check for Weil polynomials.
Include modular charpoly algorithm and improvement for mu.