Open Shiva-Konakanchi opened 3 years ago
Shiva, I can reproduce this issue. I think internally we are setting 'vpei' in params.py based on the supplied vector potential, but I'm guessing we are not taking it into account when we invoke the kernels.
Hi Shriram, thanks for taking a look. I agree, somehow '_vpei' isn't being used in the evolution kernels. Also, another hint in this direction is that if I use 'external_field' to specify the value of the external field in the GLSolver object instead of 'uniform_external_field', that value isn't being used in calculations either. Note that the documentation says that for now 'external_field' only accepts a number (although it is designed for other things in the future).
Produces this output:
However, passing A produces different output as shown below: