[vtkusers] Lineair systeme resolution

[heike hufnagel]

Hello,

there are a lot of linear system solving methods out
there which are adapted to the resepctive problem at
hand. Is your matrix sparse? or symmetric?
You could e.g. have a look at the functions available
in the gnu scientific library (gsl).

Heike


--- medamine <miladimedamine at yahoo.fr> schrieb:

> 
> 
> 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.
> 
> _______________________________________________
> This is the private VTK discussion list. 
> Please keep messages on-topic. Check the FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
> 



      __________________________________  Yahoo! Clever: Stellen Sie Fragen und finden Sie Antworten. Teilen Sie Ihr Wissen. www.yahoo.de/clever