[deal.II] Re: PETSc vs Trilinos for step41

2021-09-28 Thread 'peterrum' via deal.II User Group
Hi Reza, I don't think there is a particular reason. Trilinos was simply chosen. Making PETSc work like in spep-40 should be not an issue. There are other tutorials that use PETSc and others that support both PETSc and Trilinos. PM On Wednesday, 29 September 2021 at 03:42:04 UTC+2

[deal.II] PETSc vs Trilinos for step41

2021-09-28 Thread Mohammadreza Yaghoobi
Hi, I was wondering why Trilinos is used in Step41 and not Petsc? Thanks Reza -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google

Re: [deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Wolfgang Bangerth
On 9/28/21 4:24 PM, Morris Jowas wrote: That would be my pleasure Prof. I will endeavour to make a contribution once I test and validate my implementation. Outstanding! Feel free to ask here how to concretely submit patches for inclusion, but in short, this is how this is done:

Re: [deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Morris Jowas
That would be my pleasure Prof. I will endeavour to make a contribution once I test and validate my implementation. Best regards, MJ On Tuesday, September 28, 2021 at 9:37:36 PM UTC+1 Wolfgang Bangerth wrote: > On 9/28/21 11:59 AM, Morris Jowas wrote: > > Thank you very much Bruno. I am

Re: [deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Wolfgang Bangerth
On 9/28/21 11:59 AM, Morris Jowas wrote: Thank you very much Bruno. I am implementing your recommendation now. step-6 already deals with a case where the coefficient is spatially variable, but not tensor-valued. When you've got it figured out how to use tensor-valued coefficients, might you

Re: [deal.II] Calculating an integral over the boundary

2021-09-28 Thread Daniel Arndt
Kyle, Have a look at VectorTools::point_value( https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#a7be5c7eed52308898dfaad91c4cff204 ). Best, Daniel Am Di., 28. Sept. 2021 um 13:04 Uhr schrieb Kyle Schwiebert < kjsch...@mtu.edu>: > Thank you for the help. I think that has

[deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Morris Jowas
Thank you very much Bruno. I am implementing your recommendation now. Best regards, MJ On Tuesday, September 28, 2021 at 6:53:37 PM UTC+1 bruno.t...@gmail.com wrote: > MJ, > > fe_values.shape_grad() returns a Tensor<1,dim> so you can do something > like this > > Tensor<1,dim> k_grad; > for

[deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Morris Jowas
P.S. Here is the equation. Sorry it did not attach last time around: On Tuesday, September 28, 2021 at 6:24:38 PM UTC+1 Morris Jowas wrote: > Hello everyone, > > > > I am solving the heat equation with the Laplacian expressed in terms of > spatially varying thermal conductivities (k_x, k_y

[deal.II] Re: Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Bruno Turcksin
MJ, fe_values.shape_grad() returns a Tensor<1,dim> so you can do something like this Tensor<1,dim> k_grad; for (int d=0; d Hello everyone, > > > > I am solving the heat equation with the Laplacian expressed in terms of > spatially varying thermal conductivities (k_x, k_y and k_z) as shown

[deal.II] Assigning anisotropic thermal conductivities in the assemble loop

2021-09-28 Thread Morris Jowas
Hello everyone, I am solving the heat equation with the Laplacian expressed in terms of spatially varying thermal conductivities (k_x, k_y and k_z) as shown below: Can I please request your suggestions on how I may go about assigning the three thermal conductivities to every cell in the

Re: [deal.II] Calculating an integral over the boundary

2021-09-28 Thread Kyle Schwiebert
Thank you for the help. I think that has the pretty much sorted out for me. One additional question: On this test problem people also like to know the pressure drop across the obstacle. This basically involves calculating the pressure at two points in the mesh and finding their difference. I

[deal.II] deal.II Newsletter #184

2021-09-28 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #184. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #12779: Allow to disable fast matrix-free hanging-node algorithm

[deal.II] Dirichlet BC on block matrices

2021-09-28 Thread Мария Бронзова
Deal all, I have an issue implementing Dirichlet boundary conditions in a mixed pressure-displacement problem. In step-20 there is an example for implementing pressure boundary conditions like that: *template * *class PressureBoundaryValues : public Function