Re: [deal.II] constraints.distribute does not seem to be enforcing Dirichlet conditions

2022-08-30 Thread Wolfgang Bangerth
On 8/30/22 13:56, Lucas Myers wrote: To enforce these boundary conditions, I initialize constraints with: ``` constraints.clear(); dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints); dealii::VectorTools::     interpolate_boundary_values(dof_handler,                    

[deal.II] constraints.distribute does not seem to be enforcing Dirichlet conditions

2022-08-30 Thread Lucas Myers
Hi everyone, I'm trying to solve a nonlinear ODE with Dirichlet conditions using Newton's method. For this I give an initial guess for the solution ( y = (1/2) x ) and then update each Newton step with an update that ought to be zero at the boundaries (hence preserving the boundary conditions o

Re: [deal.II] Arpack and block systems !

2022-08-30 Thread Wolfgang Bangerth
On 8/30/22 09:30, Masoud Ahmadi wrote: I'm facing the same problem and found this topic. Is there any update on this issue or it is still the same as 7 years ago? My comment under the email you quote is still correct: "I don't think we have new suggestions other than the ones I have in the em

Re: [deal.II] Arpack and block systems !

2022-08-30 Thread Masoud Ahmadi
Dear all, I'm facing the same problem and found this topic. Is there any update on this issue or it is still the same as 7 years ago? Thanks On Friday, 3 May 2019 at 01:08:04 UTC+1 Wolfgang Bangerth wrote: > On 5/2/19 5:09 PM, Mert Can Simsek wrote: > > > > I'm facing the same problem and sinc

Re: [deal.II] Re: Solving the linear system of equations using PETSc BlockSparseMatrix

2022-08-30 Thread Masoud Ahmadi
Dear Bruno and Daniel, Thanks for your replies. I used linear_operator and then solve the problem like step-70. I guess there is no "Direct solver" for this situation, the one like we had for dealii BlockSparseMatrix. Best regards, Masoud On Tuesday, 23 August 2022 at 16:11:12 UTC+1 d.arnd...

Re: [deal.II] Gmsh reference

2022-08-30 Thread Wolfgang Bangerth
On 8/30/22 01:25, Joss G. wrote: I would like to ask if I import a Gmsh mesh, where the componeent  y coorresponds to the hight (and not z), will the everything inside Dealii follows this coordinate systems, and keep the reference (0,0) from the mesh I imported? deal.II just reads what GMSH

[deal.II] Re: Using the deal.II video lectures for finite element courses

2022-08-30 Thread blais...@gmail.com
Thank you very much for pointing this reference out. It is highly helpful and I think the community really appreciate you making all these efforts to not only move science forward, but teach it in a better fashion :) On Thursday, August 4, 2022 at 9:47:36 p.m. UTC-4 Wolfgang Bangerth wrote: >

[deal.II] Gmsh reference

2022-08-30 Thread Joss G.
Good morning, I would like to ask if I import a Gmsh mesh, where the componeent y coorresponds to the hight (and not z), will the everything inside Dealii follows this coordinate systems, and keep the reference (0,0) from the mesh I imported? Thank you -- The deal.II project is located at