Re: [deal.II] Re: A problem about cylindrical coordinates.

2023-03-23 Thread RUIJIAN LU
LU 在2023年3月24日星期五 UTC+8 00:46:33 写道: > On 3/23/23 10:32, RUIJIAN LU wrote: > > > > Thanks, you are right. Therefore, when I just create a triangulation by > > cylinder(), have the coordinates of its vertices already been in > > cylindrical coordinate(r, phi, z) sy

[deal.II] Re: A problem about cylindrical coordinates.

2023-03-23 Thread RUIJIAN LU
Hi Dr. Bangerth, Thanks for your message. I think the error is like Marc said. And if I call triangulation.refine_global() before set_manifold(), the error will go away. Best, LU 在2023年3月24日星期五 UTC+8 00:32:04 写道: > Hi Marc! > Thanks, you are right. Therefore, when I just create a triangulat

[deal.II] Re: A problem about cylindrical coordinates.

2023-03-23 Thread RUIJIAN LU
Hi Marc! Thanks, you are right. Therefore, when I just create a triangulation by cylinder(), have the coordinates of its vertices already been in cylindrical coordinate(r, phi, z) system?(Of course, the 'z' should be x.) If it is right, is fe_values.quadrature_point(q_index)(0) equal to ‘r’, fe

Re: [deal.II] A problem about cylindrical coordinates.

2023-03-23 Thread RUIJIAN LU
I think the error is like Marc said. And if I call refine_globle() before set_manifold(), the error will go away. Best, LU 在2023年3月23日星期四 UTC+8 11:48:50 写道: > On 3/22/23 01:47, RUIJIAN LU wrote: > > > > I am a beginner to dealii, now I want to solve an equation on a

[deal.II] Re: A problem about cylindrical coordinates.

2023-03-22 Thread RUIJIAN LU
Oh, the error occurs after calling triangulation.refine_global(2); I 在2023年3月22日星期三 UTC+8 15:47:47 写道: > Hi everyone, > I am a beginner to dealii, now I want to solve an equation on a > cylindrical coordinate system. I create a triangulation by cylinder() and > use CylindricalManifold for it, b

[deal.II] A problem about cylindrical coordinates.

2023-03-22 Thread RUIJIAN LU
Hi everyone, I am a beginner to dealii, now I want to solve an equation on a cylindrical coordinate system. I create a triangulation by cylinder() and use CylindricalManifold for it, but it has some error. How can I do with it? error: Number dealii::Physics::VectorRelations::signed_angle(const