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. Do you think that is a good option or do you have any other recommendations for me? Thanks, Kathrin ------------------------------------------------------------------------------ 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
