inertia-core 0.4.0

/*
 *  Copyright (C) 2021 William Youmans
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */

// TODO: op guards need work, especially *From and Assign* (DivFrom, AssignDiv etc)
//! Macros for implementing comparisons, operations, and conversions.

macro_rules! default {
    // Unary ops and From
    ($op:ident, ctx, $out_ty:ident, $in:ident) => {
        $out_ty::zero($in.context())
    };
    ($op:ident, ctx_in, $out_ty:ident, $in:ident) => {
        $out_ty::zero()
    };
    ($op:ident, matrix, $out_ty:ident, $in:ident) => {
        $out_ty::zero($in.nrows_si(), $in.ncols_si())
    };
    ($op:ident, matrix_ctx, $out_ty:ident, $in:ident) => {
        $out_ty::zero($in.nrows_si(), $in.ncols_si(), $in.context())
    };
    ($op:ident, matrix_ctx_new_only, $out_ty:ident, $in:ident) => {
        $out_ty::zero($in.nrows_si(), $in.ncols_si(), $in.context())
    };
    ($op:ident, $kw:ident, $out_ty:ident, $in:ident) => {
        $out_ty::default()
    };

    // Binary ops
    ($op:ident, ctx, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.context())
    };
    ($op:ident, ctx_lhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.context())
    };
    ($op:ident, ctx_rhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($rhs.context())
    };
    ($op:ident, matrix, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $rhs.ncols_si())
    };
    ($op:ident, matrix_ctx, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $rhs.ncols_si(), $lhs.context())
    };
    /*
    ($op:ident, matrix_ctx_lhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $rhs.ncols_si(), $lhs.context())
    };
    ($op:ident, matrix_ctx_rhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $rhs.ncols_si(), $rhs.context())
    };
    */
    ($op:ident, matrix_ctx_new_only, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $rhs.ncols_si(), $lhs.context())
    };
    ($op:ident, scalar_lhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($rhs.nrows_si(), $rhs.ncols_si())
    };
    ($op:ident, scalar_rhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::zero($lhs.nrows_si(), $lhs.ncols_si())
    };
    ($op:ident, $kw:ident, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        $out_ty::default()
    };
}

macro_rules! op_guard {
    // Unary ops and From
    (Inv, $kw:ident, $out_ty:ident, $in:ident) => {
        // TODO: check den != 0
        //$in.is_invertible()
    };
    (Pow, $kw:ident, $out_ty:ident, $in:ident) => {
        // TODO: check if exp negative or fractional for certain types
    };
    ($op:ident, $kw:ident, $out_ty:ident, $in:ident) => {
    };

    // Binary ops
    ($op:ident, ctx, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        // check contexts agree
        assert_eq!($lhs.context(), $rhs.context())
    };
    (Mul, matrix, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        assert_eq!($lhs.ncols_si(), $rhs.nrows_si())
    };
    (MulAssign, matrix, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        assert_eq!($lhs.ncols_si(), $rhs.nrows_si())
    };
    // covers all matrix ops except Mul
    ($op:ident, matrix, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        assert_eq!($lhs.nrows_si(), $rhs.nrows_si());
        assert_eq!($lhs.ncols_si(), $rhs.ncols_si())
    };
    ($op:ident, scalar_lhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
    };
    (Div, scalar_rhs, $out_ty:ident, $lhs:ident, $rhs:ident) => {
        // TODO:
        // $rhs.is_invertible()
    };
    ($op:ident, $kw:ident, $out_ty:ident, $lhs:ident, $rhs:ident) => {
    };
}

