[deal.II] Re: Q: deserialize after updating dealii

2023-11-06 Thread Najwa Alshehri
Hello Bruno, I have the solution serialized with the old version. Do you think this is the issue? I will do everything with the new version today and update you if it works. Best, Najwa On Monday, November 6, 2023 at 5:33:11 PM UTC+3 bruno.t...@gmail.com wrote: > Najwa, > > Are you

[deal.II] Re: Q: deserialize after updating dealii

2023-11-07 Thread Najwa Alshehri
UTC+3 Najwa Alshehri wrote: > Hello Bruno, > > I have the solution serialized with the old version. Do you think this is > the issue? I will do everything with the new version today and update you > if it works. > > Best, > Najwa > > On Monday, November 6, 20

[deal.II] Q: deserialize after updating dealii

2023-11-06 Thread Najwa Alshehri
Dear developers and group members, I was using an old version of Dealii where I was serializing some solutions and deserializing them where needed. In the code where I deserialize the solution I used to use the following headers: #include #include #include #include #include and I used

[deal.II] Q: Question about extracting part of a vector

2023-07-19 Thread Najwa Alshehri
Dear group members, I have one question, I am trying to calculate the L2 norm of the error of a solution. In particular, I have a *vector of the solution* u_omega and a *FeSystem* with two fes. I am interested in computing the L2 norm of the error related to the first fe using

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-20 Thread Najwa Alshehri
> for all vector components uniformly. > > Best, > Daniel > > On Wed, Jul 19, 2023 at 8:01 AM Najwa Alshehri > wrote: > >> Dear group members, >> >> I have one question, >> >> I am trying to calculate the L2 norm of the error of a sol

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-20 Thread Najwa Alshehri
s that the mask is not really picking up the first component of the solution u_omega2. Do you have any suggestions to fix this issue? Appreciate your help, Najwa On Thursday, July 20, 2023 at 11:48:04 AM UTC+3 Najwa Alshehri wrote: > Thank you Daniel for the clear quick answer. I will follow it.

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-31 Thread Najwa Alshehri
Wolfgang Bangerth wrote: > On 7/26/23 12:17, Najwa Alshehri wrote: > > > > Regarding the function "VectorFucntionFromScalarFunctionObject," I have > > observed that it works for computing the L2 norm of the error. However, > when I > > attempted to compute

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-31 Thread Najwa Alshehri
answers, >> Najwa >> >> On Wednesday, July 26, 2023 at 11:41:15 PM UTC+3 Wolfgang Bangerth wrote: >> >>> On 7/26/23 12:17, Najwa Alshehri wrote: >>> > >>> > Regarding the function "VectorFucntionFromScalarFunctionObject," I &

Re: [deal.II] Q: Question about extracting part of a vector

2023-08-03 Thread Najwa Alshehri
r. I will try to write a patch to the deal.II >>>> sources that implement the missing function. This might require some time. >>>> I will come back here if needed. >>>> >>>> >>>> Thank you all for your quick answers, >>>>

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-26 Thread Najwa Alshehri
"FeFeildFunction." Could you kindly advise me on the simplest way to accomplish this? On a side note, although "VectorFucntionFromScalarFunctionObject" works with the L2 norm of the error, it is quite computationally intensive. Best regards, Najwa On Sunday, July 23, 2023 at 4

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-23 Thread Najwa Alshehri
my_fun([&](const auto ){ return > fe_field_func.value(p);}, 0, 2); > > Luca > > Il giorno 21 lug 2023, alle ore 22:05, Najwa Alshehri > ha scritto: > > Dear Daniel, your answer makes sense, Finally it worked. > > > Dear Luca, thank you for your answer. In

Re: [deal.II] Q: Question about extracting part of a vector

2023-07-21 Thread Najwa Alshehri
ion to allocate two components: > > ExactSolution2::ExactSolution2( … ) > : Functions::Function(2) > > Best, > Luca. > > > > On 20 Jul 2023, at 12:15, Najwa Alshehri wrote: > > > > Hello again, > > > > I have a follow-up question. Does th

[deal.II] Issue with intersection of two meshes with curved boundaries

2023-11-24 Thread Najwa Alshehri
Dear developers and users, I have two meshes one is immersed in the other. I wanted to find the intersection between the two meshes, so I used the following function. NonMatching::compute_intersection(omega_grid_tools_cache,

[deal.II] Discrepancy in Computing the Inf-Sup Constants Using Matlab and deal.II

2024-05-30 Thread Najwa Alshehri
Dear all, Thank you for your always support and help. I am encountering an issue while solving an eigenvalue problem related to the computation of the discrete inf-sup constant for a saddle point problem. Specifically, I am solving the following system: Bt A^-1 B eigenvector = eigenvalue

Re: [deal.II] Discrepancy in Computing the Inf-Sup Constants Using Matlab and deal.II

2024-05-31 Thread Najwa Alshehri
; A x = \lambda M x > > > Best, > Matthias > > > On Thu, May 30, 2024, at 10:55 CDT, Najwa Alshehri > wrote: > > > Dear all, > > > > > > Thank you for your always support and help. > > > > > > I am encountering an issue while solvi

Re: [deal.II] Discrepancy in Computing the Inf-Sup Constants Using Matlab and deal.II

2024-06-01 Thread Najwa Alshehri
different results, I would greatly appreciate your input on this matter. Best, Najwa On Friday, May 31, 2024 at 11:41:27 AM UTC+3 Najwa Alshehri wrote: > Dear Wolfgang: > > Thank you for your answer. I have checked the matrices and I printed the > case for a 2 by 2 mesh from

Re: [deal.II] Discrepancy in Computing the Inf-Sup Constants Using Matlab and deal.II

2024-06-02 Thread Najwa Alshehri
Wolfgang Bangerth wrote: > On 6/1/24 14:49, Najwa Alshehri wrote: > > > > I decided to solve the exact problem directly, namely AA x = \lambda M > x. To > > achieve this, I computed the inverse of the matrix AA= Bt * A^inv * B > using a > > Conjugate Gradient (C

Re: [deal.II] Discrepancy in Computing the Inf-Sup Constants Using Matlab and deal.II

2024-06-02 Thread Najwa Alshehri
you did for the mass matrix, but replacing the inverse > operator done with umfpack with one done with CG. > > Luca > > Il giorno 2 giu 2024, alle ore 09:42, Najwa Alshehri > ha scritto: > > Wolfgang and all, > > I have a positive definite matrix M, and its inverse