[PIC] Poisson solver about plasmapy HOT 8 CLOSED

I have been thinking about it a quite a bit and my idea would be to have a function that takes every possible thing (grid, physics parameters, field values, time step, number of steps, anything else...) as a parameter and updates the fields for the desired number of steps, and finally returns the updated fields. Inside the solver, we would have cases for cartesian, cylindrical, spherical and whatever grids, and maybe a few different algorithms. A pure function without any side effects.

Then, we can have wrappers on top of it for different cases we might want.

from plasmapy.

About the EM solvers, I was wondering if we should re-invent the wheel or use an exsting EM solver. For a modular framework we would want a separate Maxwell solver for the fields.

from plasmapy.

On the one hand, use of separate library will reduce our initial workload but it will create an extra dependency. Maybe depending on the license, we could simply merge one snapshot of an existing package and take it as a starting point.

from plasmapy.

Is there a clean way to abstract the idea of an EM solver over the data structure such that you can drop-in any solver you want (provided you have a thin wrapper for some external solvers)?

from plasmapy.

Do these kinds of EM solvers go step by step, solving adjacent grid cells or does it take all the cells at once and try to come up with a global solution? If it's the former, and if the equations don't change between coordinate systems, then perhaps you can handle different grid/coordinate systems with decorators.

from plasmapy.

A reasonable object oriented structure would probably also work for code organization. I don't have strong opinions on how to organize the computation itself in a way that's parallelizable.

And with the direction we've been taking lately, I fear our code may soon approach 30% decorator applications. 😃

from plasmapy.

This issue will be closed in 90 days due to five years of inactivity unless the stale label or comment is removed.

from plasmapy.

This issue was closed because it has had no activity for five years.

from plasmapy.