macro_rules! call_unsafe {
    // Unary ops and From
    (ctx, $func:path, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $in.as_ptr(), $out.ctx_as_ptr());
        }
    };
    (matrix_ctx, $func:path, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $in.as_ptr(), $out.ctx_as_ptr());
        }
    };
    (ctx_in, $func:path, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $in.as_ptr(), $in.ctx_as_ptr());
        }
    };
    ($kw:ident, $func:path, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $in.as_ptr());
        }
    };
    (cast ctx, $func:path, $cast:ty, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), *$in as $cast, $out.ctx_as_ptr());
        }
    };
    (cast $kw:ident, $func:path, $cast:ty, $out:ident, $in:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), *$in as $cast);
        }
    };

    // Binary ops
    (ctx, $func:path, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                $lhs.as_ptr(),
                $rhs.as_ptr(),
                $lhs.ctx_as_ptr(),
            );
        }
    };
    (ctx_lhs, $func:path, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                $lhs.as_ptr(),
                $rhs.as_ptr(),
                $lhs.ctx_as_ptr(),
            );
        }
    };
    (ctx_rhs, $func:path, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                $lhs.as_ptr(),
                $rhs.as_ptr(),
                $rhs.ctx_as_ptr(),
            );
        }
    };
    (matrix_ctx, $func:path, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                $lhs.as_ptr(),
                $rhs.as_ptr(),
                $lhs.ctx_as_ptr(),
            );
        }
    };
    ($kw:ident, $func:path, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $lhs.as_ptr(), $rhs.as_ptr());
        }
    };

    // Binary ops with primitive types
    (cast_rhs ctx_lhs, $func:path, $cast:ty, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                $lhs.as_ptr(),
                *$rhs as $cast,
                $lhs.ctx_as_ptr(),
            );
        }
    };
    (cast_rhs $kw:ident, $func:path, $cast:ty, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), $lhs.as_ptr(), *$rhs as $cast);
        }
    };
    (cast_lhs ctx_rhs, $func:path, $cast:ty, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func(
                $out.as_mut_ptr(),
                *$lhs as $cast,
                $rhs.as_ptr(),
                $rhs.ctx_as_ptr(),
            );
        }
    };
    (cast_lhs $kw:ident, $func:path, $cast:ty, $out:ident, $lhs:ident, $rhs:ident) => {
        unsafe {
            $func($out.as_mut_ptr(), *$lhs as $cast, $rhs.as_ptr());
        }
    };
}

/// Macros for overloading comparison operators
#[macro_export]
macro_rules! impl_cmp {
    // a = a
    (
        eq
        $t:ident
        {
            $($code:tt)*
        }
    ) => {
        impl Eq for $t {}

        //impl Eq for &$t {}

        impl PartialEq for $t {
            #[inline]
            $($code)*
        }

        impl PartialEq<&$t> for $t {
            #[inline]
            fn eq(&self, rhs: &&$t) -> bool {
                rhs.eq(&self)
            }
        }

        impl PartialEq<$t> for &$t {
            #[inline]
            fn eq(&self, rhs: &$t) -> bool {
                self.eq(&rhs)
            }
        }
    };
    (
        partial_eq
        $t:ident
        {
            $($code:tt)*
        }
    ) => {
        impl_cmp! {
            partial_eq
            $t, $t
            {
                $($code)*
            }
        }
    };
    // a > a
    (
        ord
        $t:ident
        {
            $($code:tt)*
        }
    ) => {
        impl Ord for $t {
            #[inline]
            $($code)*
        }

        impl PartialOrd for $t {
            #[inline]
            fn partial_cmp(&self, rhs: &$t) -> Option<Ordering> {
                Some(self.cmp(rhs))
            }
        }
    };
    (
        partial_ord
        $t:ident
        {
            $($code:tt)*
        }
    ) => {
        impl_cmp! {
            partial_ord
            $t, $t
            {
                $($code)*
            }
        }
    };
    // a = b
    (
        partial_eq
        $t1:ident, $t2:ident
        {
            $($code:tt)*
        }
    ) => {
        impl PartialEq<$t2> for $t1 {
            #[inline]
            $($code)*
        }

        impl PartialEq<&$t2> for $t1 {
            #[inline]
            fn eq(&self, rhs: &&$t2) -> bool {
                (&self).eq(rhs)
            }
        }

        impl PartialEq<$t2> for &$t1 {
            #[inline]
            fn eq(&self, rhs: &$t2) -> bool {
                self.eq(&rhs)
            }
        }
    };
    // a > b
    (
        partial_ord
        $t1:ident, $t2:ident
        {
            $($code:tt)*
        }
    ) => {
        impl PartialOrd<$t2> for $t1 {
            #[inline]
            $($code)*
        }
    };

}

