use std::{
hash::{DefaultHasher, Hasher},
marker::PhantomData,
};
use std::fmt;
use crate::layouts::{
Backend, BigWord, Data, DataView, DataViewMut, DigestU64, HostDataMut, HostDataRef, VecZnxInfos, VecZnxShape, ZnxInfos,
ZnxView, ZnxViewMut, ZnxZero,
};
#[repr(C)]
pub struct VecZnxBig<D: Data, W: BigWord, B: Backend<BigWord = W>> {
pub data: D,
shape: VecZnxShape,
pub _phantom: PhantomData<(W, B)>,
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> PartialEq for VecZnxBig<D, W, B> {
fn eq(&self, other: &Self) -> bool {
self.shape == other.shape && self.data == other.data
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> Eq for VecZnxBig<D, W, B> {}
impl<D: Data + std::hash::Hash, W: BigWord, B: Backend<BigWord = W>> std::hash::Hash for VecZnxBig<D, W, B> {
fn hash<H: std::hash::Hasher>(&self, state: &mut H) {
self.shape.hash(state);
self.data.hash(state);
}
}
impl<D: HostDataRef, W: BigWord, B: Backend<BigWord = W>> DigestU64 for VecZnxBig<D, W, B> {
fn digest_u64(&self) -> u64 {
let mut h: DefaultHasher = DefaultHasher::new();
h.write(self.data.as_ref());
h.write_usize(self.n());
h.write_usize(self.cols());
h.write_usize(self.size());
h.finish()
}
}
impl<D: HostDataRef, W: BigWord, B: Backend<BigWord = W>> ZnxView for VecZnxBig<D, W, B> {
type Scalar = W;
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> ZnxInfos for VecZnxBig<D, W, B> {
fn n(&self) -> usize {
self.shape.n()
}
fn size(&self) -> usize {
self.shape.size()
}
fn poly_count(&self) -> usize {
crate::layouts::checked_product(&[self.cols(), self.size()], "polynomial count")
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> VecZnxInfos for VecZnxBig<D, W, B> {
fn cols(&self) -> usize {
self.shape.cols()
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> DataView for VecZnxBig<D, W, B> {
type D = D;
fn data(&self) -> &Self::D {
&self.data
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> DataViewMut for VecZnxBig<D, W, B> {
fn data_mut(&mut self) -> &mut Self::D {
&mut self.data
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> VecZnxBig<D, W, B> {
pub fn n(&self) -> usize {
self.shape.n()
}
pub fn cols(&self) -> usize {
self.shape.cols()
}
pub fn size(&self) -> usize {
self.shape.size()
}
pub fn shape(&self) -> VecZnxShape {
self.shape
}
}
impl<D: HostDataMut, W: BigWord, B: Backend<BigWord = W>> ZnxZero for VecZnxBig<D, W, B> {
fn zero(&mut self) {
self.raw_mut().fill(W::zero())
}
fn zero_at(&mut self, i: usize, j: usize) {
self.at_mut(i, j).fill(W::zero());
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> VecZnxBig<D, W, B> {
pub(crate) fn alloc(n: usize, cols: usize, size: usize) -> VecZnxBigOwned<B>
where
B: Backend<OwnedBuf = D>,
{
let data: <B as Backend>::OwnedBuf = B::alloc_zeroed_bytes(B::bytes_of_vec_znx_big(n, cols, size));
VecZnxBig {
data,
shape: VecZnxShape::new(n, cols, size),
_phantom: PhantomData,
}
}
pub fn from_bytes(n: usize, cols: usize, size: usize, bytes: impl Into<Vec<u8>>) -> VecZnxBigOwned<B>
where
B: Backend<OwnedBuf = D>,
{
let data: Vec<u8> = bytes.into();
assert!(data.len() == B::bytes_of_vec_znx_big(n, cols, size));
let data: <B as Backend>::OwnedBuf = B::from_host_bytes(&data);
VecZnxBig {
data,
shape: VecZnxShape::new(n, cols, size),
_phantom: PhantomData,
}
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> VecZnxBig<D, W, B> {
pub fn from_data(data: D, n: usize, cols: usize, size: usize) -> Self {
Self {
data,
shape: VecZnxShape::new(n, cols, size),
_phantom: PhantomData,
}
}
}
pub type VecZnxBigOwned<B> = VecZnxBig<<B as Backend>::OwnedBuf, <B as Backend>::BigWord, B>;
pub type VecZnxBigBackendRef<'a, B> = VecZnxBig<<B as Backend>::BufRef<'a>, <B as Backend>::BigWord, B>;
pub type VecZnxBigBackendMut<'a, B> = VecZnxBig<<B as Backend>::BufMut<'a>, <B as Backend>::BigWord, B>;
pub fn vec_znx_big_backend_ref_from_mut<'a, 'b, B: Backend + 'b>(
vec: &'a VecZnxBigBackendMut<'b, B>,
) -> VecZnxBigBackendRef<'a, B> {
VecZnxBig {
data: B::view_ref_mut(&vec.data),
shape: vec.shape,
_phantom: PhantomData,
}
}
pub trait VecZnxBigToBackendRef<B: Backend> {
fn to_backend_ref(&self) -> VecZnxBigBackendRef<'_, B>;
}
impl<B: Backend> VecZnxBigToBackendRef<B> for VecZnxBig<B::OwnedBuf, B::BigWord, B> {
fn to_backend_ref(&self) -> VecZnxBigBackendRef<'_, B> {
VecZnxBig {
data: B::view(&self.data),
shape: self.shape,
_phantom: std::marker::PhantomData,
}
}
}
impl<'b, B: Backend + 'b> VecZnxBigToBackendRef<B> for &VecZnxBig<B::BufRef<'b>, B::BigWord, B> {
fn to_backend_ref(&self) -> VecZnxBigBackendRef<'_, B> {
VecZnxBig {
data: B::view_ref(&self.data),
shape: self.shape,
_phantom: std::marker::PhantomData,
}
}
}
pub trait VecZnxBigReborrowBackendRef<B: Backend> {
fn reborrow_backend_ref(&self) -> VecZnxBigBackendRef<'_, B>;
}
impl<'b, B: Backend + 'b> VecZnxBigReborrowBackendRef<B> for VecZnxBig<B::BufMut<'b>, B::BigWord, B> {
fn reborrow_backend_ref(&self) -> VecZnxBigBackendRef<'_, B> {
vec_znx_big_backend_ref_from_mut::<B>(self)
}
}
pub trait VecZnxBigToBackendMut<B: Backend> {
fn to_backend_mut(&mut self) -> VecZnxBigBackendMut<'_, B>;
}
impl<B: Backend> VecZnxBigToBackendMut<B> for VecZnxBig<B::OwnedBuf, B::BigWord, B> {
fn to_backend_mut(&mut self) -> VecZnxBigBackendMut<'_, B> {
VecZnxBig {
data: B::view_mut(&mut self.data),
shape: self.shape,
_phantom: std::marker::PhantomData,
}
}
}
impl<'b, B: Backend + 'b> VecZnxBigToBackendMut<B> for &mut VecZnxBig<B::BufMut<'b>, B::BigWord, B> {
fn to_backend_mut(&mut self) -> VecZnxBigBackendMut<'_, B> {
VecZnxBig {
data: B::view_mut_ref(&mut self.data),
shape: self.shape,
_phantom: std::marker::PhantomData,
}
}
}
pub trait VecZnxBigReborrowBackendMut<B: Backend> {
fn reborrow_backend_mut(&mut self) -> VecZnxBigBackendMut<'_, B>;
}
impl<'b, B: Backend + 'b> VecZnxBigReborrowBackendMut<B> for VecZnxBig<B::BufMut<'b>, B::BigWord, B> {
fn reborrow_backend_mut(&mut self) -> VecZnxBigBackendMut<'_, B> {
VecZnxBig {
data: B::view_mut_ref(&mut self.data),
shape: self.shape,
_phantom: std::marker::PhantomData,
}
}
}
impl<D: HostDataRef, W: BigWord, B: Backend<BigWord = W>> fmt::Display for VecZnxBig<D, W, B> {
fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
writeln!(f, "VecZnxBig(n={}, cols={}, size={})", self.n(), self.cols(), self.size())?;
for col in 0..self.cols() {
writeln!(f, "Column {col}:")?;
for size in 0..self.size() {
let coeffs = self.at(col, size);
write!(f, " Size {size}: [")?;
let max_show = 100;
let show_count = coeffs.len().min(max_show);
for (i, &coeff) in coeffs.iter().take(show_count).enumerate() {
if i > 0 {
write!(f, ", ")?;
}
write!(f, "{coeff}")?;
}
if coeffs.len() > max_show {
write!(f, ", ... ({} more)", coeffs.len() - max_show)?;
}
writeln!(f, "]")?;
}
}
Ok(())
}
}
impl<D: Data, W: BigWord, B: Backend<BigWord = W>> VecZnxBig<D, W, B> {
pub fn into_backend<B2>(self) -> VecZnxBig<D, W, B2>
where
B2: Backend<BigWord = W>,
B: crate::layouts::VecZnxBigLayoutCompatible<B2>,
{
let shape = self.shape;
assert_eq!(
B::bytes_of_vec_znx_big(shape.n(), shape.cols(), shape.size()),
B2::bytes_of_vec_znx_big(shape.n(), shape.cols(), shape.size()),
"into_backend: byte sizes diverge despite declared layout compatibility"
);
VecZnxBig {
data: self.data,
shape,
_phantom: PhantomData,
}
}
}