[deal.II] Re: Problems in handling the userdata of Distributed Triangulation

2016-11-02 Thread Jean-Paul Pelteret
Dear Jack, You'll be happy to know that this functionality has recently been added to the developer version of deal.II via the CellDataStorage and TransferableQuadraturePointData

[deal.II] Problems in handling the userdata of Distributed Triangulation

2016-11-02 Thread Jack
I add the userdata to each cell using the cell->user_pointer. It seems that the userdata of some cells which are changed during the refining and coarsening procedure will lose. How can make sure that the refined children cells(the coarsened parent cell) inherit the userdata from the parent (t

Re: [deal.II] Re: Trouble with VisIt 2.11

2016-11-02 Thread Praveen C
So the VisIt people said they forgot to include the xdmf plugin in v2.11 and we have to wait for the next release to get it back. Meanwhile an alternative is suggested, see https://elist.ornl.gov/pipermail/visit-users/2016-November/019322.html Best praveen On Tue, Nov 1, 2016 at 4:07 PM, Pravee

Re: [deal.II] GridGeneration

2016-11-02 Thread Wolfgang Bangerth
On 11/02/2016 05:54 PM, benhour.amiria...@gmail.com wrote: I want to use half_hyper_ball for my domain. Unfortunately, I cannot refine the mesh around the center. I attached my geometry below. The code that I am using has been written as follows: Point <2> center; const double radius = 30e-9;

Re: [deal.II] Re: hp convergence for non-linear solver.

2016-11-02 Thread Wolfgang Bangerth
On 11/02/2016 07:17 PM, Jaekwang Kim wrote: I'd like check the visualized error per cell. Is there any module in deal.ii that enables this...? Yes, you just need to create a vector with as many entries as there are cell (e.g., the output of VectorTools::integrate_difference) and attach that

Re: [deal.II] Re: hp convergence for non-linear solver.

2016-11-02 Thread Jaekwang Kim
Thank you for replying me! I'd like check the visualized error per cell. Is there any module in deal.ii that enables this...? 2016년 11월 2일 수요일 오전 7시 27분 56초 UTC-5, Wolfgang Bangerth 님의 말: > > On 11/01/2016 11:06 PM, Jaekwang Kim wrote: > > > > do you have any other idea? > > > > Output the

[deal.II] Re: How can I find the neighbouring element if I'm on a given face in a given cell?

2016-11-02 Thread smetcalfephd
That looks good! Thanks for the help. Stephen On Wednesday, 2 November 2016 18:36:02 UTC+1, Bruno Turcksin wrote: > Stephen, > > On Wednesday, November 2, 2016 at 1:20:58 PM UTC-4, smetca...@gmail.com > wrote: > >> I'm currently implementing a blow-up problem with space-time adaptivity >> con

[deal.II] Re: How can I find the neighbouring element if I'm on a given face in a given cell?

2016-11-02 Thread Bruno Turcksin
Stephen, On Wednesday, November 2, 2016 at 1:20:58 PM UTC-4, smetcalfe...@gmail.com wrote: > I'm currently implementing a blow-up problem with space-time adaptivity > controlled by a derived pointwise a posteriori error bound and I'm having > trouble implementing the estimator (previously, I j

[deal.II] How can I find the neighbouring element if I'm on a given face in a given cell?

2016-11-02 Thread smetcalfephd
I'm currently implementing a blow-up problem with space-time adaptivity controlled by a derived pointwise a posteriori error bound and I'm having trouble implementing the estimator (previously, I just used meshworker and it did most of the work for me but you cannot use such an approach for po

[deal.II] Re: Extra layer around mesh

2016-11-02 Thread Daniel Arndt
Joel, Same error with your changes. I searched for subdivided in > http://www.dealii.org/developer/doxygen/deal.II/changes_between_8_3_and_8_4.html > > and found under specific improvements, that subdivided_parallelepiped > produced the wrong boundary indicators. Could this be the reason, why

Re: [deal.II] Changing geometry

2016-11-02 Thread A Banerjee
I got my mistake. Thanks Atanu -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscri

Re: [deal.II] Changing geometry

2016-11-02 Thread A Banerjee
Yes. I agree. This was stupid question. Thanks Atanu -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group

Re: [deal.II] Changing geometry

2016-11-02 Thread Wolfgang Bangerth
On 11/02/2016 09:24 AM, A Banerjee wrote: Let me tell you that I am a beginner in deal_II. I have a fair idea about finite element method whereas no idea about c++. While playing with the tutorial program step-6, and changing the domain to a square one between (0,1), I find a wrong solution. Whe

Re: [deal.II] Changing geometry

2016-11-02 Thread Deepak Gupta
Hi, I think (a guess) you forgot to change the Coefficient function. For [-1, 1], the radius of 0.5 is around (0, 0). For [0, 1], I guess it needs to be adjusted accordingly. Best, Deepak On Wed, Nov 2, 2016 at 4:24 PM, A Banerjee wrote: > Dear All, > > Let me tell you that I am a beginner in

Re: [deal.II] Re: Convergence problem arising for large number of DoFs

2016-11-02 Thread Bruno Turcksin
Hamed, 2016-11-02 10:22 GMT-04:00 Hamed Babaei : > It seems that my system_matrix is not s.p.d any more although in sequential > code it was. Then, there is a bug in your code. You can compare the matrix produced by the sequential code with the one in parallel do help you find the bug. Best, Bru

Re: [deal.II] Absolute values of deal.II vector

2016-11-02 Thread hanks0227
Dr. Bangerth, Thank you very much. It works. Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Gr

Re: [deal.II] Re: Convergence problem arising for large number of DoFs

2016-11-02 Thread Hamed Babaei
Dear Wolfgang, I wrote the sequential code first and It converged to the exact solution using SolverCG and SSOR preconditioner. However, the parallel code diverges at the first call to SolverCG but using GMRES or Bicgstab it works though up to 15 DoFs not more. First and the most important

[deal.II] Re: deal.II on youtube and twitter

2016-11-02 Thread Jean-Paul Pelteret
Awesome, thanks very much for sharing Praveen! On Wednesday, November 2, 2016 at 2:27:36 PM UTC+1, Praveen C wrote: > > Dear all > > I sometimes upload videos/pictures of solutions obtained with deal.II on > youtube and twitter. I am tagging with #dealii on twitter. > > https://twitter.com/search

[deal.II] Re: deal.II on youtube and twitter

2016-11-02 Thread Jean-Paul Pelteret
Awesome, thanks very much for sharing Praveen! On Wednesday, November 2, 2016 at 2:27:36 PM UTC+1, Praveen C wrote: > > Dear all > > I sometimes upload videos/pictures of solutions obtained with deal.II on > youtube and twitter. I am tagging with #dealii on twitter. > > https://twitter.com/search

[deal.II] deal.II on youtube and twitter

2016-11-02 Thread Praveen C
Dear all I sometimes upload videos/pictures of solutions obtained with deal.II on youtube and twitter. I am tagging with #dealii on twitter. https://twitter.com/search?q=%23dealii I dont know how to properly tag videos on youtube. I have some videos here made with deal.II https://www.youtube.co

Re: [deal.II] Absolute values of deal.II vector

2016-11-02 Thread Wolfgang Bangerth
On 11/02/2016 07:06 AM, hanks0...@gmail.com wrote: Anyway, what I'm trying to do is to get error = |exact_solution - numerical_solution| (absolute value of two Deal.II Vector) and then plot this one using Visit. But I don't know how to change all elements of error(Deal.II Vector) to absolute va

[deal.II] Absolute values of deal.II vector

2016-11-02 Thread hanks0227
Hi All, First of all, my question is very simple. But I can't find similar question. Anyway, what I'm trying to do is to get error = |exact_solution - numerical_solution| (absolute value of two Deal.II Vector) and then plot this one using Visit. But I don't know how to change all elements of e

Re: [deal.II] Re: Convergence problem arising for large number of DoFs

2016-11-02 Thread Wolfgang Bangerth
On 11/01/2016 10:31 AM, Hamed Babaei wrote: I mean non of the standard preconditioners helped to make it converged. Like Bruno mentioned, if the standard preconditioners don't work, then you will need to think about more complicated ones. Take a look at my video lectures on preconditioners,

Re: [deal.II] Re: hp convergence for non-linear solver.

2016-11-02 Thread Wolfgang Bangerth
On 11/01/2016 11:06 PM, Jaekwang Kim wrote: do you have any other idea? Output the errors per cell and visualize. Is there a pattern? Is the error localized on the boundary, for example? Best W. -- Wolfgang Bangerth

[deal.II] Re: Extra layer around mesh

2016-11-02 Thread Joel Davidsson
Dear Daniel, Same error with your changes. I searched for subdivided in http://www.dealii.org/developer/doxygen/deal.II/changes_between_8_3_and_8_4.html and found under specific improvements, that subdivided_parallelepiped produced the wrong boundary indicators. Could this be the reason, why i