@@ -436,8 +436,8 @@ pub fn sct(
436
436
. max ( dh. into_iter ( ) . sum :: < f64 > ( ) / box_size as f64 ) ;
437
437
438
438
let mut s: Mat < f64 > = Mat :: from_fn ( box_size, box_size, |i, j| {
439
- let value = ( -0.5 * ( disth. read ( i, j) / dh_mean) . powi ( 2 )
440
- - 0.5 * ( distz. read ( i, j) / args. vertical_scale ) . powi ( 2 ) )
439
+ let value = ( -0.5 * ( disth. get ( i, j) / dh_mean) . powi ( 2 )
440
+ - 0.5 * ( distz. get ( i, j) / args. vertical_scale ) . powi ( 2 ) )
441
441
. exp ( ) ;
442
442
// weight the diagonal?? (0.5 default)
443
443
if i == j {
@@ -466,14 +466,14 @@ pub fn sct(
466
466
}
467
467
468
468
let s_inv_d: Vec < f64 > = ( 0 ..box_size)
469
- . map ( |i| ( 0 ..box_size) . map ( |j| s_inv. read ( i, j) * d[ j] ) . sum ( ) )
469
+ . map ( |i| ( 0 ..box_size) . map ( |j| s_inv. get ( i, j) * d[ j] ) . sum ( ) )
470
470
. collect ( ) ;
471
471
472
472
let ares_temp: Vec < f64 > = ( 0 ..box_size)
473
- . map ( |i| ( 0 ..box_size) . map ( |j| s. read ( i, j) * s_inv_d[ j] ) . sum ( ) )
473
+ . map ( |i| ( 0 ..box_size) . map ( |j| s. get ( i, j) * s_inv_d[ j] ) . sum ( ) )
474
474
. collect ( ) ;
475
475
476
- let z_inv: Vec < f64 > = ( 0 ..box_size) . map ( |i| 1. / s_inv. read ( i, i) ) . collect ( ) ;
476
+ let z_inv: Vec < f64 > = ( 0 ..box_size) . map ( |i| 1. / s_inv. get ( i, i) ) . collect ( ) ;
477
477
478
478
let ares: Vec < f64 > = ( 0 ..box_size) . map ( |i| ares_temp[ i] - d[ i] ) . collect ( ) ;
479
479
0 commit comments