### Vectorfunction space degree and Quadrature degree

228

views

0

Hi there

Is it always necessary for the quadrature degree to be equal to be VectorFunctionSpace degree ?

I am trying to solve a nonlinear elasticity problem which gets solved when i set VectorFunctionSpace degree = 2 and Quadrature degree = 4 but the solution process gets killed when i used the degree of both VectorFunctionSpace as well as Quadrature = 4.

Have not posted code because that makes the question very long.

Thanks. (please also comment if that is not required or necessary as a rule)

Is it always necessary for the quadrature degree to be equal to be VectorFunctionSpace degree ?

I am trying to solve a nonlinear elasticity problem which gets solved when i set VectorFunctionSpace degree = 2 and Quadrature degree = 4 but the solution process gets killed when i used the degree of both VectorFunctionSpace as well as Quadrature = 4.

Have not posted code because that makes the question very long.

Thanks. (please also comment if that is not required or necessary as a rule)

Community: FEniCS Project

Please login to add an answer/comment or follow this question.

killed. Where as for VectorFuncationSpace degree = 2 and Quadrature degree =4 the code runs and solves the problem. What does killed indicate ? does it indicate an error or does it indicate excessive computations. I am using Core i 7 with 16 GB RAM.I have also added the code so that you and others interested in this post may have an idea of the complete problem.

Thanks for response.

Thanks for the response. I was able to find the answer to my querry by trying vector function space degree =1, 2 against Quadrature degree = 2 one by one and then compared the results.

Can anyone suggest how to improve the performance of mumps ?. There are 2,44,000 dof in the system of equation that i am trying to solve.For vector function space degree lesser than quadrature degree i was getting wrong or erroneous displacements where as for vector function space degree = Quadrature degree the displacement that i got was correct. I verified the solution by comparing the fem solution with solution of scipy optimize (on simplified problem). Hence my finding is that Vector function space degree for Nonlinear hyperelasticity problems related to stresses should be equal to quadrature degree. I am open to criticism on my results. If someone from community has an objection or better explanation then it would definitely help me improve my conceptual understanding too.

Now..

I have been trying figure out way to somehow use mumps for vectorfunction space degree = 4 because mumps seems to be performing better than iterative solver for vectorfunction space degree =2. so far i am using the following :-

Best regards.