Re: [deal.II] Linear operators - TrilinosWrappers

Dear Felipe, It might be that you need to set up the preconditioner operator with an exemplar matrix (since the identity operator doesn't know the size of the range and domain that it's working on). If that's not the issue then could you please try to reproduce this as a minimal example and share

Re: [deal.II] Linear operators - TrilinosWrappers

Dear Jean-Paul, Thank you for your reply, I already set up the Identity preconditioner with a block matrix, but the problem remains the same. I am not sure this is the problem because I didn't need to set up the identity operator to any matrix in the scenario where I could invert. I think my p

[deal.II] Neumann data and Kelly estimator: can't pass std::map argument

Hi everyone, I'm trying to change a bit the step-7 (Helmotz Problem ) In particular, I'm going to change the estimate function KellyErrorEstimator::estimate

Re: [deal.II] Linear operators - TrilinosWrappers

Dear Felipe, Changing this const auto op_pre = linear_operator(prec_M); to this const auto op_pre = linear_operator(M, prec_M); allowed your program "FEMDG-Dar-par_tofix" to compile for me. If you have a look at this documentation for this variant of linear_operator https://dealii.org/current/dox

Re: [deal.II] Linear operators - TrilinosWrappers

Dear Jean-Paul, Finally it is working! Thank you for your help and clear explanation! Now it is a bit more clear for me how these linear operators work here :). (Regarding the op_M operator, I was not using that variable because I declared M directly as the linear operator of the block(0,0). B

[deal.II] Re: online deal.II workshop: Friday, June 18

Hi all, This is a quick reminder that our one-day workshop will happen this Friday. We hope to see you there and celebrate the 9.3 release with us! The information on how to join the meeting via Zoom is now available at https://dealii.org/workshop-2021/ Best, Timo On Wed, Jun 2, 2021 at 1

[deal.II] Local Raviart-Thomas space

Hello, I am trying to implement a projection of a DG solution to a Raviart-Thomas space. For this, I need a local space Q_{p-1,p}(K) \times Q_{p,p-1}(K). I have been using FE_DGRaviartThomas of degree p-1, and then I use shape_value_component(i, q, 1) when I want to use the x component of the

Re: [deal.II] Re: Integrated material and spatial traction forces on boundary not equal

Hi Alex, I just got a chance to have Trilinos installed but the same error occurs regarding *EnergyFunctional*. I got the following notice when trying to test *step-31*: Error! This tutorial requires a deal.II library that was configured with the following options: DEAL_II_WITH_TRILINOS

Re: [deal.II] Re: Integrated material and spatial traction forces on boundary not equal

On 6/16/21 4:07 PM, Michael Lee wrote: Is there a way to let deal.II know Trillinos is installed? Do I have to reinstall deal.ii? Or, can I compile the code without using Trillinos with some modification? No, you have to re-install deal.II so that at compile time, it knows that Trilinos exi

Re: [deal.II] Local Raviart-Thomas space

On 6/16/21 2:00 PM, Giselle Sosa Jones wrote: I am trying to implement a projection of a DG solution to a Raviart-Thomas space. For this, I need a local space Q_{p-1,p}(K) \times Q_{p,p-1}(K). I have been using FE_DGRaviartThomas of degree p-1, and then I use shape_value_component(i, q, 1) wh

Re: [deal.II] Cahn-Hilliard Problem

On 6/16/21 12:21 PM, Kubra Karayagiz (Alumni) wrote: Could you please let me know if my understanding of the usage of this function is correct? (Note: the use of compute_nl_term() is necessary for the implicit-explicit implementation, due to the nonlinear functions) Kubra, much debugging wil

Re: [deal.II] Neumann data and Kelly estimator: can't pass std::map argument

On 6/16/21 4:29 AM, bob bill wrote: *The boundary indicator for Neumann data is 1, and I gave a pointer to a constant function as argument. **Unfortunately, the compiler says: "*no matching constructor for initialization of 'std::mapFunction<2> *>' (aka 'map *>') std::map *>(1,&neumann_func),

Re: [deal.II] Cahn-Hilliard Problem

Prof. Bangerth, Thank you for the prompt response. On Wed, Jun 16, 2021 at 7:33 PM Wolfgang Bangerth wrote: > On 6/16/21 12:21 PM, Kubra Karayagiz (Alumni) wrote: > > > > Could you please let me know if my understanding of the usage of this > function > > is correct? (Note: the use of compute_

[deal.II] Extracting a component of FESystem for use elsewhere

Hi everyone, This might be a simple question, but I haven't found a straightforward answer in the "handling vector-valued problems" documentation. If I have a Vector which stores the solution to an FESystem, say 3 FE_DGQ fields for u, v, p, and I want to extract a new Vector which corresponds

[deal.II] Re: Extracting a component of FESystem for use elsewhere

I've found one way to do this, but I imagine it's not optimal: - call DoFRenumbering::component_wise(fesys_dofh) - Do the finite element solve - then copy from one vector to another using the number of dofs in that component and an offset - This isn't great, however, because I'd ra

Re: [deal.II] Re: Extracting a component of FESystem for use elsewhere

Hi Corbin, The cell-wise solution corresponding to the p field could be determined as follows: ``` dealii::FEValuesExtractorScalar pressureExtractor(2); //assuming you have three scalar fields and pressure is the last one dealii::Vector solVecPerCell(nDOFsPerCell); feValues[pressureExtractor].get_