[deal.II] How can we get the factorized system matrix when using SparseDirectMUMPS?

2019-06-26 Thread Phạm Ngọc Kiên
Hi colleagues, I am trying to build codes using PETScWrappers::SparseDirectMUMPS solver. Step-62 shows a very good example for using it. However, I really want to know if we can choose the factorization method in MUMPS by deal.II codes. The second question is: do we have a way to get the factor

[deal.II] use compute_no_normal_flux_constraints within step-35

2019-06-26 Thread Giovanni Formica
Dear All, I'm trying to slightly modify the step-35 code in order to better understand how boundary conditions work in this case study. The boundary id 3 is the right-hand line of the domain ("nsbench2.inp" 2D example), and I'd like to impose a zero tangential velocity in an alternative way the

Re: [deal.II] Re: NedelecSZ with adaptive mesh refinement

2019-06-26 Thread Daniel Garcia-Sanchez
Wolfgang, On Wednesday, June 26, 2019 at 3:50:22 PM UTC+2, Wolfgang Bangerth wrote: > > > OK, attached is some code that builds the constraints that correspond to > this > rather esoteric paper (but also one I enjoyed working on the most because > it's > just so weird!): > > https://www.

Re: [deal.II] Re: NedelecSZ with adaptive mesh refinement

2019-06-26 Thread Wolfgang Bangerth
On 6/26/19 3:59 AM, Alexander wrote: > indeed, this should be possible to do it that way. Would you mind putting an > example to the mentioned github issue? I might give it a thought in the > background if I see more details. OK, attached is some code that builds the constraints that correspond

Re: [deal.II] Re: NedelecSZ with adaptive mesh refinement

2019-06-26 Thread Alexander
Hi Daniel, > Let me know if your are going to take care of this or need help! > If you take care of this that would be great otherwise I can give it try > Please feel free to take a leading role here. I'm happy to help and look at your code, but I won't be able to lead this myself in the next

Re: [deal.II] Re: NedelecSZ with adaptive mesh refinement

2019-06-26 Thread Daniel Garcia-Sanchez
Hi Wolfang, On Tuesday, June 25, 2019 at 9:49:15 PM UTC+2, Wolfgang Bangerth wrote: > > > But one doesn't have to go through this function. One could just as well > write > a function that is specific to the FE_NedelecSZ element and that computes > these kinds of constraints and puts them into

Re: [deal.II] Re: NedelecSZ with adaptive mesh refinement

2019-06-26 Thread Alexander
Wolfgang, indeed, this should be possible to do it that way. Would you mind putting an example to the mentioned github issue? I might give it a thought in the background if I see more details. Best Alexander On Tuesday, June 25, 2019 at 9:49:15 PM UTC+2, Wolfgang Bangerth wrote: > > On 6/25/19

Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-26 Thread Ramprasad R
Hi Daniel, Thank you so much for your detailed Explanation. As you said earlier, I have to do the eigen value analysis. You explanation has given me an approach to start with. I can not thank you enough. Warm regards Ramprasad On Tue, Jun 25, 2019 at 6:27 PM Daniel Garcia-Sanchez < danielgarsa