From: Vijay M <unknownreference@gm...>  20071214 18:31:06

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 matrixfree 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 1D 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 