Re: [deal.II] Re: dealii::SolverControl::NoConvergence with Bicgstab

2021-04-15 Thread Bruno Turcksin
Le jeu. 15 avr. 2021 à 11:57, bunel...@gmail.com a écrit : > Hi, > > thank you very much for your answer. > I have actually tried to solve using GMRES like this : > > SolverControl > solver_control(phi_system_rhs.size()*2,1e-10); > PreconditionJacobi<> preconditioner; >

[deal.II] deal.II Newsletter #162

2021-04-15 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #162. 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: #12046: Remove calls to create_triangulation_compatibility in

[deal.II] Re: dealii::SolverControl::NoConvergence with Bicgstab

2021-04-15 Thread bunel...@gmail.com
Hi, thank you very much for your answer. I have actually tried to solve using GMRES like this : SolverControl solver_control(phi_system_rhs.size()*2,1e-10); PreconditionJacobi<> preconditioner; preconditioner.initialize(phi_system_matrix, 1.0);

[deal.II] Re: dealii::SolverControl::NoConvergence with Bicgstab

2021-04-15 Thread Bruno Turcksin
Hi, Bicgstab is not guaranteed to converge especially when using a bad preconditioner. You are using Jacobi as preconditioner which is a simple but not very good preconditioner. As you can see from the error message, you performed 1959 bicgstab iterations. That's a lot and it probably

[deal.II] Reading a list from parameter file

2021-04-15 Thread Paras Kumar
Dear All, I am trying to read a list of integers from the parameter file using the dealii::Patterns::List class. However, I get a runtime error which can be reproduced using the attached MWE. I also tried the approach described (for tensor) in

[deal.II] Re: LineMinimization and p-Laplace

2021-04-15 Thread Bruno Turcksin
Julie, There is a tutorial in progress https://github.com/dealii/dealii/pull/11953 that is based on step-15 but uses KINSOL . KINSOL is a package dedicated to non-linear solve. It's probably worth using it in

[deal.II] LineMinimization and p-Laplace

2021-04-15 Thread J.Y. Merten
Dear all, I am trying to solve a p-Laplace Dirichlet problem with zero boundary and used step-15 as the base. The code works well for p=1.1-2.9 with Newton and standard linesearch procedure. But when p>=3 the program diverges, continuation method brought me to p=2.99 but not beyond. Recently I

Re: [deal.II] dealii/spack on apple m1

2021-04-15 Thread Praveen C
Thank you for sharing your experience. So it looks like we cannot yet install spack natively, right ? > On 15-Apr-2021, at 12:08 PM, Alberto Salvadori > wrote: > > 1 - After opening a terminal rosetta, installing gcc, installing deal.ii via > spack and using gfortran from gcc installation,

[deal.II] dealii::SolverControl::NoConvergence with Bicgstab

2021-04-15 Thread bunel...@gmail.com
Hi, I'm having a problem for some values of parameters in my code. I get an error "dealii::SolverControl::NoConvergence", almost instantly after the start of the solving process. The status at the end is this : "Iterative method reported convergence failure in step 1959. The residual in the

Re: [deal.II] dealii/spack on apple m1

2021-04-15 Thread Alberto Salvadori
Hi Praveen, here is my experience on apple M1 and deal.ii. 1 - After opening a terminal rosetta, installing gcc, installing deal.ii via spack and using gfortran from gcc installation, it works just fine 2 - Even in a rosetta terminal, the current dmg package does not seem to work 3 - After