> > 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.
I tried, but I don't understand what's happening inside the cube_indices() function. 51 // The only way to make any sense of this 52 // is to look at the mgflo/mg2/mgf documentation Where can I find this documentation for the cube_indices() function? 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. I'm using MKL FEAST for calculating the eigenvalues. It doesn't use an iterative method and first calculates the number of eigenvalues in an interval. I get a 0 at that step for an order 4 Hierarchic basis in HEX27 geometry. On Thu, Jan 8, 2015 at 6:32 PM, John Peterson <[email protected]> wrote: > 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
