Re: [deal.II] Re: How to get the coordinates of a given degree of freedom on an edge?

Hi, The attachment is my small program that I created. My fe system have 24 dofs for both real and imaginary of 12 edges of a cell. No matter how I change the parameter "order" the number of dofs is unchanged. I think I have problem when constructing the object. I would like to thank you very

[deal.II] Re: How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Dear Prof.Arndt, Thank you very much! I understand what you say. It helps me a lot! Best, Chucui 在 2019年2月21日星期四 UTC+8下午7:31:20，Daniel Arndt写道： > > Chucui, > > > For Question 1, I write a code to compute the L2 norm of (solution_1 - >> solution_2): >> [...] >> >> Is that right? >> > You can

[deal.II] deal.II at SIAM CSE (Spokane, WA)

Dear all, I had a quick chat with Wolfgang today and we were wondering who else is attending the SIAM CSE conference in Spokane next week? Also (in very German tradition) we decided that we absolutely have to get an informal get-together of users and developers going! Therefore, let's meet on

[deal.II] deal.II Newsletter #68

Hello everyone! This is deal.II newsletter #68. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #7745: Fix FEFaceEvaluation with faces in non-standard orientation

Re: [deal.II] Re: How to get the coordinates of a given degree of freedom on an edge?

On 2/20/19 7:01 PM, Phạm Ngọc Kiên wrote: > I tested in my codes using > > template > CSEM::CSEM() > ://mapping(1), > dof_handler(triangulation), > fe(FE_NedelecSZ(order),1,//(order), multiplicity > FE_NedelecSZ(order),1) {} > > When I changed the parameter order, I saw only

Re: [deal.II] Accuracy of Dirichlet condition for p in step-20

> I am trying to solve a system of equations that do this: > Stokes to solve for v_r and p_r for one fluid (viscous rock), I use these > solutions on the RHS of a Darcy type equation solved like step-20 for the > pressure p_f in the fluid in the domain. Using the 3 solutions, I update >

Re: [deal.II] Parallel Nonlinear Poisson code in DEALII

Dear Jean-Paul, Thank you very much for the quick response. Yes, it gives an error during execution. I will also look at the Step-40 tutorial also. Yours sincerely Boron On Thu, Feb 21, 2019 at 6:03 PM Jean-Paul Pelteret wrote: > Dear Boron, > > As a new user, I’d like to welcome to the forum

Re: [deal.II] Parallel Nonlinear Poisson code in DEALII

Dear Boron, As a new user, I’d like to welcome to the forum and deal.II! Am I correct that you only see this issue when executing your program in parallel? I think that the problem results from the way that you initialise the solution vector. That is, this line: >

[deal.II] Re: How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Chucui, For Question 1, I write a code to compute the L2 norm of (solution_1 - > solution_2): > [...] > > Is that right? > You can just take the difference of the two variables and use VectorTools::integrate_difference with a Functions::ZeroFunction as "exact solution". > > For question

[deal.II] Re: Accuracy of Dirichlet condition for p in step-20

Jane, can you write down the discretized weak formulation you using? Is there any paper showing well-posedness and suitable error estimates for that discretization? Best, Daniel Am Mittwoch, 6. Februar 2019 17:22:28 UTC+1 schrieb jane...@jandj-ltd.com: > > Hi all, > > I am trying to solve a

[deal.II] Re: Accuracy of Dirichlet condition for p in step-20

Hi all, Does anyone have any suggestions on this? I'm still struggling to get the values to equal more precisely. Thanks On Wednesday, February 6, 2019 at 7:22:28 PM UTC+3, jane...@jandj-ltd.com wrote: > > Hi all, > > I am trying to solve a system of equations that do this: > Stokes to solve

[deal.II] Parallel Nonlinear Poisson code in DEALII

Dear all, I am Boron. I am a new to DEALII. I am currently trying to write a parallel code in DEALII for solving nonlinear Poisson's equation. The file is also attahed below. My doubt is "How do we pass history variable while constructing the cell_matrix?" A code snippet is (Line No 211-225) :

Re: [deal.II] How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Dear Jean-Paul, Thank you very much for your quick reply! I will read it now! Best, Chucui 在 2019年2月21日星期四 UTC+8下午5:17:17，Jean-Paul Pelteret写道： > > Hi Chucui, > > I can only offer you a very brief reply right now: You might want to look > into the ConvergenceTable >

[deal.II] Re: How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Hi, all, If there is no function to compute convergence above, I need to write the code by myself. For Question 1, I write a code to compute the L2 norm of (solution_1 - solution_2): template double StokesProblem::norm_compute (const Vector solu_1,

Re: [deal.II] How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Hi Chucui, I can only offer you a very brief reply right now: You might want to look into the ConvergenceTable class to help you with (1). The introduction to the class directs you to some tutorials that indicate its

[deal.II] How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Hi all, I have 2 questions of computing convergence rate of L2 norm of error: 1. If I don't have exact solution, does deal.ii have functions to compute convergence rate directly? (Like VectorTools::integrate_difference, which deals with the problem with exact solution ) 2. If I want to