@@ -271,14 +271,12 @@ where
271
271
where
272
272
Sb : DataMut < Elem = A > ,
273
273
{
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
+ ) ?;
282
280
Ok ( rhs)
283
281
}
284
282
fn solve_t_tridiagonal_inplace < ' a , Sb > (
@@ -288,14 +286,12 @@ where
288
286
where
289
287
Sb : DataMut < Elem = A > ,
290
288
{
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
+ ) ?;
299
295
Ok ( rhs)
300
296
}
301
297
fn solve_h_tridiagonal_inplace < ' a , Sb > (
@@ -305,14 +301,12 @@ where
305
301
where
306
302
Sb : DataMut < Elem = A > ,
307
303
{
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
+ ) ?;
316
310
Ok ( rhs)
317
311
}
318
312
}
@@ -566,7 +560,7 @@ where
566
560
A : Scalar + Lapack ,
567
561
{
568
562
fn factorize_tridiagonal_into ( self ) -> Result < LUFactorizedTridiagonal < A > > {
569
- Ok ( unsafe { A :: lu_tridiagonal ( self ) ? } )
563
+ Ok ( A :: lu_tridiagonal ( self ) ?)
570
564
}
571
565
}
572
566
@@ -576,7 +570,7 @@ where
576
570
{
577
571
fn factorize_tridiagonal ( & self ) -> Result < LUFactorizedTridiagonal < A > > {
578
572
let a = self . clone ( ) ;
579
- Ok ( unsafe { A :: lu_tridiagonal ( a) ? } )
573
+ Ok ( A :: lu_tridiagonal ( a) ?)
580
574
}
581
575
}
582
576
@@ -587,7 +581,7 @@ where
587
581
{
588
582
fn factorize_tridiagonal ( & self ) -> Result < LUFactorizedTridiagonal < A > > {
589
583
let a = self . extract_tridiagonal ( ) ?;
590
- Ok ( unsafe { A :: lu_tridiagonal ( a) ? } )
584
+ Ok ( A :: lu_tridiagonal ( a) ?)
591
585
}
592
586
}
593
587
@@ -677,7 +671,7 @@ where
677
671
A : Scalar + Lapack ,
678
672
{
679
673
fn rcond_tridiagonal ( & self ) -> Result < A :: Real > {
680
- unsafe { Ok ( A :: rcond_tridiagonal ( & self ) ?) }
674
+ Ok ( A :: rcond_tridiagonal ( & self ) ?)
681
675
}
682
676
}
683
677
0 commit comments