Once #893 is merged, the hyper-resitivity will spatially vary due to 2 factors:
1- it will depend on the level
2- it will depend on local B/n
therefore it won't be obvious to locally know what was the hyper-resistivity.
This issue is for computing the hyper-resistive term in a run getter run.GetHyperResistive() that would essentially calculate the hyper-resistive term locally as a VectorField.
The hyper-resistive coef depends on the Simulation parameter hyper_mode. If equal to constant then the hyper-resistive coef is just the $\nu_0$ provided in Simulation. If set to spatial then it is $\nu_0 B/n$. Warning: this coef is numerically multiplied to the x,y,z components of the laplacian and therefore has a different centering for each component (we project B components on E ones, etc.)
Once #893 is merged, the hyper-resitivity will spatially vary due to 2 factors: 1- it will depend on the level 2- it will depend on local B/n
therefore it won't be obvious to locally know what was the hyper-resistivity.
This issue is for computing the hyper-resistive term in a run getter
run.GetHyperResistive()that would essentially calculate the hyper-resistive term locally as aVectorField.The hyper-resistive coef depends on the
Simulationparameterhyper_mode. If equal toconstantthen the hyper-resistive coef is just the $\nu_0$ provided inSimulation. If set tospatialthen it is $\nu_0 B/n$. Warning: this coef is numerically multiplied to the x,y,z components of the laplacian and therefore has a different centering for each component (we project B components on E ones, etc.)