Re: [deal.II] Re: Process finished with exit code 9 ?

2019-05-03 Thread Phạm Ngọc Kiên
Thank you very much. I will take by myself. I don't think that I am running out of memory because my computer can solve larger problem previously. Or somewhere I get the memory leakage. Best Pham Ngoc Kien Vào Th 5, 2 thg 5, 2019 vào lúc 19:40 đã viết: > Hi, > > Is it UMFPACK that returns the

Re: [deal.II] SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Wolfgang Bangerth
On 5/3/19 3:34 PM, Hamed Babaei wrote: > > I believe it will be less time-consuming if I simply avoid using the > SymmetricTensor class and use the Tensor class instead and insert the > symmetric elements manually. For sure. And for once I'm not even going to argue that this might be a nice

Re: [deal.II] Non-tensor polynomials with bilinear elements

2019-05-03 Thread Wolfgang Bangerth
On 5/3/19 4:14 PM, Doug wrote: > Ok, so my initial understanding of deal.II was accurate. It all makes sense > now. Although FE_DGP is spanning {1,x,y} in reference space, it is not > spanning those bases in physical space. Instead, it is spanning some > combination of {1,x,y,xy} with 3 degrees of

[deal.II] FESystem for non-mixed DG finite element

2019-05-03 Thread Doug
Hello again, I am now looking into adding vector-valued equations to my current framework. The goal will be to do some hp-adaptive Navier-Stokes flow simulation using Discontinuous Galerkin. I have read "Handling vector valued problems", went through most of the vector-valued examples, and saw

Re: [deal.II] Non-tensor polynomials with bilinear elements

2019-05-03 Thread Doug
Prof. Bangerth, Ok, so my initial understanding of deal.II was accurate. It all makes sense now. Although FE_DGP is spanning {1,x,y} in reference space, it is not spanning those bases in physical space. Instead, it is spanning some combination of {1,x,y,xy} with 3 degrees of freedom, that only

Re: [deal.II] SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Hamed Babaei
Dear Wolfgang, I believe it will be less time-consuming if I simply avoid using the SymmetricTensor class and use the Tensor class instead and insert the symmetric elements manually. Thanks anyway for your incredible support! -- The deal.II project is located at http://www.dealii.org/ For

Re: [deal.II] SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Wolfgang Bangerth
On 5/3/19 2:07 PM, Hamed Babaei wrote: > > $ *SymmetricTensor<6, dim> C; * > > > The error is: > > /home/hbabaei/deal.ii-candi/deal.II-v9.0.1/include/deal.II/base/symmetric_tensor.h:563:29: > > error: incomplete type > ‘dealii::internal::SymmetricTensorAccessors::StorageType<6, 3, double>’

[deal.II] Re: SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Hamed Babaei
Dear Bruno, Thank you very much for your hint. I added the extra include but it did not help and I am still getting the same error! I am wondering what other things I may want to try? Thanks -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Hamed Babaei
Dear Bruno, Thank you very much for your hint. I added the extra include but it did not help and I am still getting the same error! I am wondering what other things I may want to try? Thanks On Friday, May 3, 2019 at 3:24:04 PM UTC-5, bruno@gmail.com wrote: > > Hi, > > Try to add #include

[deal.II] Re: SymmetricTensor rank 6 implementation issue

2019-05-03 Thread bruno . turcksin
Hi, Try to add #include We cannot do explicit instantiation for every thing so if you have something that is less common you need the extra include. Best Bruno On Friday, May 3, 2019 at 4:07:56 PM UTC-4, Hamed Babaei wrote: > > Hello, > > I am planning to use a symmetric tensor rank-6 for

[deal.II] SymmetricTensor rank 6 implementation issue

2019-05-03 Thread Hamed Babaei
Hello, I am planning to use a symmetric tensor rank-6 for third-order elastic constants in a nonlinear elasticity code. The tensor is supposed to be symmetric within each three pair indices as well as all possible orders of pair permutations. Namely, C_ijklmn = C_jiklmn = C_ijlkmn = C_ijklnm

Re: [deal.II] Non-tensor polynomials with bilinear elements

2019-05-03 Thread Wolfgang Bangerth
On 5/2/19 7:26 PM, Doug wrote: > > Based on this, I assume that FEValues, which takes a FiniteElement and a > Mapping, uses them to evaluate the Jacobian field instead of using a bilinear > mapping as the documentation of Mapping suggests. So FE_DGP would only be > able > to parametrize the