On Thu, Dec 8, 2016 at 11:03 AM, Kong, Fande <[email protected]> wrote:
> On Thu, Dec 8, 2016 at 10:44 AM, Kathrin Smetana <[email protected]> wrote: > > > Dear libmesh users/developers, > > > > I have to solve a linear system of equations (system size approximately > > 10^6) very often (about 5000 times). The linear system of equations is > the > > result of a FEM discretization of 3D linear elasticity. > > > > I thought about using a sparse Cholesky decomposition as the matrix is > > symmetric or a sparse LU decomposition, depending on availability. I have > > had a look at the Eigen package and their direct LU factorization for > > instance with Umfpack or SuperLU support looks very promising to me. > > > > I was wondering why you are going to use a direct solver? Because she wants to solve the same system many times, and therefore saving the factors would be worthwhile. -- John ------------------------------------------------------------------------------ Developer Access Program for Intel Xeon Phi Processors Access to Intel Xeon Phi processor-based developer platforms. With one year of Intel Parallel Studio XE. Training and support from Colfax. Order your platform today.http://sdm.link/xeonphi _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
