[deal.II] Re: Help with the weak form

Hi Dragan, Did you do a literature review to see what researchers are doing in your field? Most likely you need to linearize it w.r.t. to "c" so that you will be solving for \Delta c around some points c_0 employing something like Newton-Raphson method. Regards, Denis. On Thursday, September

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

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[

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

I have tried to solve this problem including what you suggested. But,I'm sorry to tell you, Everything didn't work out well... It seems that csc(theta) or sec(theta) can't be used in cylindrical coord system But, I solve GradShaf Eq that I was originally trying to solve using step-6 Anyway, T

Re: [deal.II] Problem creating a pull request on code-gallery

Got it! I also forgot to commit but I think now I was able to make the pull request. Thanks! :) On Thursday, September 8, 2016 at 4:44:28 PM UTC-5, Wolfgang Bangerth wrote: > > > Hi Manuel, > > > I want to contribute with an small project to the code-gallery. > > Fantastic! > > > > I followed

Re: [deal.II] Problem creating a pull request on code-gallery

Hi Manuel, I want to contribute with an small project to the code-gallery. Fantastic! > I followed the instructions at (http://www.dealii.org/code-gallery.html) but I think I might have had a problem on the very last step; i.e., creating a "pull request". I went to my github account and cl

[deal.II] Problem creating a pull request on code-gallery

Hi, I want to contribute with an small project to the code-gallery. I followed the instructions at (http://www.dealii.org/code-gallery.html) but I think I might have had a problem on the very last step; i.e., creating a "pull request". I went to my github account and clicked on "New pull reque

[deal.II] Help with the weak form

Hello guys, I would need a small help with deriving the weak form. The equation without the unnecessary terms is the following: log(c/(1-c)) - k*Laplacian(c) = 0 I am not sure how to treat the term log(c/(1-c)). Thank you. Regards, Dragan -- The deal.II project is located at http://www.dea

Re: [deal.II] FE_TraceQ using GaussLobatto nodes

Hello Martin > You are right, we should allow for interpolation of initial condition. We > also allow this for FE_FaceQ. I will soon open a PR. On the other hand, we > will probably not be able to project onto FE_TraceQ with our tools because > this element does not define basis functions inside

Re: [deal.II] FE_TraceQ using GaussLobatto nodes

Dear Praveen, > > When I use FE_TraceQ and am setting initial condition with > VectorTools::interpolate, I get this error > > > > An error occurred in line <127> of file > > in function > > void dealii::VectorTools::interpolate(const Ma

Re: [deal.II] FE_TraceQ using GaussLobatto nodes

Hello Martin When I use FE_TraceQ and am setting initial condition with VectorTools::interpolate, I get this error An error occurred in line <127> of file in function void dealii::VectorTools::interpolate(const Mapping &, const DoFHa