There are hooks in the Vlasov solver that allows for running the advection solve across nodes - add these 'back' and deal with boundary conditions. This will also require new treatment for the poisson solve. For now it may be easiest to broadcast all densities to one node, solve Poisson, and send the E field back to the nodes. In future we can do this with Trilinos.
There are hooks in the Vlasov solver that allows for running the advection solve across nodes - add these 'back' and deal with boundary conditions. This will also require new treatment for the poisson solve. For now it may be easiest to broadcast all densities to one node, solve Poisson, and send the E field back to the nodes. In future we can do this with Trilinos.