[deal.II] Installation problem using spack

2022-08-09 Thread Marco Nawijn
Dear All, I run into an issue when trying to install dealii using spack. The error is as follows: Could not find the hdf5 library! Insufficient hdf5 installation found! hdf5 has to be configured with MPI support. Please ensure that a suitable hdf5 library is installed on your computer

[deal.II] Component dof index

2022-08-09 Thread Uclus Heis
Dear all, Is it possible to obtain the index of the total dof distribution for each component? I work with two components, having a 2*N size system and I need to know where in my solution vector is each component. I found functions like *count_dofs_per_component* or* ComponentMask *however,

[deal.II] Stiffness Matrix

2022-08-09 Thread Prashanth
Hey Everone, I have recently started studying FEM and I have some basic doubts. Hope I can get some clarity here. 1. Is the stiffness matrix defined for element or node? In a simple truss with 2 linear elements, we calculate the stiffness matrix for K(1,1), K (1,2), K(2,1), K(2,2). In K(e,n):

Re: [deal.II] FEInterfaceValues and hp::FECollection

2022-08-09 Thread Simon Sticko
Hi, Unfortunately no-one has implemented FEInterfaceValues for hp yet (I think it would be great if someone did). The only workaround I know is to NOT use FEInterfaceValues and instead: 1. Use 2 hp::FEFaceValues objects, one initialized with cell 0 and one initialized with cell 1. 2. Split t

Re: [deal.II] Installation problem using spack

2022-08-09 Thread Daniel Arndt
Marco, It appears that -- Insufficient hdf5 installation found: hdf5 has to be configured with MPI support. is the problem. HDF5 has not been compiled with MPI support. Thus, you can either try to make sure that HDF5 is built with MPI support or, in case you are not going to use HDF5 anyway, jus

Re: [deal.II] Installation problem using spack

2022-08-09 Thread Timo Heister
Your issue might be related to the bug https://github.com/dealii/dealii/issues/14065 (see the discussion and the linked PR that should solve the problem for deal.II master, but not 9.4) On Tuesday, August 9, 2022 at 1:23:55 PM UTC-4 d.arnd...@gmail.com wrote: > Marco, > > It appears that > > --

Re: [deal.II] Component dof index

2022-08-09 Thread Wolfgang Bangerth
On 8/9/22 06:22, Uclus Heis wrote: Is it possible to obtain the index of the total dof distribution for each component? I work with two components, having a 2*N size system and I need to know where in my solution vector is each component. You want to look at DoFTools::extract_dofs(). It re

Re: [deal.II] Stiffness Matrix

2022-08-09 Thread Wolfgang Bangerth
On 8/9/22 06:53, Prashanth wrote: 1. Is the stiffness matrix defined for element or node? In a simple truss with 2 linear elements, we calculate the stiffness matrix for K(1,1), K (1,2), K(2,1), K(2,2). In K(e,n): e is the element number and n is the node number. So is the stiffness matrix defi

Re: [deal.II] fourth-order referential deviatoric tensor Dev_P: why is 'S' used in the definition?

2022-08-09 Thread Wolfgang Bangerth
On 8/8/22 15:11, Simon Wiesheier wrote: grafik.png grafik.png The only place I use Dev_P is to compute the above double contractions ('P' in the above is what dealii returns as Dev_P) As you can see, the fourth order tangent C_bar consists of outer products of second order symmetric tensors: bo

Re: [deal.II] Matrix-free DG with h-adaptive mesh refinement

2022-08-09 Thread Wolfgang Bangerth
On 8/8/22 02:44, M. Bänsch wrote: *** Caution: EXTERNAL Sender *** Dear community, I am currently trying to implement h-adaptive mesh refinement into step-67 (by adding a refine_grid() functionality). As references I looked into step-31