linearsolveexample.cpp

finite element mesh 参数化有限元网格划分

#include "stdafx.h"
#include "TaucsSolver.h"
#include "Renderer.h"

void solve_linear_system() {
	TaucsSolver solver;

	taucs_logfile("taucs_logs.txt"); 

	unsigned int N = 3;
	std::vector<double>   b(N), x(N);

	// setup matrix and rhs
	for (unsigned int i=0; i<N; ++i)
	{ 
		// Start a new matrix row
		solver.begin_row();

		// set up rhs
		b[i] = 1;
		// Add a value for a specific column in the row
		solver.add_value(i, 2*i+1);
		
		solver.end_row();
	}
	
	if (!solver.factorize() || !solver.solve(b, x))
	{
		MessageBox(NULL,"Solving failed","Error",MB_OK);
	}
	else 
	{
		char a[256];
		sprintf(a," %g %g %g", x[0], x[1], x[2]);
		renderer->displayText(a);
	}
}