tmp

src/optimization/lcg_solver_ex.cpp

@@ -29,7 +29,6 @@
 #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
@@ -110,15 +109,12 @@ 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);
 
 	ex1 test;
+	test.set_blas_backend("accelerate");
 
 	// 计算共轭梯度B项
 	gctl::array<double> B(N);