use crate::bindings as ll_bindings;
use crate::metadata;
use crate::sys;
use crate::Position;
use crate::SizeType;
use crate::Time;
use crate::{tsk_id_t, TskitError};
use crate::{MigrationId, NodeId, PopulationId};
#[derive(Debug)]
pub struct MigrationTableRow {
pub id: MigrationId,
pub left: Position,
pub right: Position,
pub node: NodeId,
pub source: PopulationId,
pub dest: PopulationId,
pub time: Time,
pub metadata: Option<Vec<u8>>,
}
impl PartialEq for MigrationTableRow {
fn eq(&self, other: &Self) -> bool {
self.id == other.id
&& self.node == other.node
&& self.source == other.source
&& self.dest == other.dest
&& crate::util::partial_cmp_equal(&self.left, &other.left)
&& crate::util::partial_cmp_equal(&self.right, &other.right)
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& self.metadata == other.metadata
}
}
fn make_migration_table_row(table: &MigrationTable, pos: tsk_id_t) -> Option<MigrationTableRow> {
Some(MigrationTableRow {
id: pos.into(),
left: table.left(pos)?,
right: table.right(pos)?,
node: table.node(pos)?,
source: table.source(pos)?,
dest: table.dest(pos)?,
time: table.time(pos)?,
metadata: table.raw_metadata(pos).map(|m| m.to_vec()),
})
}
pub(crate) type MigrationTableRefIterator<'a> =
crate::table_iterator::TableIterator<&'a MigrationTable>;
pub(crate) type MigrationTableIterator = crate::table_iterator::TableIterator<MigrationTable>;
impl<'a> Iterator for MigrationTableRefIterator<'a> {
type Item = MigrationTableRow;
fn next(&mut self) -> Option<Self::Item> {
let rv = make_migration_table_row(self.table, self.pos);
self.pos += 1;
rv
}
}
impl Iterator for MigrationTableIterator {
type Item = crate::migration_table::MigrationTableRow;
fn next(&mut self) -> Option<Self::Item> {
let rv = make_migration_table_row(&self.table, self.pos);
self.pos += 1;
rv
}
}
#[derive(Debug)]
pub struct MigrationTableRowView<'a> {
table: &'a MigrationTable,
pub id: MigrationId,
pub left: Position,
pub right: Position,
pub node: NodeId,
pub source: PopulationId,
pub dest: PopulationId,
pub time: Time,
pub metadata: Option<&'a [u8]>,
}
impl<'a> MigrationTableRowView<'a> {
fn new(table: &'a MigrationTable) -> Self {
Self {
table,
id: MigrationId::NULL,
left: Position::from(f64::NAN),
right: Position::from(f64::NAN),
node: NodeId::NULL,
source: PopulationId::NULL,
dest: PopulationId::NULL,
time: Time::from(f64::NAN),
metadata: None,
}
}
}
impl<'a> PartialEq for MigrationTableRowView<'a> {
fn eq(&self, other: &Self) -> bool {
self.id == other.id
&& self.node == other.node
&& self.source == other.source
&& self.dest == other.dest
&& crate::util::partial_cmp_equal(&self.left, &other.left)
&& crate::util::partial_cmp_equal(&self.right, &other.right)
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& self.metadata == other.metadata
}
}
impl<'a> Eq for MigrationTableRowView<'a> {}
impl<'a> PartialEq<MigrationTableRow> for MigrationTableRowView<'a> {
fn eq(&self, other: &MigrationTableRow) -> bool {
self.id == other.id
&& self.node == other.node
&& self.source == other.source
&& self.dest == other.dest
&& crate::util::partial_cmp_equal(&self.left, &other.left)
&& crate::util::partial_cmp_equal(&self.right, &other.right)
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& optional_container_comparison!(self.metadata, other.metadata)
}
}
impl PartialEq<MigrationTableRowView<'_>> for MigrationTableRow {
fn eq(&self, other: &MigrationTableRowView) -> bool {
self.id == other.id
&& self.node == other.node
&& self.source == other.source
&& self.dest == other.dest
&& crate::util::partial_cmp_equal(&self.left, &other.left)
&& crate::util::partial_cmp_equal(&self.right, &other.right)
&& crate::util::partial_cmp_equal(&self.time, &other.time)
&& optional_container_comparison!(self.metadata, other.metadata)
}
}
impl<'a> streaming_iterator::StreamingIterator for MigrationTableRowView<'a> {
type Item = Self;
row_lending_iterator_get!();
fn advance(&mut self) {
self.id = (i32::from(self.id) + 1).into();
self.left = self.table.left(self.id).unwrap_or_else(|| f64::NAN.into());
self.right = self.table.right(self.id).unwrap_or_else(|| f64::NAN.into());
self.node = self.table.node(self.id).unwrap_or(NodeId::NULL);
self.source = self.table.source(self.id).unwrap_or(PopulationId::NULL);
self.dest = self.table.dest(self.id).unwrap_or(PopulationId::NULL);
self.time = self.table.time(self.id).unwrap_or_else(|| f64::NAN.into());
self.metadata = self.table.raw_metadata(self.id);
}
}
#[derive(Debug)]
#[repr(transparent)]
pub struct MigrationTable {
table_: sys::LLMigrationTableRef,
}
impl MigrationTable {
pub(crate) fn new_from_table(
migrations: *mut ll_bindings::tsk_migration_table_t,
) -> Result<Self, TskitError> {
let table_ = sys::LLMigrationTableRef::new_from_table(migrations)?;
Ok(MigrationTable { table_ })
}
pub(crate) fn as_ref(&self) -> &ll_bindings::tsk_migration_table_t {
self.table_.as_ref()
}
pub fn num_rows(&self) -> SizeType {
self.as_ref().num_rows.into()
}
raw_metadata_getter_for_tables!(MigrationId);
pub fn left<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<Position> {
sys::tsk_column_access::<Position, _, _, _>(row.into(), self.as_ref().left, self.num_rows())
}
pub fn right<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<Position> {
sys::tsk_column_access::<Position, _, _, _>(
row.into(),
self.as_ref().right,
self.num_rows(),
)
}
pub fn node<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<NodeId> {
sys::tsk_column_access::<NodeId, _, _, _>(row.into(), self.as_ref().node, self.num_rows())
}
pub fn source<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<PopulationId> {
sys::tsk_column_access::<PopulationId, _, _, _>(
row.into(),
self.as_ref().source,
self.num_rows(),
)
}
pub fn dest<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<PopulationId> {
sys::tsk_column_access::<PopulationId, _, _, _>(
row.into(),
self.as_ref().dest,
self.num_rows(),
)
}
pub fn time<M: Into<MigrationId> + Copy>(&self, row: M) -> Option<Time> {
sys::tsk_column_access::<Time, _, _, _>(row.into(), self.as_ref().time, self.num_rows())
}
pub fn metadata<T: metadata::MigrationMetadata>(
&self,
row: MigrationId,
) -> Option<Result<T, TskitError>> {
let buffer = self.raw_metadata(row)?;
Some(decode_metadata_row!(T, buffer).map_err(|e| e.into()))
}
pub fn iter(&self) -> impl Iterator<Item = MigrationTableRow> + '_ {
crate::table_iterator::make_table_iterator::<&MigrationTable>(self)
}
pub fn lending_iter(&self) -> MigrationTableRowView {
MigrationTableRowView::new(self)
}
pub fn row<M: Into<MigrationId> + Copy>(&self, r: M) -> Option<MigrationTableRow> {
let ri = r.into().into();
table_row_access!(ri, self, make_migration_table_row)
}
pub fn row_view<M: Into<MigrationId> + Copy>(&self, r: M) -> Option<MigrationTableRowView> {
let view = MigrationTableRowView {
table: self,
id: r.into(),
left: self.left(r)?,
right: self.right(r)?,
node: self.node(r)?,
source: self.source(r)?,
dest: self.dest(r)?,
time: self.time(r)?,
metadata: self.raw_metadata(r.into()),
};
Some(view)
}
build_table_column_slice_getter!(
=> left, left_slice, Position);
build_table_column_slice_getter!(
=> left, left_slice_raw, f64);
build_table_column_slice_getter!(
=> right, right_slice, Position);
build_table_column_slice_getter!(
=> right, right_slice_raw, f64);
build_table_column_slice_getter!(
=> time, time_slice, Time);
build_table_column_slice_getter!(
=> time, time_slice_raw, f64);
build_table_column_slice_getter!(
=> node, node_slice, NodeId);
build_table_column_slice_getter!(
=> node, node_slice_raw, ll_bindings::tsk_id_t);
build_table_column_slice_getter!(
=> source, source_slice, PopulationId);
build_table_column_slice_getter!(
=> source, source_slice_raw, ll_bindings::tsk_id_t);
build_table_column_slice_getter!(
=> dest, dest_slice, PopulationId);
build_table_column_slice_getter!(
=> dest, dest_slice_raw, ll_bindings::tsk_id_t);
}
build_owned_table_type!(
=> OwningMigrationTable,
MigrationTable,
crate::sys::LLOwningMigrationTable,
crate::bindings::tsk_migration_table_t
);
impl OwningMigrationTable {
migration_table_add_row!(=> add_row, self, self.as_mut_ptr());
migration_table_add_row_with_metadata!(=> add_row_with_metadata, self, self.as_mut_ptr());
}