Repository navigation
Expand file tree
/
Copy pathdensity_wrapper.cpp
More file actions
128 lines (116 loc) · 5.5 KB
/
Copy pathdensity_wrapper.cpp
File metadata and controls
128 lines (116 loc) · 5.5 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
/*
* Licensed to the Apache Software Foundation (ASF) under one
* or more contributor license agreements. See the NOTICE file
* distributed with this work for additional information
* regarding copyright ownership. The ASF licenses this file
* to you under the Apache License, Version 2.0 (the
* "License"); you may not use this file except in compliance
* with the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing,
* software distributed under the License is distributed on an
* "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
* KIND, either express or implied. See the License for the
* specific language governing permissions and limitations
* under the License.
*/
#include <vector>
#include <memory>
#include <nanobind/nanobind.h>
#include <nanobind/intrusive/counter.h>
#include <nanobind/stl/string.h>
#include <nanobind/stl/vector.h>
#include <nanobind/make_iterator.h>
#include <nanobind/eval.h>
#include <numpy/arrayobject.h>
#include "kernel_function.hpp"
#include "density_sketch.hpp"
namespace nb = nanobind;
template<typename T, typename K>
void bind_density_sketch(nb::module_ &m, const char* name) {
using namespace datasketches;
nb::class_<density_sketch<T, K>>(m, name)
.def("__init__", [](density_sketch<T, K>* sk, uint16_t k, uint32_t dim, kernel_function* kernel)
{ K holder(kernel);
new (sk) density_sketch<T, K>(k, dim, holder);
},
nb::arg("k"), nb::arg("dim"), nb::arg("kernel"),
"Creates a new density sketch\n\n"
":param k: controls the size and error of the sketch\n:type k: int\n"
":param dim: dimension of the input data\n:type dim: int\n"
":param kernel: instance of a kernel\n:type kernel: KernelFunction\n"
)
.def("__copy__", [](const density_sketch<T,K>& sk){ return density_sketch<T,K>(sk); })
.def("update", static_cast<void (density_sketch<T, K>::*)(const std::vector<T>&)>(&density_sketch<T, K>::update), nb::arg("vector"),
"Updates the sketch with the given vector")
.def("merge", static_cast<void (density_sketch<T, K>::*)(const density_sketch<T, K>&)>(&density_sketch<T, K>::merge), nb::arg("sketch"),
"Merges the provided sketch into this one")
.def("is_empty", &density_sketch<T, K>::is_empty,
"Returns True if the sketch is empty, otherwise False")
.def_prop_ro("k", &density_sketch<T, K>::get_k,
"The configured parameter k")
.def_prop_ro("dim", &density_sketch<T, K>::get_dim,
"The configured parameter dim")
.def_prop_ro("n", &density_sketch<T, K>::get_n,
"The length of the input stream")
.def_prop_ro("num_retained", &density_sketch<T, K>::get_num_retained,
"The number of retained items (samples) in the sketch")
.def("is_estimation_mode", &density_sketch<T, K>::is_estimation_mode,
"Returns True if the sketch is in estimation mode, otherwise False")
.def("get_estimate", &density_sketch<T, K>::get_estimate, nb::arg("point"),
"Returns an approximate density at the given point")
.def("__str__", [](const density_sketch<T, K>& sk) { return sk.to_string(); },
"Produces a string summary of the sketch")
.def("to_string", &density_sketch<T, K>::to_string, nb::arg("print_levels")=false, nb::arg("print_items")=false,
"Produces a string summary of the sketch")
.def("__iter__", [](const density_sketch<T, K> &sk) {
return nb::make_iterator(nb::type<density_sketch<T,K> >(),
"density_iterator",
sk.begin(),
sk.end());
},
nb::keep_alive<0,1>())
.def("serialize",
[](const density_sketch<T, K>& sk) {
auto bytes = sk.serialize();
return nb::bytes(reinterpret_cast<const char*>(bytes.data()), bytes.size());
},
"Serializes the sketch into a bytes object"
)
.def_static(
"deserialize",
[](const nb::bytes& bytes, kernel_function* kernel) {
K holder(kernel);
return density_sketch<T, K>::deserialize(bytes.c_str(), bytes.size(), holder);
},
nb::arg("bytes"), nb::arg("kernel"),
"Reads a bytes object and returns the corresponding density_sketch"
);
}
int prepare_numpy() {
import_array1(0);
return 0;
}
void init_density(nb::module_ &m) {
using namespace datasketches;
prepare_numpy();
// generic kernel function
nb::class_<kernel_function, KernelFunction>(m, "KernelFunction",
nb::intrusive_ptr<kernel_function>(
[](kernel_function *kf, PyObject *po) noexcept { kf->set_self_py(po); }),
"A generic base class from which user-defined kernels must inherit.")
.def(nb::init())
.def("__call__", &kernel_function::operator(), nb::arg("a"), nb::arg("b"),
"A method to evaluate a kernel with given inputs a and b.\n\n"
":param a: An input vector\n:type a: numpy array\n"
":param b: An input vector\n:type b: numpy array\n"
":return: A vector similarity score\n:rtype: float"
)
;
// the old sketch names can almost be defined, but the kernel_function_holder won't work in init()
//bind_density_sketch<float, gaussian_kernel<float>>(m, "density_floats_sketch");
//bind_density_sketch<double, gaussian_kernel<double>>(m, "density_doubles_sketch");
bind_density_sketch<double, kernel_function_holder>(m, "density_sketch");
}