Re: [deal.II] Getting RHS values at nodes with DBC

_owned_dofs.is_element(d) && constraints.is_constrained(d) > && !hanging_node_constraints.is_constrained(d)) > const double value = unconstrained_rhs_vector(d); > > Best, > Jean-Paul > > On 01 Aug 2018, at 05:33, RAJAT ARORA > > wrote: > > Hi

Re: [deal.II] Shape function derivative wrt to a different cell

angerth wrote: > > On 03/07/2018 12:19 PM, RAJAT ARORA wrote: > > > > I want to do the following thing: I have a scalar quantity available > > only at the Gauss points of a cell. (cell = 2d rectangular element) > > I want to get the derivative of this quantity at the cente

[deal.II] Shape function derivative wrt to a different cell

Hello all, I am using deal.ii to solve a 2d solid mechanics problem. I want to do the following thing: I have a scalar quantity available only at the Gauss points of a cell. (cell = 2d rectangular element) I want to get the derivative of this quantity at the center of the cell. One way to do thi

Re: [deal.II] Deal.ii installation problem:- Step 40 runtime error

. On Monday, January 22, 2018 at 11:09:41 PM UTC-5, Wolfgang Bangerth wrote: > > On 01/22/2018 08:48 AM, RAJAT ARORA wrote: > > > > Running with PETSc on 2 MPI rank(s)... > > Cycle 0: > > Number of active cells: 1024 > > Number of degrees of

[deal.II] Deal.ii installation problem:- Step 40 runtime error

Hello all, I recently installed deal.ii using candi but after everything successfully finishes, I am not able to run step-40 on more than 1 processor. For, 1 processor it runs fine. To install everything, I used candi and and loaded following modules. 1) psc_path/1.1 2) slurm/17.02.5

[deal.II] Documenting code for Code gallery

Hello all, I have written a code for a small strain elasto-plastic solid mechanics problem in 2d using deal.ii. I wanted to contribute it to code gallery and have gone through the instructions given here: http://www.dealii.org/code-gallery.html#instructions I wanted to know how can I document

[deal.II] Re: Moving Mesh when using Fe_Q(p) p>1

node and then calculate the Jacobian and other information? Thanks. On Monday, December 11, 2017 at 8:39:20 AM UTC-5, Bruno Turcksin wrote: > > Hi, > > On Sunday, December 10, 2017 at 9:55:49 PM UTC-5, RAJAT ARORA wrote: > >> My question is: >> Is the mesh movement inde

[deal.II] Moving Mesh when using Fe_Q(p) p>1

