On Thu, Jan 8, 2015 at 12:38 PM, Harshad Sahasrabudhe <[email protected]> wrote: > Do I need to do something different from the normal loop over elements and > DOFs for integration when using Hierarchical elements of order >3? I'm > evaluating a matrix of which I find eigenvalues. This works for order <=3 > for Hierarchical and Lagrange (<=2) on HEX27 geometry. When I switch to > order 4 Hierarchical on the same geometry, I stop getting any eigenvalues > for the matrix in the same range.
Are you sure it converged? I know as you increase p, the condition number gets higher/problem gets harder to solve. Run your code with -ksp_monitor, etc. It's also possible there's a bug with 4th-order hierarchics in 3D, the cube_indices() function in fe_hierarchic_shape_3D.C scares me. -- John ------------------------------------------------------------------------------ Dive into the World of Parallel Programming! The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
