Function l2_normalize

vexdb_pg/src/floatvector.cpp:561–588

PG_FUNCTION_INFO_V1(l2_normalize);
Datum l2_normalize(PG_FUNCTION_ARGS)
{
    FloatVector *a = PG_GETARG_FLOATVECTOR_P(0);
    const float *ax = a->x;
    double norm = -NEGATIVE_INNER_PRODUCT_DIST(ax, ax, a->dim);
    norm = sqrt(norm);

    FloatVector *result = InitFloatVector(a->dim);
    if (norm <= 0) { /* Return zero vector for zero norm */
        PG_RETURN_POINTER(result);
    }

    float *rx = result->x;
    if (norm == 1) {
        memcpy(rx, ax, a->dim * sizeof(float));
        PG_RETURN_POINTER(result);
    }

    for (int16 i = 0; i < a->dim; ++i) {
        rx[i] = ax[i] / norm;
    }
    for (int16 i = 0; i < a->dim; ++i) {
        if (isinf(rx[i])) {
            float_overflow_error();
        }
    }
    PG_RETURN_POINTER(result);
}

void floatvector_normalize(float *ax, int16 dim)
{