Hello all, I am using the following function to move a p:d:triangulation. template void FEM:: move_mesh (const TrilinosWrappers::MPI::Vector &displacement) const { std::vector vertex_touched(triangulation.n_vertices(), false); for (typename DoFHandler::active_cell_iterator cell = dof_h

[deal.II] Re: Moving vertices of parallel triangulation breaks periodic face pair match

Hello Sambit, Can you try doing this? Also move the vertices of the ghost cell and avoid calling dftPtr->triangulation.communicate_locally_moved_vertices(locally_owned_vertices); When I tried to use triangulation.communicate_locally_moved_vertices(locally_owned_vertices) last time, somethin

Re: [deal.II] New Assert class (Discussion)

I_exceptions::internals::abort_on_exception,\ > __FILE__, __LINE__, __PRETTY_FUNCTION__, #cond, #exc, > exc); \ > } > #else > #define ModerateAssert(cond, exc) > \ > {} > #endif > > > > On Sun, Dec 3, 2017 at 5:17 PM, RAJAT A

[deal.II] New Assert class (Discussion)

Hello all, This is more of a discussion than a question. I was wondering if there is a need for a new class which does something like Assert (Lets call it myAssert for now). Its function is as defined below. While developing an application code, once I am confident that I am not making a mist

[deal.II] Re: compile master without examples using Candi

. -D DEAL_II_COMPILE_EXAMPLES=OFF" This turns off automatic compilation of examples when building master. On Friday, December 1, 2017 at 8:48:14 PM UTC-5, RAJAT ARORA wrote: > > Hello, > > From one of the posts, I came to know that > DEAL_II_COMPILE_EXAMPLES=ON by default

[deal.II] compile master without examples using Candi

Hello, >From one of the posts, I came to know that DEAL_II_COMPILE_EXAMPLES=ON by default for master. I am installing deal.ii master branch using candi. I do the build myself, I can pass that flag = OFF during the "cmake ." step but can anyone please tell me, how can I switch it off if I am us

[deal.II] Re: Questions related to FE_Nothing

o use triangulation.save() or tria.load() for serialization and deserialization. How can I use checkpoint/restart with parallel:shared:tria. Thanks. On Friday, November 17, 2017 at 6:38:25 PM UTC-5, RAJAT ARORA wrote: > > Hello all, > > I am using deal.ii to solve a system of Hamilton Ja

[deal.II] Questions related to FE_Nothing

Hello all, I am using deal.ii to solve a system of Hamilton Jacobi equations. I need to use FE_Nothing element in a part of my domain. I read a couple of posts in this forum regarding this but I have some questions. I will appreciate if anyone can answer these. Since I need to use MPI, and I

Re: [deal.II] adaptive mesh refinement doubts

Hello Professor, Thanks a lot for your earlier answers. I implemented all the details and it works like a charm. I still have a few questions that I want to ask. The questions also have solutions/suggestions. Sorry for being too verbose. Also, I am not sure if this should be a separate post so

[deal.II] adaptive mesh refinement doubts

Hello everyone, I am using deal.ii to solve a 3D solid mechanics problem. My code uses Parallel::Distributed::Triangulation and Petsc Wrappers. Until now, the code does not support AMR, however, I want to implement it in the code. For the last couple of days, I am looking at the example codes (

[deal.II] Line Integral of function

Hello everyone, I am using deal.ii to solve a 3D solid mechanics problem. My code uses PetSc and P4est. The geometry is a cubical box [0 1] X [0 1] X [0 1] with a uniform rectilinear mesh in the domain. I have obtained the scalar solution vector u in the domain. I need to find the line integr

Re: [deal.II] fixing one component of solution to the same value

Hello Daniel, Thank you so much for pointing this out. I have been stuck here for a long time and could not figure out. I understood that the sparsity pattern class was some how failing to know that space needs to be allocated for coupling with *the_dof even though it was not locally relevant

Re: [deal.II] fixing one component of solution to the same value

Hello Daniel, Thanks for the reply. I am actually doing what you suggested. Please look at the updated code that I had posted earlier and works with deal.ii 8.5. These lines are already there in the code. Surprisingly, such an error is still showing up. I am not sure what is still causing thi

Re: [deal.II] fixing one component of solution to the same value

Hello Professor, I am now using deal.ii 8.5 and I am attaching the updated code. The error message still remains the same. On Monday, April 17, 2017 at 4:43:35 PM UTC-4, RAJAT ARORA wrote: > > > Hello Professor, > > I have created a test case to show this wierd behaviour (b

Re: [deal.II] fixing one component of solution to the same value

: An error with error number 63 occurred while calling a PETSc function On Wednesday, April 12, 2017 at 7:09:17 PM UTC-4, Wolfgang Bangerth wrote: > > On 04/12/2017 05:03 PM, RAJAT ARORA wrote: > > > > > > I am

Re: [deal.II] Deal.ii installation Error

modern MPI and compiler you can try? We > have had many problems with intel compilers in the past, so I would > try a gcc with mpich or openmpi unless you are sure in what you are > doing. > > On Sat, Apr 15, 2017 at 5:19 PM, RAJAT ARORA > wrote: > > Hello Everyone, >

Re: [deal.II] Deal.ii installation Error

modern MPI and compiler you can try? We > have had many problems with intel compilers in the past, so I would > try a gcc with mpich or openmpi unless you are sure in what you are > doing. > > On Sat, Apr 15, 2017 at 5:19 PM, RAJAT ARORA > wrote: > > Hello Everyone, >

Re: [deal.II] fixing one component of solution to the same value

ement. And the dof 6 is for the x component and dof 2249 is for z component. Can you please guide me as to what is wrong and how should I proceed? Thanks. On Thursday, March 23, 2017 at 2:20:58 PM UTC-4, RAJAT ARORA wrote: > > Thanks Professor. > > > On Wednesday, March 15, 2017

Re: [deal.II] Fe_values->shape_grad() wrt to reference mesh

Hello Professor, The error is resolved when I used the template parameter "PETScWrappers::Vector". Thanks a lot for the help. I understand that this was not the most elegant way. To clarify, I had put the declaration in the header file of the class so that I am sure that the vector is not de

Re: [deal.II] Fe_values->shape_grad() wrt to reference mesh

his scope On Friday, March 31, 2017 at 11:15:54 PM UTC-4, RAJAT ARORA wrote: > > Hello, > > Since I am moving the mesh physically, the current coordinates are of the > current configuration. I am trying to make a new fe_values object with > mapping as shown below but I am getti

Re: [deal.II] Fe_values->shape_grad() wrt to reference mesh

Hello, Since I am moving the mesh physically, the current coordinates are of the current configuration. I am trying to make a new fe_values object with mapping as shown below but I am getting an error An error occurred in line <114> of file in function virtual dealii::Subscriptor::~Subsc

[deal.II] Fe_values->shape_grad() wrt to reference mesh

Hello all, I am using deal.ii to solve a 3D solid mechanics problem. I am using p:d:triangulation and Petsc in the code. I am updating my mesh after every time-increment. So, cell->vertex(v) contains the current coordinates of the mesh. I know the coordinates of the cell (and whole mesh) at t

Re: [deal.II] fixing one component of solution to the same value

Thanks Professor. On Wednesday, March 15, 2017 at 3:35:42 PM UTC-4, Wolfgang Bangerth wrote: > > > > I am using deal.ii to solve a 3D solid mechanics problem which uses > > p::d::triangulation. > > > > I am solving equilibrium equations to solve for the displacement in the > domain. > > The

[deal.II] fixing one component of solution to the same value

Hello all, I am using deal.ii to solve a 3D solid mechanics problem which uses p::d::triangulation. I am solving equilibrium equations to solve for the displacement in the domain. The domain of the problem is a cylinder with the z-axis aligned along the axis of the cylinder. To avoid the war

[deal.II] Re: Using constraints for already assembled RHS

system_rhs); } On Monday, September 12, 2016 at 1:35:52 PM UTC-4, RAJAT ARORA wrote: > > Hi Daniel, > > Thanks for the reply. > > This indeed is one way of doing it. > > I was wondering if there is any other way in which this can be done. > > Than

[deal.II] Re: Question regarding P:D:Tria

Hello Daniel, Thanks for the prompt reply. Yes, I thought about this as deal.ii uses P4est which does not support anisotropic refinement. I wanted to confirm that this was indeed the case and I am not missing something. Thanks again for the clarification. On Sunday, March 5, 2017 at 7:48:56

[deal.II] Question regarding P:D:Tria

Hello all, I am using deal.ii along with PetSc and P4est to solve a 3D solid mechanics problem. I am using the GridGenerator::hyper_rectangle function to create a rectangular grid with 1000 X 500 X 1 (0.5 million) elements. My question is: - since the coarsest mesh contains 0.5 million element

[deal.II] Re: How to output a single scalar in a parallel code

Hello, Please look at the solution by Timo. In your case, you can do something like. if (Utilities::MPI::this_mpi_process(mpi_communicator) == 0) { std::ofstream myfile; myfile.open ("resultant_strass.txt"); myfile<< resultant_stress< > Dear Rajat, > > >> ... which I then write in a file on

[deal.II] Re: How to output a single scalar in a parallel code

Hello, I am solving a similar problem. What I do is to sum the contributions from all the processors, and then use Utilities::MPI::Sum() to get the overall force which I then write in a file on a master process, For the contribution of a single cell, you have to integrate stress X normal X are

[deal.II] Query using hdf5 with deal.ii

hat contains the data from all the gauss points spread over multiple processors ? Or can this data be included while doing triangulation.save() which saves triangulation and data that is associated with it using Parallel solution transfer class ? Thanks for the help. -Rajat Arora PhD. Stud

Re: [deal.II] Problem with MUMPS Solver

Hello Timo, Thanks for your detailed reply giving many possible reasons. There are no constraints on the equation. It is a first order time dependent problem, so there are just initial conditions that are taken care of separately. The code is for 3D brick elements with 9 dofs per node. I just u

[deal.II] Problem with MUMPS Solver

Hello all, I am using deal.ii to solve a 3D solid mechanics problem along with Petsc and P4est. The problem has 9 Dofs per node. I run my code on stampede. This error may or may not be directly related to an issue with deal.ii (may be sparsity pattern ?) or its usage. But, given a diverse comm

Re: [deal.II] Doubt with set_boundary_id()

Hello Professor, The problem was because of the mesh. I regenerated the mesh and this time there is no such issue that arises. I think the options used to export and make the mesh earlier were wrong. Thanks for your help again. On Thursday, December 8, 2016 at 11:26:49 PM UTC-5, RAJAT ARORA

Re: [deal.II] Doubt with set_boundary_id()

oFHandler, 3>(dealii::DoFHandler<3, 3> const&, unsigned int, std::set, std::allocator > const&, dealii::ConstraintMatrix&, dealii::Mapping<3, 3> const&) #2 ./fdm: BugFinder<3>::BugFinder(int) #3 ./fdm: main ---

[deal.II] Doubt with set_boundary_id()

Hello all, I am using deal.ii to solve a 3D solid mechanics problem which uses Parallel::Distributed::Triangulation. I have a question related to setting the boundary indicators on the triangulation when no automatic partitioning is allowed. To set the boundary ids, 1) Should I only loop ove

Re: [deal.II] Help with triangulation::face_iterator

Thank you Professor. That worked but I am not sure why a reference didn't work but a const. reference did. On Tuesday, November 22, 2016 at 7:31:47 PM UTC-5, Wolfgang Bangerth wrote: > > > Rajat, > I don't know the exact cause of the problem, but... > > > Can we pass Triangulation::face_iter

[deal.II] Re: Move vertices in P::D::Triangulation

Thanks Daniel, I implemented this and it solved the issue that I was having. Just a small question, instead of using Triangulation::communicate_locally_moved_vertices , there is no harm in moving all locallyOwnedCells and ghost cells by one self. Right ? On Thursday, November 17, 2016 at 7:15

[deal.II] Help with triangulation::face_iterator

Hello, Can we pass Triangulation::face_iterator by reference ? I am trying to call the function template types::boundary_id Fdm::get_face_boundary_id(typename Triangulation ::face_iterator &f_it) { } by using typename DofHandler::cell_iterator cell = dof_handler.begin(); for( face_num =

[deal.II] Move vertices in P::D::Triangulation

Hello all, I am using parallel::distributed::triangulation, P4est, and Petsc to solve a 3D solid mechanics problem using deal.ii I am moving my mesh after every time step using the below written code. I have 2 questions regarding this. template void PlasticityContactProblem:: move_mesh (c

Re: [deal.II] Multiplie mpi instances of code

, Wolfgang Bangerth wrote: > > On 10/28/2016 11:49 AM, RAJAT ARORA wrote: > > > > I think that openmpi and Mpich are clashing now. (I didn't even know > > that both were installed upto now). > > > > I think, if I can disable openmpi somehow, I will be in goo

Re: [deal.II] Multiplie mpi instances of code

I am surprised why it is not working. I can't recall what has changed. I don't remember installing any new libraries. It was working till Monday. Also, I have installed petsc with --download mpich tag but donot have any other mpi installation. and more importantly, it was working until monday.

Re: [deal.II] Multiplie mpi instances of code

, 2016 at 11:48:51 AM UTC-4, Wolfgang Bangerth wrote: > > On 10/28/2016 09:47 AM, RAJAT ARORA wrote: > > > > Yes, you are right :). > > When running with 4 mpi processes, all threads print their rank as 0. > > > > Why is it broken all of a sudden. Wh

Re: [deal.II] Multiplie mpi instances of code

10/28/2016 09:34 AM, RAJAT ARORA wrote: > > > > To run the code on n processes, I used to run the command > > > > mpirun -np ./ > > > > But, now, when I run it using this command, it runs n programs with n > > mpi processes each. > > > &g

[deal.II] Multiplie mpi instances of code

Hello all, I am working on a 3D solid mechanics problem using deal.ii. To run the code on n processes, I used to run the command mpirun -np ./ But, now, when I run it using this command, it runs n programs with n mpi processes each. It is like I have executed mpirun -np ./ command n times

[deal.II] Finite Difference with deal.ii and Parallel Distributed Triangulation

Hello all, I am curious to know if it is feasible to use the P4est mesh structure with deal.ii to do a finite difference simulation. For a 2D structured grid, in order to do finite differences, we need to know the index (or dof index) of vertices to the left as well as to the right of a give

[deal.II] Doubt with data output in parallel

Hello all, I have a small question related to data output in parallel. My program uses Parallel::distributed::triangulation. If I want to output a data array associated with cell, which one of the following is the correct way of doing it ? And why ? (This cell data type vector is to be viewed i

[deal.II] Re: Using constraints for already assembled RHS

Hi Daniel, Thanks for the reply. This indeed is one way of doing it. I was wondering if there is any other way in which this can be done. Thanks again. On Monday, September 12, 2016 at 5:04:47 AM UTC-4, Daniel Arndt wrote: > > Rajat, > > The first equation is the usual equilibrium equation.

[deal.II] Using constraints for already assembled RHS

equation. Thanks, - Rajat Arora -- 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"

Re: [deal.II] Applying Multi Point Constraints

Respected Prof. Wolfgang, I really appreciate your taking the time to reply. Thank you very much. Yes, by middle element I mean the central layer of the mesh. Yes you are right. It turns out this is the problem. I was creating a hyper_rectangle using grid generation class. Then 4 and 5 correspo

Re: [deal.II] Applying Multi Point Constraints

Respected Prof. Wolfgang, I really appreciate your taking the time to reply. Thank you very much. Yes, by middle element I mean the central layer of the mesh. Yes you are right. It turns out this is the problem. I was creating a hyper_rectangle using grid generation class. Then 4 and 5 correspo

[deal.II] Applying Multi Point Constraints

Hello all, I am trying to apply multi point constraints to a 3D solid mechanics problem. My code used P4est, and PetSc. Lets say I have a mesh with 10 X 10 X 3 elements. I want to apply constraints such that all the dofs on the +z surface of the middle element are constrained to be equal to th

[deal.II] Doubt regarding constraints

Hello all, I have a question regarding use of constraints when using parallel programming. Suppose I want to constraint dof1 owned by processor 1 to be equal to dof2 owned by processor 2. I can do that by calling ConstraintMatrix constraints; constraints.reinit(locally_relevant_dofs); cons

[deal.II] Re: PETScWrappers::SparseDirectMUMPS solver error

this post, a new version of MUMPS, 5.0.2, has been released on 17th July 2016. On Thursday, July 14, 2016 at 4:25:13 AM UTC-4, Denis Davydov wrote: > > > On Wednesday, July 13, 2016 at 10:56:27 PM UTC+2, RAJAT ARORA wrote: >> >> Hello Denis, >> >> Thanks for your

[deal.II] Re: PETScWrappers::SparseDirectMUMPS solver error

, Denis Davydov wrote: > > Hi Rajat, > > On Tuesday, July 12, 2016 at 8:37:11 AM UTC+2, RAJAT ARORA wrote: >> >> Hello Jean, >> >> Thanks for the reply. >> I understand your concern for the lack of proper information but this is >> all I have. >&

[deal.II] Re: PETScWrappers::SparseDirectMUMPS solver error

27;re using > quite an old version of deal.II, and the root of the problem may have > already been fixed in a later version. > > Regards, > Jean-Paul > > On Tuesday, July 12, 2016 at 3:38:54 AM UTC+2, RAJAT ARORA wrote: >> >> Hello all, >> >> I am r

[deal.II] PETScWrappers::SparseDirectMUMPS solver error

Hello all, I am recently encountering an issue with MUMPS solver. The error is coming from the dealii::PETScWrappers::SparseDirectMUMPS class in these lines. dealii::PETScWrappers::SparseDirectMUMPS solver(solver_control, mpi_communicator); solver.solve (system_matrix,

[deal.II] PETScWrappers::SparseDirectMUMPS

Hello all, I am recently encountering an issue with MUMPS solver. The error is coming from the dealii::PETScWrappers::SparseDirectMUMPS class in these lines. dealii::PETScWrappers::SparseDirectMUMPS solver(solver_control, mpi_communicator); solver.solve (system_matrix,

[deal.II] Re: Triangulation.save() not working as expected

e > <https://github.com/dealii/dealii/issues/2661>, but we should wait for > confirmation from someone with more knowledge on how the triangulation > class works. > > Regards, > J-P > > On Friday, July 8, 2016 at 9:56:07 AM UTC+2, RAJAT ARORA wrote: >> >

Re: [deal.II] Triangulation.save() not working as expected

Hey Timo, Thanks for your prompt reply. Here is the documentation where in the serializing section it says that serializing stores vertices and its location. I am confused if it does store vertex locations now. Also, I

[deal.II] Triangulation.save() not working as expected

Hello all, I am observing some strange behavior with Parallel::distributed::triangulation::save(const * char) function. I make a triangulation and then move its vertices and then save it to a file. After saving the triangulation and I load it in a different function. However, this loaded t

[deal.II] Doubt in Exceptions Documention

Hi, In the Exception and Assertions documentation, the documentation for the function ExcNumberNotFinite

[deal.II] Re: Error durind checkpoint / restart using parallel distributed solution transfer

that you > found the issue, and this will be helpful to know. Perhaps we can > mention/reinforce it in the documentation. > > Regards, > J-P > > On Sunday, July 3, 2016 at 9:22:06 PM UTC+2, RAJAT ARORA wrote: >> >> Hello, >> >> This is because of the v

[deal.II] Re: Error durind checkpoint / restart using parallel distributed solution transfer

Hello, This is because of the vector locally_owned_x. As per the documentation, during the serialization, the vector should be ghosted. If I make the change, the code works fine. On Sunday, July 3, 2016 at 2:52:57 PM UTC-4, RAJAT ARORA wrote: > > Hello all, > > I am tryin

[deal.II] Error durind checkpoint / restart using parallel distributed solution transfer

Hello all, I am trying to use the checkpoint / restart code with the help of the example given in the documentation here I am getting a long error while saving a parallel vector and parallel dis

[deal.II] Restoring Gauss Point History during checkpoint/restart

ead it during a checkpoint/restart. I am using parallel distributed triangulation. The Quadrature point history is a std::vector<> of a struct PointHistory just like in step 18 <https://dealii.org/8.4.1/doxygen/deal.II/step_18.html> . Is there any way to do this ? Regards, Rajat Aror

Re: [deal.II] About MPI vectors, ghosted vectors

Hi Ehsan, On Wednesday, June 29, 2016 at 3:28:40 PM UTC-4, Ehsan Esfahani wrote: > > Dear users, > > I'm asking my question here because I think it's somehow related to this > post, and I don't think it would be necessary to open a new topic. Could > you please tell me what are differences betw

[deal.II] Linking Deal with Atom or QT (Text Editor)

Hello all, I have a simple question. I have been using Eclipse to write code and get auto complete suggestions from deal.ii functions and classes. However, Eclipse uses a lot of memory. Is it possible to include Deal.ii in the atom text editor so that we can get auto-complete suggestions from

[deal.II] Re: Doubt in FEvalues

t data is completely reevaluated. However, I > can't be certain as I'm not familiar with the inner working of these > classes. > > On Wednesday, June 1, 2016 at 2:20:35 PM UTC+2, RAJAT ARORA wrote: >> >> Hello Jean, >> >> Thanks for your help. It wou

[deal.II] Re: Doubt in FEvalues

up a bug report on your behalf here > <https://github.com/dealii/dealii/issues/2661>. So, in summary, you'll > just have to wait for a bit until the problem is fixed, or use at least one > level of global refinement :-) > J-P > > On Wednesday, June 1, 2016 at 2:04:19 A

[deal.II] Re: Doubt in FEvalues

idTools: 0.75 Calculated: 1 > > However, once you perform some refinement then this error would be picked > up in debug mode (are you running your tests in debug mode?). > > Does that help? If you solve the issue, then I'd be interested to know > what the problem w

[deal.II] Re: Doubt in FEvalues

(i.e. > the spatial gradient). However, you will note that in step-18, only the > incremental displacement gradient is computed as the stresses are > systematically incremented using this result. > > On Tuesday, May 31, 2016 at 2:54:32 AM UTC+2, RAJAT ARORA wrote: >> >

[deal.II] Doubt in FEvalues

Hello all, I have a question about FEValues object that we use in the computation of element stiffness and other places. If I declare fe_values object and get then move the mesh and then try to call fe_values.get_function_gradients() will this give me gradients on the deformed mesh or the mesh

[deal.II] Re: Doubt in Mesh Movement Step 42

, > Daniel > > Am Samstag, 28. Mai 2016 22:34:41 UTC+2 schrieb RAJAT ARORA: >> >> Hi Daniel, >> Thanks for your prompt reply. >> >> While in case of mesh movement, the vertices are independent. However, if >> want to add something to the solution vector at dofs

[deal.II] Re: Doubt in Mesh Movement Step 42

Hi Daniel, Thanks for your prompt reply. While in case of mesh movement, the vertices are independent. However, if want to add something to the solution vector at dofs of a vertex shared between processors, This mesh movement code wont be able to do that, It will add that value twice for the d

[deal.II] Assembling RHS coming from different equation with same triangulation but different dofhandler

Hello all, I am using deal.ii to solve two vector valued equations by FEM:- one for the velocity field "*v*" and the other for a vector field "*f*". I created different classes for solving these two equations separately. I am using same parallel::distributed::triangulation for both the equation

[deal.II] Doubt in Mesh Movement Step 42

Hi all, I have a doubt regarding the mesh movement in step42.cc Below is the code taken from the example. Lets say cell 1 and cell 2 are two adjacent cells of a fem mesh (share 2 nodes / vertex) and they live on processor 1 and 2 respectively. Processor 1 only keeps track of vertices that it h