Restore unsafe-passthrough

Author: jnsiemer

src/integer_mod_q/mat_polynomial_ring_zq.rs

@@ -31,6 +31,7 @@ mod sort;
 mod tensor;
 mod to_string;
 mod transpose;
+mod unsafe_functions;
 mod vector;
 
 /// [`MatPolynomialRingZq`] is a matrix with entries of type [`PolynomialRingZq`](crate::integer_mod_q::PolynomialRingZq).

src/integer_mod_q/mat_polynomial_ring_zq/reduce.rs

@@ -40,7 +40,7 @@ impl MatPolynomialRingZq {
     pub(crate) fn reduce(&mut self) {
         for row_num in 0..self.matrix.get_num_rows() {
             for column_num in 0..self.matrix.get_num_columns() {
-                self.reduce_entry(row_num, column_num);
+                unsafe { self.reduce_entry(row_num, column_num) };
             }
         }
     }
@@ -62,9 +62,9 @@ impl MatPolynomialRingZq {
     /// let poly_mat = MatPolyOverZ::from_str("[[4  -1 0 1 1, 1  42],[0, 2  1 2]]").unwrap();
     /// let mut poly_ring_mat = MatPolynomialRingZq::from((&poly_mat, &modulus));
     ///
-    /// poly_ring_mat.reduce_entry(0, 0)
+    /// unsafe { poly_ring_mat.reduce_entry(0, 0) };
     /// ```
-    pub(crate) fn reduce_entry(&mut self, row: i64, column: i64) {
+    pub(crate) unsafe fn reduce_entry(&mut self, row: i64, column: i64) {
         let entry = unsafe { fmpz_poly_mat_entry(&self.matrix.matrix, row, column) };
         if (unsafe { *entry }).length > 0 {
             unsafe {

src/integer_mod_q/mat_polynomial_ring_zq/tensor.rs

@@ -154,12 +154,12 @@ unsafe fn set_matrix_window_mul(
                     ),
                     &scalar.poly,
                     fmpz_poly_mat_entry(&matrix.matrix.matrix, i_other, j_other),
-                )
+                );
+                out.reduce_entry(
+                    row_left * rows_other + i_other,
+                    column_upper * columns_other + j_other,
+                );
             }
-            out.reduce_entry(
-                row_left * rows_other + i_other,
-                column_upper * columns_other + j_other,
-            );
         }
     }
 }

src/integer_mod_q/mat_polynomial_ring_zq/unsafe_functions.rs

@@ -0,0 +1,27 @@
+// Copyright © 2025 Niklas Siemer
+//
+// This file is part of qFALL-math.
+//
+// qFALL-math is free software: you can redistribute it and/or modify it under
+// the terms of the Mozilla Public License Version 2.0 as published by the
+// Mozilla Foundation. See <https://mozilla.org/en-US/MPL/2.0/>.
+
+//! This module contains public functions that enable access to underlying
+//! [FLINT](https://flintlib.org/) structs. Therefore, they require to be unsafe.
+use super::MatPolynomialRingZq;
+use crate::macros::unsafe_passthrough::{unsafe_getter_indirect, unsafe_setter_indirect};
+use flint_sys::fmpz_poly_mat::fmpz_poly_mat_struct;
+
+unsafe_getter_indirect!(
+    MatPolynomialRingZq,
+    matrix,
+    get_fmpz_poly_mat_struct,
+    fmpz_poly_mat_struct
+);
+
+unsafe_setter_indirect!(
+    MatPolynomialRingZq,
+    matrix,
+    set_fmpz_poly_mat_struct,
+    fmpz_poly_mat_struct
+);

src/integer_mod_q/polynomial_ring_zq.rs

@@ -33,6 +33,7 @@ mod reduce;
 mod sample;
 mod set;
 mod to_string;
+mod unsafe_functions;
 
 /// [`PolynomialRingZq`] represents polynomials over the finite field
 /// [`PolyOverZq`](crate::integer_mod_q::PolyOverZq)/f(X) where f(X) is a polynomial over [`Zq`](super::Zq).

src/integer_mod_q/polynomial_ring_zq/unsafe_functions.rs

@@ -0,0 +1,27 @@
+// Copyright © 2025 Niklas Siemer
+//
+// This file is part of qFALL-math.
+//
+// qFALL-math is free software: you can redistribute it and/or modify it under
+// the terms of the Mozilla Public License Version 2.0 as published by the
+// Mozilla Foundation. See <https://mozilla.org/en-US/MPL/2.0/>.
+
+//! This module contains public functions that enable access to underlying
+//! [FLINT](https://flintlib.org/) structs. Therefore, they require to be unsafe.
+use super::PolynomialRingZq;
+use crate::macros::unsafe_passthrough::{unsafe_getter_indirect, unsafe_setter_indirect};
+use flint_sys::fmpz_poly::fmpz_poly_struct;
+
+unsafe_getter_indirect!(
+    PolynomialRingZq,
+    poly,
+    get_fmpz_poly_struct,
+    fmpz_poly_struct
+);
+
+unsafe_setter_indirect!(
+    PolynomialRingZq,
+    poly,
+    set_fmpz_poly_struct,
+    fmpz_poly_struct
+);

src/utils/reduce.rs

@@ -6,7 +6,7 @@
 // the terms of the Mozilla Public License Version 2.0 as published by the
 // Mozilla Foundation. See <https://mozilla.org/en-US/MPL/2.0/>.
 
-//! Implementations to reduce a [`MatPolynomialRingZq`] with the
+//! Implementations to reduce a [`MatPolynomialRingZq`](crate::integer_mod_q::MatPolynomialRingZq) with the
 //! [`ModulusPolynomialRingZq`](crate::integer_mod_q::ModulusPolynomialRingZq).
 //! Only contains internal functions for reductions.