[deal.II] installation issue on quarks

2018-08-27 Thread Peng Wei
Dear all, I was trying to install dealii version 9.0.0 on department cluster using candi, git clone https://github.com/dealii/candi cd candi ./candi.sh --platform=./deal.II-toolchain/platforms/supported/linux_cluster.platform but I got the following error (it worked well for the installati

Re: [deal.II] how to remove rows from system matrix

2018-08-27 Thread Loylick
Thank you, Wolfgang, for your suggestion. I've read step-11 tutorial but I still can not figure out what kind of constraints I can impose in my case. Here is a picture of a void with enriched nodes. Blue triangls are enriched with Heaviside function, red rectangles are zero degrees of freedom. Wh

Re: [deal.II] Error in running Quasi_Static Compressible example on Ubuntu

2018-08-27 Thread Wolfgang Bangerth
On 08/27/2018 08:29 AM, Lucky Qin wrote: I am trying to run the example in code gallery named 'Quasi_static_Finite_strain_Compressible_Elasticity'. However, the errors appear as shown in the picture. It seems to remind me that UMFPack is not installed. From the guidance in the installation pa

Re: [deal.II] Create some kind of anisotropic FE_Q element

2018-08-27 Thread Wolfgang Bangerth
On 08/27/2018 02:06 AM, Dustin Kumor wrote: I like to create a three dimensional Lagrange finite element that uses ansatz functions of the following form: N_i (x,y,z) = q_i(x,y) p_i (z), where q_i (x,y) are the basis functions of a two dimensional Lagrange finite element of arbitrary degree,

Re: [deal.II] how to remove rows from system matrix

2018-08-27 Thread Wolfgang Bangerth
On 08/27/2018 08:36 AM, Loylick wrote: So I get a singular matrix with zero degrees of freedom inside voids. For this system to be solved I should eliminate rows corresponding to zero degrees of freedom from my system matrix. I need your advise how to implement it in dealii in an optimal way. I

[deal.II] DFT-FE Open Source Code

2018-08-27 Thread vikramg
Dear Colleagues, I am writing to bring to your attention the open source beta release of DFT-FE. DFT-FE is a real-space density functional theory code based on adaptive finite-element basis (using the deal.II framework) for ground-state pseudopotential and all-electron calculations (including

Re: [deal.II] Re: How do I choose my spaces to enforce gradient continuity across an interface?

2018-08-27 Thread Wolfgang Bangerth
On 08/27/2018 07:52 AM, georgios.sialou...@gmail.com wrote: What I meant is for the interface condition to hold weakly on the boundary and that, as you correctly pointed out, the jump in the gradient should equal the pressure times the normal vector.  As such, I specified a pressure which I kno

[deal.II] how to remove rows from system matrix

2018-08-27 Thread Loylick
Hi! I'm trying to implement XFEM method for elasticity problem with voids. I use FE_Collection: fe_collection.push_back (FESystem (FE_Q(1), 1, FE_Nothing(), 1)); fe_collection.push_back (FESystem (FE_Q(1), 1, FE_Q(1), 1)); Triangulation remains always unchanged, I do not cut voids in the mesh.

Re: [deal.II] Re: How do I choose my spaces to enforce gradient continuity across an interface?

2018-08-27 Thread georgios . sialounas
Dear Prof. Bangerth, Thank you very much for your reply! Apologies for the ambiguity. What I meant is for the interface condition to hold weakly on the boundary and that, as you correctly pointed out, the jump in the gradient should equal the pressure times the normal vector. As such, I spec

[deal.II] Create some kind of anisotropic FE_Q element

2018-08-27 Thread Dustin Kumor
Dear all, I like to create a three dimensional Lagrange finite element that uses ansatz functions of the following form: N_i (x,y,z) = q_i(x,y) p_i (z), where q_i (x,y) are the basis functions of a two dimensional Lagrange finite element of arbitrary degree, basically FE_Q<2> (degree) and p_i a