Check if R is being passed correctly in `find_d1d2`

[dbrakenhoff]

The value R is multiplied by some lambda twice. It is first passed to bessellsv2 as R = self.rzero * self.aq.lababs[i, j] first, then in find_d1d2() it is multiplied by lab[0], where lab[0] is the first entry in aq.lab2[i, j, :] with dimensions (naq, nint, npint).

[mbakker7]

Can you point to the file where this happens? I cannot find it in besselnumba.py

[dbrakenhoff]

Multiplication of rzero with lambda in potinf of linesink:

ttim/ttim/linesink.py

Line 948 in 1420f99

self.rzero * self.aq.lababs[i, j],

Multiplication of R with the first (and largest lambda):

ttim/ttim/besselnumba.py

Line 477 in 1420f99

d1, d2 = find_d1d2(z1, z2, complex(x, y), R * np.abs(lab[0]))

[mbakker7]

Hmmm. This does look fishy indeed. Good thing is that if it is a bug at least it is conservative. Let's investigate.