diff --git a/src/optimization/lcg_solver_ex.cpp b/src/optimization/lcg_solver_ex.cpp
index acdd40f..81bfcd4 100644
--- a/src/optimization/lcg_solver_ex.cpp
+++ b/src/optimization/lcg_solver_ex.cpp
@@ -29,6 +29,7 @@
 #include "gctl/core/matrix_algorithm.h"
 #include "gctl/optimization/lcg.h"
 #include "gctl/graphic/gnuplot.h"
+#include "geoblas.h"
 
 #define M 1000
 #define N 800
@@ -109,6 +110,10 @@ void ex1::LCG_Mx(const gctl::array<double> &x, gctl::array<double> &mx)
 
 int main(int argc, char const *argv[])
 {
+	// 选择BLAS后端
+	geoblas_backend_set("accelerate");
+	std::cout << "BLAS backend: " << geoblas_backend_name() << "\n";
+
 	// 生成一组正演解
 	gctl::array<double> fm(N);
 	random_float(fm, 1.0, 2.0, gctl::RdUniform);