RE: how do i solve sparse linear equation like A*x=b

Hi, 

We use MTL matrices in our GPL - Finite Element Library FELyX (http://felyx.sourceforge.net/). Therefore we implemented a sparse solver that works with MTL envelope matrices, and which is reasonably fast, especially if the provided ATLAS interface is used for the most inner dot product.

If you want to have a really FAST solver for linear equation systems, use PARDISO from the University of Basel. It is available within the latest MKL Release (Math Libs from Intel) or directly from the developers (http://www.computational.unibas.ch/cs/scicomp/software/pardiso/). In the source code of FELyX, you will also find an interface to use this solver with MTL sparse matrices. With PARDISO, we solve a Finite Element model with app. 300'000 degrees of freedom ( = number of equations) in about 75s ( 1 Xeon 2.8GHz processor)...

Greetings
OK

----------------------------------------------------
Oliver Koenig

Center of Structure Technologies
Institute of Mechanical Systems
Leonhardstr. 27, LEO C4
ETH Zurich
CH-8092 Zurich

Phone : +41 1 632 23 28
Fax : +41 1 633 11 25
Mail : [email protected]
WWW : http://www.imes.ethz.ch/st/
----------------------------------------------------

________________________________________
From: [email protected] [mailto:[email protected]] On Behalf Of yanyaqin
Sent: Mittwoch, 1. September 2004 14:16
To: [email protected]
Subject: MTL: how do i solve sparse linear equation like A*x=b

Hi folks,
it seems that MTL can only solve linear equation like A*x=b for a dense or trianguler sparse  matrix A.
how do i solve for a normal sparse or symmetric sparse matrix A ?
 
If MTL is not very convenience, would you sugguest another one?
 
Thanks and best regards
 
Yan Yaqin 
2004-09-01
________________________________________
The mail have filted by the McAfee anti-virus gateway 

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/