optimize AD derivatives and add Aqua.jl tests

[fgerick]

this gives about 2x speedup for matrix assembly by reducing massively the number of evaluations of Jacobi polynomials.

[codecov]

Codecov Report

✅ All modified and coverable lines are covered by tests.
✅ Project coverage is 82.83%. Comparing base (a8912a5) to head (8c4a2ef).
⚠️ Report is 6 commits behind head on main.

Additional details and impacted files

@@            Coverage Diff             @@
##             main      #28      +/-   ##
==========================================
+ Coverage   82.47%   82.83%   +0.35%     
==========================================
  Files          32       32              
  Lines        3481     3554      +73     
==========================================
+ Hits         2871     2944      +73     
  Misses        610      610              

☔ View full report in Codecov by Harness.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:

• ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.