Hi all, I am currently looking for a library that can work well with PETSc and can provide me an FEM framework to handle a set of coupled nonlinear PDEs in 2 and 3 dimensions.
I hope to compare the usability of LibMesh and Deal II for this purpose. I have already installed LibMesh with the current version of PETSc on an OCS cluster. Since the number of DOFs for my system can be in the order of millions, it becomes important to ask two primary questions: 1) Are there any benchmark results for LibMesh while using with PETSc ? To rephrase that question, how does LibMesh scale with the problem ? (Weak scaling) 2) What additional code changes would a user be required to make in order to get matrix-free solution algorithms to work with LibMesh ? For example, only local assemblies need be performed and global assemblies that require matrix to be stored are not needed. For an iterative linear Krylov solve, this would be done at every matvec product request. If somebody can provide me with some pointers for the scalability results and possibly a sample code for even a simple 1-D diffusion problem, that would be fantastic. Thanks in advance for all your help. Vijay No virus found in this outgoing message. Checked by AVG Free Edition. Version: 7.5.503 / Virus Database: 269.17.2/1184 - Release Date: 12/14/2007 11:29 AM ------------------------------------------------------------------------- SF.Net email is sponsored by: Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marketplace _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
