31 template<class Interface>
39     return "lu_solve_"+Interface::name();  in name()
51     typename Interface::stl_matrix A_stl;  in calculate()
52     typename Interface::stl_vector B_stl;  in calculate()
53     typename Interface::stl_vector X_stl;  in calculate()
61     typename Interface::gene_matrix A;  in calculate()
62     typename Interface::gene_vector B;  in calculate()
63     typename Interface::gene_vector X;  in calculate()
65     typename Interface::gene_matrix LU;   in calculate()
67     Interface::matrix_from_stl(A,A_stl);  in calculate()
68     Interface::vector_from_stl(B,B_stl);  in calculate()
69     Interface::vector_from_stl(X,X_stl);  in calculate()
70     Interface::matrix_from_stl(LU,A_stl);  in calculate()
74     typename Interface::Pivot_Vector pivot; // pivot vector  in calculate()
75     Interface::new_Pivot_Vector(pivot,size);  in calculate()
88       Interface::copy_matrix(A,LU,size);  in calculate()
89       Interface::LU_factor(LU,pivot,size);  in calculate()
93       Interface::LU_solve(LU,pivot,B,X,size);  in calculate()
105     typename Interface::stl_vector B_new_stl(size);  in calculate()
106     Interface::vector_to_stl(X,X_stl);  in calculate()
108 …  STL_interface<typename Interface::real_type>::matrix_vector_product(A_stl,X_stl,B_new_stl,size);   in calculate()
110     typename Interface::real_type error=  in calculate()
111       STL_interface<typename Interface::real_type>::norm_diff(B_stl,B_new_stl);  in calculate()
115       STL_interface<typename Interface::real_type>::display_vector(B_stl);  in calculate()
116       STL_interface<typename Interface::real_type>::display_vector(B_new_stl);  in calculate()
122     Interface::free_matrix(A,size);  in calculate()
123     Interface::free_vector(B);  in calculate()
124     Interface::free_vector(X);  in calculate()
125     Interface::free_Pivot_Vector(pivot);  in calculate()