Re: [deal.II] Using PetscWrappers::TimeStepper with LinearAlgebra::distributed::Vector

2024-10-02 Thread Wolfgang Bangerth
On 10/2/24 08:17, Praveen C wrote: The documentation says https://dealii.org/current/doxygen/deal.II/classPETScWrappers_1_1TimeStepper.html VectorType must provide a petsc_vector() function but the LA distributed vector does not seem to have this, right ? In C++20 we actually check which

Re: [deal.II] Using PetscWrappers::TimeStepper with LinearAlgebra::distributed::Vector

2024-10-02 Thread Praveen C
Hello Marc The documentation says https://dealii.org/current/doxygen/deal.II/classPETScWrappers_1_1TimeStepper.html VectorType must provide a petsc_vector() function but the LA distributed vector does not seem to have this, right ? best praveen > On 2 Oct 2024, at 7:26 PM, Marc Fehling wrote

[deal.II] Re: Using PetscWrappers::TimeStepper with LinearAlgebra::distributed::Vector

2024-10-02 Thread Marc Fehling
Hi Praveen, PETScWrappers::TimeStepper has only been explicitly instantiated for the PETScWrappers Vectors: https://github.com/dealii/dealii/blob/25495718f4bb6aa45a8f6c8cc4f2aedbe8d62231/source/lac/petsc_ts.cc#L82-L89 Will deal.II compile if you add the corresponding instantiations for the Line

Re: [deal.II] Evaluate FEM solution at any point in domain

2024-10-02 Thread Peter Munch
I guess what you need is something like VectorTools::point_values() (https://www.dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#a4638b2ccdbf190a9b1f6bbbee98a2978) and VectorTools::point_gradients(https://www.dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#aa96d29e

Re: [deal.II] Evaluate FEM solution at any point in domain

2024-10-02 Thread Praveen C
May be this does what you are asking for https://dealii.org/current/doxygen/deal.II/classFunctions_1_1FEFieldFunction.html best praveen > On 2 Oct 2024, at 6:29 PM, 'Xuefeng Li' via deal.II User Group > wrote: > > I apologize for hitting the Send button before completing my questions. > > I

Re: [deal.II] Evaluate FEM solution at any point in domain

2024-10-02 Thread 'Xuefeng Li' via deal.II User Group
I apologize for hitting the Send button before completing my questions. I realize that the evaluation of a FEM function is an expensive process due to the fact that we'll need to determine the cell in the mesh where a point belongs before interpolation takes place. That's why I'm just asking for t