C++ 数学计算库 TooN

TooN 是一个C++ 数学计算库，其目的是有效运作的大量小型矩阵，并提供方便一些算法包括矩阵分解和优化。

clear all 
out={};
symmetric=1

for S=5:5
	T = num2str（S）;

	m=sym（zeros（SS））;
	for r=1:S
		for c=1:S
			if r >= c | symmetric == 0
				v = [ ‘A‘ num2str（r-1） num2str（c-1） ];
			else
				v = [ ‘A‘ num2str（c-1） num2str（r-1） ];
			end
			m（r c） = sym（v）;
		end
	end

	m
	pause

	inverse =simplify（inv（m） * det（m））;
	idet = ccode（det（m））;


	out{end+1} = [ ‘template void solve_direct（const Matrix<‘ T ‘>& A const Matrix<‘ T ‘ Cols>& b Matrix<‘ T ‘Cols>& x）‘ ];
	out{end+1} = [ ‘{‘ ];

	idet = [ idet   ‘                                      ‘];
	q=[];
	generated = zeros（1100）;

	for i=1:length（idet）
		if idet（i） >= ‘0‘ & idet（i） <= ‘9‘  & （idet（i-1） == ‘A‘ | idet（i-2） == ‘A‘）
			q = [q ‘[‘ idet（i） ‘]‘ ];
		elseif idet（i） == ‘M‘ & idet（i+13） == ‘=‘ | idet（i） == ‘t‘
			genvar_num  = str2num（idet（i+11））;

			if generated（genvar_num） == 1
				q = [q idet（i）];
			else
				q = [q ‘double ‘ idet（i）];
				generated（genvar_num）=1;
			end

		elseif idet（i） == ‘;‘
			out{end+1} = [q ‘;‘];
			q=[];
		else
			q = [q idet（i）];
		end

	end

	out{end+1} = ‘	double idet = 1/t0;‘;

	for r=1:S
		for c=1:S
			
			p = ccode（inverse（rc））;
			p = [p ‘                                            ‘];

			for i=1:length（p）
				if p（i） >= ‘0‘ & p（i） <= ‘9‘  & （p（i-1） == ‘A‘ | p（i-2） == ‘A‘）
					q = [q ‘[‘ p（i） ‘]‘ ];
				elseif p（i） == ‘;‘
					out{end+1} = [q ‘;‘];
					q=[];
				else
					q = [q p（i）];
				end

			end


			out{end+1} = [ ‘	double i‘ num2str（r-1） num2str（c-1） ‘ = t0*idet;‘];
		end
	end
	
	out{end+1} = ‘	for（unsigned int i=0; i < x.num_cols（）; i++）‘;
	out{end+1} = ‘	{‘;
	for r=1:S
		row = [‘		x[‘ num2str（r-1） ‘][i] = ‘];
		for c=1:S
			row = [row ‘i‘ num2str（r-1） num2str（c-1） ‘*b[‘ num2str（c-1） ‘][i] + ‘ ];
		end
		out{end+1} = [row（1:end-3） ‘;‘];
	end
	out{end+1} = ‘	}‘;
	out{end+1} = ‘}‘;
end

out{end+1} = [ ‘static const int highest_solver= ‘ T ‘;‘];

slv = fopen（‘solvers.cc‘ ‘w‘）;
for i=1:length（out）
	fprintf（slv ‘%s\n‘ out{i}）;
end
fclose（slv）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-02-08 20:52  TooN-2.0\
     文件         171  2012-02-06 17:53  TooN-2.0\TooN.pc
     文件         160  2012-02-06 17:53  TooN-2.0\TooN.pc.in
     文件        2080  2012-02-06 17:53  TooN-2.0\generated.h
     目录           0  2012-02-06 17:53  TooN-2.0\regressions\
     文件        1003  2012-02-06 17:53  TooN-2.0\regressions\eigen-sqrt.cc
     文件         877  2012-02-06 17:53  TooN-2.0\regressions\qr.txt
     文件         593  2012-02-06 17:53  TooN-2.0\regressions\eigen-sqrt.txt
     文件          52  2012-02-06 17:53  TooN-2.0\regressions\slice.txt
     文件        1230  2012-02-06 17:53  TooN-2.0\regressions\chol_lapack.txt
     文件        6559  2012-02-06 17:53  TooN-2.0\regressions\determinant.txt
     文件        1554  2012-02-06 17:53  TooN-2.0\regressions\gauss_jordan.cc
     文件        3226  2012-02-06 17:53  TooN-2.0\regressions\lu.txt
     文件        1715  2012-02-06 17:53  TooN-2.0\regressions\chol_toon.cc
     文件         250  2012-02-06 17:53  TooN-2.0\regressions\so3.cc
     文件          85  2012-02-06 17:53  TooN-2.0\regressions\so3.txt
     文件         822  2012-02-06 17:53  TooN-2.0\regressions\fill.cc
     文件         571  2012-02-06 17:53  TooN-2.0\regressions\gr_svd.txt
     文件        1202  2012-02-06 17:53  TooN-2.0\regressions\determinant.m
     文件       21297  2012-02-06 17:53  TooN-2.0\regressions\sym_eigen.cc
     文件         101  2012-02-06 17:53  TooN-2.0\regressions\simplex.txt
     文件         212  2012-02-06 17:53  TooN-2.0\regressions\regression.h
     文件        1301  2012-02-06 17:53  TooN-2.0\regressions\chol_toon.txt
     文件         518  2012-02-06 17:53  TooN-2.0\regressions\complex.txt
     文件        1553  2012-02-06 17:53  TooN-2.0\regressions\chol_lapack.cc
     文件         456  2012-02-06 17:53  TooN-2.0\regressions\slice.cc
     文件         104  2012-02-06 17:53  TooN-2.0\regressions\fill.txt
     文件         603  2012-02-06 17:53  TooN-2.0\regressions\qr.cc
     文件        3156  2012-02-06 17:53  TooN-2.0\regressions\gauss_jordan.txt
     文件          58  2012-02-06 17:53  TooN-2.0\regressions\vector_resize.txt
     文件        1238  2012-02-06 17:53  TooN-2.0\regressions\simplex.cc
............此处省略796个文件信息