diff --git a/Cargo.toml b/Cargo.toml
index 4ce8e42..c003a95 100644
--- a/Cargo.toml
+++ b/Cargo.toml
@@ -18,5 +18,5 @@ chrono = "^0.4.10"
 criterion = { version = "0.4", features = ["html_reports"] }
 
 [[bench]]
-name = "benchmark"
+name = "benchmarks"
 harness = false
\ No newline at end of file
diff --git a/benches/benchmarks.rs b/benches/benchmarks.rs
new file mode 100644
index 0000000..cd4cf9f
--- /dev/null
+++ b/benches/benchmarks.rs
@@ -0,0 +1,139 @@
+// Combined benchmarks for rustframe
+use criterion::{criterion_group, criterion_main, Criterion};
+use rustframe::{
+    frame::{Frame, RowIndex},
+    matrix::{BoolMatrix, Matrix},
+    utils::{BDateFreq, BDatesList},
+};
+use chrono::NaiveDate;
+
+fn bool_matrix_operations_benchmark(c: &mut Criterion) {
+    let sizes = [1, 100, 1000];
+
+    for &size in &sizes {
+        let data1: Vec<bool> = (0..size * size).map(|x| x % 2 == 0).collect();
+        let data2: Vec<bool> = (0..size * size).map(|x| x % 3 == 0).collect();
+        let bm1 = BoolMatrix::from_vec(data1.clone(), size, size);
+        let bm2 = BoolMatrix::from_vec(data2.clone(), size, size);
+
+        c.bench_function(&format!("bool_matrix_and ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 & &bm2;
+            });
+        });
+
+        c.bench_function(&format!("bool_matrix_or ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 | &bm2;
+            });
+        });
+
+        c.bench_function(&format!("bool_matrix_xor ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 ^ &bm2;
+            });
+        });
+
+        c.bench_function(&format!("bool_matrix_not ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = !&bm1;
+            });
+        });
+    }
+}
+
+fn matrix_boolean_operations_benchmark(c: &mut Criterion) {
+    let sizes = [1, 100, 1000];
+
+    for &size in &sizes {
+        let data1: Vec<bool> = (0..size * size).map(|x| x % 2 == 0).collect();
+        let data2: Vec<bool> = (0..size * size).map(|x| x % 3 == 0).collect();
+        let bm1 = BoolMatrix::from_vec(data1.clone(), size, size);
+        let bm2 = BoolMatrix::from_vec(data2.clone(), size, size);
+
+        c.bench_function(&format!("boolean AND ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 & &bm2;
+            });
+        });
+
+        c.bench_function(&format!("boolean OR ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 | &bm2;
+            });
+        });
+
+        c.bench_function(&format!("boolean XOR ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &bm1 ^ &bm2;
+            });
+        });
+
+        c.bench_function(&format!("boolean NOT ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = !&bm1;
+            });
+        });
+    }
+}
+
+fn matrix_operations_benchmark(c: &mut Criterion) {
+    let n_periods = 4;
+    let dates: Vec<NaiveDate> =
+        BDatesList::from_n_periods("2024-01-02".to_string(), BDateFreq::Daily, n_periods)
+            .unwrap()
+            .list()
+            .unwrap();
+
+    let col_names: Vec<String> = vec!["a".to_string(), "b".to_string()];
+
+    let ma = Matrix::from_cols(vec![vec![1.0, 2.0, 3.0, 4.0], vec![5.0, 6.0, 7.0, 8.0]]);
+    let mb = Matrix::from_cols(vec![vec![4.0, 3.0, 2.0, 1.0], vec![8.0, 7.0, 6.0, 5.0]]);
+
+    let fa = Frame::new(
+        ma.clone(),
+        col_names.clone(),
+        Some(RowIndex::Date(dates.clone())),
+    );
+    let fb = Frame::new(mb, col_names, Some(RowIndex::Date(dates)));
+
+    c.bench_function("element-wise multiply", |b| {
+        b.iter(|| {
+            let _result = &fa * &fb;
+        });
+    });
+
+    let sizes = [1, 100, 1000];
+
+    for &size in &sizes {
+        let data: Vec<f64> = (0..size * size).map(|x| x as f64).collect();
+        let ma = Matrix::from_vec(data.clone(), size, size);
+
+        c.bench_function(&format!("scalar add ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &ma + 1.0;
+            });
+        });
+
+        c.bench_function(&format!("scalar subtract ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &ma - 1.0;
+            });
+        });
+
+        c.bench_function(&format!("scalar multiply ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &ma * 2.0;
+            });
+        });
+
+        c.bench_function(&format!("scalar divide ({}x{})", size, size), |b| {
+            b.iter(|| {
+                let _result = &ma / 2.0;
+            });
+        });
+    }
+}
+
+criterion_group!(combined_benches, bool_matrix_operations_benchmark, matrix_boolean_operations_benchmark, matrix_operations_benchmark);
+criterion_main!(combined_benches);
\ No newline at end of file