/// Macros for overloading comparison operators with unsafe functions.
macro_rules! impl_cmp_unsafe {
    (
        eq
        $t:ident
        $func:path
    ) => {
        impl_cmp! {
            eq
            $t
            {
                fn eq(&self, rhs: &$t) -> bool {
                    unsafe { $func(self.as_ptr(), rhs.as_ptr()) != 0 }
                }
            }
        }
    };
    (
        partial_eq
        $t1:ident, $t2:ident
        $func:path
    ) => {
        impl_cmp! {
            partial_eq
            $t1, $t2
            {
                fn eq(&self, rhs: &$t2) -> bool {
                    unsafe { $func(self.as_ptr(), rhs.as_ptr()) != 0 }
                }
            }
        }
        impl_cmp! {
            partial_eq
            $t2, $t1
            {
                fn eq(&self, rhs: &$t1) -> bool {
                    unsafe { $func(rhs.as_ptr(), self.as_ptr()) != 0 }
                }
            }
        }
    };
    (
        partial_eq
        $t1:ident, $cast:ident {$($t2:ident)+}
        $func:path
    ) => ($(
        impl_cmp! {
            partial_eq
            $t1, $t2
            {
                fn eq(&self, rhs: &$t2) -> bool {
                    unsafe { $func(self.as_ptr(), *rhs as $cast) != 0 }
                }
            }
        }
        impl_cmp! {
            partial_eq
            $t2, $t1
            {
                fn eq(&self, rhs: &$t1) -> bool {
                    unsafe { $func(rhs.as_ptr(), *self as $cast) != 0 }
                }
            }
        }
    )+);
    (
        ord
        $t:ident
        $func:path
    ) => {
        impl_cmp! {
            ord
            $t
            {
                fn cmp(&self, rhs: &$t) -> Ordering {
                    let cmp = unsafe { $func(self.as_ptr(), rhs.as_ptr()) };
                    if cmp == 0 {
                        Equal
                    } else if cmp < 0 {
                        Less
                    } else {
                        Greater
                    }
                }
            }
        }
    };
    (
        partial_ord
        $t1:ident, $t2:ident
        $func:path
    ) => {
        impl_cmp! {
            partial_ord
            $t1, $t2
            {
                fn partial_cmp(&self, rhs: &$t2) -> Option<Ordering> {
                    let cmp = unsafe { $func(self.as_ptr(), rhs.as_ptr()) };
                    if cmp == 0 {
                        Some(Equal)
                    } else if cmp < 0 {
                        Some(Less)
                    } else {
                        Some(Greater)
                    }
                }
            }
        }
        impl_cmp! {
            partial_ord
            $t2, $t1
            {
                fn partial_cmp(&self, rhs: &$t1) -> Option<Ordering> {
                    let cmp = unsafe { $func(rhs.as_ptr(), self.as_ptr()) };
                    if cmp == 0 {
                        Some(Equal)
                    } else if cmp > 0 {
                        Some(Less)
                    } else {
                        Some(Greater)
                    }
                }
            }
        }
    };
    (
        partial_ord
        $t1:ident, $cast:ident {$($t2:ident)+}
        $func:path
    ) => ($(
        impl_cmp! {
            partial_ord
            $t1, $t2
            {
                fn partial_cmp(&self, rhs: &$t2) -> Option<Ordering> {
                    let cmp = unsafe { $func(self.as_ptr(), *rhs as $cast) };
                    if cmp == 0 {
                        Some(Equal)
                    } else if cmp < 0 {
                        Some(Less)
                    } else {
                        Some(Greater)
                    }
                }
            }
        }
        impl_cmp! {
            partial_ord
            $t2, $t1
            {
                fn partial_cmp(&self, rhs: &$t1) -> Option<Ordering> {
                    let cmp = unsafe { $func(rhs.as_ptr(), *self as $cast) };
                    if cmp == 0 {
                        Some(Equal)
                    } else if cmp > 0 {
                        Some(Less)
                    } else {
                        Some(Greater)
                    }
                }
            }
        }
    )+)
}

