[vtkusers] Lineair systeme resolution

[medamine]

Hi ;
i want to solve the linear systeme AX=b;A is a matrix n*n,X and b are 2
vector.
i writed this code to get solution.the problem that it  get more then 1h to
return  a result because n represent the numer of point in the mesh.
Have any one a solution or idea to resolve this problem.
Thanks


float* Resolve(float** A,float* B,int taille)
{  	float  *X=new float[taille];
	for(int i=0;i<taille;i++)
        {X[i]=B[i];}
	cout<<"aa"<<endl;
	vtkMath::SolveLinearSystem((double**)A,(double*)X,taille);
	cout<<"aa"<<endl;
	for(int i=0;i<taille;i++)
		cout<<X[i]<<endl;
    return X;
}
-- 
View this message in context: http://www.nabble.com/Lineair-systeme-resolution-tf3840234.html#a10873256
Sent from the VTK - Users mailing list archive at Nabble.com.