[deal.II] Re: shift-invert transformation problem with slepc

2021-08-23 Thread simon...@gmail.com
Hi, I think you need to set both which eigen pair and the target value. For example, if you look for the eigenvalue closest in magnitude to 0, you need to set the following const double target_eigenvalue = 0; eigensolver.set_which_eigenpairs(EPS_TARGET_MAGNITUDE); eigensolver.set_target_eigenva

[deal.II] shift-invert transformation problem with slepc

2021-08-23 Thread Anton Ermakov
Dear All, I am trying to use a shift-invert transformation in step-36. I did minimal changes in the solve function adding the following code that describes what spectral transformation I need to do: SolverControl solver_control (dof_handler.n_dofs(), 1e-9); SLEPcWrappers::SolverKrylovSchur eige

Re: [deal.II] Re: Subdivided cylinder and boundary ids

2021-08-23 Thread blais...@gmail.com
This issue is now fixed by a PR merged yesterday. Best Bruno On Monday, June 28, 2021 at 6:40:04 a.m. UTC-4 alexanderc...@gmail.com wrote: > Hi Bruno, > > Sure, I can try making a PR this week. > > Best, > Alex > > On Sunday, June 27, 2021 at 9:36:50 p.m. UTC+2 blais...@gmail.com wrote: > >> I