/// Macros for overloading unary operators.
#[macro_export]
macro_rules! impl_unop {
    (
        // assign
        $t:ident
        $op:ident {$meth:ident}
        {
            $($code:tt)*
        }
        $op_assign:ident {$meth_assign:ident}
        {
            $($code_assign:tt)*
        }
    ) => {
        impl $op for $t {
            type Output = $t;
            #[inline]
            fn $meth(mut self) -> $t {
                self.$meth_assign();
                self
            }
        }

        impl $op for &$t {
            type Output = $t;
            #[inline]
            $($code)*
        }

        impl $op_assign for $t {
            #[inline]
            $($code_assign)*
        }
    };
    (
        // no assign
        $t:ident, $out:ident
        $op:ident {$meth:ident}
        {
            $($code:tt)*
        }
    ) => {
        impl $op for $t {
            type Output = $out;
            #[inline]
            $($code)*
        }

        impl $op for &$t {
            type Output = $out;
            #[inline]
            $($code)*
        }
    };
    (
        // no assign
        $t:ident, Option<$out:ident>
        $op:ident {$meth:ident}
        {
            $($code:tt)*
        }
    ) => {
        impl $op for $t {
            type Output = Option<$out>;
            #[inline]
            $($code)*
        }

        impl $op for &$t {
            type Output = Option<$out>;
            #[inline]
            $($code)*
        }
    };
}

/// Macros for overloading unary operators with unsafe functions.
macro_rules! impl_unop_unsafe {
    (
        $kw:ident
        $t:ident
        $op:ident {$meth:ident}
        $op_assign:ident {$meth_assign:ident}
        $func:path
    ) => {
        impl_unop! {
            $t
            $op {$meth}
            {
                fn $meth(self) -> $t {
                    let mut res = default!($op, $kw, $t, self);
                    call_unsafe!($kw, $func, res, self);
                    res
                }
            }
            $op_assign {$meth_assign}
            {
                fn $meth_assign(&mut self) {
                    call_unsafe!($kw, $func, self, self);
                }
            }
        }
    };
    (
        $kw:ident
        $t:ident, $out:ident
        $op:ident {$meth:ident}
        $func:path
    ) => {
        impl_unop! {
            $t, $out
            $op {$meth}
            {
                fn $meth(self) -> $out {
                    let mut res = default!($op, $kw, $out, self);
                    //unsafe { $func(res.as_mut_ptr(), self.as_ptr()); }
                    call_unsafe!($kw, $func, res, self);
                    res
                }
            }
        }
    };
}

