On Mon, 27 Aug 2012, Lei Shi wrote: > So there is no class support Lagrange based discontinuous polynomial > space? If the answer is no, I think it is pretty easy to add one. > Just need the shape function and the DOF locations.
I'd forgotten that David had added this, but in any case we'd need a new codepath in System::project_vector() to handle non-hierarchic projections on coarsening or refinement. For adaptivity on steady problems those projections are necessary to get an efficient result; for transient adaptivity those projections are necessary to get a correct result. > About the efficiency, is there any paper or presentation talking > about the efficiency? Or is there any people compare the efficiency > of solver based on libmesh or other lib like deal.ii etc? It is > really a big issue which should be considered before people decided > to put a lot of effort on a lib. Thanks a lot. Last I checked (long ago) the deal.II built-in algebra was significantly faster than PETSc for supported non-MPI solver options, but do any problems on distributed systems or with advanced solvers you'd have to switch to PETSc anyway. Element residual/jacobian efficiency can be nearly as important as solver efficiency, but it's hard to make apples-to-apples comparisons there since mostly the efficiency depends on the user code. FEM shape value initialization, mesh adaptivity, etc. is typically a small enough part of runtime that I don't know if anyone's done comparisons. --- Roy ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
