[deal.II] Re: Function to output stress and strain

Hi Lisa, To project from quadrature points to (continuous) FE space you can use this function in current master branch of GitHub: https://www.dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#a2fdbf33bdf8531ce0e4febf656cd92da you would need to do for each component of stress and st

[deal.II] Re: Help in value_list

Thanks a lot for the reference. I got the point. Please help in one more thing. While defining the virtual member functions, e.g., RightHandSide::value, at the end of the function declaration "const" is added. Why is it so? Can you please give a reference? Thanks Atanu On Thursday, November

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

> > Dear Prof. Bangerth, > Many thanks for your suggestions and help! > I will try this with your reference of Aspect. I hope I can figure out this > smoothly. > Thanks again! > Regards, > Jack > 在 2016年11月3日星期四 UTC+8下午10:54:44，Wolfgang Bangerth写道： > > On 11/03/2016 03:32 AM, Denis Davydov wrote

Re: [deal.II] DWARF version?

> > /work/00434/eijkhout/trilinos/installation-12.2.1/lib/cmake/STK/STKConfig.cmake > says that this file needs tot exist > /work/00434/eijkhout/trilinos/installation-12.2.1/lib/cmake/STK/../STKDoc_tests/STKDoc_testsConfig.cmake > Can you check if it is there? Ok, I fixed this one by editing o

[deal.II] Re: Function to output stress and strain

Lisa, I solved a problem using deal ii and now in the out put step I have problem > exporting nodal stress and strain values. > in each cell I can calculate stress and strains at Gauss points but I > don't know how to calculate the stress and strains at nodes because they > are not continues. >

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

Thanks for replying me ! I moved to squared domain to avoid other possible error when I imply boundary condition, but I still have same manufactured solution e^(-x^2-y^2) I have to say that I have some changed error curve after moving to new mesh domain. Which is slightly better, but still uns

[deal.II] Function to output stress and strain

Hello, I solved a problem using deal ii and now in the out put step I have problem exporting nodal stress and strain values. in each cell I can calculate stress and strains at Gauss points but I don't know how to calculate the stress and strains at nodes because they are not continues. Can you

[deal.II] Re: Help in value_list

You're welcome. Since you say that you're a C++ beginner, I'm not sure how much you know about polymorphism: just in case you need some further clarification, here's an example that highlights what's actually going on here with the base and

[deal.II] Re: Help in value_list

Ok... Thank you Jean. -- 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 unsubscribe from thi

[deal.II] Re: Help in value_list

Dear Atanu, The RightHandSide class in step-7 derives from the Function class. This in turn defines the value_list

[deal.II] Help in value_list

Hi, This may be a trivial issue, but as a beginner in c++, I am unable to to understand the following. Someone please explain. In step-7, there is no member function named "value_list" is defined for the "right_hand_side" object. Then how is it used in the assemble_system. Thanks Atanu --

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

On 11/03/2016 03:32 AM, Denis Davydov wrote: I see. In serial mode all is doable because u can go from current cell to neighbors to quadrature data and do whatever. For the case of MPI you would still need to somehow ship the data to ghost cells. Because at some point neighbors of locally owned

Re: [deal.II] Postprocessing output of TensorFunction<1,dim>

On 11/03/2016 07:54 AM, Julian Dorn wrote: template class AdvectionField : public TensorFunction<1,dim> { // define 2D velocity vector field } How to do its Paraview/GMV-output? The easiest way is probably to create an FESystem(FE_Q, dim) (i.e., a finite element with as many c

[deal.II] Postprocessing output of TensorFunction<1,dim>

hello everybody, in (taken from the example 9): template class AdvectionField : public TensorFunction<1,dim> { // define 2D velocity vector field } How to do its Paraview/GMV-output? Thank you! -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum o

[deal.II] Re: Visit error...?

Oh, I see I'm sorry I thought about it in the wrong way. Thank you very much your kind reply. 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

[deal.II] Re: Visit error...?

Kyusik, On Thursday, November 3, 2016 at 5:51:20 AM UTC-4, hanks0...@gmail.com wrote: > > > And, I also calculated the {l2_norm of (exact_solution - numeric_solution) > / sqrt(#cells)} (I divided l2norm by sqrt(#cell) Because I want to know > kind of average value of difference between two solu

[deal.II] Re: Visit error...?

Arndt, First of all, Thank you very much for your reply. > What do you mean by "average value" here? You probably just wnat to > compute the error in the L2-norm, isn't it? > I'm sorry this is not average, it is RMS value. In deal.ii, from step7, L2 norm is calculated like this... L2_n

[deal.II] Re: Visit error...?

Kyusik, Am Donnerstag, 3. November 2016 10:51:20 UTC+1 schrieb hank...@gmail.com: > > Hi all, > > I'm trying to calculate absolute values of (exact_solution - > numeric_solution) and visualize it using Visit. > > So, I added the following in my code. > > error = solution; > error -= exact_solut

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

Hi Denis, Thanks. Yes, that’s right. For distributed MPI parallel computing, this is a big problem. I have just learned deal.II for less than three months, so I’m still new. And also I’m not working on numerical mathematics, so it is difficult for me to solve such problem but I will try.

[deal.II] Re: Visit error...?

I'm sorry, I missed the sqrt(#cell) in the calculation of L2norm I did std::cout<<"l2norm(error)="

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

Hi Jack, One more approach: you first loop over all locally owned cells and determine the local flag (tru/false). Then you still need to ship this to ghost cells on other MPI processors. On doing that you can do your assembly as usual. Kind regards, Denis > On 3 Nov 2016, at 10:12, Jack wrot

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

I see. In serial mode all is doable because u can go from current cell to neighbors to quadrature data and do whatever. For the case of MPI you would still need to somehow ship the data to ghost cells. Because at some point neighbors of locally owned cells will be ghost cells. I think it is do

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

Dear Denis, Thanks for your replies. I do not have some data attached to ghost cells. But I have a variable that changes with time and its change depends on the data of neighbor cells. For example, const double value1 = 1, value2 =-1; for (cell = tr.begin_active(); cell != tr.end();

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

Hi Jack, Could you please elaborate why you need to read the data on ghost cells? The CellDataStorage and TransferableQuadraturePointData classes are designed to do assembly, where we always loop over *locally owned cells* and calculate local matrices to be assembled. If the ownership of cells

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

Dear Jean， Thanks so much for your suggestions. I’m really very happy for such new development. I understand that using the CellDataStorage and TransferableQuadraturePointData classes is better than user_pointer when the mesh requires coarsening and refining. *Yet I do not understan