Re: [deal.II] Error reading a comsol mphtxt file

2022-10-24 Thread 'yy.wayne' via deal.II User Group
That's great. Thanks! :) 在2022年10月24日星期一 UTC+8 23:36:44 写道: > On 10/16/22 01:45, 'yy.wayne' via deal.II User Group wrote: > > *** Caution: EXTERNAL Sender *** > > > > I'm currently using deal.II 9.4.0 version (virtual box 1.22). When I > > read a comsol mphtxt mesh(attached below, a 2D plane in

Re: [deal.II] Elasticity probelm: same matrices - different solution vector.

2022-10-24 Thread Wolfgang Bangerth
On 10/24/22 13:16, Мария Бронзова wrote: *** Caution: EXTERNAL Sender *** Dear all, I have been trying to implement orthotropic elasticity (similar as in step-18) and in order to verify my results I decided to set all E and G modulus values as in an isotropic case and to compare the outcome w

[deal.II] Elasticity probelm: same matrices - different solution vector.

2022-10-24 Thread Мария Бронзова
Dear all, I have been trying to implement orthotropic elasticity (similar as in step-18) and in order to verify my results I decided to set all E and G modulus values as in an isotropic case and to compare the outcome with the one for the same isotropic elastic problem but with another formulat

Re: [deal.II] large condition number after solving linear systems

2022-10-24 Thread Simon Wiesheier
Sorry that I come back to this again, but I would appreciate your feedback to the issue I am confronted with: " You would generally choose the (normalizes) eigenvectors of the J^T J you have as the unit directions and the new set of parameters are then multipliers for these directions. (So far you

Re: [deal.II] Error reading a comsol mphtxt file

2022-10-24 Thread Wolfgang Bangerth
On 10/16/22 01:45, 'yy.wayne' via deal.II User Group wrote: *** Caution: EXTERNAL Sender *** I'm currently using deal.II 9.4.0 version (virtual box 1.22). When I read a comsol mphtxt mesh(attached below, a 2D plane in 3D space), it throws an error on grid_in.cc line 1557 that s!="4 Mesh". I

Re: [deal.II] Hanging node constraints in FE_NedelecSZ

2022-10-24 Thread Vinta Yinda Wang
Thank you Sebastian and thank you for your new contribution. I will keep on focusing on your paper for better understanding of this work! So if the 3-D cases are consisted of structed meshes, the oritention is, all-appearance, not conflicted, the computation of such mesh works as well as the ada

Re: [deal.II] Hanging node constraints in FE_NedelecSZ

2022-10-24 Thread Sebastian Kinnewig
Yes, the FE_Nedelec elements will work for not-so-complicated meshes, for example, a simple cube will work just fine. Best, Sebastian Vinta Yinda Wang schrieb am Montag, 24. Oktober 2022 um 10:52:19 UTC+2: > Thank you Wolfgang and thank you Alexander! > In the Todo list of FE_Nedelec, 3-D cases

Re: [deal.II] Hanging node constraints in FE_NedelecSZ

2022-10-24 Thread Vinta Yinda Wang
Thank you Wolfgang and thank you Alexander! In the Todo list of FE_Nedelec, 3-D cases ought to be carefully treated, because (from my understanding) some complex geometries in 3-D with non-standard grids (maybe unstructed or distorted initially) will introduce sign conflict issues present in tra

Re: [deal.II] Hanging node constraints in FE_NedelecSZ

2022-10-24 Thread Sebastian Kinnewig
Hello everyone, I recently finished my implementation of hanging nodes for the FE_NedelecSZ class with arbitrary polynomial degrees for the 2-D and 3-D cases. I will submit my implementation soon to deal.II on GitHub, but as the implementation was quite cumbersome and shall be part of my disse

Re: [deal.II] Hanging node constraints in FE_NedelecSZ

2022-10-24 Thread Alexander
Dear Vinta for 2-D, you could consider resorting to FE_Nedelec, which supports hanging nodes. In contrast to the 3-D case, FE_Nedelec should work even for (some) non-standard grids in 2-D. Best Alexander On Monday, October 24, 2022 at 4:46:32 AM UTC+2 Wolfgang Bangerth wrote: > On 10/23/22 0