In [1]:
%matplotlib inline
import os
# os.environ['DD_USE_LATEX'] = '0'
import lib_dd.apps.ccd_single as ccd_single
  • type of kernel
  • lambda
  • c value
  • max number of iterations
  • nr of terms per decade
  • range of tau values
  • data weighting scheme
  • plots (do we need this? probably for the l-curve)

functionality todo

  • data weighting
  • lambda selection: also add possiblity to use automatic lambda

In [2]:
app = ccd_single.ccd_single_app('frequencies2.dat', 'data2.dat')
app.show()


Data summary

Number of frequencies: 25
Minimum frequency: 0.001 Hz
Maximum frequency: 1000.0 Hz

number of SIP spectra: 2


INFO:root:running CCD
INFO:root:single processing
INFO:root:Fitting spectrum 1 of 2
INFO:NDimInv.main:Running inversion
INFO:root:Iteration: 1
INFO:root:Iteration: 2
INFO:root:Iteration: 3
INFO:root:Iteration: 4
INFO:root:Iteration: 5
INFO:root:RMS update below threshold: 0.0001945739028149796 - 0.00019458787689415391 < 1e-05
INFO:root:Fitting spectrum 2 of 2
INFO:NDimInv.main:Running inversion
INFO:root:Iteration: 1
INFO:root:Iteration: 2
INFO:root:Iteration: 3
INFO:root:Iteration: 4
INFO:root:Iteration: 5
INFO:root:RMS update below threshold: 0.0001945739028149796 - 0.00019458787689415391 < 1e-05
finished
INFO:root:running CCD
INFO:root:single processing
INFO:root:Fitting spectrum 1 of 2
INFO:NDimInv.main:Running inversion
INFO:root:Iteration: 1
INFO:root:Iteration: 2
INFO:root:Iteration: 3
INFO:root:Iteration: 4
INFO:root:Iteration: 5
INFO:root:RMS update below threshold: 0.0001945739028149796 - 0.00019458787689415391 < 1e-05
INFO:root:Fitting spectrum 2 of 2
INFO:NDimInv.main:Running inversion
INFO:root:Iteration: 1
INFO:root:Iteration: 2
INFO:root:Iteration: 3
INFO:root:Iteration: 4
INFO:root:Iteration: 5
INFO:root:RMS update below threshold: 0.0001945739028149796 - 0.00019458787689415391 < 1e-05
INFO:root:plotting ... this may take a while
INFO:root:Plotting l-curve for iteration 4
INFO:root:Plotting l-curve for iteration 4
/home/mweigand/.virtualenvs/ccd/lib/python3.5/site-packages/matplotlib/pyplot.py:524: RuntimeWarning: More than 20 figures have been opened. Figures created through the pyplot interface (`matplotlib.pyplot.figure`) are retained until explicitly closed and may consume too much memory. (To control this warning, see the rcParam `figure.max_open_warning`).
  max_open_warning, RuntimeWarning)
finished