ddemidov

#include <vector>
#include <iostream>

#include <amgcl/backend/builtin.hpp>
#include <amgcl/value_type/static_matrix.hpp>
#include <amgcl/adapter/crs_tuple.hpp>
#include <amgcl/adapter/block_matrix.hpp>

#include <amgcl/mpi/distributed_matrix.hpp>
#include <amgcl/mpi/make_solver.hpp>

ddemidov

#include <vector>
#include <array>    
#include <iostream>
#include <chrono>
#include <cassert>
#include <cmath>
#include <omp.h>

template<
        class TContainerType,

ddemidov

#include <iostream>
#include <vector>
#include <atomic>
#include <chrono>
#include <cmath>
#include <omp.h>

double now() {
    return omp_get_wtime();
}

ddemidov

VEX_FUNCTION(float, l2_distance, (size_t, idx)(uint32_t, num_dim)(float*, query)(float*, candidates),
              float d = 0;
              for (uint i = 0; i < num_dim; ++i) {
                  d += pow(query[i] - candidates[idx * num_dim + i], 2);
              }
              return d;
            );
            
Knn::ValueVec Knn::search_h(const ValueVec& query, size_t query_offset) const {
  prof_.tic_cl("distance");

ddemidov

#include <vexcl/vexcl.hpp>

vex::vector<cl_float2>
hilbert(vex::vector<float> &in, std::size_t rows, std::size_t cols) {
    // No need to pass context in:
    const auto &ctx = in.queue_list();

    // Use function instead of vector. Should save global memory IO.
    VEX_FUNCTION(h, float(size_t/*prm1 = idx*/, size_t/*prm2 = rows*/),
            VEX_STRINGIZE_SOURCE(

ddemidov

/*
 * third.cpp
 * Date: 2012-11-01
 * Author: Karsten Ahnert (karsten.ahnert@gmx.de)
 */

#include <iostream>
#include <boost/utility/enable_if.hpp>

#define tab "\t"