round of all floats to 2 decimal places

Author: WardDeb

src/bamcoverage.rs

@@ -208,7 +208,7 @@ pub fn r_bamcoverage(
                             (fields[0].to_string(), Value {
                                 start: fields[1].parse::<u32>().unwrap(),
                                 end: fields[2].parse::<u32>().unwrap(),
-                                value: fields[3].parse::<f32>().unwrap() * sf
+                                value: (fields[3].parse::<f32>().unwrap() * sf * 100.0).round() / 100.0,
                             })
                         )
                     }

src/calc.rs

@@ -134,33 +134,39 @@ pub fn calc_ratio(
         "log2" => {
             let num: f32 = (cov1 * *sf1) + *pseudocount1;
             let den: f32 = (cov2 * *sf2) + *pseudocount2;
-            return (num / den).log2();
+            let fcov: f32 = (num / den).log2();
+            return (fcov * 100.0).round() / 100.0;
         }
         "ratio" => {
             let num: f32 = (cov1 * *sf1) + *pseudocount1;
             let den: f32 = (cov2 * *sf2) + *pseudocount2;
-            return num / den;
+            let fcov: f32 = num / den;
+            return (fcov * 100.0).round() / 100.0;
         }
         "reciprocal_ratio" => {
             let num: f32 = (cov1 * *sf1) + *pseudocount1;
             let den: f32 = (cov2 * *sf2) + *pseudocount2;
             let ratio: f32 = num / den;
             if ratio >= 1.0 {
-                return den / num;
+                let fcov: f32 = den / num;
+                return (fcov * 100.0).round() / 100.0;
             } else {
-                return -num / den;
+                let fcov: f32 = -num / den;
+                return (fcov * 100.0).round() / 100.0;
             }
         }
         "subtract" => {
             let num: f32 = (cov1 * *sf1) + *pseudocount1;
             let den: f32 = (cov2 * *sf2) + *pseudocount2;
-            return num - den;
+            let fcov: f32 = num - den;
+            return (fcov * 100.0).round() / 100.0;
         }
         _ => {
             // No operation is never allowed (on the py arg level, so just default to log2)
             let num: f32 = (cov1 * *sf1) + *pseudocount1;
             let den: f32 = (cov2 * *sf2) + *pseudocount2;
-            return (num / den).log2();
+            let fcov: f32 = (num / den).log2();
+            return (fcov * 100.0).round() / 100.0;
         }
     }
 }

src/filehandler.rs

@@ -645,7 +645,12 @@ pub fn write_matrix(
             region.strand,               // Strand field persisted from bedfile
         );
         writerow.push_str(
-            &row.iter().map(|x| (scale_regions.scale * x).to_string()).collect::<Vec<String>>().join("\t")
+            &row
+                .iter()
+                .map(
+                    |x| ((scale_regions.scale * x * 100.0).round() / 100.0).to_string()
+                )
+                .collect::<Vec<String>>().join("\t")
         );
         writerow.push_str("\n");
         encoder.write_all(writerow.as_bytes()).unwrap();

src/multibamsummary.rs

@@ -240,7 +240,10 @@ pub fn r_mbams(
                         .map(|x| x.split('\t').collect())
                         .map(|x: Vec<&str> | ( x[0].to_string(), x[1].to_string(), x[2].to_string(), x[3].parse::<f32>().unwrap() ) )
                         .collect();
-                    let counts = lines.par_iter().map(|x| x.3).collect::<Vec<_>>();
+                    let counts = lines
+                        .par_iter()
+                        .map(|x| (x.3 * 100.0).round() / 100.0)
+                        .collect::<Vec<_>>();
                     let regions: (String, String, String) = (lines[0].0.clone(), lines[0].1.clone(), lines[0].2.clone());
                     _matvec.push(counts);
                     _regions.push(regions);