Re: [deal.II] Standard way of applying a function to a finite element vector

2021-07-12 Thread Wolfgang Bangerth
On 7/12/21 5:56 PM, Lucas Myers wrote: I'd like to apply a method to a finite element function (that is, something that looks like the solution vectors) in order to get another finite element function. For example, suppose my finite element function is a symmetric rank-2 tensor in type, and I

[deal.II] Standard way of applying a function to a finite element vector

2021-07-12 Thread Lucas Myers
Hi all, I'd like to apply a method to a finite element function (that is, something that looks like the solution vectors) in order to get another finite element function. For example, suppose my finite element function is a symmetric rank-2 tensor in type, and I would like a finite element func

[deal.II] Re: Step-77 not working using Sundials 5.7.0

2021-07-12 Thread blais...@gmail.com
In addition, Step-77 does not compile with the Sundials 3.1.0 distributed with candi. The following error is obtained at run time in debug mode: An error occurred in line <427> of file in function unsigned int dealii::SUNDIALS::KINSOL::

[deal.II] Re: Step-77 not working using Sundials 5.7.0

2021-07-12 Thread Laura Prieto
Follow up: just tested with Sundials 3.1 and I get the same error with similar values for the norms El lunes, 12 de julio de 2021 a las 12:18:52 UTC-4, Laura Prieto escribió: > Dear deal.ii community, > > I am testing step-77 using the latest deal.ii release with sundials on > (version 5.7.0).

[deal.II] Step-77 not working using Sundials 5.7.0

2021-07-12 Thread Laura Prieto
Dear deal.ii community, I am testing step-77 using the latest deal.ii release with sundials on (version 5.7.0). The example compiles and runs but not as expected, as I get the following output: Mesh refinement step 0 Target_tolerance: 0.001 Computing residual vector... norm=0.231202 Com

Re: [deal.II] How to add 2nd iterator to MeshWorker::mesh_loop ?

2021-07-12 Thread Wolfgang Bangerth
On 7/12/21 4:18 AM, Sylvain Mathonnière wrote: _*Question:*_ How can I modify the cells worker loop (on dof1) initialisation to have a second iterator linked with a second dof_handler looping at the same time ? or how can I add a concurrentely running for loop to the cell_worker with a differ

Re: [deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

2021-07-12 Thread Wolfgang Bangerth
On 7/12/21 7:28 AM, Matthew Rich wrote: At one time did dataout have an active_cell_iterator property. The code gallery stuff may not get the same upkeep as the tutorials and maybe the library updated and broke this piece of the code. Yes, that seems to be the case. Would you be interested i

[deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

2021-07-12 Thread Bruno Turcksin
Yes, it looks like DataOut used to have DataOut::active_cell_iterator. I looked back at an old version of deal.II and DataOut::active_cell_iterator was just a typedef for Triangulation::cell_iterator. Bruno On Monday, July 12, 2021 at 9:28:50 AM UTC-4 mjri...@gmail.com wrote: > OK, > > At one

[deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

2021-07-12 Thread Matthew Rich
OK, At one time did dataout have an active_cell_iterator property. The code gallery stuff may not get the same upkeep as the tutorials and maybe the library updated and broke this piece of the code. I see both typename DoFHandler::active_cell_iterator and DataOut::active_cell_iterator, whi

Re: [deal.II] Contraction between Tensor and SymmetricTensor

2021-07-12 Thread Wolfgang Bangerth
On 7/12/21 1:08 AM, Lixing Zhu wrote: For the case (grad *w*, grad^s *u*), how should I carry out the contraction between a tensor and a symmetric tensor? It seems that operator * does not work here. The scalar_product() function taking two Tensor arguments is what you are looking for! Be

Re: [deal.II] step-35 in modified

2021-07-12 Thread Wolfgang Bangerth
On 7/12/21 12:59 AM, andrew jack urombo wrote: The error I am trying to have suggestions to fix is: //home/dealii/dealii/examples/step-35a/step-35a.cc: In member function ‘void Step35a::NavierStokesProjection::refine_mesh(unsigned int, unsigned int)’:// ///home/dealii/dealii/examples/step-35a/s

[deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

2021-07-12 Thread Bruno Turcksin
Hi, On Monday, July 12, 2021 at 5:40:40 AM UTC-4 mjri...@gmail.com wrote: > > Should I just change active_cell_iterator to cell_iterator? I am worried > that they wont work the same creating problems downstream. > Only cell_iterator is accessible not active_cell_iterator. So yes, active_cell_i

[deal.II] How to add 2nd iterator to MeshWorker::mesh_loop ?

2021-07-12 Thread Sylvain Mathonnière
Dear deal.II users, I am trying to mix step-12 and step-31 for my project. I have two coupled equations (named eq1 and eq2) that I am solving in parallel just like in step-31 but for one of them I am using the discontinuous Galerkin method like in step-12 (eq1) and for the other one I use stand

[deal.II] Compile Error: no type named ‘active_cell_iterator’

2021-07-12 Thread Matthew Rich
Hi all, Once again perplexed by the c++ compiler. I am getting the following error when trying to compile a program from the code gallery. prog1.cc:In instantiation of ‘typename dealii::DataOut::cell_iterator Prog1::FilteredDataOut::first_cell() [with int dim = 3; DH = dealii::DoFHandler<3, 3>

[deal.II] Contraction between Tensor and SymmetricTensor

2021-07-12 Thread Lixing Zhu
Dear Community, Regarding the term (grad *w*, grad *u*) where *u *is the unknown vector trial solutions and *w *is the corresponding weighting function, In the context of linear elasticity, we may have (grad^s *w*, grad^s *u*), and we utilize symgrad_phi_u[i] * symgrad_phi_u[j] as shown in Ste