Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

Hi I am solving a problem using Trilinos and parallel::distributed triangulation. It works well. I transfer solution data and old solution data during mesh

Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

> > > Thank you. Is this class available in the release version 8.3.0? > I set separate dofhandler (history_dof_handler) for FE_DGQ. I just > confusing because the history_field vector has different size respect to > the main solution vector and it also has different DoFHandler. I attached >

[deal.II] Re: How to solve each block separately when working with TrilinosWrappers::BlockSparseMatrix and AMG precondition.

Dear Arndt Thank you for answer. I want to solve the system during a time step using a split method. First assembling block(0,0) and solving the displacement field, then using the updated displacement, the second field must be solved and so block(1,1) of system matrix is determined. I

Re: [deal.II] Re: extract global node index (or DOF index) which belong to a specific boundary id

Hi I am trying to extract all degrees of freedom which are at the boundary and belong to specified components using DoFTools::extract_boundary_dofs() function and IndexSet for a parallel distributed triangulation. Now, how I can iterate over the IndexSet and how I can access to the Dof index

[deal.II] Re: Extract boundary DoFs using DoFTools::extract_boundary_dofs() and Indexset for parallel distributed triangulation

Please, ignore that this method will apply a value correspond to a common node between two adjacent element twice. I will correct that. I just want to be sure what I am calculating is correct or not. On Friday, July 1, 2016 at 12:05:24 PM UTC+4:30, dealii...@gmail.com wrote: > > Dear Daniel >

[deal.II] Re: Extract boundary DoFs using DoFTools::extract_boundary_dofs() and Indexset for parallel distributed triangulation

Dear Daniel My name is Pasha. Thank you for your response. I want to extract sum of values from system RHS ( which is a dealii:TrilinosWrappers::MPI::BlockVector) correspond to all degrees of freedom which are at the boundary. I am working with a parallel::distributed::Triangulation. I

[deal.II] Re: Extract boundary DoFs using DoFTools::extract_boundary_dofs() and Indexset for parallel distributed triangulation

Dear Daniel Thank you. You really help me. Best, Pasha On Friday, July 1, 2016 at 5:21:08 PM UTC+4:30, Daniel Arndt wrote: > > Pasha, > > First of all: The sum of values in the right-hand side vector only > corresponds to a sum of values of the represented right-hand side function, > if you

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Dear all Regarding using MappingQEulerian class in parallel, I did what you said as bellow: dealii::TrilinosWrappers::MPI::BlockVector rel_solution2(partition_relevant); rel_solution2 = solution; rel_solution2 *= -1.0; MappingQEulerian

[deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Thank you for the answer. I am working on a solid mechanical model with large deformations. The first RHS is the residual load vector (external load – internal load). The second RHS is the global external load vector. Since I use update Lagrangian formulations, the mesh is changed during

[deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Sorry, the local_rhs is available for the second global RHS during assemble process of the glbal matrix and the first global RHS. On Tuesday, January 17, 2017 at 1:26:29 PM UTC+3:30, dealii...@gmail.com wrote: > > Thank you for the answer. > > I am working on a solid mechanical model with large

[deal.II] Re: Nonlinear Material

Dear Lisa Take a look at to step-18, step-42 (A solver for an elasto-plastic contact problem, running on parallel machines) and Step-44 (the subject of this tutorial is nonlinear solid mechanics. Finite strain hyperelasticity based on a three-field formulation.), Best, Pasha On Wednesday,

[deal.II] Wall time vs number of cores for MPI run

Dear all I wrote a code using deal.II. Now, I want to compare effect of the number of cores on the simulation time. The code is developed based on step-42 from deal.II tutorial. The code is executed in Ubuntu 14.04, Intel core i5 and 4GB ram as, mpirun –n ./code parameter_file.prm

[deal.II] compute sign of the determinant of the global tangent stiffness matrix

Dear all I am working on a arc-length-Newton-Raphson solution algorithm. To derive the arc-length method, I allow incremental load factor ∆λ to become a variable and redefine the residual equation. I used cylindrical arc-length method. In this method the iterative load factor is normally

[deal.II] compute sign of the determinant of the global tangent stiffness matrix

Dear all I am working on a arc-length-Newton-Raphson solution algorithm. To derive the arc-length method, I allow incremental load factor ∆λ to become a variable and redefine the residual equation. I used cylindrical arc-length method. In this method the iterative load factor is normally

[deal.II] Re: compute sign of the determinant of the global tangent stiffness matrix

Hi Denis. I also don't know how the first method can be implemented. I read the above models in the book By Neto, Peric and Owen “Computational method for plasticity”. They have been computed the determinant of the stiffness matrix during solving the linear system using the classical frontal