/// Macros for overloading binary operators.
#[macro_export]
macro_rules! impl_binop {
    (
        // a + a = a
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            {
                $($code:tt)*
            }
            $op_assign:ident {$meth_assign:ident}
            {
                $($code_assign:tt)*
            }
            $op_from:ident {$meth_from:ident}
            {
                $($code_from:tt)*
            }
            $assign_op:ident {$assign_meth:ident}
            {
                $($assign_code:tt)*
            }
        )*
    ) => ($(

        impl $op<&$t2> for &$t1 {
            type Output = $out;
            #[inline]
            $($code)*
        }
       
        impl $op<$t2> for &$t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, mut rhs: $t2) -> $out {
                rhs.$meth_from(self);
                rhs
            }
        }

        impl $op<&$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(mut self, rhs: &$t2) -> $out {
                self.$meth_assign(rhs);
                self
            }
        }

        impl $op<$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(mut self, rhs: $t2) -> $out {
                self.$meth_assign(&rhs);
                self
            }
        }

        impl_binop! {@op_assign
            $t1, $t2, $out
            $op_assign {$meth_assign}
            {
                $($code_assign)*
            }
        }

        impl_binop! {@op_from
            $t1, $t2, $out
            $op_from {$meth_from}
            {
                $($code_from)*
            }
        }

        impl_binop! {@assign_op
            $t1, $t2, $out
            $assign_op {$assign_meth}
            {
                $($assign_code)*
            }
        }
    )*);
    (
        // a + b = a
        op_assign
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            {
                $($code:tt)*
            }
            $op_assign:ident {$meth_assign:ident}
            {
                $($code_assign:tt)*
            }
            $assign_op:ident {$assign_meth:ident}
            {
                $($assign_code:tt)*
            }
        )*
    ) => ($(

        impl $op<&$t2> for &$t1 {
            type Output = $out;
            #[inline]
            $($code)*
        }

        impl $op<$t2> for &$t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, rhs: $t2) -> $out {
                self.$meth(&rhs)
            }
        }

        impl $op<&$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(mut self, rhs: &$t2) -> $out {
                self.$meth_assign(rhs);
                self
            }
        }

        impl $op<$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(mut self, rhs: $t2) -> $out {
                self.$meth_assign(&rhs);
                self
            }
        }

        impl_binop! {@op_assign
            $t1, $t2, $out
            $op_assign {$meth_assign}
            {
                $($code_assign)*
            }
        }

        impl_binop! {@assign_op
            $t1, $t2, $out
            $assign_op {$assign_meth}
            {
                $($assign_code)*
            }
        }
    )*);
    (
        // a + b = b
        op_from
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            {
                $($code:tt)*
            }
            $op_from:ident {$meth_from:ident}
            {
                $($code_from:tt)*
            }
            $assign_op:ident {$assign_meth:ident}
            {
                $($assign_code:tt)*
            }
        )*
    ) => ($(

        impl $op<&$t2> for &$t1 {
            type Output = $out;
            #[inline]
            $($code)*
        }

        impl $op<$t2> for &$t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, mut rhs: $t2) -> $out {
                rhs.$meth_from(self);
                rhs
            }
        }

        impl $op<&$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, rhs: &$t2) -> $out {
                (&self).$meth(rhs)
            }
        }

        impl $op<$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, mut rhs: $t2) -> $out {
                rhs.$meth_from(self);
                rhs
            }
        }

        impl_binop! {@op_from
            $t1, $t2, $out
            $op_from {$meth_from}
            {
                $($code_from)*
            }
        }
        impl_binop! {@assign_op
            $t1, $t2, $out
            $assign_op {$assign_meth}
            {
                $($assign_code)*
            }
        }
    )*);
    (
        // a + b = c
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            {
                $($code:tt)*
            }
            $assign_op:ident {$assign_meth:ident}
            {
                $($assign_code:tt)*
            }
        )*
    ) => ($(
        impl $op<&$t2> for &$t1 {
            type Output = $out;
            #[inline]
            $($code)*
        }

        impl $op<$t2> for &$t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, rhs: $t2) -> $out {
                self.$meth(&rhs)
            }
        }

        impl $op<&$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, rhs: &$t2) -> $out {
                (&self).$meth(rhs)
            }
        }

        impl $op<$t2> for $t1 {
            type Output = $out;
            #[inline]
            fn $meth(self, rhs: $t2) -> $out {
                (&self).$meth(&rhs)
            }
        }

        impl_binop! {@assign_op
            $t1, $t2, $out
            $assign_op {$assign_meth}
            {
                $($assign_code)*
            }
        }
    )*);
    (
        @op_assign
        $t1:ident, $t2:ident, $out:ident
        $op_assign:ident {$meth_assign:ident}
        {
            $($code_assign:tt)*
        }
    ) => {
        impl $op_assign<&$t2> for $t1 {
            #[inline]
            $($code_assign)*
        }

        impl $op_assign<$t2> for $t1 {
            #[inline]
            fn $meth_assign(&mut self, rhs: $t2) {
                self.$meth_assign(&rhs);
            }
        }
    };
    (
        @op_from
        $t1:ident, $t2:ident, $out:ident
        $op_from:ident {$meth_from:ident}
        {
            $($code_from:tt)*
        }
    ) => {
        impl $op_from<&$t1> for $t2 {
            #[inline]
            $($code_from)*
        }

        impl $op_from<$t1> for $t2 {
            #[inline]
            fn $meth_from(&mut self, lhs: $t1) {
                self.$meth_from(&lhs);
            }
        }
    };
    (
        @assign_op
        $t1:ident, $t2:ident, $out:ident
        $assign_op:ident {$assign_meth:ident}
        {
            $($assign_code:tt)*
        }
    ) => {
        impl $assign_op<&$t1, &$t2> for $out {
            #[inline]
            $($assign_code)*
        }

        impl $assign_op<$t1, &$t2> for $out {
            #[inline]
            fn $assign_meth(&mut self, lhs: $t1, rhs: &$t2) {
                self.$assign_meth(&lhs, rhs);
            }
        }

        impl $assign_op<&$t1, $t2> for $out {
            #[inline]
            fn $assign_meth(&mut self, lhs: &$t1, rhs: $t2) {
                self.$assign_meth(lhs, &rhs);
            }
        }

        impl $assign_op<$t1, $t2> for $out {
            #[inline]
            fn $assign_meth(&mut self, lhs: $t1, rhs: $t2) {
                self.$assign_meth(&lhs, &rhs);
            }
        }
    };
}

