wavelet_operator.cc File Reference

#include <chrono>
#include <benchmark/benchmark.h>
#include "benchmarks/utilities.h"
#include "purify/wavelet_operator_factory.h"

Include dependency graph for wavelet_operator.cc:

Go to the source code of this file.

Classes
class	WaveletOperatorFixture
class	WaveletOperatorAdjointFixture

Functions
void	wavelet_operator_constructor (benchmark::State &state)
	RangeMultiplier (2) -> Range(1024, 1024<< 10) ->UseManualTime() ->Repetitions(10) ->ReportAggregatesOnly(true) ->Unit(benchmark::kMillisecond)
	BENCHMARK_DEFINE_F (WaveletOperatorFixture, Apply)(benchmark
	BENCHMARK_DEFINE_F (WaveletOperatorAdjointFixture, Apply)(benchmark
	BENCHMARK_MAIN ()

Function Documentation

BENCHMARK_DEFINE_F() [1/2]

BENCHMARK_DEFINE_F	(	WaveletOperatorAdjointFixture	,
		Apply
	)

Definition at line 91 of file wavelet_operator.cc.

                                                                                {
  t_uint m_imsizex = state.range(0);
  t_uint m_imsizey = state.range(0);
  sopt::wavelets::SARA m_sara{
      std::make_tuple("Dirac", 3u), std::make_tuple("DB1", 3u), std::make_tuple("DB2", 3u),
      std::make_tuple("DB3", 3u),   std::make_tuple("DB4", 3u), std::make_tuple("DB5", 3u),
      std::make_tuple("DB6", 3u),   std::make_tuple("DB7", 3u), std::make_tuple("DB8", 3u)};
 
  sopt::LinearTransform<Vector<t_complex>> m_Psi =
      sopt::linear_transform<t_complex>(m_sara, m_imsizey, m_imsizex);
 
  // Get the number of wavelet coefs
  t_uint const n_wave_coeff = m_sara.size() * m_imsizex * m_imsizey;
 
  // Apply Psi to a temporary vector
  Vector<t_complex> temp;
  Vector<t_complex> const x = Vector<t_complex>::Random(m_imsizex * m_imsizey);
 
  while (state.KeepRunning()) {
    auto start = std::chrono::high_resolution_clock::now();
    temp = m_Psi.adjoint() * x;
    auto end = std::chrono::high_resolution_clock::now();
    state.SetIterationTime(b_utilities::duration(start, end));
  }
}

References b_utilities::duration().

BENCHMARK_DEFINE_F() [2/2]

BENCHMARK_DEFINE_F	(	WaveletOperatorFixture	,
		Apply
	)

Definition at line 65 of file wavelet_operator.cc.

                                                                         {
  t_uint m_imsizex = state.range(0);
  t_uint m_imsizey = state.range(0);
  sopt::wavelets::SARA m_sara{
      std::make_tuple("Dirac", 3u), std::make_tuple("DB1", 3u), std::make_tuple("DB2", 3u),
      std::make_tuple("DB3", 3u),   std::make_tuple("DB4", 3u), std::make_tuple("DB5", 3u),
      std::make_tuple("DB6", 3u),   std::make_tuple("DB7", 3u), std::make_tuple("DB8", 3u)};
 
  sopt::LinearTransform<Vector<t_complex>> m_Psi =
      sopt::linear_transform<t_complex>(m_sara, m_imsizey, m_imsizex);
 
  // Get the number of wavelet coefs
  t_uint const n_wave_coeff = m_sara.size() * m_imsizex * m_imsizey;
 
  // Apply Psi to a temporary vector
  Vector<t_complex> temp;
  Vector<t_complex> const x = Vector<t_complex>::Random(n_wave_coeff);
 
  while (state.KeepRunning()) {
    auto start = std::chrono::high_resolution_clock::now();
    temp = m_Psi * x;
    auto end = std::chrono::high_resolution_clock::now();
    state.SetIterationTime(b_utilities::duration(start, end));
  }
}

References b_utilities::duration().

BENCHMARK_MAIN()

BENCHMARK_MAIN

(

)

RangeMultiplier()

RangeMultiplier

(

2

)

-> Range(1024, 1024<< 10) ->UseManualTime() ->Repetitions(10) ->ReportAggregatesOnly(true) ->Unit(benchmark::kMillisecond)

wavelet_operator_constructor()

void wavelet_operator_constructor

(

benchmark::State &

state

)

Definition at line 10 of file wavelet_operator.cc.

                                                         {
  // Image size
  t_uint m_imsizex = state.range(0);
  t_uint m_imsizey = state.range(0);
 
  // benchmark the creation of measurement operator
  while (state.KeepRunning()) {
    auto start = std::chrono::high_resolution_clock::now();
 
    const sopt::wavelets::SARA m_sara{
        std::make_tuple("Dirac", 3u), std::make_tuple("DB1", 3u), std::make_tuple("DB2", 3u),
        std::make_tuple("DB3", 3u),   std::make_tuple("DB4", 3u), std::make_tuple("DB5", 3u),
        std::make_tuple("DB6", 3u),   std::make_tuple("DB7", 3u), std::make_tuple("DB8", 3u)};
 
    sopt::LinearTransform<Vector<t_complex>> Psi =
        sopt::linear_transform<t_complex>(m_sara, m_imsizey, m_imsizex);
 
    auto end = std::chrono::high_resolution_clock::now();
 
    state.SetIterationTime(b_utilities::duration(start, end));
  }
}

References b_utilities::duration().