doris/be/src/exprs/function/array/function_array_distance.cpp at 4ef2fbb7a434c6ee67fa6cf550c3ffc0116f1f45 · apache/doris · GitHub

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
// 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 "exprs/function/array/function_array_distance.h"

#include "exprs/function/simple_function_factory.h"

namespace doris {

FAISS_PRAGMA_IMPRECISE_FUNCTION_BEGIN
float CosineDistance::distance(const float* x, const float* y, size_t d) {
    float dot_prod = 0;
    float squared_x = 0;
    float squared_y = 0;
    for (size_t i = 0; i < d; ++i) {
        dot_prod += x[i] * y[i];
        squared_x += x[i] * x[i];
        squared_y += y[i] * y[i];
    }
    if (squared_x == 0 or squared_y == 0) {
        return 2.0f;
    }
    return 1 - dot_prod / sqrt(squared_x * squared_y);
}
FAISS_PRAGMA_IMPRECISE_FUNCTION_END

FAISS_PRAGMA_IMPRECISE_FUNCTION_BEGIN
float CosineSimilarity::distance(const float* x, const float* y, size_t d) {
    float dot_prod = 0;
    float squared_x = 0;
    float squared_y = 0;
    for (size_t i = 0; i < d; ++i) {
        dot_prod += x[i] * y[i];
        squared_x += x[i] * x[i];
        squared_y += y[i] * y[i];
    }
    if (squared_x == 0 or squared_y == 0) {
        return 0.0f;
    }
    return dot_prod / sqrt(squared_x * squared_y);
}
FAISS_PRAGMA_IMPRECISE_FUNCTION_END

void register_function_array_distance(SimpleFunctionFactory& factory) {
    factory.register_function<FunctionArrayDistance<L1Distance>>();
    factory.register_function<FunctionArrayDistance<L2Distance>>();
    factory.register_function<FunctionArrayDistance<CosineDistance>>();
    factory.register_function<FunctionArrayDistance<CosineSimilarity>>();
    factory.register_function<FunctionArrayDistance<InnerProduct>>();
    factory.register_function<FunctionArrayDistance<L2DistanceApproximate>>();
    factory.register_function<FunctionArrayDistance<InnerProductApproximate>>();
    factory.register_function<FunctionArrayDistance<CosineDistanceApproximate>>();
    factory.register_function<FunctionArrayDistance<CosineSimilarityApproximate>>();
    factory.register_function<FunctionArrayDistance<HammingDistance, TYPE_BIGINT>>();
    factory.register_function<FunctionArrayDistance<HammingDistanceApproximate, TYPE_BIGINT>>();
}

} // namespace doris