Now, probes check all particles in the domain.
This is done only when probes are outputed to save CPU time.
However, still some issues and distribution functions are not properly
being calculated.
The code nows offer the possibility to obtain the distribution function
for a specific species in a 3D velocity grid at a determined position.
This is a simple method that just scatter the particles in one cell into
the velocity grid.
Possibility to input initial species distributions (density, velocity
and temperature) via an input file for each species.
New moduleRandom includes function to generate random numbers in
different ways (still uses) the implicit RANDOM_NUMBER().
Maxwellian and Diract Delta distributions have been implemented.
The input for injection of particles should be rewritten to allow more
clear input file.
Now, the solver needs to be an input parameter of the case, to select if
it is for charged or neutral particles.
Resolution of Poisson equation with Dirichlet boundary conditions is
possible. The source vector is the charge density. This resolution is
done in two steps to save computational time:
1. When reading the mesh, the PLU factorization of the K matrix is
computed.
2. In each iteration, the system K*u = f is solved, in which f is the
source vector (charge density) and u is the solution (potential) in
each node.
No case has been added to the repository. This will be done in next
commit.
The 'non-analog' scheme has been commented. It still needs to split
the particle to avoid 'overweight' particles.
