Re: [deal.II] Re: Anistropic refinement DG - saddle point problem

2020-08-17 Thread Juan Felipe Giraldo
Thank you Wolfgang, I just realized the new version of step 12. I was working with the older one, which made it complicated, especially on the inner faces. Regards, Juan On Tue, Aug 18, 2020 at 7:28 AM Wolfgang Bangerth wrote: > On 8/17/20 6:02 AM, jfgir...@gmail.com wrote: > > > > I would lik

[deal.II] METIS Issue during Installation Process

2020-08-17 Thread Edward Read
Hello, I am trying to install Deal.ii on my new cluster, from the build screen output it looks there is a problem in finding METIS (which is installed). I am not exactly sure what is going on? Is there a way to hard code the METIS Path during the make process? I have attached the screen out fr

Re: [deal.II] PETSc iteration does not converge

2020-08-17 Thread yuesu jin
OK, I see. I will do it, thank you! I will let you know if I get any result. On Mon, Aug 17, 2020 at 7:06 PM Wolfgang Bangerth wrote: > On 8/17/20 6:01 PM, yuesu jin wrote: > > I did nothing to verify those properties. because the single thread > CG > > solver converged well. I used differ

Re: [deal.II] PETSc iteration does not converge

2020-08-17 Thread Wolfgang Bangerth
On 8/17/20 6:01 PM, yuesu jin wrote:     I did nothing to verify those properties. because the single thread CG solver converged well. I used  different preconditioners in parallel version and single thread version. In the parallel version I used block Jacobi and in the single thread version I

Re: [deal.II] PETSc iteration does not converge

2020-08-17 Thread yuesu jin
Dear Dr.Bangerth, I did nothing to verify those properties. because the single thread CG solver converged well. I used different preconditioners in parallel version and single thread version. In the parallel version I used block Jacobi and in the single thread version I used Jacobi. How can I

Re: [deal.II] PETSc iteration does not converge

2020-08-17 Thread Wolfgang Bangerth
   I tried both, first I tried 1e-4*system_rhs.l2_norm(), it failed. Then either the matrix is not symmetric/positive definite/whatever other property your iterative solver requires, or your preconditioner is unsuitable. What have you don to verify that your matrix has the necessary proper

Re: [deal.II] mmult memory leak with petsc

2020-08-17 Thread Wolfgang Bangerth
Richard, I am working on incompressible flow problems and stumbled upon an issue when calling PETScWrappers::SparseMatrix::mmult(), but before I describe the problem in more detail, let me comment on the basic building blocks of the MWE: (i) parallel::distributed::Triangulation & either PETS

Re: [deal.II] Re: Anistropic refinement DG - saddle point problem

2020-08-17 Thread Wolfgang Bangerth
On 8/17/20 6:02 AM, jfgir...@gmail.com wrote: I would like to open again the topic with another question. Is there any way to use the Meshworker to solve the DG formulation but using block matrix and vectors?   I couldn't find a proper way to do it with the Meshworker because I have to choose

[deal.II] deal.II Newsletter #130

2020-08-17 Thread Rene Gassmoeller
Hello everyone! This is deal.II newsletter #130. 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: #10833: Create 1D Advection-Diffusion Equation (proposed by syedtahirbu

[deal.II] mmult memory leak with petsc

2020-08-17 Thread Richard Schussnig
Hi everyone, I am working on incompressible flow problems and stumbled upon an issue when calling PETScWrappers::SparseMatrix::mmult(), but before I describe the problem in more detail, let me comment on the basic building blocks of the MWE: (i) parallel::distributed::Triangulation & either PET

[deal.II] Re: Anistropic refinement DG - saddle point problem

2020-08-17 Thread jfgir...@gmail.com
Dear Comunity, I would like to open again the topic with another question. Is there any way to use the Meshworker to solve the DG formulation but using block matrix and vectors? I couldn't find a proper way to do it with the Meshworker because I have to choose the component of the shape func