search

leon-vv / Traceon

Electromagnetic solver and electron tracer
https://traceon.org/
GNU Affero General Public License v3.0
2 stars 4 forks source link

Explore better integrator technique than simpsons rule #6

Closed leon-vv closed 1 year ago

leon-vv commented 1 year ago

Currently we use Simpson integration rule (see solver/util.py). It might be beneficent to instead use a Gaussian quadrature. We need to explore the performance/accuracy benefits. Hopefully using the better tests of issue #4.

leon-vv commented 1 year ago

I tested Gaussian quadrature, but could not easily get better results. Finally settled on Simpson 3/8 rule which is slightly better than what we had before. More research is necessary..

leon-vv commented 1 year ago

Related to #8