[deal.II] A querry about the inversion of M in step 20

2023-11-21 Thread Hussein Shaito
Hello I am working on step 20 in deal. I believe that the schur inverse calls op_M_inv multiple times is that correct? How can I know how many times the inversion of M is done in that case? I tried to write a wrapper but couldn't do it myself. Help is appreciated. Thanks in advance. Shaito --

[deal.II] Re: A querry about the inversion of M in step 20

2023-11-24 Thread Hussein Shaito
Hello I moved on to step 22 so I am using classes to construct my preconditioner now and can create a counter within the class. This works. On Tuesday, November 21, 2023 at 9:48:10 PM UTC+2 Hussein Shaito wrote: > Hello > I am working on step 20 in deal. > I believe that the schur inve

Re: [deal.II] A querry about the inversion of M in step 20

2023-11-28 Thread Hussein Shaito
ctor that has these elements as ghost entries. I have built from the source code. I pulled dealii from github and built again but I still got the same error. I am using trilinos-release-13-2-0. Are you getting the same Issue? Hussein Shaito On Saturday, November 25, 2023 at 1:57:38 AM U

[deal.II] Seeking Advice on Implementing Solver with Trillinos and PETSc Vectors

2025-02-21 Thread Hussein Shaito
Dear community, I am trying to write a solver and I want to support both trillinos and petsc vectors at runtime. I am doing this by having the vector type as a class template for my solver. Since you guys are more experienced than I am, is there a disadvantage for this approach? Will this bloat

[deal.II] Re: Seeking Advice on Implementing Solver with Trillinos and PETSc Vectors

2025-02-22 Thread Hussein Shaito
lop our matrix > free solver. > > TLDR: Use a template parameter, compile both versions, you won't suffer > from it. > > Feel free to reach out if you have any questions. > > On Friday, February 21, 2025 at 6:48:21 a.m. UTC-5 Hussein Shaito wrote: > >> Dear comm