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

On 11/06/2016 11:15 AM, Hamed Babaei wrote: I am attaching the sequential and parallel codes as well as parameters.prm file. It would be appreciated if you could compare tangent_matrix (the system_matrix of the elasticity part of the problem) of parallel and sequential codes to find out where I

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

Hi, everyone I’m not familiar to ASPECT, could someone tell me how to get the data of ghost cells? In the ASPECT, I find that it is implemented through the following functions to obtain ghost particles in ~/aspect/source/particle/world.cc template std::map World::get_subdomain_id_to_neighbor

Re: [deal.II] How to use a function avoiding to pass some of its arguments

On 11/08/2016 05:42 PM, Hamed Babaei wrote: Since I use Newton method and I want to apply boundary displacement as constraint, I should apply constraint only in zero Newton iteration. So I can not do so in setup_system but after assembly. I followed step-44 for that. Then keep two ConstraintM

Re: [deal.II] How to use a function avoiding to pass some of its arguments

Dear Wolfgang, Since I use Newton method and I want to apply boundary displacement as constraint, I should apply constraint only in zero Newton iteration. So I can not do so in setup_system but after assembly. I followed step-44 for that. Thanks, On Tuesday, November 8, 2016 at 6:28:17 PM UT

Re: [deal.II] How to use a function avoiding to pass some of its arguments

On 11/08/2016 05:19 PM, Hamed Babaei wrote: DoFTools::make_sparsity_pattern (dof_handler, sp, constraints, false, Utilities::MPI::this_mpi_process

[deal.II] How to use a function avoiding to pass some of its arguments

Hi all, I want to use the following function: DoFTools::make_sparsity_pattern (dof_handler, sp, constraints, false, Utilities::MPI::this_mpi_process

Re: [deal.II] Re: using .vtk file as initial condition

It seems that using block_write() and block_read would be easier to store data and reuse them as initial conditions on identical geometry and mesh compared to using boost:serialization. In sequential computations, that is probably true. (i) Suppose I want to store my 'Solution' vector to the z

Re: [deal.II] Re: Learning Finite Element Method

Hi Divya, The first resource that I would recommend is the online lectures of Prof. Bangerth. The second resource is the book "A First Course in the Finite Element Method" by William B. Bickford. Oded -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options,

Re: [deal.II] Re: using .vtk file as initial condition

Hello all, It seems that using block_write() and block_read would be easier to store data and reuse them as initial conditions on identical geometry and mesh compared to using boost:serialization. However, I have some issues and quires (I do not use parallelization in my code). I shall be than

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

Dear Wolfgang, below is the code, where the function template void HeatEquation::tensor1_output() { ... } performs the gmv-output of template class AdvectionField : public TensorFunction<1,dim> { // define 2D velocity vector field } Best wishes, JD ***

Re: [deal.II] Re: Question about Step 44 shape function gradient computations

On 11/08/2016 01:36 AM, Jean-Paul Pelteret wrote: You're welcome. If you go through the formulation and still don't have full clarity, then let me know. I do have this derivation documented somewhere... Claire -- If, on the other hand, you do figure it out, write it up into a one-page latex d

Re: [deal.II] Re: Learning Finite Element Method

The following course introduces FE and deal.II https://www.youtube.com/playlist?list=PLJhG_d-Sp_JHKVRhfTgDqbic_4MHpltXZ A > On 8 Nov 2016, at 12:41, Jean-Paul Pelteret wrote: > > Dear Divya, > > Its incredibly diffic

Re: [deal.II] Meshworker

Sudarshan, there is a Triangulation::cell_iterator and a DoFHandler::cell_iterator and dinfo1.cell returns the first one while you are asking for the first one. Maybe have a look at [1] for the difference between the two. What exactly do you need the cell_iterator for? Best, Daniel [1] https:/

[deal.II] Re: Learning Finite Element Method

Dear Divya, Its incredibly difficult to answer your question because the answer is a deeply personal one. Instead of embarking on some long diatribe, I will simply encourage you to do your own research to find a set of resources that work for you. Understandably the lecture videos provided b

[deal.II] Re: Question about Step 44 shape function gradient computations

Dear Claire, You're welcome. If you go through the formulation and still don't have full clarity, then let me know. I do have this derivation documented somewhere... Regards, Jean-Paul On Tuesday, November 8, 2016 at 9:29:23 AM UTC+1, Claire wrote: > > Dear Jean-Paul, > > Thank you very much fo

[deal.II] Re: Question about Step 44 shape function gradient computations

Dear Jean-Paul, Thank you very much for your answer. I still have to figure out some of the calculations you suggested me to do, but I got the concepts and the reason why the gradients are defined on the deformed configuration rather than on the initial one. Best, Claire Le lundi 7 novembre