PICfun

This repository implements the PIC algorithm described by Chen et al. With this project I'm testing how to accelerate iterative linear solvers in simple PIC algorithms with ML.

This code has been used for arXiv:2110.12444

PIC equations

The code solves the equations

on a 1-dimensional grid. The new electric field is determined by

. Here

$$SM[Q] = (Q_{i-1} + 2Q_i + Q_{i+1}) / 4$$

is a binomial smoothing operator and angular brackets denote spatial averaging.

A good point to start the simulation is by running

$ julia implicit_gmres.jl