use std::fmt;
use bytemuck::{Pod, Zeroable};
use crate::distance::Distance;
use crate::node::Item;
use crate::spaces::simple::dot_product;
use crate::unaligned_vector::UnalignedVector;
#[derive(Debug, Clone)]
pub enum Cosine {}
#[repr(C)]
#[derive(Pod, Zeroable, Clone, Copy)]
pub struct NodeHeaderCosine {
norm: f32,
}
impl fmt::Debug for NodeHeaderCosine {
fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
f.debug_struct("NodeHeaderCosine").field("norm", &format!("{:.4}", self.norm)).finish()
}
}
impl Distance for Cosine {
type Header = NodeHeaderCosine;
type VectorCodec = f32;
fn name() -> &'static str {
"cosine"
}
fn new_header(vector: &UnalignedVector<Self::VectorCodec>) -> Self::Header {
NodeHeaderCosine { norm: Self::norm_no_header(vector) }
}
fn distance(p: &Item<Self>, q: &Item<Self>) -> f32 {
let pn = p.header.norm;
let qn = q.header.norm;
let pq = dot_product(&p.vector, &q.vector);
let pnqn = pn * qn;
if pnqn > f32::EPSILON {
let cos = pq / pnqn;
let cos = cos.clamp(-1.0, 1.0);
(1.0 - cos) / 2.0
} else {
0.0
}
}
fn norm_no_header(v: &UnalignedVector<Self::VectorCodec>) -> f32 {
dot_product(v, v).sqrt()
}
}