@@ -13,14 +13,22 @@ fn cholesky() {
1313 println!( "a = \n {:?}" , a_orig) ;
1414
1515 let upper = a_orig. cholesky( UPLO :: Upper ) . unwrap( ) ;
16- assert_close_l2!( & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) , & a_orig, $rtol) ;
16+ assert_close_l2!(
17+ & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) ,
18+ & a_orig,
19+ $rtol
20+ ) ;
1721
1822 let lower = a_orig. cholesky( UPLO :: Lower ) . unwrap( ) ;
1923 assert_close_l2!( & lower. dot( & lower. t( ) . mapv( |elem| elem. conj( ) ) ) , & a_orig, $rtol) ;
2024
2125 let a: Array2 <$elem> = replicate( & a_orig) ;
2226 let upper = a. cholesky_into( UPLO :: Upper ) . unwrap( ) ;
23- assert_close_l2!( & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) , & a_orig, $rtol) ;
27+ assert_close_l2!(
28+ & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) ,
29+ & a_orig,
30+ $rtol
31+ ) ;
2432
2533 let a: Array2 <$elem> = replicate( & a_orig) ;
2634 let lower = a. cholesky_into( UPLO :: Lower ) . unwrap( ) ;
@@ -29,7 +37,11 @@ fn cholesky() {
2937 let mut a: Array2 <$elem> = replicate( & a_orig) ;
3038 {
3139 let upper = a. cholesky_inplace( UPLO :: Upper ) . unwrap( ) ;
32- assert_close_l2!( & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) , & a_orig, $rtol) ;
40+ assert_close_l2!(
41+ & upper. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) ,
42+ & a_orig,
43+ $rtol
44+ ) ;
3345 }
3446 assert_close_l2!( & a. t( ) . mapv( |elem| elem. conj( ) ) . dot( & upper. view( ) ) , & a_orig, $rtol) ;
3547
@@ -103,7 +115,11 @@ fn cholesky_det() {
103115 ( $elem: ty, $atol: expr) => {
104116 let a: Array2 <$elem> = random_hpd( 3 ) ;
105117 println!( "a = \n {:?}" , a) ;
106- let ln_det = a. eigvalsh( UPLO :: Upper ) . unwrap( ) . mapv( |elem| elem. ln( ) ) . scalar_sum( ) ;
118+ let ln_det = a
119+ . eigvalsh( UPLO :: Upper )
120+ . unwrap( )
121+ . mapv( |elem| elem. ln( ) )
122+ . scalar_sum( ) ;
107123 let det = ln_det. exp( ) ;
108124 assert_aclose!( a. factorizec( UPLO :: Upper ) . unwrap( ) . detc( ) , det, $atol) ;
109125 assert_aclose!( a. factorizec( UPLO :: Upper ) . unwrap( ) . ln_detc( ) , ln_det, $atol) ;
0 commit comments