use crate::metadata;
use crate::metadata::MutationMetadata;
use crate::sys;
use crate::SizeType;
use crate::Time;
use crate::TskitError;
use crate::{MutationId, NodeId, SiteId};
use ll_bindings::tsk_id_t;
use sys::bindings as ll_bindings;
#[derive(Debug)]
pub struct MutationTableRow {
pub id: MutationId,
pub site: SiteId,
pub node: NodeId,
pub parent: MutationId,
pub time: Time,
pub derived_state: Option<Vec<u8>>,
pub metadata: Option<Vec<u8>>,
}
impl PartialEq for MutationTableRow {
fn eq(&self, other: &Self) -> bool {
self.id == other.id
&& self.site == other.site
&& self.node == other.node
&& self.parent == other.parent
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& self.derived_state == other.derived_state
&& self.metadata == other.metadata
}
}
fn make_mutation_table_row(table: &MutationTable, pos: tsk_id_t) -> Option<MutationTableRow> {
let index = ll_bindings::tsk_size_t::try_from(pos).ok()?;
match index {
i if i < table.num_rows() => {
let derived_state = table.derived_state(pos).map(|s| s.to_vec());
Some(MutationTableRow {
id: pos.into(),
site: table.site(pos)?,
node: table.node(pos)?,
parent: table.parent(pos)?,
time: table.time(pos)?,
derived_state,
metadata: table.table_.raw_metadata(pos).map(|m| m.to_vec()),
})
}
_ => None,
}
}
pub(crate) type MutationTableRefIterator<'a> =
crate::table_iterator::TableIterator<&'a MutationTable>;
pub(crate) type MutationTableIterator = crate::table_iterator::TableIterator<MutationTable>;
impl Iterator for MutationTableRefIterator<'_> {
type Item = MutationTableRow;
fn next(&mut self) -> Option<Self::Item> {
let rv = make_mutation_table_row(self.table, self.pos);
self.pos += 1;
rv
}
}
impl Iterator for MutationTableIterator {
type Item = MutationTableRow;
fn next(&mut self) -> Option<Self::Item> {
let rv = make_mutation_table_row(&self.table, self.pos);
self.pos += 1;
rv
}
}
#[derive(Debug)]
pub struct MutationTableRowView<'a> {
table: &'a MutationTable,
pub id: MutationId,
pub site: SiteId,
pub node: NodeId,
pub parent: MutationId,
pub time: Time,
pub derived_state: Option<&'a [u8]>,
pub metadata: Option<&'a [u8]>,
}
impl<'a> MutationTableRowView<'a> {
fn new(table: &'a MutationTable) -> Self {
Self {
table,
id: MutationId::NULL,
site: SiteId::NULL,
node: NodeId::NULL,
parent: MutationId::NULL,
time: f64::NAN.into(),
derived_state: None,
metadata: None,
}
}
}
impl PartialEq for MutationTableRowView<'_> {
fn eq(&self, other: &Self) -> bool {
self.id == other.id
&& self.site == other.site
&& self.node == other.node
&& self.parent == other.parent
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& self.derived_state == other.derived_state
&& self.metadata == other.metadata
}
}
impl Eq for MutationTableRowView<'_> {}
impl PartialEq<MutationTableRow> for MutationTableRowView<'_> {
fn eq(&self, other: &MutationTableRow) -> bool {
self.id == other.id
&& self.site == other.site
&& self.node == other.node
&& self.parent == other.parent
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& optional_container_comparison!(self.derived_state, other.derived_state)
&& optional_container_comparison!(self.metadata, other.metadata)
}
}
impl PartialEq<MutationTableRowView<'_>> for MutationTableRow {
fn eq(&self, other: &MutationTableRowView) -> bool {
self.id == other.id
&& self.site == other.site
&& self.node == other.node
&& self.parent == other.parent
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& optional_container_comparison!(self.derived_state, other.derived_state)
&& optional_container_comparison!(self.metadata, other.metadata)
}
}
impl crate::StreamingIterator for MutationTableRowView<'_> {
type Item = Self;
row_lending_iterator_get!();
fn advance(&mut self) {
self.id = (i32::from(self.id) + 1).into();
self.site = self.table.site(self.id).unwrap_or(SiteId::NULL);
self.node = self.table.node(self.id).unwrap_or(NodeId::NULL);
self.parent = self.table.parent(self.id).unwrap_or(MutationId::NULL);
self.time = self.table.time(self.id).unwrap_or_else(|| f64::NAN.into());
self.derived_state = self.table.derived_state(self.id);
self.metadata = self.table.table_.raw_metadata(self.id);
}
}
#[derive(Debug, Default)]
#[repr(transparent)]
pub struct MutationTable {
table_: sys::MutationTable,
}
impl MutationTable {
pub(crate) unsafe fn new_from_table(
mutations: *mut ll_bindings::tsk_mutation_table_t,
) -> Result<Self, TskitError> {
let ptr = std::ptr::NonNull::new(mutations).unwrap();
let table_ = unsafe { sys::MutationTable::new_borrowed(ptr) };
Ok(MutationTable { table_ })
}
pub(crate) fn as_ref(&self) -> &ll_bindings::tsk_mutation_table_t {
self.table_.as_ref()
}
pub fn num_rows(&self) -> SizeType {
self.as_ref().num_rows.into()
}
pub fn site<M: Into<MutationId> + Copy>(&self, row: M) -> Option<SiteId> {
self.table_.site(row.into())
}
pub fn node<M: Into<MutationId> + Copy>(&self, row: M) -> Option<NodeId> {
self.table_.node(row.into())
}
pub fn parent<M: Into<MutationId> + Copy>(&self, row: M) -> Option<MutationId> {
self.table_.parent(row.into())
}
pub fn time<M: Into<MutationId> + Copy>(&self, row: M) -> Option<Time> {
self.table_.time(row.into())
}
pub fn derived_state<M: Into<MutationId>>(&self, row: M) -> Option<&[u8]> {
self.table_.derived_state(row.into())
}
pub fn metadata<T: metadata::MutationMetadata>(
&self,
row: impl Into<MutationId>,
) -> Option<Result<T, TskitError>> {
let buffer = self.table_.raw_metadata(row)?;
Some(decode_metadata_row!(T, buffer).map_err(|e| e.into()))
}
pub fn iter(&self) -> impl Iterator<Item = MutationTableRow> + '_ {
crate::table_iterator::make_table_iterator::<&MutationTable>(self)
}
pub fn lending_iter(&'_ self) -> MutationTableRowView<'_> {
MutationTableRowView::new(self)
}
pub fn row<M: Into<MutationId> + Copy>(&self, r: M) -> Option<MutationTableRow> {
let ri = r.into().into();
table_row_access!(ri, self, make_mutation_table_row)
}
pub fn row_view<M: Into<MutationId> + Copy>(
&'_ self,
r: M,
) -> Option<MutationTableRowView<'_>> {
let view = MutationTableRowView {
table: self,
id: r.into(),
site: self.site(r)?,
node: self.node(r)?,
parent: self.parent(r)?,
time: self.time(r)?,
derived_state: self.derived_state(r),
metadata: self.table_.raw_metadata(r.into()),
};
Some(view)
}
build_table_column_slice_getter!(
=> node, node_slice, NodeId);
build_table_column_slice_getter!(
=> node, node_slice_raw, crate::sys::bindings::tsk_id_t);
build_table_column_slice_getter!(
=> site, site_slice, SiteId);
build_table_column_slice_getter!(
=> site, site_slice_raw, crate::sys::bindings::tsk_id_t);
build_table_column_slice_getter!(
=> time, time_slice, Time);
build_table_column_slice_getter!(
=> time, time_slice_raw, f64);
build_table_column_slice_getter!(
=> parent, parent_slice, MutationId);
build_table_column_slice_getter!(
=> parent, parent_slice_raw, crate::sys::bindings::tsk_id_t);
pub fn clear(&mut self) -> Result<i32, TskitError> {
handle_tsk_return_value!(self.table_.clear())
}
pub fn add_row<S, N, P, T>(
&mut self,
site: S,
node: N,
parent: P,
time: T,
derived_state: Option<&[u8]>,
) -> Result<MutationId, TskitError>
where
S: Into<SiteId>,
N: Into<NodeId>,
P: Into<MutationId>,
T: Into<Time>,
{
let rv = self.table_.add_row(
site.into().into(),
node.into().into(),
parent.into().into(),
time.into().into(),
derived_state,
)?;
handle_tsk_return_value!(rv, rv.into())
}
pub fn add_row_with_metadata<S, N, P, T, M>(
&mut self,
site: S,
node: N,
parent: P,
time: T,
derived_state: Option<&[u8]>,
metadata: &M,
) -> Result<MutationId, TskitError>
where
S: Into<SiteId>,
N: Into<NodeId>,
P: Into<MutationId>,
T: Into<Time>,
M: MutationMetadata,
{
let md = crate::metadata::EncodedMetadata::new(metadata)?;
let rv = self.table_.add_row_with_metadata(
site.into().into(),
node.into().into(),
parent.into().into(),
time.into().into(),
derived_state,
md.as_slice(),
)?;
handle_tsk_return_value!(rv, rv.into())
}
}