/// Macros for overloading binary operators with unsafe functions.
macro_rules! impl_binop_unsafe {
    (
        // a + a = a
        $kw:ident
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $op_assign:ident {$meth_assign:ident}
            $op_from:ident {$meth_from:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )+
    ) => ($(
        impl_binop! {
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    op_guard!($op, $kw, $out, self, rhs);
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!($kw, $func, res, self, rhs);
                    res
                }
            }
            $op_assign {$meth_assign}
            {
                fn $meth_assign(&mut self, rhs: &$t2) {
                    op_guard!($op, $kw, $out, self, rhs);
                    call_unsafe!($kw, $func, self, self, rhs);
                }
            }
            $op_from {$meth_from}
            {
                fn $meth_from(&mut self, lhs: &$t2) {
                    op_guard!($op, $kw, $out, lhs, self);
                    call_unsafe!($kw, $func, self, lhs, self);
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    op_guard!($op, $kw, $out, lhs, rhs);
                    call_unsafe!($kw, $func, self, lhs, rhs);
                }
            }
        }
    )+);
    (
        // a + b = a
        $kw:ident
        op_assign
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $op_assign:ident {$meth_assign:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )+
    ) => ($(
        impl_binop! {
            op_assign
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!($kw, $func, res, self, rhs);
                    res
                }
            }
            $op_assign {$meth_assign}
            {
                fn $meth_assign(&mut self, rhs: &$t2) {
                    call_unsafe!($kw, $func, self, self, rhs);
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!($kw, $func, self, lhs, rhs);
                }
            }
        }
    )+);
    (
        // a + b = a, b primitive
        $kw:ident
        op_assign
        $t1:ident, $cast:ty {$($t2:ident)+}, $out:ident

        $op:ident {$meth:ident}
        $op_assign:ident {$meth_assign:ident}
        $assign_op:ident {$assign_meth:ident}
        $func:path;

        $($next:tt)*
    ) => ($(
        impl_binop_unsafe! {@inner
            $kw
            op_assign
            $t1, $cast {$t2}, $out

            $op {$meth}
            $op_assign {$meth_assign}
            $assign_op {$assign_meth}
            $func;
        })+

        impl_binop_unsafe! {
            $kw
            op_assign
            $t1, $cast {$($t2)+}, $out
            $($next)*
        }
    );
    (@inner
        $kw:ident
        op_assign
        $t1:ident, $cast:ty {$t2:ident}, $out:ident
        $(
            $op:ident {$meth:ident}
            $op_assign:ident {$meth_assign:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )*
    ) => ($(
        impl_binop! {
            op_assign
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!(cast_rhs $kw, $func, $cast, res, self, rhs);
                    res
                }
            }
            $op_assign {$meth_assign}
            {
                fn $meth_assign(&mut self, rhs: &$t2) {
                    call_unsafe!(cast_rhs $kw, $func, $cast, self, self, rhs);
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!(cast_rhs $kw, $func, $cast, self, lhs, rhs);
                }
            }
        }
    )*);
    (
        $kw:ident
        op_assign
        $t1:ident, $cast:ty {$($t2:ident)+}, $out:ident
    ) => {};
    (
        // a + b = b
        $kw:ident
        op_from
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $op_from:ident {$meth_from:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )+
    ) => ($(
        impl_binop! {
            op_from
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!($kw, $func, res, self, rhs);
                    res
                }
            }
            $op_from {$meth_from}
            {
                fn $meth_from(&mut self, lhs: &$t1) {
                    call_unsafe!($kw, $func, self, lhs, self);
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!($kw, $func, self, lhs, rhs);
                }
            }
        }
    )+);
    (
        // a + b = b, a primitive
        $kw:ident
        op_from
        $cast:ty {$($t1:ident)+}, $t2:ident, $out:ident

        $op:ident {$meth:ident}
        $op_from:ident {$meth_from:ident}
        $assign_op:ident {$assign_meth:ident}
        $func:path;

        $($next:tt)*
    ) => ($(
        impl_binop_unsafe! {@inner
            $kw
            op_from
            $cast {$t1}, $t2, $out

            $op {$meth}
            $op_from {$meth_from}
            $assign_op {$assign_meth}
            $func;
        })+

        impl_binop_unsafe! {
            $kw
            op_from
            $cast {$($t1)+}, $t2, $out
            $($next)*
        }
    );
    (@inner
        $kw:ident
        op_from
        $cast:ty {$t1:ident}, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $op_from:ident {$meth_from:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )*
    ) => ($(
        impl_binop! {
            op_from
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!(cast_lhs $kw, $func, $cast, res, self, rhs);
                    res
                }
            }
            $op_from {$meth_from}
            {
                fn $meth_from(&mut self, lhs: &$t1) {
                    call_unsafe!(cast_lhs $kw, $func, $cast, self, lhs, self);
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!(cast_lhs $kw, $func, $cast, self, lhs, rhs);
                }
            }
        }
    )*);
    (
        $kw:ident
        op_from
        $cast:ty {$($t1:ident)+}, $t2:ident, $out:ident
    ) => {};
    (
        // a + b = c
        $kw:ident
        $t1:ident, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )+
    ) => ($(
        impl_binop! {
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!($kw, $func, res, self, rhs);
                    res
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!($kw, $func, self, lhs, rhs);
                }
            }
        }
    )+);
    (
        // a + b = c, a primitive
        $kw:ident
        $cast:ty {$($t1:ident)+}, $t2:ident, $out:ident

        $op:ident {$meth:ident}
        $assign_op:ident {$assign_meth:ident}
        $func:path;

        $($next:tt)*
    ) => ($(
        impl_binop_unsafe! {@inner
            $kw
            $cast {$t1}, $t2, $out

            $op {$meth}
            $assign_op {$assign_meth}
            $func;
        })+

        impl_binop_unsafe! {
            $kw
            $cast {$($t1)+}, $t2, $out
            $($next)*
        }
    );
    (@inner
        $kw:ident
        $cast:ty {$t1:ident}, $t2:ident, $out:ident
        $(
            $op:ident {$meth:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )*
    ) => ($(
        impl_binop! {
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!(cast_lhs $kw, $func, $cast, res, self, rhs);
                    res
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!(cast_lhs $kw, $func, $cast, self, lhs, rhs);
                }
            }
        }
    )*);
    (
        $kw:ident
        $cast:ty {$($t1:ident)+}, $t2:ident, $out:ident
    ) => {};
    (
        // a + b = c, b primitive
        $kw:ident
        $t1:ident, $cast:ty {$($t2:ident)+}, $out:ident

        $op:ident {$meth:ident}
        $assign_op:ident {$assign_meth:ident}
        $func:path;

        $($next:tt)*
    ) => ($(
        impl_binop_unsafe! {@inner
            $kw
            $t1, $cast {$t2}, $out

            $op {$meth}
            $assign_op {$assign_meth}
            $func;
        })+

        impl_binop_unsafe! {
            $kw
            $t1, $cast {$($t2)+}, $out
            $($next)*
        }
    );
    (@inner
        $kw:ident
        $t1:ident, $cast:ty {$t2:ident}, $out:ident
        $(
            $op:ident {$meth:ident}
            $assign_op:ident {$assign_meth:ident}
            $func:path;
        )*
    ) => ($(
        impl_binop! {
            $t1, $t2, $out
            $op {$meth}
            {
                fn $meth(self, rhs: &$t2) -> $out {
                    let mut res = default!($op, $kw, $out, self, rhs);
                    call_unsafe!(cast_rhs $kw, $func, $cast, res, self, rhs);
                    res
                }
            }
            $assign_op {$assign_meth}
            {
                fn $assign_meth(&mut self, lhs: &$t1, rhs: &$t2) {
                    call_unsafe!(cast_rhs $kw, $func, $cast, self, lhs, rhs);
                }
            }
        }
    )*);
    (
        $kw:ident
        $t1:ident, $cast:ty {$($t2:ident)+}, $out:ident
    ) => {};
}

