Whamp (matlab): No solution found for array of wavevectors

[ikitcheng]

Step 1: Latest code?

• This issue is present in the most recent code (i.e. after running git pull this is still a problem).

Step 2: Describe your environment

• irfu-matlab branch (e.g. master, devel): master
• Matlab version used (e.g. Matlab R2020b): R2020a
• Operating system (e.g. Windows 10, Mac OS X 10.15.6, Linux Ubuntu 20.04): Windows 10

• A 64-bit installation of Matlab was used.

Step 3: Describe the problem

For the example run (as in the whamp documentation): Solutions were found for certain arrays of wavevectors e.g.kperp=0, kpar = [0.0008 0.0002 0.0025]. But, if kperp=0, kpar = [0.0006 0.0002 0.0025], then no solutions were found, even though this set of wavevectors include the kperp and kpar points which did have solution before.

N.B. the numbers in the square bracket mean [k_start k_step k_end].

Expected behavior

For example run with kpar and kperp modified to: kperp=0, kpar = [0.0006 0.0002 0.0025]
Output.flagSolutionFound = [0,1,1,1,1,1,1,1,1,1]

Actual behavior

For example run with k modified to: kperp=0, kpar = [0.0006 0.0002 0.0025]
Output.flagSolutionFound = [0,0,0,0,0,0,0,0,0,0]

Steps to reproduce:

See "Relevant code" section below.

Relevant code:

To reproduce unexpected behaviour:

Oxygen = struct('m',16,'n',1,'t',10,'a',5,'vd',1);
Electrons = struct('m',0,'n',1,'t',100);
PlasmaModel = struct('B',100);
PlasmaModel.Species = {Oxygen,Electrons};
InputParameters = struct('fstart',0.1,'kperp',0,'kpar',[0.0006 0.0002 0.0025],'useLog',0);
Output = whamp.run(PlasmaModel,InputParameters)