mpi_wavelet_factory.cc

Go to the documentation of this file.

#include "catch2/catch_all.hpp"
 
#include <vector>
#include "purify/distribute.h"
#include "purify/logging.h"
#include "purify/mpi_utilities.h"
#include "purify/utilities.h"
#include "purify/wavelet_operator_factory.h"
#include <sopt/mpi/communicator.h>
using namespace purify;
 
TEST_CASE("Serial vs Distributed Operator") {
  purify::logging::set_level("debug");
  auto const world = sopt::mpi::Communicator::World();
  const std::vector<std::tuple<std::string, t_uint>> wavelets{
      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)};
  CAPTURE(world.size());
  CAPTURE(wavelets.size());
  auto const imsizey = 128;
  auto const imsizex = 128;
  const auto sara = sopt::wavelets::SARA(wavelets.begin(), wavelets.end());
  const auto distribute_sara = sopt::wavelets::distribute_sara(sara, world);
  auto op = std::make_shared<sopt::LinearTransform<Vector<t_complex>>>(
      sopt::linear_transform<t_complex>(distribute_sara, imsizey, imsizex, world));
  auto factory_op = factory::wavelet_operator_factory<Vector<t_complex>>(
      factory::distributed_wavelet_operator::mpi_sara, wavelets, imsizey, imsizex);
 
  SECTION("forward") {
    const Vector<t_complex> input =
        Vector<t_complex>::Random(distribute_sara.size() * imsizex * imsizey);
    CAPTURE(input.size());
    REQUIRE((*op * input).isApprox(*factory_op * input));
  }
  SECTION("backward") {
    const Vector<t_complex> input = Vector<t_complex>::Random(imsizex * imsizey);
    REQUIRE((op->adjoint() * input).isApprox(factory_op->adjoint() * input));
  }
}