two sided spectra

[apatlpo]

I would need for my applications to keep two-sided spectra.
This doesn't seem to be implemented.
Do you think it would be possible to modify the library to do that ?

[krischer]

What do you mean by two sided spectra? Also the negative frequencies? The Python wrapper currently only wraps the single and double precision floating point variants of the spectral computation.

But a single precision complex version does exist on the fortran side:

mtspec/mtspec/src/src/mtspec.f90

Line 1520 in c34b648

subroutine mtspec_c (npts,dt,x,tbp,kspec,nf,freq,spec, &

It should be fairly simple to add to the existing wrappers. I'm a bit short on time recently so I cannot tackle it but I'm happy to merge a PR and help along the way!

[apatlpo]

I mean the spectrum of a complex variable, for which we need to retain negative frequencies (as they will differ from positive ones).

What would the steps to wrap mtspec_c ?

Would you have another wrapping that could be used as a template?

[krischer]

I think it should just be two steps or so:

1. Add a complex64 variant here:

   mtspec/mtspec/multitaper.py

   Line 752 in c34b648

   class _MtspecType(object):

   (Basically maps numpy dtype to the corresponding fortran subroutine).

2. Maybe that is it already. If not: Fiddle a bit here:

   mtspec/mtspec/multitaper.py

   Line 255 in c34b648

   # Verbose mode on or off.

Testing should also be fairly simple: Passing in data with all complex values set to zero should result in a symmetric spectrum which can be tested against the already existing wrappers.

[apatlpo]

started working on it: #30
it is not as simple as I hoped for