@@ -101,19 +101,19 @@ fn cholesky_inverse() {
101
101
#[ test]
102
102
fn cholesky_det ( ) {
103
103
macro_rules! cholesky_det {
104
- ( $elem: ty, $rtol : expr) => {
104
+ ( $elem: ty, $atol : expr) => {
105
105
let a: Array2 <$elem> = random_hpd( 3 ) ;
106
106
println!( "a = \n {:?}" , a) ;
107
107
let det = a. eigvalsh( UPLO :: Upper ) . unwrap( ) . mapv( |elem| elem. ln( ) ) . scalar_sum( ) . exp( ) ;
108
- assert_rclose !( a. detc( ) . unwrap( ) , det, $rtol ) ;
109
- assert_rclose !( a. factorizec( UPLO :: Upper ) . unwrap( ) . detc( ) , det, $rtol ) ;
110
- assert_rclose !( a. factorizec( UPLO :: Lower ) . unwrap( ) . detc( ) , det, $rtol ) ;
108
+ assert_aclose !( a. detc( ) . unwrap( ) , det, $atol ) ;
109
+ assert_aclose !( a. factorizec( UPLO :: Upper ) . unwrap( ) . detc( ) , det, $atol ) ;
110
+ assert_aclose !( a. factorizec( UPLO :: Lower ) . unwrap( ) . detc( ) , det, $atol ) ;
111
111
}
112
112
}
113
113
cholesky_det ! ( f64 , 1e-9 ) ;
114
- cholesky_det ! ( f32 , 1e-4 ) ;
114
+ cholesky_det ! ( f32 , 1e-3 ) ;
115
115
cholesky_det ! ( c64, 1e-9 ) ;
116
- cholesky_det ! ( c32, 1e-4 ) ;
116
+ cholesky_det ! ( c32, 1e-3 ) ;
117
117
}
118
118
119
119
#[ test]
0 commit comments