Re: [deal.II] singularity error due to 1/0

2016-09-09 Thread Wolfgang Bangerth
On 09/08/2016 10:17 PM, hanks0...@gmail.com wrote: Thanks for your answer. But, I still don't know how I can deal with this error. Then you need to describe more clearly what you want to do. As I mentioned here... > K_Inv[0][0]=2*r/(cos(theta)*cos(theta)), > K_Inv[1][1]=2*r/(sin(the

Re: [deal.II] singularity error due to 1/0

2016-09-09 Thread Daniel Arndt
Kyusik, Can you summarize what you have tried and which errors/problems you got with these approaches? What exactly do you mean with "csc(theta) or sec(theta) can't be used in cylindrical coord system"? Best, Daniel Am Freitag, 9. September 2016 06:17:44 UTC+2 schrieb hank...@gmail.com: > > Th

Re: [deal.II] singularity error due to 1/0

2016-09-08 Thread hanks0227
Thanks for your answer. But, I still don't know how I can deal with this error. Anyway, Thank you very much. Kyusik. 2016년 9월 5일 월요일 오후 9시 47분 21초 UTC+9, Wolfgang Bangerth 님의 말: > > On 09/05/2016 06:04 AM, hank...@gmail.com wrote: > > > > K_Inv[0][0]=2*r/(cos(theta)*cos(theta)), > > K_Inv[

Re: [deal.II] singularity error due to 1/0

2016-09-05 Thread Wolfgang Bangerth
On 09/05/2016 06:04 AM, hanks0...@gmail.com wrote: K_Inv[0][0]=2*r/(cos(theta)*cos(theta)), K_Inv[1][1]=2*r/(sin(theta)*sin(theta)), K_Inv[2][2]=r r is calculated by "r=sqrt(x^2+y^2) ", and theta is calculated by "theta=x/y" But, as you can expect that ,on the points where |x|<0.001 or |y|