Re: [deal.II] Linear operators - TrilinosWrappers

2021-06-16 Thread Juan Felipe Giraldo
it, as op_M then acts like a preconditioner for M, which > I don’t think is what you want.) > > I didn’t try to run the code after checking that it compiled, but I hope > that this now sorts things out for you. If not, just let us know. > > Best, > Jean-Paul > > On 16. Jun 2021

Re: [deal.II] Linear operators - TrilinosWrappers

2021-06-16 Thread Juan Felipe Giraldo
working > scenario and the problematic one. That way we could investigate further and > try to figure out what's going wrong here. > > Best, > Jean-Paul > > Sent from my mobile device. Please excuse my brevity and any typos. > > On Wed, 16 Jun 2021, 08:50 Juan Felipe Girald

Re: [deal.II] Linear operators - TrilinosWrappers

2021-06-15 Thread Juan Felipe Giraldo
change you’d be able to get your > program running. If not, then please do let us know so that we can try to > help further. > > Best, > Jean-Paul > > > On 14. Jun 2021, at 19:18, Juan Felipe Giraldo wrote: > > Hello everyone, > > I have implemented a resid

[deal.II] Linear operators - TrilinosWrappers

2021-06-14 Thread Juan Felipe Giraldo
Hello everyone, I have implemented a residual-minimization framework that somehow is similar to DPG. I want to extend my results by using parallelization using MPI with PETSc or Trilinos. So far, I have solved the saddle point problem using the Schur complement exactly how it is described in s

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] Re: Anistropic refinement DG - saddle point problem

2020-08-11 Thread Juan Felipe Giraldo
y bad consequences. So commenting the assert in > /include/deal.II/meshworker/loop.h line 357 and recompiling deal.II should > be safe and fix your problem. > > Best, > > Bruno > > On Tuesday, August 11, 2020 at 6:57:24 AM UTC-4, Juan Felipe Giraldo wrote: >> >> Dea

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

2020-08-11 Thread Juan Felipe Giraldo
Dear community, I am working on an adaptive stabilized finite element method which consist in a saddle point problem to obtain: - A continuous solution - A discontinuous error estimator (which I use as an adaptive refinement). I successfully implemented the method using the sample of the st

Re: [deal.II] get_dof_indices() only works on active cells - blocks vector

2019-12-09 Thread Juan Felipe Giraldo
ion of deal.II), or > in the tests/meshworker/mesh_loop_*.cc files and > tests/meshworker/scratch_data_*.cc > > Best, > Luca. > > > On 9 Dec 2019, at 4:55, Juan Felipe Giraldo > wrote: > > > > Dear David, > > > > Thank you for your quick

Re: [deal.II] get_dof_indices() only works on active cells - blocks vector

2019-12-08 Thread Juan Felipe Giraldo
ghbor() is inactive. > > This is only a guess. Could you show us the part of your code that > triggers this exception? That will help us give you a better answer :) > > Thanks, > David > > On Fri, Dec 6, 2019 at 3:53 AM Juan Felipe Giraldo > wrote: > &g

[deal.II] get_dof_indices() only works on active cells - blocks vector

2019-12-06 Thread Juan Felipe Giraldo
e of years ago, but I have not seen a proper solution to solve it. Anyone knows how can I solve it? Thank you so much, Regards, Juan Felipe Giraldo -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=e