Re: [deal.II] EigenValue Problem using dealii

2024-06-22 Thread Wolfgang Bangerth
On 6/20/24 22:09, kai huang wrote:    I am working on some problems involving eigenvalues and eigenvectors.    I am having some issues when assembling the stiffness matrix. Generally, there are three methods to handle Dirichlet boundary conditions, such as the penalty method (multiplying by a

[deal.II] EigenValue Problem using dealii

2024-06-21 Thread kai huang
Dear all, I am working on some problems involving eigenvalues and eigenvectors. I am having some issues when assembling the stiffness matrix. Generally, there are three methods to handle Dirichlet boundary conditions, such as the penalty method (multiplying by a large number); setting the

[deal.II] EigenValue Problem

2017-09-19 Thread Denis Davydov
Jean-Paul is right. Most likely you don't adjust the quadrature order which breaks LU factorization of the mass matrix, used by default for Generalized Hermitian Eigenvalue Problem. Essentially your mass matrix is not positive definite anymore. Regards, Denis -- The deal.II project is locate

Re: [deal.II] EigenValue Problem

2017-09-19 Thread Jean-Paul Pelteret
Dear Nitish, I’m not familiar with this tutorial and I don’t have much to suggest, but I have one basic question: When increasing the polynomial order of the finite element, do you increase the corresponding quadrature order as well? Regards, Jean-Paul > On 19 Sep 2017, at 14:53, Nitish Anand

[deal.II] EigenValue Problem

2017-09-19 Thread Nitish Anand
Hey, I wish to solve an eigenvalue problem for a rotating disk and to understand more about the implementation I looked into the step-36 example. Unfortunately, when I increase the degrees of freedom in the code (changing line 70 to fe(3) or fe(2)) I get the following error. Number of acti