[deal.II] Evaluating the gradient of the solution at the centroid of the element

I have gone through some of the tutorials previously and thought it was about time I actually tried to implement something on my own in the library. I am attempting to implement a finite strain code using the F-Bar method. Essentially this is a single field formulation which requires the calc

[deal.II] Re: Postprocessing on Vector Valued Problem

You were right Oh I just got the result! I didn't know that dof_handler can automatically matching the sequence of the cell. Thank you for your time Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/foru

[deal.II] Re: Postprocessing on Vector Valued Problem

Jaekwang, On Wednesday, December 21, 2016 at 4:05:54 PM UTC-5, Jaekwang Kim wrote: > > > but How can I relate my 'cellwise_shear_rate' with dof handler so that I > finally able to see the values in 'vtk' format? > for now, it just seems that no more than a list of arrays > You mean a list o

Re: [deal.II] deal.II FE structure

Also another thing which gives me a headache is that why is the jacobian value twice the value it should be. For a rectangular bilinear lagrange element we can compute simply J by a/2 and b/2, but we would receive 1 and 1 for the above case where a=b=2. This doesn't correspond to what deal.II ou

Re: [deal.II] Postprocessing on Vector Valued Problem

On 12/21/2016 02:05 PM, Jaekwang Kim wrote: I tried to code following data post processing for my goal. template void ShearProblem::post_processing () { // std::cout << " Post Processing..." << std::endl; constMappingQ mapping (degree); QGa

[deal.II] Postprocessing on Vector Valued Problem

Hi, all. I have a question on post processing of solution values. My solution is block vector solution that is consist of two parts. First part is velocities(u_x,u_y) and pressure. My purpose is to calculate shear-rate at each dof (or at each cell) which only deals with dof of velocities par

[deal.II] Re: deal.II FE structure

Dear Seyed, You may also be interested in looking at the module on the interplay between finite elements and mappings , and the description to the mapping class. Note that th

[deal.II] Re: deal.II FE structure

Hi Seyed, On Wednesday, December 21, 2016 at 7:16:32 PM UTC+1, mohseng...@gmail.com wrote: > > Hi, > > I am trying to understand and learn the main structure of deal.II's > design, because I cannot work with a program until I do not understand it > fully. > Hence, I am checking a simple 2D cube

[deal.II] deal.II FE structure

Hi, I am trying to understand and learn the main structure of deal.II's design, because I cannot work with a program until I do not understand it fully. Hence, I am checking a simple 2D cube with hand-written FE solutions to see how things are stored in deal.II. The model I use is created direc

[deal.II] restart/checkpoint strategies in user applications/libraries

Hi all, This is a bit off-topic, but i would like to ask how do you prefer to implement restart/checkpoints in your user applications? I see two options: 1. Make restarting usable with the same input parameter file so that calculations will continue exactly from the point of restart. For exampl