/// Macros for implementing `From` for conversions.
#[macro_export]
macro_rules! impl_from {
    (
        $t1:ident, $t2:ident
        {
            $($code:tt)*
        }
    ) => {
        impl From<$t2> for $t1 {
            #[inline]
            fn from(src: $t2) -> $t1 {
                <$t1>::from(&src)
            }
        }

        impl From<&$t2> for $t1 {
            #[inline]
            $($code)*
        }
    };
}

/// Macros for implementing `From` for conversions with unsafe functions.
macro_rules! impl_from_unsafe {
    (
        $kw:ident
        $t1:ident, $t2:ident
        $func:path
    ) => (
        impl_from! {
            $t1, $t2
            {
                fn from(src: &$t2) -> $t1 {
                    let mut res = default!(From, $kw, $t1, src);
                    call_unsafe!($kw, $func, res, src);
                    res
                }
            }
        }
    );
    (
        // a -> b, a primitive
        $kw:ident
        $t1:ident, $cast:ident {$($t2:ident)*}
        $func:path
    ) => ($(
        impl_from! {
            $t1, $t2
            {
                fn from(src: &$t2) -> $t1 {
                    let mut res = default!(From, $kw, $t1, src);
                    call_unsafe!(cast $kw, $func, $cast, res, src);
                    res
                }
            }
        }
    )*);
    /*
    (
        // a -> b, with third argument (precision, etc)
        $kw:ident
        $t1:ident, $t2:ident, $arg:expr;
        $func:path
    ) => (
        impl_from! {
            $t1, $t2
            {
                fn from(src: &$t2) -> $t1 {
                    let mut res = default!(From, $kw, $t1, src);
                    unsafe { $func(res.as_mut_ptr(), src.as_ptr(), $arg); }
                    res
                }
            }
        }
    );
    */
}

