[vtkusers] Lineair systeme resolution
heike hufnagel
hufiheike at yahoo.de
Wed May 30 12:08:41 EDT 2007
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
More information about the vtkusers
mailing list