'nl3_1919_lessvisc' ,
'nl3_135135_lessvisc' ,
'tm1_135135_lessvisc' ,
'tma_135135_lessvisc' ,
'sfhx_135135_lessvisc',
'sfho_135135_lessvisc']
npcs = len(waveform_names)
#
# Build Catalogues
#
fshift_center=1000
print "building catalogues"
(freqaxis, low_cat, high_cat, shift_cat, original_cat, fpeaks, low_sigmas,
high_sigmas) = pmns_pca_utils.build_catalogues(waveform_names,
fshift_center)
delta_f = np.diff(freqaxis)[0]
# Convert to magnitude/phase
shift_mag, shift_phase = pmns_pca_utils.complex_to_polar(shift_cat)
original_mag, original_phase = pmns_pca_utils.complex_to_polar(original_cat)
#
# PCA
#
print "Performing PCA"
shift_pca = pmns_pca_utils.pca_magphase(shift_cat, freqaxis, flow=10)
#
# Compute idealised minimal matches
#
'apr_135135_lessvisc', 'shen_135135_lessvisc', 'dd2_135135_lessvisc',
'dd2_165165_lessvisc', 'nl3_1919_lessvisc', 'nl3_135135_lessvisc',
'tm1_135135_lessvisc', 'tma_135135_lessvisc', 'sfhx_135135_lessvisc',
'sfho_135135_lessvisc'
]
npcs = len(waveform_names)
#
# Build Catalogues
#
fshift_center = 1000
print "building catalogues"
(freqaxis, low_cat, high_cat, shift_cat, original_cat, fpeaks, low_sigmas,
high_sigmas,
false_freqs_shift) = pmns_pca_utils.build_catalogues(waveform_names,
fshift_center)
delta_f = np.diff(freqaxis)[0]
# Convert to magnitude/phase
full_mag, full_phase = pmns_pca_utils.complex_to_polar(original_cat)
shift_mag, shift_phase = pmns_pca_utils.complex_to_polar(shift_cat)
low_mag, low_phase = pmns_pca_utils.complex_to_polar(low_cat)
high_mag, high_phase = pmns_pca_utils.complex_to_polar(high_cat)
waveform = pmns_utils.Waveform(waveform_names[2])
waveform.reproject_waveform()
# ******** #
# Plotting #
# ******** #
imageformats = ['png', 'eps', 'pdf']