[deal.II] Problem with user_index()

typename DoFHandler::active_cell_iterator cell = dof_handler.begin_active(), endc = dof_handler.end(); for (; cell!=endc; ++cell) { std::cout

Re: [deal.II] problem with user_index()

Dear Sudharshan, Probably you are looking for cell->index(). Am I right? best Deepak On Sun, Nov 13, 2016 at 7:03 PM, Sudarshan Kumar wrote: > > typename DoFHandler::active_cell_ > iterator > cell = dof_handler.begin_active(), > endc = dof_handler.end(); > >

[deal.II] Re: storing sparse matrix in .txt and scan in Matlab

Dear Anup, > 1. if the way of writing the sparse matrix on a .txt file is appropriate > for post-processing in matlab, if not, please suggest > a way; > There are three print functions in the deal.II SparseMatrix. As explained in the documentation

[deal.II] Re: Very basic question on setting boundary id

Jaekwang, It seems that you are missing constraints.distribute(solution); after solving the linear system. The constrained dofs are set to zero during assembly and this call is necessary to get the correct values for these in the solution as well. You might want to have a look at

[deal.II] Re: storing sparse matrix in .txt and scan in Matlab

Anup, > >> 1. if the way of writing the sparse matrix on a .txt file is appropriate >> for post-processing in matlab, if not, please suggest >> a way; >> > > There are three print functions in the deal.II SparseMatrix. As explained > in the documentation >

Re: [deal.II] problem with user_index()

Thanks a lot Gupta, It is right, working fine now. However I am wondering why the user_index() is not giving the answer. On Sunday, November 13, 2016 at 3:16:05 PM UTC-3, Deepak Gupta wrote: > > Dear Sudharshan, > > Probably you are looking for cell->index(). Am I right? > > best >

[deal.II] Re: Very basic question on setting boundary id

Yes, well this is going to give strange results because you're setting the boundary ID of internal faces (x=0.5 is the centreline of your geometry). Before you set any boundary ID's, you should always first check to see that you're actually on a boundary. The first few tutorials do go through

[deal.II] Very basic question on setting boundary id

Hi all, I am learning grid generator how to set boundary condition on given domain. It is very basic, but I couldn't even set boundary id for squared mesh I want to generate squared computational domain and I want to give boundary id of 1 to all the outer boundary of the square. So I used

[deal.II] Re: Very basic question on setting boundary id

Hi Jaekwang, With what you've shown I'm going to guess that your problem might simply be related to numerical / round-off error associated with floating point calculations. if ( cell->face(f)->center()[0]==1 ) // The LHS value is a computed double, > and the RHS value is an integer which

[deal.II] Re: Very basic question on setting boundary id

Thank you for quick reply!!! I tried your suggestions but I couldn't resolve the problem. The more weird thing is if I commands as.. GridGenerator::hyper_cube (triangulation, 0, 1); triangulation.refine_global (); for (typename Triangulation::active_cell_iterator

[deal.II] Re: Very basic question on setting boundary id

Thank you for being kindness! I followed your suggestions but I couldn't resolve the problem yet. Or did I make other mistake in another places? >From a slight modification from STEP-4 tutorial , I have been trying solve Poisson Equation. Actually, I have manufactured solution to see my

[deal.II] Re: Very basic question on setting boundary id

Problem Solved! Thank you!! 2016년 11월 13일 일요일 오후 5시 25분 20초 UTC-6, Daniel Arndt 님의 말: > > Jaekwang, > > It seems that you are missing > constraints.distribute(solution); > after solving the linear system. The constrained dofs are set to zero > during assembly and > this call is necessary to

Re: [dealii-developers] Re: automated bulk github issue creation

> Agreed. > > What do people think about a couple of days in the first half of > January? Maybe in the first week of January, and then bring out a > release right after? I'm in. -- You received this message because you are subscribed to the Google Groups "deal.II developers" group. To