Worked on the inaccurate results some. Problem boils down to bad
numerical integration methods. Rewriting the wrong matrix entries gives
a good result. Going to attempt to implement Sauter-Schwab quadrature
rules.
Prettied up the kernels a bit.
Implemented the Telles and Piessens quadrature methods. Working super
well! However, there's still some sort of bug in the integration for the
G matrices.
Forgot to commit assembler.py