@@ -271,14 +271,12 @@ where
271271 where
272272 Sb : DataMut < Elem = A > ,
273273 {
274- unsafe {
275- A :: solve_tridiagonal (
276- & self ,
277- rhs. layout ( ) ?,
278- Transpose :: No ,
279- rhs. as_slice_mut ( ) . unwrap ( ) ,
280- ) ?
281- } ;
274+ A :: solve_tridiagonal (
275+ & self ,
276+ rhs. layout ( ) ?,
277+ Transpose :: No ,
278+ rhs. as_slice_mut ( ) . unwrap ( ) ,
279+ ) ?;
282280 Ok ( rhs)
283281 }
284282 fn solve_t_tridiagonal_inplace < ' a , Sb > (
@@ -288,14 +286,12 @@ where
288286 where
289287 Sb : DataMut < Elem = A > ,
290288 {
291- unsafe {
292- A :: solve_tridiagonal (
293- & self ,
294- rhs. layout ( ) ?,
295- Transpose :: Transpose ,
296- rhs. as_slice_mut ( ) . unwrap ( ) ,
297- ) ?
298- } ;
289+ A :: solve_tridiagonal (
290+ & self ,
291+ rhs. layout ( ) ?,
292+ Transpose :: Transpose ,
293+ rhs. as_slice_mut ( ) . unwrap ( ) ,
294+ ) ?;
299295 Ok ( rhs)
300296 }
301297 fn solve_h_tridiagonal_inplace < ' a , Sb > (
@@ -305,14 +301,12 @@ where
305301 where
306302 Sb : DataMut < Elem = A > ,
307303 {
308- unsafe {
309- A :: solve_tridiagonal (
310- & self ,
311- rhs. layout ( ) ?,
312- Transpose :: Hermite ,
313- rhs. as_slice_mut ( ) . unwrap ( ) ,
314- ) ?
315- } ;
304+ A :: solve_tridiagonal (
305+ & self ,
306+ rhs. layout ( ) ?,
307+ Transpose :: Hermite ,
308+ rhs. as_slice_mut ( ) . unwrap ( ) ,
309+ ) ?;
316310 Ok ( rhs)
317311 }
318312}
@@ -566,7 +560,7 @@ where
566560 A : Scalar + Lapack ,
567561{
568562 fn factorize_tridiagonal_into ( self ) -> Result < LUFactorizedTridiagonal < A > > {
569- Ok ( unsafe { A :: lu_tridiagonal ( self ) ? } )
563+ Ok ( A :: lu_tridiagonal ( self ) ?)
570564 }
571565}
572566
@@ -576,7 +570,7 @@ where
576570{
577571 fn factorize_tridiagonal ( & self ) -> Result < LUFactorizedTridiagonal < A > > {
578572 let a = self . clone ( ) ;
579- Ok ( unsafe { A :: lu_tridiagonal ( a) ? } )
573+ Ok ( A :: lu_tridiagonal ( a) ?)
580574 }
581575}
582576
@@ -587,7 +581,7 @@ where
587581{
588582 fn factorize_tridiagonal ( & self ) -> Result < LUFactorizedTridiagonal < A > > {
589583 let a = self . extract_tridiagonal ( ) ?;
590- Ok ( unsafe { A :: lu_tridiagonal ( a) ? } )
584+ Ok ( A :: lu_tridiagonal ( a) ?)
591585 }
592586}
593587
@@ -677,7 +671,7 @@ where
677671 A : Scalar + Lapack ,
678672{
679673 fn rcond_tridiagonal ( & self ) -> Result < A :: Real > {
680- unsafe { Ok ( A :: rcond_tridiagonal ( & self ) ?) }
674+ Ok ( A :: rcond_tridiagonal ( & self ) ?)
681675 }
682676}
683677
0 commit comments