/// Macros for implementing `TryFrom` for conversions.
#[macro_export]
macro_rules! impl_tryfrom {
    (
        $t1:ident, $t2:ident
        {
            $($code:tt)*
        }
    ) => {
        impl TryFrom<$t2> for $t1 {
            type Error = &'static str;
            #[inline]
            fn try_from(src: $t2) -> Result<Self,Self::Error> {
                <$t1>::try_from(&src)
            }
        }

        impl TryFrom<&$t2> for $t1 {
            type Error = &'static str;
            #[inline]
            $($code)*
        }
    };
}

/// Macros for implementing `From` for conversions with unsafe functions.
macro_rules! impl_tryfrom_unsafe {
    (
        $kw:ident
        $t1:ident, $t2:ident
        $func:path
    ) => (
        impl_tryfrom! {
            $t1, $t2
            {
                fn from(src: &$t2) -> $t1 {
                    let mut res = default!(From, $kw, $t1, src);
                    call_unsafe!($kw, $func, res, src);
                    res
                }
            }
        }
    );
    (
        // a -> b, a primitive
        $kw:ident
        $t1:ident, $cast:ident {$($t2:ident)*}
        $func:path
    ) => ($(
        impl_tryfrom! {
            $t1, $t2
            {
                fn from(src: &$t2) -> $t1 {
                    let mut res = default!(From, $kw, $t1, src);
                    call_unsafe!(cast $kw, $func, $cast, res, src);
                    res
                }
            }
        }
    )*);
}

/// Macros for implementing `Assign` for in-place assignment.
#[macro_export]
macro_rules! impl_assign {
    (
        $t1:ident, $t2:ident
        {
            $($code:tt)*
        }
    ) => {
        impl Assign<$t2> for $t1 {
            #[inline]
            fn assign(&mut self, src: $t2) {
                self.assign(&src);
            }
        }

        impl Assign<&$t2> for $t1 {
            #[inline]
            $($code)*
        }
    };
}

macro_rules! impl_assign_unsafe {
    (
        $kw:ident
        $t1:ident, $t2:ident
        $func:path
    ) => (
        impl_assign! {
            $t1, $t2
            {
                fn assign(&mut self, src: &$t2) {
                    call_unsafe!($kw, $func, self, src);
                }
            }
        }
    );
    (
        // a -> b, a primitive
        $kw:ident
        $t1:ident, $cast:ident {$($t2:ident)*}
        $func:path
    ) => ($(
        impl_assign! {
            $t1, $t2
            {
                fn assign(&mut self, src: &$t2) {
                    call_unsafe!(cast $kw, $func, $cast, self, src);
                }
            }
        }
    )*);
}