1) move from 3D to 1D
-> done, but we see that initial amplitude (phase?) depends on the number of processors: 1 versus 8
random number initial cond -> const phase -> still a problem

-> check if without conductivity we have a problem with processors -> No problem!

Calculation only on one processor for conductivity?
Global current (mean-field/Tomo) is coded only for one processor. Not trustable for multiple procs! only for 1D, not 2D and 3D due to a sum to find Erms, Brms.
Local current is?

2) Test the case with local conductivity that we coded.
Local conductivity works!

3) Include causality arguments - 
https://arxiv.org/abs/1110.2891

4) Numerical instability/artifacts (oscillations in phi and phidot) after the end of inflation, not dependent on resolution 512, 1024, 2048. Second-order derivatives have infinite errors.

5) Vacuum=T we do not compute sigma
Benefit: compute sigma even if it is small, but not take into account.

6) Why E(rms) and B(rms) blow up at the end?
time step related?

7) check smaller values of alpha=10?

----------------------------
8) check By

8') Should we include helicity?

9) BD modes and Schwinger scale?

10) charge and its running

11) keep in mind that phi0 is not being computed but set in the input by hand!
dotphi0=-Sqrt[1/12Pi]mphi Mpl
t0 is computed

11) Final point: How it affects magnetogenesis?

12) Check that local and non local cases work with the new charge implementation.
Implement in new framework.

13) another average for <sigmaE E^2>, <sigmaE> <E^2>?

14) test that everything works correctly on multiple processors.

New things to do:

1) prescription for E, B->0
2) integration of the current
3) Broken power-law /|\ for phi and all the fields?
(kpeak/k1)^3



For the draft:

0a) BD conditions

0b) Make sure the local case works with the new implementation of echarge.

1) Linear regime compare with backreation regime/ with Schwinger effect
 
2)comparison with rho_chi, echarge

3) Consequences for magnetogenesis