''' These are the unittests for the pycbc PSD module. ''' import sys import os import tempfile import pycbc import pycbc.psd from pycbc.types import TimeSeries, FrequencySeries from pycbc.fft import ifft from pycbc.fft.fftw import set_measure_level import unittest import numpy from utils import parse_args_all_schemes, simple_exit set_measure_level(0) _scheme, _context = parse_args_all_schemes("PSD") class TestPSD(unittest.TestCase): def setUp(self): self.scheme = _scheme self.context = _context self.psd_len = 1024 self.psd_delta_f = 0.1 self.psd_low_freq_cutoff = 10. # generate 1/f noise for testing PSD estimation noise_size = 524288 sample_freq = 4096. delta_f = sample_freq / noise_size numpy.random.seed(132435)
parser.add_option('--num-threads', type=int, help='set the number of threads to use', default=1) parser.add_option('--import-float-wisdom', type=str, help='import an FFTW float wisdom file') parser.add_option('--export-float-wisdom', type=str, help='export an FFTW float wisdom file') (options, args) = parser.parse_args() #Changing the optvalues to a dict makes them easier to read _options = vars(options) if _options['scheme'] == 'cpu': ctx = CPUScheme(num_threads=options.num_threads) if options.backend == 'fftw': from pycbc.fft.fftw import set_measure_level, set_threads_backend with ctx: set_measure_level(options.measure_level) if options.num_threads != 1: set_threads_backend('openmp') if _options['scheme'] == 'cuda': ctx = CUDAScheme(device_num=_options['devicenum']) if _options['scheme'] == 'opencl': ctx = OpenCLScheme(device_num=_options['devicenum']) niter = options.iterations if type(ctx) is CUDAScheme: print "RUNNING ON ", ctx.device.name() else: print "RUNNING ON CPU"
# ''' These are the unittests for the pycbc PSD module. ''' from __future__ import division import os import tempfile import pycbc import pycbc.psd from pycbc.types import TimeSeries, FrequencySeries from pycbc.fft import ifft from pycbc.fft.fftw import set_measure_level import unittest import numpy from utils import parse_args_all_schemes, simple_exit set_measure_level(0) _scheme, _context = parse_args_all_schemes("PSD") class TestPSD(unittest.TestCase): def setUp(self): self.scheme = _scheme self.context = _context self.psd_len = 1024 self.psd_delta_f = 0.1 self.psd_low_freq_cutoff = 10. # generate 1/f noise for testing PSD estimation noise_size = 524288 sample_freq = 4096. delta_f = sample_freq / noise_size
help='import an FFTW float wisdom file') parser.add_option('--export-float-wisdom', type=str, help='export an FFTW float wisdom file') (options, args) = parser.parse_args() #Changing the optvalues to a dict makes them easier to read _options = vars(options) if _options['scheme'] == 'cpu': ctx = CPUScheme(num_threads=options.num_threads) if options.backend == 'fftw': from pycbc.fft.fftw import set_measure_level, set_threads_backend with ctx: set_measure_level(options.measure_level) if options.num_threads != 1: set_threads_backend('openmp') if _options['scheme'] == 'cuda': ctx = CUDAScheme(device_num=_options['devicenum']) if _options['scheme'] == 'opencl': ctx = OpenCLScheme(device_num=_options['devicenum']) niter = options.iterations if type(ctx) is CUDAScheme: print("RUNNING ON ", ctx.device.name()) else: print("RUNNING ON CPU") print(type(ctx))