Signal Subspace

1-D and ensemble signal subspace analysis with methods such as Esprit and RootMusic in Fortran, C, and Python

based in part upon the Spectral Analysis Lib public domain code

The core subspace code is written in Fortran 2008 and is called from other languages (Python, C). Since the programs are Fortran / Python based, they should compile and run in virtually any platform from embedded to supercomputer.

cmake -B build

cmake --build build

If Lapack is not available, it is built automatically.

In 2023-2024, F2PY and Numpy are going through a transition of build systems. CMake script build the f2py bindings are used to build the f2py targets that allow Python use of this library.

cmake -B build -Dpython=yes

cmake --build build

In the examples below, observe the frequency estimates printed along with their corresponding eigenvalues. A larger eigenvalue is increased confidence in that particular frequency estimate.

ESPRIT examples

All example use a noisy sinusoid. Some are "complex" using complex numbers, while others use real numbers.

There are two versions of the complex program, one a full accuracy using double complex numbers, and the other using single real numbers as input. The single real (4 bytes/number) runs about 4 times faster than the double complex (16 bytes/number) program.

build/f_esprit_cmpl

build/f_esprit_real

Call Fortran Esprit from C using real single precision float:

build/c_esprit

Call Fortran Esprit from C++ using real single precision float:

build/cpp_esprit

Plots comparing real vs. complex ESPRIT:

python BasicEspritExample.py

pip install issue Procedure ‘random_init’ called at (1) is not explicitly declared

                 from build/src.linux-x86_64-3.6/build/src.linux-x86_64-3.6/fortranobject.c:2:
/home/lulu/.local/lib/python3.6/site-packages/numpy/core/include/numpy/npy_1_7_deprecated_api.h:17:2: warning: #warning "Using deprecated NumPy API, disable it with " "#define NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION" [-Wcpp]
 #warning "Using deprecated NumPy API, disable it with " \
  ^~~~~~~
compiling Fortran 90 module sources
Fortran f77 compiler: /usr/bin/gfortran -Wall -g -ffixed-form -fno-second-underscore -fPIC -O3 -funroll-loops
Fortran f90 compiler: /usr/bin/gfortran -Wall -g -fno-second-underscore -fPIC -O3 -funroll-loops
Fortran fix compiler: /usr/bin/gfortran -Wall -g -ffixed-form -fno-second-underscore -Wall -g -fno-second-underscore -fPIC -O3 -funroll-loops
creating build/temp.linux-x86_64-3.6/src
compile options: '-Ibuild/src.linux-x86_64-3.6/build/src.linux-x86_64-3.6 -I/home/lulu/.local/lib/python3.6/site-packages/numpy/core/include -I/usr/include/python3.6m -c'
extra options: '-Jbuild/temp.linux-x86_64-3.6/ -Ibuild/temp.linux-x86_64-3.6/'
gfortran:f90: src/comm_f2py.f90
src/comm_f2py.f90:18:44:

 call random_init(repeatable, image_distinct)
                                            1
Error: Procedure ‘random_init’ called at (1) is not explicitly declared
error: Command "/usr/bin/gfortran -Wall -g -fno-second-underscore -fPIC -O3 -funroll-loops -Ibuild/src.linux-x86_64-3.6/build/src.linux-x86_64-3.6 -I/home/lulu/.local/lib/python3.6/site-packages/numpy/core/include -I/usr/include/python3.6m -c -c src/comm_f2py.f90 -o build/temp.linux-x86_64-3.6/src/comm_f2py.o -Jbuild/temp.linux-x86_64-3.6/ -Ibuild/temp.linux-x86_64-3.6/" failed with exit status 1

----------------------------------------

Command "/usr/bin/python3 -u -c "import setuptools, tokenize;file='/tmp/pip-build-6fdl9jiu/signal-subspace/setup.py';f=getattr(tokenize, 'open', open)(file);code=f.read().replace('\r\n', '\n');f.close();exec(compile(code, file, 'exec'))" install --record /tmp/pip-jk2fl2tn-record/install-record.txt --single-version-externally-managed --compile --user --prefix=" failed with error code 1 in /tmp/pip-build-6fdl9jiu/signal-subspace/

scivision / signal_subspace Goto Github PK

signal_subspace's Introduction

Signal Subspace

ESPRIT examples

signal_subspace's People

Contributors

Stargazers

Watchers

Forkers

signal_subspace's Issues

no setup.cmake file

pip install issue Procedure ‘random_init’ called at (1) is not explicitly declared

Recommend Projects

React

Vue.js

Typescript

TensorFlow

Django

Laravel

D3

Recommend Topics

javascript

web

server

Machine learning

Visualization

Game

Recommend Org

Facebook

Microsoft

Google

Alibaba

D3

Tencent