Skip to content

MATLAB hold-overs in Fourier Bessel zero finding method #609

@chris-langfield

Description

@chris-langfield

The following should be modified to be a list of arrays (possibly of variable lengths, for each ell), rather than a zero padded NumPy array.

cf. #599

max_num_zeros = max(len(z) for z in zeros)
for i, z in enumerate(zeros):
zeros[i] = np.hstack(
(z, np.zeros(max_num_zeros - len(z), dtype=self.dtype))
)
self.r0 = m_reshape(np.hstack(zeros), (-1, self.ell_max + 1)).astype(self.dtype)

Metadata

Metadata

Labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions