Re: [deal.II] Error in the second derivative of solution...?

2017-05-17 Thread Jean-Paul Pelteret
Thats great! I'm glad that you've managed to work it out :-) Best, Jean-Paul On Wednesday, May 17, 2017 at 10:06:18 AM UTC+2, hanks0...@gmail.com wrote: > > Thank you very much!! > > It seems that I can solve my problem. > > Thank you! > > Kyusik. > -- The deal.II project is located at

Re: [deal.II] Error in the second derivative of solution...?

2017-05-17 Thread hanks0227
Thank you very much!! It seems that I can solve my problem. Thank you! 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

Re: [deal.II] Error in the second derivative of solution...?

2017-05-16 Thread Jean-Paul Pelteret
Hi Kyusik, What makes me confused here is that, according to your explanation, two > different DoFHandler could be used simultaneously. > Yes, that it correct. You can convert between cell iterators. So I do the following: typename Triangulation::active_cell_iterator cell

Re: [deal.II] Error in the second derivative of solution...?

2017-05-16 Thread hanks0227
Dear Jean-Paul, Thank you again for a detailed explanation. However, please excuse my ignorance... > > 1. Create a vector-valued finite element system (maybe an > FESystem(FE_Q(1),dim) ) and use it to distribute DoFs onto a new > DoFHandler. > 2. Create a mass matrix "M" with the new

Re: [deal.II] Error in the second derivative of solution...?

2017-05-15 Thread Jean-Paul Pelteret
Dear Kyusik, > I'm sorry to tell you that actually I don't know how to project > grad_solution into the vector-valued FE space > > Do you happen to know some example tutorial code that includes it? > Off the top of my head, no I don't know of an example tutorial that demonstrates it.

Re: [deal.II] Error in the second derivative of solution...?

2017-05-15 Thread hanks0227
Dear Jean-Paul Pelteret, Thank you very much for your help. It is very helpful. > However, you'll note that this function only makes sense (and its > therefore only defined) for a vector field. Your solution field is scalar > though, which presents a bit of an issue. On the face of it you'd

Re: [deal.II] Error in the second derivative of solution...?

2017-05-15 Thread Jean-Paul Pelteret
Dear Kyusik, So my new question is , in deal.II, could I get local divergence of certain > function to get surface integral? > >> Yes it is. You can extract the divergence of a function using, for example, the FEValuesView function get_function_divergences

Re: [deal.II] Error in the second derivative of solution...?

2017-05-14 Thread hanks0227
Dr. Wolfgang Bangerth, Thank you very much for your kind reply. But, could I have one more question? > > I don't know whether the problem really comes from that, but you can't > take > the second derivative of a finite element solution without problem. The > issue > is that it exists,

Re: [deal.II] Error in the second derivative of solution...?

2017-05-14 Thread Wolfgang Bangerth
On 05/14/2017 01:57 AM, hanks0...@gmail.com wrote: I guess this problem comes from the second derivatives of solution. I don't know whether the problem really comes from that, but you can't take the second derivative of a finite element solution without problem. The issue is that it exists,